朱錦新 鄭永濤 孔維賓 畢奧燃

摘? 要：針對(duì)預(yù)修正快速傅里葉變換（P-FFT）方法的近場矩陣填充速度和效率不高的問題，利用快速填充技術(shù)減少近場矩陣的計(jì)算時(shí)間，分析電大目標(biāo)的電磁特性。采用預(yù)修正快速傅里葉變換方法，不僅能夠減少存儲(chǔ)需求，而且能夠加快遠(yuǎn)場的矩陣矢量積。重點(diǎn)研究了基于三角形配對(duì)的快速填充方法，通過提取基函數(shù)測試配對(duì)積分過程中的共享積分部分，對(duì)近場矩陣元素進(jìn)行填充，提高了計(jì)算效率。數(shù)值結(jié)果表明，基于三角形配對(duì)的預(yù)修正快速傅里葉變換方法可以減少近場部分的計(jì)算時(shí)間，從而提高算法的效率。

關(guān)鍵詞：矩量法;預(yù)修正快速傅里葉變換方法;電磁散射;近場矩陣

中圖分類號(hào)：TP311? ? ?文獻(xiàn)標(biāo)識(shí)碼：A

Research on Analyzing Electromagnetic Scattering Problems of Targets?by P-FFT based on Matrix Fast Filling

ZHU Jinxin， ZHENG Yongtao， KONG Weibin， BI Aoran

（School of Information Technology， Yancheng Institute of Technology， Yancheng 224051， China）

Abstract： Aiming at the low filling speed efficiency of the near-field matrix in the P-FFT （Pre-corrected Fast Fourier Transform） method， this paper proposes to use fast filling technique to reduce calculation time of the near-field matrix， which is used to analyze the electromagnetic characteristics of the electrically large target. The P-FFT method not only reduces storage requirements， but also speeds up the far-field matrix-vector product. The fast filling method based on triangle pairing is mainly studied. By extracting the shared integral part of the paired integration process of the basis function test， the near-field matrix elements are filled and the calculation efficiency is improved. Numerical results show that the P-FFT method based on triangle pairing can reduce the computation time of the near-field part， thereby improving the efficiency of the algorithm.

Keywords： method of moment; pre-corrected fast Fourier transform method; electromagnetic scattering; near-field matrix

1? ?引言（Introduction）

矩量法作為一種精確的數(shù)值方法，被廣泛應(yīng)用于分析各種復(fù)雜的電磁問題[1-2]。然而，矩量法的存儲(chǔ)量和計(jì)算量都很大，限制了求解電磁問題的規(guī)模。為了解決矩量法的局限性，許多人開始研究基于矩量法的快速算法?？焖偎惴ù篌w上可以分為兩類：一類是與積分核相關(guān)的算法，如快速多級(jí)子算法（FMM）[3]、多層快速多極子（MLFMA）[4-5]和多層格林函數(shù)插值方法（MLGFIM）[6]，以及基于快速傅立葉變換的方法，如自適應(yīng)積分方法（AIM）、預(yù)修正快速傅里葉變換方法（P-FFT）、積分方程快速傅里葉變換（IE-FFT）等[7-11]。另一類是矩陣壓縮方法，如自適應(yīng)交叉逼近（ACA）算法[12]、骨架化[13]等。

本文提出了一種結(jié)合近場矩陣快速填充的預(yù)修正快速傅里葉變換（P-FFT）方法來分析電大目標(biāo)的電磁特性。采用P-FFT方法可加快矩陣矢量乘法的求解速度。此外，阻抗矩陣存儲(chǔ)稀疏，減少了內(nèi)存需求。同時(shí)，采用基于三角形配對(duì)的快速填充方法生成近場矩陣，能夠提高計(jì)算速度。

2? ?理論基礎(chǔ)（Theoretical basis）

2.1? ?矩量法

為了用矩量法數(shù)值求解電場積分方程（EFIE）和磁場積分方程（MFIE），首先將導(dǎo)體表面劃分為三角形，然后用RWG基函數(shù)展開式表示的電流密度代入：

3? ?數(shù)值結(jié)果（Numerical results）

在本部分中給出了兩個(gè)數(shù)值算例，以驗(yàn)證P-FFT方法的效率和準(zhǔn)確性。在基于快速傅里葉方法的快速算法中選擇笛卡爾網(wǎng)格間距相等，即hx=hy=hz=h，展開的階數(shù)為M。

3.1? ?球體

如圖2所示為一個(gè)半徑為2λh的PEC球。球的雷達(dá)散射截面使用Mie級(jí)數(shù)和P-FFT計(jì)算。球面使用1/10波長離散，共產(chǎn)生18，468 個(gè)RWG基函數(shù)。計(jì)算所得雙站RCS曲線被顯示在圖3上。從圖3中可以看到，三條曲線是一致的，這說明了基于三角形填充的P-FFT的正確性。從表1可以看出，本文提出的方法比直接P-FFT在填充近場矩陣速度上具有很大的優(yōu)勢(shì)。

3.2? ?導(dǎo)彈模型

如圖4所示，給出了一個(gè)PEC導(dǎo)彈模型。導(dǎo)彈被剖分為98，475 個(gè)RWG基函數(shù)。從圖5中可以看出，基于基函數(shù)填充的P-FFT、AIM和基于三角形填充的P-FFT的計(jì)算結(jié)果的曲線是吻合的，說明基于三角形填充的P-FFT在計(jì)算復(fù)雜結(jié)構(gòu)的物體電磁散射問題上是正確的。在表1中，本文方法與AIM、基于基函數(shù)填充的P-FFT在近場矩陣填充的速度上進(jìn)行了對(duì)比，結(jié)果表明本文方法的計(jì)算效率得到了提高。

4? ?結(jié)論（Conclusion）

本文采用P-FFT方法結(jié)合快速填充方法求解電磁場積分方程，通過引入一種三角形配對(duì)的方式，提取其中的積分共享部分填充矩陣元素。從本文所給的數(shù)值算例可以看出，三角形配對(duì)的采用，減少了近場矩陣元素計(jì)算的時(shí)間，從而提高了預(yù)修正所需要的計(jì)算效率。

參考文獻(xiàn)（References）

[1] HARRINGTON R F. Field computation by moment methods[M]. New York， NY， USA： MacMillian， 1968：20-26.

[2] 李弘祖，郭立新，董春雷，等.基于八叉樹優(yōu)化的MoM-PO/PTD混合算法分析目標(biāo)電磁散射及輻射問題[J/OL].系統(tǒng)工程與電子技術(shù).[2021-8-17].https：//kns.cnki.net/kcms/detail/11.2422.TN.20210817.0947.002.html.

[3] ENGHETA N， MURPHY W D， ROKHLIN V， et al. The fast multipole method （FMM） for electromagnetic scattering problems[J]. IEEE Transactions on Antennas and Propagation， 1992， 40（6）：634-641.

[4] SONG J M， LU C C， CHEW W C. Multilevel fast multipole algorithm for electromagnetic scattering by large complex objects[J]. IEEE Transactions on Antennas and Propagation， 1997， 45（10）：1488-1493.

[5] DELGADO C， CATEDRA F. Fast monostatic RCS computation using the near-field sparse approximate inverse and the multilevel fast multipole algorithm[J]. Applied Computational Electromagnetics Society， 2020， 35（7）：735-741.

[6] LI L， WANG H G， CHAN C H. An improved multilevel Green's function interpolation method with adaptive phase compensation[J]. IEEE Transactions on Antennas and Propagation， 2008， 56（6）：1381-1393.

[7] BLESZYBSKI E， BLESZYNSKI M， JAROSZEWICZ T. AIM： Adaptive integral method for solving large-scale electromagnetic scattering and radiation problems[J]. Radio Science， 1996（31）：1225-1251.

[8] PHILLIPS J R， WHITE J K. A Precorrected-FFT method for electrostatic analysis of complicated 3-D structures[J]. IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems， 1997（16）：1059-1072.

[9]? SEO S M， LEE J F. A fast IE-FFT algorithm for solving PEC scattering problems[J]. IEEE Transactions on Magnetics， 2005， 41（5）：1476-1479.

[10] 韓曉冰，張瀟，王露潔，等.采用分部外推BCGS-FFT方法快速求解電磁散射問題[J/OL].電波科學(xué)學(xué)報(bào).[2021-8-17].https：//doi.org/10.13443/j.cjors.20 20 081801.

[11] LI M， CHEN R S， WANG H， et al.? A multilevel FFT method for the 3-D capacitance extraction[J]. IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems， 2013， 32（2）：318-322.

[12] LI M， SU T， CHEN R S. Equivalence principle algorithm with body of revolution equivalence surface for the modeling of large multiscale structures[J]. IEEE Transactions on Antennas and Propagation， 2016， 64（5）：1818-1828.

[13] RASOOL H， CHEN J， PAN X M， et al. Skeletonization accelerated solution of Crank-Nicolson method for solving three-dimensional parabolic equation[J]. Applied Computational Electromagnetics Society， 2020， 35（9）：1006-1011.

作者簡介：

朱錦新（1978-），男，碩士，講師.研究領(lǐng)域：計(jì)算智能，計(jì)算機(jī)并行計(jì)算.

鄭永濤（1998-），男，碩士生.研究領(lǐng)域：系統(tǒng)與自動(dòng)控制.

孔維賓（1982-），男，博士，副教授.研究領(lǐng)域：計(jì)算電磁學(xué)，智能算法，電子信息.

畢奧燃（1999-），男，本科生.研究領(lǐng)域：電子信息，智能控制.