基于体面积分方程的预条件技术.docVIP

南京理工大学 硕士学位论文 基于体面积分方程的预条件技术 姓名：曹清清 申请学位级别：硕士 专业：电子与通信工程 指导教师：陈如山；曹宏海 201203  基于体面积分方程的预条件技术 矩量法是电磁仿真中精度较高的方法，本文研究的体面积分方程方法是基于矩量法 的，可以用于计算混合目标的电磁散射特性。体积分算子即使在网格剖分很密的情况下 仍然具有良好的条件数，但面积分算子条件数较差。在计算复杂的金属、介质混合目标 时迭代收敛较慢，甚至不能收敛。基于这种情况，我们采用下列两种预条件来改善矩阵 的性态。 其次，本文研究了基于网格细分的多分辨基函数，并在此基础上构造了多分辨预条  ???? 硕士论文 ??????? ??? ??? ?? ????????? ??? ??????????? ?????? ??? ? ????????? ????瑃? ????? ????? ??? ? ????產? ?．??? ????? ????．?? ??? ?????．? ?????????????? ????? ?? ?????? ??????? ?????? ?????????? ?????? ????????． ??????? ? ???瑃? ??? ??? ????： ??琕??琾??????? ??????? ??????????琈??????????????琑?? Ⅱ  基于体面积分方程的预条件技术 ??本文的研究背景 随着雷达技术、隐身与反隐身技术、微波遥感和目标识别等技术的不断发展和应用 要求的不断提高，迫切需要对三维复杂目标的电磁散射特性进行分析。雷达散射截面是 描述目标电磁波散射强度的一个非常重要的量化指标。虽然电磁散射特性数据可以经过 实际测量获得，但是需要昂贵的测量仪器、理想的测量条件，同时还要耗费大量的人力、 物力和财力。某些情况下还需要较长的测量周期，总而言之，需要付出较高的代价。当 然通过理论建模和数值模拟同样可以获得结果。数值方法的优点是能够在不受外界条件 阵，但是未知量分布在整个传播空间中，即未知量的数目较大，必须对传播空间进行剖 分，容易产生网格色散，还必须辅助截断边界，通常这也是解的误差来源。以矩量法为 代表的积分方程数值方法需要在散射体区域中进行网格划分，未知量较少并且很容易控 制，并自然满足截断边界条件，但是缺点是产生的阻抗矩阵为满阵。矩量法出现于? 世纪?年代，虽然能够给出精确的结果，但是受客观条件的影响，只能计算电尺寸很 小目标的电磁散射特性。随着?年代快速多极子?儿。?的出现，数值方法的计算能力得到 了极大的提高，矩量法也得到了快速的发展。 至今为止，预条件技术已经大力发展并取得了显著的成果。基于稀疏矩阵的常见方 小迭代收敛步数不明显。对称超松弛预条件???技术相比对角预条件含有较多的矩阵  硕士论文 的信息，效果较对角预条件有所改善，但是仍然不是很明显且系数矩阵的对角线元素不 本文引入的体面积分方程方法【??【?】是基于矩量法的，此方法可以用来求解导体和 介质混合目标的电磁散射特性。因为矩量法是一种严格的数值方法，结果精确，所以经 体面积分运算的结果也较为精确。体积分方程作用在介质部分，面积分方程作用在金属 部分，所以体面积分方程对任何导体介质混合结构都适用，对模型没有特殊要求。对于 非均匀介质，体面积分方程也是适用的。我们假设介质体离散的每个小四面体单元内的 介电常数近似为一常数。但是我们都知道体面积分方程存在一定的弊端，就是需要对介 质体进行剖分，使得未知量较大、内存消耗较多，但是这可以通过引入多层快速多极子 但是这些快速方法都只是减少每一步迭代求解的计算复杂度和内存需求，并不能减少迭 代求解过程中的迭代求解的总步数。当阻抗矩阵的未知量较大或者网格离散密度趋近于 零时，不可避免地存在缺陷，使得积分方程的特征值趋近于零或者是无穷大，产生的阻 抗矩阵的条件数较大，矩阵方程迭代收敛较慢甚至不能收敛，所以出现了预条件技术。 但是传统的预条件技术都是单纯考虑阻抗矩阵的信息而不是考虑整个矩阵方程，从阻抗 矩阵中提出某些信息作为预条件矩阵，使得这些方法只对某些算例有效。基于这种考虑， ??????? 矩阵的性态。与单纯基于阻抗矩阵信息构造的上述预条件技术相比，本文中使用的这两 种预条件技术适用性更广泛，尤其是在网格剖分密度很大时预条件技术效果更加显著。 通过这两种方法，可以实现对阻抗矩阵方程快速、高效的求解。 ??????热嗽谖南譡?】将????乘式预条件技术应用到电场积分 方程???中，可以用于分析开放和闭合结构的目标，并且与现有的矩量法很容易结合。 ??????热嗽谖南譡?】又将此方法应用到磁场积分方程中。??年 面积分方程中面积分算子性态差的问题，西? ????乘式预