三维大纵横比目标散射的快速精确求解.pdfVIP

三维大纵横比目标散射的快速精确求解 胡俊聂在平王军姚海英王浩刚 电子科技大学微波工程系，成都610054 [摘要】本文采用积分方程法严格求解三维大纵横比目标的电磁散 射．在积分方程法的迭代求解中用快速多极子法(FMM)加速矩阵与矢 量的相乘计算，同时运用快速傅立叶变换(FFT)进一步提高快速多极 子方法中的转换项计算。数值结果表明：这种快速多极子法一快速傅 立叶变换方法(FMM-FFT)特别适合于三维大纵横比目标的散射求解． 该法也适合于介质体、随机粗糙表面的散射分析． 美键词：大纵横比积分方程法快速多极子法快速傅立叶变换 1 引言 近年来，由于工程与应用的强烈需要，如何高效、精确求解三维电大尺寸目标的电磁 散射问题十分令人关注。积分方程法(IEM)直接满足辐射条件、网格色散误差极小，加之该 法仅在散射体表面或内部离散未知量较少，因此被广泛采用。得益于计算机技术和许多数 值方法的发展，许多基于积分方程法的传统方法都拓宽了应用领域。如矩量法040M)现已 method)便是用在矩量法的迭代求解中加速矩阵与矢量计算的高效方法。对于一 multipole Ⅳ阶矩阵，直接计算矩阵与矢量相乘的工作量为OIN!j量级，采用快速多极子方法计算其 5】量级。对于三维大纵横比目标的电磁散射，采用快速多极子法一快速 工作量则降为O(N1 傅立叶变换法(FMM．FFT)Ⅲ9能得到更高的计算效率。 粗糙表面的电磁散射。快速多极予方法用于加速非附近组间的互耦计算，通过聚合、转换、 配置三步实现，其中的转换项计算为一线性卷积，可用快速傅立叶变换计算：而附近组问 的互耦则仍采用直接计算。对于随机粗糙表面目标散射，该法比快速多极子法需要更少的 计算时间和计算机存储，因此效率很高。该法已被我们成功用于三维电大尺寸平面型目标 的电磁散射求解【3l。本文则将该法用于三维大纵横比目标的散射分析，取得了很好的计算 效果。 2快速多极子法一快速傅立叶变换法(FMM．FFT) 在三维导体目标的矢量散射问题中，电场积分方程可写为 ·2l· ㈣ im，n≯加w]譬册百41ri㈣，一 其中r，r’分别为表面S上场点、源点矢径，f为S上任意单位切向矢量，E’为外来入射 场，k为自由空间中波数，R=j，一，’『。方程(1)可用并矢格林函数C(r，，’J表达为 i，f瓢，r_肌邶7=等他¨J (2) 其中 r-d 一 r一 1 ]。“1 Gr一7—2i卜古V矿i哥T ④ 快速多极子法的基本数学工具是加法定理。首先将导电目标用一合适的长方体罩住， 在该长方体的Ⅳ．，r．Z三个方向上均匀划分为Nj格，N、格和Ⅳ：格，得到N；-N—N：个 立方体。判断每一基函数的中心与各立方体的中心的距离，得到待求的非空立方体。对于 附近的非空立方体的场源互耦．采用直接计算：非附近的非空立方体间的场源互耦则用加 法定理计算。运用加法定理，得丑J141 d2狮一瓢“n”“’。帆‘¨I，mo,r＼r扩‰＼㈣ 己frj，r户去l 嘭。。．r‘，i／=∑f?f，2l+1J^j1。r帆，．，鼻，‘。．石／ (5) 其中，，，ri分别为场点、源点矢径，匕。；‰一r。．，～，r，，分别为场立方体、源立方体 的中心矢径。^∥和只分别为第一类球汉克尔函数和勒让德函数， 瞄2￡是谱空间单位球 D为立方体对角线长度。 选用合适的基函数-，。与权函数，，进行内积，得到线性代数方程组为 三 乞以口。=‘ ／=1,2，…Ⅳ