论文部分内容阅读
随机粗糙表面、电大尺寸复杂目标,以及目标与粗糙面复合模型的散射计算,对雷达制导与截获技术、目标识别与特征提取、浅层地下目标勘探等具有十分重要的研究价值,日趋频繁的工程应用需求对数值模型的求解速度和效率提出了较高要求。本文主要研究随机粗糙面上目标电磁散射的计算方法,提出FBM-CG和KA-MoM等混合算法,在保证计算精度的前提下,大大节省了计算时间,为复杂环境中的目标散射计算提供了快速有效的解决方法。首先研究粗糙面散射的锥形入射波问题。在粗糙面散射的数值模拟中,为消除粗糙面的有限截断产生的边缘效应,往往引入锥形波入射。锥形波的宽度参数g和粗糙面的截断长度L直接影响数值模拟的存储量和计算量,以及数值结果的有效性。但是,锥形波参数的选择没有一个明确而有效的标准。本文分析锥形波满足的Helmholtz波动方程,从误差最小的角度推导出宽度参数g与入射角度θ_i的定量关系。结合随机粗糙面的相关特性,以及入射功率的有效截断等条件,提出粗糙面长度L的选择依据。数值结果验证了参数选择的有效性。其次研究体目标电磁散射的快速多极子算法(Fast Multipole Method,FMM)和多层快速多极子算法(Multi-Level Fast Multipole Algorithm,MLFMA)。多极子方法的出现极大地加快了矩量法求解的运算速度,使复杂电大尺寸目标的电磁散射分析成为可能。但是,多极子方法中的数值积分计算存在很大的误差,尤其是对节点单元分组较大的情况。本文分析了多极子算法中平面波展开的角谱积分,从采样定理的角度考虑,提出了更合适的数值积分标准,提高了多极子算法的计算精度。接着研究目标与粗糙面复合散射模型的高效建模与数值仿真。基于半空间Green函数的性质,推导出关于目标感应电流和粗糙面的差值感应电流的一组新的耦合积分方程,相比Johnson关于差场散射的数值模型,本文直接计算差值散射场,而无需分别求解有无目标两种数值模型。由于差场散射只包含目标自身的体散射以及目标与粗糙面之间的相互作用,而未计入粗糙面自身的面散射,从而消除了面-体联合散射中存在的量纲不一致问题,并且与入射锥形波的照明宽度无关。本文提出了目标的共轭梯度(Conjugate Gradient,CG)求解与粗糙面的前后向迭代计算(Forward Backward Method,FBM)相结合的快速互耦迭代算法,迭代过程通过更新目标和下垫粗糙面上的激励场来考虑目标与粗糙面之间的相互作用。指出粗糙面对目标的散射贡献主要来自对准目标的镜面方向上的一小段粗糙面,并据此加速计算。讨论了锥形波的宽度参数g以及粗糙面的截断长度L的取值问题,给出了相应的选择依据。结合Monte-Carlo方法,将这一方法应用于P-M(Pierson-Moskowitz)谱粗糙海面上方二维导体目标的散射计算和数值结果分析。本文进一步将这一算法推广到介质目标和多目标的情况,对多目标的情况提出选择每个目标的镜面方向上一小块粗糙面的散射贡献,以及切割各自合适的粗糙面长度用于迭代计算。提出多目标散射的并行计算方法,对每个目标采用不同的处理器并行处理,目标之间的耦合通过多个处理器之间交换数据实现。最后,针对粗糙面数值求解的存储量和计算量大的问题,提出将解析和数值两种方法有效结合的混合迭代算法,对粗糙面采用基尔霍夫近似(Kirchhoff Approximation,KA)解析计算,对目标采用矩量法(Method of Moment,MoM)和共轭梯度法(CG)数值求解,并数值验证了该混合算法的有效性和收敛性。由于粗糙面的KA计算只需对目标表面的一次积分,节省了大量的运算量和存储量,为目标与粗糙面的复合散射计算提供了新的可行方案。结合Monte-Carlo方法,数值分析了二维和三维随机粗糙面上不同形状和取向目标的散射特性。本文的研究工作为目标与随机粗糙面的复合建模与数值仿真,以及复杂环境中的目标识别和特征提取提供了一套快速有效的理论和计算方法。