李 超，何思远＊，朱国强，邓方顺，陶 华，肖柏勋

1 武汉大学电子信息学院，武汉 430079

2 中国船舶重工集团第722研究所，武汉 430079

3 中国人民解放军92474部队

4 长江大学工程地球物理研究中心，湖北荆州 434023

1 引 言

目标与粗糙面复合电磁散射研究在雷达目标探测、目标识别以及微波遥感等领域有着广泛的应用[1-20］.求解粗糙面上方或下方目标的电磁散射问题必须充分考虑目标与粗糙面之间的相互作用.Li等[1］提出用广义前后向迭代算法与谱加速算法结合（GFBM／SAA）计算海洋粗糙面上方目标的电磁散射场.Deng等[2］借助传统UV算法，结合分层思想，发展得到适用于三维矢量波散射问题的MLUV算法，并成功将其应用到目标位于粗糙面上方和目标位于粗糙面上时的复合散射问题.对于埋地问题，相关文献[3-4］用分层媒质并矢格林函数分析了平面分层大地中埋入目标的散射特性.在此基础上，Zhang等[5］引入离散复镜像技术（DCIM），加快了分层媒质并矢格林函数的计算，提高了计算效率.徐利明等[6］则利用积分方程以及半空间并矢格林函数的快速算法对埋地目标矢量电磁散射进行了计算.在一些埋入体散射问题中，需要计及大地表面起伏特性对散射场的影响.对于呈复杂曲面的地表面，埋入体的散射问题将变得更为复杂.为求解该类问题，一些文献采用基于矩量法的快速算法计算粗糙地面埋入体的电磁散射问题，如El-Shenawee等[7-8］采用最陡下降的快速多级子方法（SDFMM），Zhang等[9］则采用了稀疏矩阵正则网格（SMCG）算……