基于混合位積分方程的矩量法的微帶天線(xiàn)RCS分析

在采用基于混合位積分方程的矩量法分析微帶天線(xiàn)RCS時(shí)，首先采用二級離散復鏡像法求解空域格林函數，從而大大提高了矩量法的計算效率，然后利用三角形網(wǎng)格剖分計算目標，也使得矩量法更適合分析復雜結構。采用本文方法的計算結果也和相關(guān)文獻的計算結果相一致，從而驗證了本文方法的高效性和正確性。
　　
　　關(guān)鍵詞：矩量法；離散復鏡像法；微帶天線(xiàn)；RCS
　　
　　Calculations of Microstrip Antenna RCS
　　
　　SHI Lei,DONG Xu,JIAO Yi (Xi′an Communication Institute,Xi′an,710106,China)
　　オ
　　Abstract：In analyzing the RCS of microstrip antennas by MPIE-MoM,the two-level approximation discrete complex image method is applied to obtain spatial-domain Green′s function,which enhance computing efficiency.Triangular facets are used to approximate the surface of object.Hence the MoM is suitable to analyze complex structure.The result is as the same as the document given,so the method in this text is fast and right.
　　
　　Keywords：MoM;discrete complex image method;microstrip antenna;RCS
　　オ
　　微帶天線(xiàn)具有尺寸小、重量輕、成本低、易于廣角掃描，并能與飛行器表面共形等許多優(yōu)點(diǎn)，有著(zhù)廣泛的應用前景，因此對于微帶天線(xiàn)的研究也越來(lái)越受到人們的重視。目前研究微帶天線(xiàn)輻射和散射的主要方法之一就是采用基于混合位積分方程的空域矩量法。該方法利用基于廣義函數束方法(GPOF) [1]的二級離散復鏡像法[2]計算了微帶結構下的空域格林函數，避免了對表面波項和準動(dòng)態(tài)項的提取，大大提高了矩量法的計算效率，并采用RWG基[3]作為MoM的基函數，從而使該方法能夠很好地對大型微帶陣列以及復雜散射體的表面進(jìn)行剖分。本文將以此方法計算微帶天線(xiàn)RCS，并驗證其正確性。
　　
　　1 二級離散復鏡像方法
　　
　　本文采用基于GPOF方法的二級DCIM方法，克服了Prony方法[4]受噪聲影響大、計算精度差的缺點(diǎn)，從而提高了DCIM方法的指數擬合精度，也使DCIM方法的適用范圍增大。在DCIM方法的應用中，首先由譜域格林函數進(jìn)行漢克爾反變換得到空域格林函數：
　　
　　其中G和分別表示空域和譜域格林函數，H(2)0表示第二類(lèi)的零階漢克函數，SIP表示Sommerfeld積分路徑。本文選取一條二級積分路徑：
　　
　　從原始的函數(kρ)中減去近似函數f(kρ)所得的余函數，沿Cap2對其進(jìn)行均勻采樣。對C┆ap2而言無(wú)需很多采樣點(diǎn)，就可以很好地拾取余函數的特性。再利用GPOF方法擬合，于是有：
　　
　　2 混合位積分方程
　　
　　空域混合位積分方程(MPIE)是將混合積分方程離散化,該方法的關(guān)鍵問(wèn)題在于將譜域格林函數轉變?yōu)榭沼蚋窳趾瘮?。設E[WTBX]i為入射場(chǎng)，則混合位積分方程為：
　　
　　其中N是貼片上的基函數個(gè)數，即對微帶陣列進(jìn)行剖分時(shí)產(chǎn)生的公共邊個(gè)數。于是混合積分方程就轉化為矩陣方程：
　　
　　其中[WTHX]I[WTBX]為長(cháng)度N的電流列矢量，他包含了式(9)中的未知系數In，[WTHX]Z[WTBX]為N×N元阻抗矩陣，[WTHX]V[WTBX]是長(cháng)度為N的激勵列矢量。阻抗矩陣[WTHX]Z[WTBX]中的元素為：
　　
　　其中f[WTBX]i和f[WTBX]j分別指權函數和基函數，Ti和T′j分別指f[WTBX]i和f[WTBX]j的有效面元。通過(guò)高斯積分公式可以求得表面電流J[WTBX]，再通過(guò)互易定理求得遠場(chǎng)的方向圖。

3 計算結果及分析
　　
　　應用本文方法計算單片微帶天線(xiàn)的雷達散射截面。假設基帶厚度h=158 mm，介電常數εr=217，貼片大小在x和y方向分別為L(cháng)=366mm和W=26mm。計算其在頻率范圍為f=2~5GHz，入射角為θi=60°,φi=45°的入射波照射下，單站RCS隨頻率變化的趨勢，如圖1所示。
　　將該結果與文獻[5]中的結果進(jìn)行對比，如圖2所示，可以發(fā)現頻率在2~5 GHz時(shí)本文計算的單站RCS隨頻率的變化趨勢和文獻中的結果還是比較吻合的，從而證明了本文方法的正確性。