本發(fā)明涉及一種目標(biāo)電磁散射特性快速計(jì)算方法，尤其涉及一種基于自適應(yīng)交叉近似加速的不連續(xù)伽遼金積分方程計(jì)算方法。
背景技術(shù)：
：目標(biāo)電磁散射特性快速精確計(jì)算是分析目標(biāo)隱身性能及反隱身設(shè)計(jì)的重要手段之一，隨著進(jìn)氣道、機(jī)翼前沿、雷達(dá)艙、頻率選擇表面等復(fù)雜部件在軍事目標(biāo)上的大量應(yīng)用，急需開展快速精確計(jì)算電大復(fù)雜目標(biāo)的電磁散射特性，為武器系統(tǒng)參數(shù)設(shè)計(jì)和電磁波感應(yīng)機(jī)理統(tǒng)計(jì)提供依據(jù)。如包含耶路撒冷十字頻率選擇(FSS)表面的某導(dǎo)彈雷達(dá)罩，若根據(jù)基函數(shù)尺寸網(wǎng)格剖分，則可得到96327個(gè)三角形面元，但基函數(shù)整體網(wǎng)格不足以描述頻率選擇表面，必須進(jìn)一步細(xì)分剖分網(wǎng)格尺度，一個(gè)FSS單元至少需12個(gè)三角形貼片，此時(shí)雷達(dá)罩整體產(chǎn)生835748個(gè)三角形貼片，極大的增加了計(jì)算量。另外一種可取的剖分策略為：FSS上選用小貼片單元，雷達(dá)罩主體采用大貼片單元，此時(shí)三角形貼片個(gè)數(shù)為261853個(gè)，可在保證計(jì)算精度的情況下極大的提高計(jì)算效率。呂志清，安翔等在2009年發(fā)表在全國微波毫米波會議優(yōu)秀論文集上的非專利文獻(xiàn)“基于自適應(yīng)交叉近似的非匹配網(wǎng)格區(qū)域分解快速算法”給出了一種基于自適應(yīng)交叉近似加速的非匹配網(wǎng)格的區(qū)域分解算法，并利用超松弛迭代法加速收斂過程，在貼片天線陣列、頻率選擇表面等周期結(jié)構(gòu)電磁計(jì)算方面展示了優(yōu)越的性能。但該方法是基于有限元方法給出的，計(jì)算精度不足且計(jì)算規(guī)模不大。ZhenPeng，Kheng-HweeLim等在2013年發(fā)表在IEEETransactionsonAntennasandPropagation上的非專利文獻(xiàn)“ADiscontinuousGalerkinSurfaceIntegralEquationMethodforElectromagneticWaveScatteringFromNonpenetrableTargets”公開了一種不連續(xù)伽遼金表面積分方程，詳細(xì)給出了計(jì)算非共形網(wǎng)格電大尺寸目標(biāo)的迭代求解方法，相比于傳統(tǒng)的電磁散射積分方程求解方程，該方法利用L2型基函數(shù)處理不同部件交界處的連續(xù)性問題，保證了目標(biāo)表面電磁流的連續(xù)性。但L2型基函數(shù)增加了阻抗矩陣規(guī)模，限制了目標(biāo)的求解能力。YuleiZhao，CHuangmingTong等在2011年發(fā)表在CrossStraitQuad-RegionalRadioScienceandWirelessTechnologyConference上的非專利文獻(xiàn)“TheACAAlgorithmintheElectromagneticScatterofMultipleTarget”提出了一種計(jì)算多尺度目標(biāo)的快速電磁散射方法，該方法利用純代數(shù)壓縮手段極大的提升了積分方程求解目標(biāo)電磁散射問題的能力，具有很好的擴(kuò)展性和實(shí)用性。但該方法針對目標(biāo)整體計(jì)算，尚未考慮目標(biāo)局部細(xì)節(jié)與目標(biāo)主體間的對應(yīng)關(guān)系。技術(shù)實(shí)現(xiàn)要素：本發(fā)明提供一種目標(biāo)電磁散射特性快速計(jì)算方法，解決了現(xiàn)有技術(shù)中多尺度電大復(fù)雜目標(biāo)的電子散射快速仿真計(jì)算的難題，為多尺度電大尺寸目標(biāo)電磁散射特性獲取提供技術(shù)手段，適應(yīng)性強(qiáng)、計(jì)算精度高，極大的拓展了積分方程的計(jì)算效率和計(jì)算能力，減少了剖分難度，易于工程實(shí)現(xiàn)。為了達(dá)到上述目的，本發(fā)明提供一種目標(biāo)電磁散射特性快速計(jì)算方法，包含以下步驟：步驟S1、通過測繪獲取多尺度目標(biāo)的幾何外形數(shù)據(jù)及部件間的相互位置關(guān)系，建立包含目標(biāo)精細(xì)部件的幾何外形非共形網(wǎng)格模型；步驟S2、根據(jù)目標(biāo)幾何外形及局部細(xì)節(jié)，不同部件采用不同尺寸的剖分網(wǎng)格，對目標(biāo)進(jìn)行整體剖分；步驟S3、在目標(biāo)剖分網(wǎng)格上選用電磁基函數(shù)模擬目標(biāo)表面的感應(yīng)電磁流，利用不連續(xù)伽遼金積分方程和邊界條件構(gòu)建阻抗矩陣元素，建立待測目標(biāo)感應(yīng)電磁流與照射場間的散射關(guān)系；電磁基函數(shù)為RWG基函數(shù)，在部件連接處選用半個(gè)基函數(shù)；步驟S4、采用樹形分塊分組策略對目標(biāo)進(jìn)行分組，形成便于計(jì)算的阻抗矩陣方程；步驟S5、利用自適應(yīng)交叉近似算法加速矩陣-矢量相乘速度；對樹形結(jié)構(gòu)，非相鄰組包含的基函數(shù)的互耦矩陣Zm×n進(jìn)行壓縮，形成兩個(gè)子矩陣Um×r,Vr×n，其中m,n為兩個(gè)盒子包含的基函數(shù)個(gè)數(shù)，r為壓縮后得到的矩陣維數(shù)；步驟S6、采用廣義最小殘量法迭代計(jì)算照射源下的目標(biāo)表面電磁流感應(yīng)系數(shù)；S7、根據(jù)目標(biāo)表面的感應(yīng)電磁流系數(shù)，開展目標(biāo)表面的感應(yīng)電磁流的輻射場計(jì)算，獲取二次散射場特性數(shù)據(jù)。所述的步驟S2中，電磁波頻率范圍為3.0GHz，照射波激勵(lì)方向?yàn)?俯仰，方位)＝(0°,0°)，散射場方向?yàn)?俯仰，方位)＝(0～180°，180°)，間隔1°，大部件剖分網(wǎng)格為1/5波長，小部件細(xì)節(jié)剖分網(wǎng)格為1/8波長，大小部件連接處的網(wǎng)格是非共形的。所述的步驟S3中，在部件內(nèi)部的RWG基函數(shù)定義在兩個(gè)三角形貼片的公共邊上，在部件邊界的貼片上定義的RWG基函數(shù)僅在一個(gè)三角形貼片上，RWG的定義形式為：fn(r)=Ln2An+ρn+(r)rinTn+Ln2An-ρn-(r)rinTn-0others---(1)]]>其中，r為三角形貼片內(nèi)部的位置矢量，ρ為三角形的邊對應(yīng)的頂點(diǎn)到r的位置矢量，表示與基函數(shù)相關(guān)聯(lián)的兩個(gè)三角形貼片，±為RWG基函數(shù)對應(yīng)的兩個(gè)三角形貼片，如果RWG基函數(shù)對應(yīng)的三角形在部件的邊界上，則沒有對應(yīng)的三角形對形成互耦，此時(shí)RWG僅定義在一個(gè)貼片上；利用構(gòu)造的RWG基函數(shù)組構(gòu)造阻抗矩陣元素，根據(jù)電磁場理論，任意金屬目標(biāo)的邊界積分方程為：Eti(r)=[jk0A(r)-1jk0▿▿·A(r)]t---(2)]]>其中，為目標(biāo)表面的激勵(lì)電場切向分量，k0為波數(shù)，A(r)為自由空間的矢量位函數(shù)；根據(jù)步驟S2中的剖分策略，目標(biāo)整體分為多個(gè)部件，目標(biāo)邊界可描述為：為目標(biāo)原始邊界，Sn為每個(gè)剖分部件的邊界，利用不連續(xù)伽遼金方法，每個(gè)剖分部件的感應(yīng)電磁流的殘差可描述為：RSmEFIE=Et,mi(r)-Σk=1N[jk0Ak(r)-1jk0▿▿·Ak(r)]t,m=1,2,...N---(3)]]>對上式兩邊利用測試函數(shù)組tn做內(nèi)積運(yùn)算，可得到參數(shù)的弱形式為：<tn,RSmEFIE>Sm=<tn,Et,mi(r)>Sm-<tn,Σk=1N[jk0Ak(r)-1jk0▿▿·Ak(r)]t>Sm,n=1,2,...N---(4)]]>同樣的，可以得到每個(gè)部件邊界上的磁場殘差余量的弱形式為：<tn,RSmMFIE>Sm=<tn,Ht,mi(r)>Sm-<tn,12J(r)>Sm+<tn,Σk=1NFk(r)>Sm,n=1,2,...N---(5)]]>其中，為第m個(gè)剖分部件的切向磁場激勵(lì)分量，J(r)為目標(biāo)感應(yīng)電流，F(xiàn)k(r)為自由空間的磁矢量位函數(shù)；為保證互相連接的部件間的電磁流連續(xù)性，在兩個(gè)相連部件邊界輪廓上強(qiáng)加邊界條件：<t,RCcurrent>=<t,1ik0ΣCtmn·Jm+tmn·Jn><t,RCpotential>=<t,1ik0ΣC∫Cmn(tmn·Jm(r)+tmn·Jn(r))Gdr>---(6)]]>整個(gè)矩陣運(yùn)算的誤差為式(4)、(5)、(6)三個(gè)部分的線性累積，通過迭代運(yùn)算求解使得整個(gè)殘差最小的電磁流系數(shù)即可相應(yīng)的感應(yīng)電磁流值。所述的步驟S4中，采用樹形結(jié)構(gòu)對目標(biāo)進(jìn)行分組，首先找到能包圍整個(gè)目標(biāo)的最小長方體盒子，然后在x、y、z三個(gè)方向?qū)Φ确纸M成八個(gè)子盒子，依次將每個(gè)子盒子細(xì)分，直至最細(xì)層盒子尺寸接近半個(gè)波長停止，分組過程中，不包含目標(biāo)網(wǎng)格的盒子舍棄，僅保留有效盒子組。所述的步驟S5具體包含以下步驟：步驟S5.1、初始化第一個(gè)行索引I1，此時(shí)近似矩陣步驟S5.2、初始化近似誤差矩陣的第一行：步驟S5.3、在第一行中找尋最大值從而確定第一個(gè)列索引J1:步驟S5.4、得到V矩陣的第一行：步驟S5.5、初始化近似誤差矩陣的第一列：步驟S5.6、計(jì)算U矩陣的第一列：步驟S5.7、計(jì)算：||Z~(1)||2=||Z~(0)||2+||u1||2||v1||2(||Z~(k)||2=||Z~(k-1)||2+2Σj=1k-1|ujTuk|·|vjTvk|||u1||2||v1||2);]]>其中，||·||為矩陣的二范數(shù)；步驟S5.8、判斷收斂誤差：如果迭代終止，否則繼續(xù)；步驟S5.9、找到第一列中的最大值作為第二個(gè)索引I2，且I2≠I1；重復(fù)上述步驟，至到循環(huán)終止跳出，最終兩個(gè)非近鄰盒子形成的組行元素壓縮為：Zm×n≈Um×r×Vr×n(7)其中，r為迭代步數(shù)；在迭代求解運(yùn)算中，矩陣和矢量的乘積表示為：Zm×n·I＝Um×r·(Vr×n·I)(8)。自適應(yīng)交叉近似算法將矩陣計(jì)算的存儲量和復(fù)雜度由O(m×n)降到O(r×(m+n))。所述的步驟S7中，根據(jù)感應(yīng)電磁流系數(shù)求解輻射場信息，計(jì)算遠(yuǎn)區(qū)電磁場值及近區(qū)相位變化，在保證求解精度的前提下極大的提升求解效率。本發(fā)明具有以下優(yōu)點(diǎn)：采用非共形網(wǎng)格對目標(biāo)幾何模型進(jìn)行剖分，極大的簡化了剖分模擬難度，對于復(fù)雜的目標(biāo)可采用更加靈活的方式處理，采用RWG基函數(shù)模擬目標(biāo)表面的感應(yīng)電磁流，在不同部件連接處采用半個(gè)RWG基函數(shù)的形式極大的節(jié)省了計(jì)算機(jī)資源，采用自適應(yīng)交叉近似算法加速迭代求解過程，不需通過復(fù)雜的場源關(guān)系及微分變化過程即可改善求解能力，計(jì)算簡單，代碼移植性高。附圖說明圖1為本發(fā)明提供的一種目標(biāo)電磁散射特性快速計(jì)算方法的流程圖。圖2為實(shí)施例中掃描得到某導(dǎo)彈幾何模型。圖3為實(shí)施例中剖分網(wǎng)格與現(xiàn)有網(wǎng)格技術(shù)分布對比。圖4為實(shí)施例采用的RWG基函數(shù)結(jié)構(gòu)。圖5為實(shí)施例采用的樹形分組結(jié)構(gòu)。圖6為實(shí)施例中目標(biāo)感應(yīng)電流分布與現(xiàn)有方法計(jì)算結(jié)果對比。圖7為實(shí)施例中目標(biāo)遠(yuǎn)場RCS分布與商業(yè)軟件計(jì)算對比結(jié)果。具體實(shí)施方式以下根據(jù)圖1～圖7，具體說明本發(fā)明的較佳實(shí)施例。如圖1所示，本發(fā)明提供一種目標(biāo)電磁散射特性快速計(jì)算方法，包含以下步驟：步驟S1、通過測繪獲取多尺度目標(biāo)的幾何外形數(shù)據(jù)及部件間的相互位置關(guān)系，建立包含目標(biāo)精細(xì)部件的幾何外形非共形網(wǎng)格模型；如圖2所示，本實(shí)施例中，測繪目標(biāo)為某導(dǎo)彈模型，目標(biāo)原始尺寸(m)為：5.58×2.49×1.06，測繪儀器選用激光掃描儀，掃描精度為1mm，構(gòu)建某導(dǎo)彈目標(biāo)模型，目標(biāo)幾何外形建模精度為99％；步驟S2、根據(jù)目標(biāo)幾何外形及局部細(xì)節(jié)，不同部件采用不同尺寸的剖分網(wǎng)格，對目標(biāo)進(jìn)行整體剖分；將步驟S1中得到的多尺度目標(biāo)幾何外形進(jìn)行分類區(qū)分，依據(jù)照射電磁波頻率，選擇足以描述目標(biāo)幾何、電磁模型的剖分網(wǎng)格尺寸對目標(biāo)進(jìn)行整體擬合，為降低對目標(biāo)整體剖分帶來的壓力；采用局部細(xì)分的方法減輕計(jì)算量，即針對不同的部件，采用不同的三角形網(wǎng)格進(jìn)行擬合，且不同部件間連接處不必是連續(xù)的，即非共形的；電磁波頻率范圍為3.0GHz，照射波激勵(lì)方向?yàn)?俯仰，方位)＝(0°,0°)，散射場方向?yàn)?俯仰，方位)＝(0～180°，180°)，間隔1°，大部件剖分網(wǎng)格為1/5波長，小部件細(xì)節(jié)剖分網(wǎng)格為1/8波長，大小部件連接處的網(wǎng)格是非共形的，非共形網(wǎng)格處采用特殊的邊界條件保證電磁流連續(xù)性，進(jìn)而保證目標(biāo)整體感應(yīng)電磁流模擬的準(zhǔn)確性和穩(wěn)定性；如圖3所示，本實(shí)施例中，原始剖分網(wǎng)格得到261853個(gè)三角形貼片，采用本發(fā)明的方法，照射波頻率為3.0GHz，目標(biāo)主體剖分網(wǎng)格尺寸為20mm，機(jī)翼部分剖分網(wǎng)格為12.8mm，進(jìn)氣道剖分網(wǎng)格為10mm，機(jī)翼、進(jìn)氣道與主體連接處的網(wǎng)格不必是連續(xù)的，得到109638個(gè)三角形貼片，大大減少了貼片數(shù)量；步驟S3、在目標(biāo)剖分網(wǎng)格上選用合適的電磁基函數(shù)模擬目標(biāo)表面的感應(yīng)電磁流，利用不連續(xù)伽遼金積分方程和邊界條件構(gòu)建阻抗矩陣元素，建立待測目標(biāo)感應(yīng)電磁流與照射場間的散射關(guān)系；本實(shí)施例中，電磁基函數(shù)為RWG基函數(shù)，在部件連接處選用半個(gè)基函數(shù)；如圖4所示，在部件內(nèi)部的RWG基函數(shù)定義在兩個(gè)三角形貼片的公共邊上，在部件邊界的貼片上定義的RWG基函數(shù)僅在一個(gè)三角形貼片上，RWG的定義形式為：fn(r)=Ln2An+ρn+(r)rinTn+Ln2An-ρn-(r)rinTn-0others---(1)]]>其中，r為三角形貼片內(nèi)部的位置矢量，ρ為三角形的邊對應(yīng)的頂點(diǎn)到r的位置矢量，表示與基函數(shù)相關(guān)聯(lián)的兩個(gè)三角形貼片，±為RWG基函數(shù)對應(yīng)的兩個(gè)三角形貼片，如果RWG基函數(shù)對應(yīng)的三角形在部件的邊界上，則沒有對應(yīng)的三角形對形成互耦，此時(shí)RWG僅定義在一個(gè)貼片上；利用構(gòu)造的RWG基函數(shù)組構(gòu)造阻抗矩陣元素，根據(jù)電磁場理論，任意金屬目標(biāo)的邊界積分方程為：Eti(r)=[jk0A(r)-1jk0▿▿·A(r)]t---(2)]]>其中，為目標(biāo)表面的激勵(lì)電場切向分量，k0為波數(shù)，A(r)為自由空間的矢量位函數(shù)；根據(jù)步驟S2中的剖分策略，目標(biāo)整體分為多個(gè)部件，目標(biāo)邊界可描述為：為目標(biāo)原始邊界，Sn為每個(gè)剖分部件的邊界，利用不連續(xù)伽遼金方法，每個(gè)剖分部件的感應(yīng)電磁流的殘差可描述為：RSmEFIE=Et,mi(r)-Σk=1N[jk0Ak(r)-1jk0▿▿·Ak(r)]t,m=1,2,...N---(3)]]>對上式兩邊利用測試函數(shù)組tn(本發(fā)明中測試函數(shù)與基函數(shù)表達(dá)形式一致)做內(nèi)積運(yùn)算，可得到參數(shù)的弱形式為：<tn,RSmEFIE>Sm=<tn,Et,mi(r)>Sm-<tn,Σk=1N[jk0Ak(r)-1jk0▿▿·Ak(r)]t>Sm,n=1,2,...N---(4)]]>同樣的，可以得到每個(gè)部件邊界上的磁場殘差余量的弱形式為：<tn,RSmMFIE>Sm=<tn,Ht,mi(r)>Sm-<tn,12J(r)>Sm+<tn,Σk=1NFk(r)>Sm,n=1,2,...N---(5)]]>其中，為第m個(gè)剖分部件的切向磁場激勵(lì)分量，J(r)為目標(biāo)感應(yīng)電流，F(xiàn)k(r)為自由空間的磁矢量位函數(shù)；為保證互相連接的部件間的電磁流連續(xù)性，在兩個(gè)相連部件邊界輪廓上強(qiáng)加邊界條件：<t,RCcurrent>=<t,1ik0ΣCtmn·Jm+tmn·Jn><t,RCpotential>=<t,1ik0ΣC∫Cmn(tmn·Jm(r)+tmn·Jn(r))Gdr>---(6)]]>整個(gè)矩陣運(yùn)算的誤差為式(4)、(5)、(6)三個(gè)部分的線性累積，通過迭代運(yùn)算求解使得整個(gè)殘差最小的電磁流系數(shù)即可相應(yīng)的感應(yīng)電磁流值；步驟S4、采用樹形分塊分組策略對目標(biāo)進(jìn)行分組，形成便于計(jì)算的阻抗矩陣方程；在目標(biāo)剖分過程中，模擬目標(biāo)幾何模型的網(wǎng)格數(shù)是巨大的，即便當(dāng)前的超級計(jì)算機(jī)也無法在短時(shí)間內(nèi)完成對式(4)、(5)、(6)的迭代運(yùn)算，例如對一個(gè)網(wǎng)格貼片數(shù)為100000的目標(biāo)，其所需的存儲空間為200GHz的內(nèi)存空間，且隨著貼片網(wǎng)格的增加，其所需內(nèi)存空間按照平方量級增加，在分析電大尺寸目標(biāo)時(shí)是非常耗時(shí)的；如圖5所示，本發(fā)明采用樹形結(jié)構(gòu)對目標(biāo)進(jìn)行分組，首先找到能包圍整個(gè)目標(biāo)的最小長方體盒子，然后在x、y、z三個(gè)方向?qū)Φ确纸M成八個(gè)子盒子，依次將每個(gè)子盒子細(xì)分，直至最細(xì)層盒子尺寸接近半個(gè)波長停止，分組過程中，不包含目標(biāo)網(wǎng)格的盒子舍棄，僅保留有效盒子組；本實(shí)施例中，采用的樹形結(jié)構(gòu)為八叉樹形式，最底層八叉樹的邊長為40mm，通過迭代形式最終形成矩陣塊；步驟S5、利用自適應(yīng)交叉近似算法加速矩陣-矢量相乘速度；自適應(yīng)交叉近似算法是一種高效的純代數(shù)矩陣壓縮方法，利用壓縮計(jì)算減少阻抗矩陣的存儲量，本實(shí)施例中，自適應(yīng)交叉近似算法采用的壓縮門限為0.01，迭代收斂門限為0.003，在目標(biāo)電磁場計(jì)算中，目標(biāo)不同部件間的電磁波互耦隨距離變化劇烈，對步驟S4中的樹形結(jié)構(gòu)，非相鄰組包含的基函數(shù)的互耦矩陣Zm×n進(jìn)行壓縮，形成兩個(gè)子矩陣Um×r,Vr×n，其中m,n為兩個(gè)盒子包含的基函數(shù)個(gè)數(shù)，r為壓縮后得到的矩陣維數(shù)；所述的自適應(yīng)交叉近似算法具體包含以下步驟：步驟S5.1、初始化第一個(gè)行索引I1，此時(shí)近似矩陣步驟S5.2、初始化近似誤差矩陣的第一行：步驟S5.3、在第一行中找尋最大值從而確定第一個(gè)列索引J1:步驟S5.4、得到V矩陣的第一行：步驟S5.5、初始化近似誤差矩陣的第一列：步驟S5.6、計(jì)算U矩陣的第一列：步驟S5.7、計(jì)算：||Z~(1)||2=||Z~(0)||2+||u1||2||v1||2(||Z~(k)||2=||Z~(k-1)||2+2Σj=1k-1|ujTuk|·|vjTvk|||u1||2||v1||2);]]>其中，||·||為矩陣的二范數(shù)；步驟S5.8、判斷收斂誤差：如果迭代終止，否則繼續(xù)步驟S5.9、找到第一列中的最大值作為第二個(gè)索引I2，且I2≠I1；重復(fù)上述步驟，至到循環(huán)終止跳出，最終兩個(gè)非近鄰盒子形成的組行元素壓縮為：Zm×n≈Um×r×Vr×n(7)其中，r為迭代步數(shù)；在迭代求解運(yùn)算中，矩陣和矢量的乘積表示為：Zm×n·I＝Um×r·(Vr×n·I)(8)顯然，自適應(yīng)交叉近似算法將矩陣計(jì)算的存儲量和復(fù)雜度由O(m×n)降到O(r×(m+n))，特別是r往往遠(yuǎn)低于m,n，從而達(dá)到節(jié)省計(jì)算量的目的。步驟S6、采用廣義最小殘量法(Gmres)迭代計(jì)算照射源下的目標(biāo)表面電磁流感應(yīng)系數(shù)；采用Gmres迭代算法僅利用一次矩陣矢量相乘即可完成迭代運(yùn)算，利用矩陣求解空間向殘量小的方向開展，從而逐次逼近真實(shí)解；S7、根據(jù)目標(biāo)表面的感應(yīng)電磁流系數(shù)，開展目標(biāo)表面的感應(yīng)電磁流的輻射場計(jì)算，獲取二次散射場特性數(shù)據(jù)；得到諸如RCS等近遠(yuǎn)場特性數(shù)據(jù)；根據(jù)感應(yīng)電磁流系數(shù)求解輻射場信息，計(jì)算遠(yuǎn)區(qū)電磁場值及近區(qū)相位變化，在保證求解精度的前提下極大的提升求解效率。圖6和圖7給出了利用本發(fā)明的計(jì)算結(jié)果與商業(yè)軟件計(jì)算結(jié)果的對比，圖6中左側(cè)為商業(yè)軟件電流結(jié)果，右側(cè)為本發(fā)明的電流計(jì)算結(jié)果，表明本發(fā)明計(jì)算的感應(yīng)電流分布與商業(yè)軟件一致，驗(yàn)證了本發(fā)明的正確性，圖7中給出了通過感應(yīng)電流計(jì)算的遠(yuǎn)場散射RCS值，驗(yàn)證了計(jì)算遠(yuǎn)區(qū)輻射場的正確性。本發(fā)明具有以下優(yōu)點(diǎn)：采用非共形網(wǎng)格對目標(biāo)幾何模型進(jìn)行剖分，極大的簡化了剖分模擬難度，對于復(fù)雜的目標(biāo)可采用更加靈活的方式處理，采用RWG基函數(shù)模擬目標(biāo)表面的感應(yīng)電磁流，在不同部件連接處采用半個(gè)RWG基函數(shù)的形式極大的節(jié)省了計(jì)算機(jī)資源，采用自適應(yīng)交叉近似算法加速迭代求解過程，不需通過復(fù)雜的場源關(guān)系及微分變化過程即可改善求解能力，計(jì)算簡單，代碼移植性高。盡管本發(fā)明的內(nèi)容已經(jīng)通過上述優(yōu)選實(shí)施例作了詳細(xì)介紹，但應(yīng)當(dāng)認(rèn)識到上述的描述不應(yīng)被認(rèn)為是對本發(fā)明的限制。在本領(lǐng)域技術(shù)人員閱讀了上述內(nèi)容后，對于本發(fā)明的多種修改和替代都將是顯而易見的。因此，本發(fā)明的保護(hù)范圍應(yīng)由所附的權(quán)利要求來限定。當(dāng)前第1頁1 2 3