基于线性相位RWG函数的矩量法矩阵元素的奇异性处理.pdfVIP

基于线性相位RWG函数的矩量法矩阵元素的奇异性处理 ． 常得江周后型洪伟 东南大学信息科学与工程学院毫米波国家重点实验室，南京210096 Email：djchang@emfield．org 摘要：在矩量法中，将事先正确预测的相位信息加入RWG函数构成线性相位RWG函数 (LP．RUG)，能在保证精度的前提下大大减少未知量的个数，从而降低数据的存储需求和计 算量，但在应用Galerkin法进行测试时需计算二重奇异积分，在应用LP～RwG函数作为展 开和测试函数时，积分核更是含有复杂的复指数函数，必需进行特殊地处理。本文在文献中 已有的关于RWG函数积分的奇异性提取技术的基础上提出了关于线性相位RWG函数奇异 性积分的一种高效处理方法，很好的解决了这一问题。 关键词：矩量法线性相位RWG基函数奇异值提取技术 1 引言： 函数之后【2J，矩量法在电磁辐射和散射问题中得到了更加广泛的应用。经典矩量法的数据存 储量和计算复杂度都很高，大大限制了矩量法的应用范围。因此，研究和发展基于矩量法的 快速算法成了当今计算电磁学领域里的一个热点。研究快速矩量法可以从不同的角度着手， 著名的多层快速多极子算法(MLFMA)、自适应积分法(AIM)pJ等都着眼于加速矩阵向量积和 只存储矩阵的一个稀疏部分。另一个途径是在保证精度的前提下，用尽可能少的基函数来展 开未知函判41，这也意味着矩量法矩阵的阶数减小，存储量随之大大降低。这点可以通过在 基函数中加入正确预测的相位信息得以实现。James首先将这一思想用于解决电磁问题【4J。 Altman等例证了在散射体的光滑部分的感应电流可以由一组复指数函数高效表示[51。文献【61 进一步将其应用于RWG函数(即线性相位RWG基函数LP-RWG)，并得到了理想的结果。 应用Galerkin法进行测试时需计算二重奇异积分，以LP．RWG函数作为展开和测试函 数时，积分核更是含有复杂的复指数函数，必需进行特殊处理。因此矩量法矩阵元素的奇异 性和几乎奇异性的精确处理成为这种矩量法模型实现的一个重要基础。在文献【6】中，作者 在局部坐标系下针对两个完全重合的三角形单元提出了一种奇异性处理方法，但这个方法不 能直接推广到包括几乎奇异性的情形。我们在文献【7】提出的关于RWG函数积分的奇异性处 理技术的基础上，推导了关于线性相位RWG函数积分的一般奇异性处理公式。文中提供了 数值算例，验证了所提出公式正确性。 、2 矩量法矩阵元素的公式 考虑三维理想导体的散射问题。引入等效电磁流，并应用等效原理，我们得电场积分 方程(EFIE)： ，-L(J)lf勰=E。l协 (1) r 厂 1 、 其中，三(乃=豇矽上凼’l歹+吉V(V7·．7)}G扩，尸)，．7p)为散射目标表面等效电流密度， 。一肚J，一，1 G扩，尸)。南为自由空间Green函数。为了应用矩量法求解方程(1)’将散射体表面 45 ^， ． 用三角形网稽离散，用线性相位RwG函数{爵)‘展开表面等效电流．7(尹)=∑％磊(尹)。这 n=1 个展开式与我们通常使用RWG函数来展开在形式上是一样的。线性相位RWG基函数的定 义如下‘引： 五c芦，={爻；{；；三：羔：茎兰乏；茎霉芦∈巧 (2) 【人-(芦)P叫％∽”一“’， 其中，人：(，)为通常的RWG矢量基函数，定义为乜1： 根r)_±去砩尹∈节 厶为公共边的长度，钟为三角形矸的面积，乃：为由自由顶点指向