基于快速多极子算法的电磁散射计算加速性能研究.pdfVIP

现代制造技术与装备 2015第1期总第224期 6 基于陕速多极子算法的电磁散射计算加速性能研究 马俊成 张佳芬 (电子科技大学成都学院 微电子技术系，成都 611731) 摘 要：计算电磁散射问题 时常用的方法是矩量法，但计算电大尺寸 目标及密网格 问题 时，计算量很大 ，对 计算机的存储要求很高。快速 多极子算法作为矩量法的改进算法 ，现 已广泛应用于此类问题的求解。本文就是 在验证快速多极子算法正确性的基础上，研 究快速多极子算法的加速性能。 关键词 ：矩量法 快速 多极子算法 引言 积分在角谱空间中的积分点数 K。 矩量法 (MOM)是求解积分 方程 的主要 方法 ，有精度 K=2L2 (4) 高 、所使用 的格林函数可直接满足辐射条件和无需设置 设 目标总的未知数为 N，未知数 N包含在 M个组内， 吸收边界等优点 ，被广泛应用于求解 电磁散射 问题。然而 这样每个组所含未知数 的平均个数为 A 一N／M，总平面波 传统的矩量法(MOM)仅限于求解低频区和谐振区 目标的 数为 K矩阵 向量的乘积过程分为三步 ，具体如下 ： 电磁散射问题。因为应用矩量法所得到 的是一个满阵，对 (1)聚合 于高频 区电大尺寸 目标的求解 ，很多情况下因其所需要 ， (忘)=∑ ，(，)，，(t) r、 的存储量极大而无法实现。例如 ：在求解具有 10万未知 数的问题时 ，要生成而且并存储一个含 10m个单精度复数 其中，m =1,2…．，M ，k=1,2，…， 。 类型的元素 的矩阵，虽然现在 的主流计算机能够满足 (2)转移 80GB的存储要求 ，但是其计算量也是相 当惊人的，应用迭 ()= (嬲 ， ) ( ) 代法时其计算复杂度为 0(N 10m)，应用直接求逆法其计 咄 (6) 其中，m=l…2 ．．，M,k=1，2，…，Ko 代表第龃 的所有远场组的集 算复杂度则将达到 0(10 )。 为了提高传统矩量法 的求解效率 ，提出了很多改进方 (3)配置 法 ，其 中以快速多极子算法 (FMM)最具代表性 ，因为快速 上 = ( )= Z ， 多极子算法 (FMM)能使 目标的积分方程法求解复杂度变 为只有 0(10 )。 薹z + gni自 (7) 1快速多极子算法的数学描述 快速多极子方法的数学基础是矢量加法定理 。首先利 薹z，+争K 。()gni() 用加法定理对积分方程 中的标量格林 函数进行处理 ，然 其 中，=l，2，KNo 是向 第n个元素 后在角谱空间中展开，利用平面波进行算子对角化 ，最终 式 (7)中，方程右边第一项表示的是近场作用 ，可直接 将稠密矩阵与向量的乘积运算转化为几个稀疏矩阵与向 利用稀疏矩阵 z一计算。式 (5)中的向量表示所有源点组 量的乘积运算。 的辐射方向图或外向平面波展开 。将每个源点组的外 向平 展开形式如下 ： 面波展开 S，乘以转移 因子 TIl，并且对远场组 D 中所有源 (rIrL1，+VVI，r，)