路東偉　闕肖峰　綦　鑫　聶在平

(電子科技大學電子工程學院　成都　610054)

米波頻段飛機目標全極化散射仿真與特性分析

路東偉*闕肖峰綦鑫聶在平

(電子科技大學電子工程學院成都610054)

該文基于全極化雷達對空探測和目標識別的應用需求，研究了飛機目標全極化散射特性，特別是交叉極化散射的產(chǎn)生原理和可利用性。在米波頻段對隱身和非隱身飛機目標開展了數(shù)值建模和仿真計算；分析對比了不同飛機目標的極化散射空域分布特性，分析了不同的外形結(jié)構(gòu)產(chǎn)生強交叉極化的條件。結(jié)果表明：非隱身飛機具有較豐富的極化散射特性，有利于雷達目標的極化識別。

全極化雷達；電磁散射；交叉極化；隱身目標

引用格式：路東偉, 闕肖峰, 綦鑫, 等.米波頻段飛機目標全極化散射仿真與特性分析[J].雷達學報, 2016, 5(2): 182–189.DOI: 10.12000/JR16030.

Reference format: Lu Dongwei, Que Xiaofeng, Qi Xin, et al..Simulation and analysis of the fully polarimetric scattering characteristics of aircraft in UHF band[J].Journal of Radars, 2016, 5(2): 182–189.DOI: 10.12000/JR16030.

1　引言

雷達目標極化散射特性是近年來研究的熱點。由于目標全極化散射回波中有著與目標相關的豐富信息，例如目標的幾何形狀、材料、介電常數(shù)等[1]，因此隨著電磁波極化特性的加入，使得處理接收到的目標回波信號時擴展了信息維度，能夠?qū)δ繕诉M行更加精細的識別和分析[2]。目前在電磁波極化的科學研究和工程應用中，全極化SAR研究取得了較大進展，技術較為成熟[3]。對于全極化雷達對空探測的研究則主要集中在極化信息處理，例如極化濾波[4]，相干斑去噪[5]，極化目標特征提取等[6]。相關研究仍然面臨著較大的挑戰(zhàn)，如通常情況下，雷達目標散射回波中的交叉極化分量電平遠小于主極化電平，這使得雷達難以接收到全極化信息；另外，對于目標極化散射機理研究不夠深刻，傳統(tǒng)的研究通常將電磁散射歸結(jié)為平面、二面角、球標等標準體目標的散射[7,8]，而人造目標的復雜性導致的變極化效應難以進行準確的描述。根據(jù)電磁散射理論，對于常規(guī)飛行目標，米波頻段為電磁散射的諧振區(qū)，目標的散射截面積(RCS)劇烈震蕩，回波得以增強[9]。因此在諧振頻段目標的極化散射特性是需要進一步研究的內(nèi)容。

本文首先分析了典型結(jié)構(gòu)的極化散射機理；然后采用基于積分方程––多層快速多極子算法(MultiLevel Fast Multipole Algorithm, MLFMA)對隱身和非隱身飛機目標在米波頻段開展了極化散射數(shù)值建模計算；分析對比了隱身和非隱身飛機目標的結(jié)構(gòu)和遠場極化散射特性的差異，揭示全極化回波特別是強交叉極化回波的產(chǎn)生機理，為目標全極化信息利用和特征識別提供了一定的參考。

2　典型結(jié)構(gòu)全極化散射機理

從雷達極化學的角度，電磁波與目標的作用可以表示為[10]：

式(1)中，S表示全極化散射矩陣，為目標的固有屬性并且是與目標姿態(tài)角有關的參量，其中下標x, y代表任意的正交極化基，矩陣中的每個元素刻畫了目標對不同極化電磁波的散射特性；主對角線元素的模值反映了主極化散射分量電平強度，次對角線元素模值反映了交叉極化散射分量電平強度。Ei為入射波，Es為散射波。給定散射矩陣，入射波方向改變，RCS也會變化。同樣目標姿態(tài)變化也會引起RCS的波動，這種波動會在某些角度上達到最大值，實現(xiàn)雷達目標極化自適應匹配照射[11]。此時如果接收天線極化方向與散射電場極化方向相同，可以使得接收到的目標散射回波功率最大，滿足極化自適應匹配接收條件[12]。極化信息利用的一大難點在于目標全極化散射回波中，通常情況下主極化散射分量遠大于交叉極化散射分量，這使得全極化信息利用遇到了極大的限制。但是由于目標運動時具有不同的姿態(tài)，雷達波照射角度也不盡相同，交叉極化散射回波在某些照射條件下較大，可以實現(xiàn)極化信息的有效應用。

針對交叉極化的產(chǎn)生機理，考察了典型結(jié)構(gòu)，如金屬二面角和偶極子的極化散射特性。如圖1所示的垂直放置二面角散射回波中，主極化分量占絕對優(yōu)勢，交叉極化為零。將二面角旋轉(zhuǎn)一定角度后，入射波方向不變，交叉極化散射回波可以得到顯著增強[13]。圖1(a)到圖1(b)中二面角散射矩陣變換關系為：

這里的R(θ)為旋轉(zhuǎn)矩陣，旋轉(zhuǎn)角θ=p/4時，散射回波中只有交叉極化分量，主極化分量為零。對金屬偶極子也能得到類似的結(jié)論。圖2(a)到圖2(b)中金屬偶極子散射矩陣變換過程為[14]：

圖1　二面角Fig.1　Dihedrals

圖2　金屬偶極子Fig.2　Mental dipole

3　基于積分方程的極化散射數(shù)值分析

為準確分析目標全極化散射特性，采用了基于積分方程MLFMA的數(shù)值分析方法。對于3維導體目標，首先建立電場積分方程[15]:

式(5)中，S表示目標表面，r和r'表示場點和源點位置矢量；t表示目標表面S上切向矢量，Ei表示入射電場，R=|r–r'|。并矢格林函數(shù)為：

其中，ri和rj表示場點和源點矢量，rjm, rim'表示場點和源點所在組的組中心矢量，αmm'為轉(zhuǎn)移因子。在采用迭代方法求解以上的積分方程中，涉及大量的矩矢相乘運算。MLFMA將矩矢運算按照一定互耦關系分組分層，層間嵌套，逐層遞推實現(xiàn)快速運算。MLFMA將傳統(tǒng)矩量法的計算復雜度從O(N2)量級降低至O(NlogN)[16]。

圖3　飛機目標幾何示意圖Fig.3　Airplanes

4　飛機目標全極化散射分析

如圖3所示隱身和非隱身飛機目標幾何示意圖，目標均為全尺寸模型，飛機機頭正對x軸正方向。頻率450 MHz雷達波照射角度為：(θ, φ)[45°, 135°]×[0°, 180°]。其中θ和φ為球坐標系下的坐標變量。

4.1非隱身飛機全極化散射分析

圖4給出給定角度域內(nèi)不同極化雷達波照射時的單站RCS?？梢钥闯?，主極化散射較強的區(qū)域?qū)谀繕饲跋路健?cè)向和后向側(cè)照射等，如圖4(a)和圖4(d)中φ[70°, 110°]––側(cè)面照射時RCS達到峰值；相對而言交叉極化散射分量小于主極化分量。但是對應于飛機前下方和側(cè)下方照射，在(θ, φ)[95°, 130°]×[0°, 40°](區(qū)域1)和(θ, φ)[100°, 135°]×[100°, 140°](區(qū)域2)范圍內(nèi)，對應于飛機側(cè)下方照射時，交叉極化散射分量也達到了15 dBsm；而且φθ和θφ極化回波RCS分布相同。如圖5所示，選取θ=100°, φ[0°, 15°](區(qū)域1內(nèi))和θ=114°, φ[110°, 125°](區(qū)域2內(nèi))作圖，可以看出在該角度范圍內(nèi)，交叉極化分量與主極化大小相當。

非隱身飛機具有較多的典型強散射源：進氣道、掛彈、翼面邊緣等，飛機的整體散射由各結(jié)構(gòu)散射共同產(chǎn)生。以下對其極化散射特性進行分析[17]：

(1)邊緣散射

目標機翼和垂尾對不同極化的電磁波照射非常敏感。入射波電場矢量與機翼前緣平行時會產(chǎn)生強散射，電場矢量與翼面垂直時在機翼后緣的散射較強。電磁波斜入射時，電場可以矢量分解為平行和垂直于翼面的兩個分量。散射回波經(jīng)過矢量疊加后在與入射波垂直的極化方向上有一定的投影，產(chǎn)生交叉極化。

(2)面結(jié)構(gòu)散射

目標表面法向與入射波方向平行時會產(chǎn)生極大的面散射。目標翼面、機身部分都會有面散射而產(chǎn)生強主極化分量。如圖7所示，外掛物上的翼面會產(chǎn)生角反射器，外掛物與機翼也會形成多次反射，同樣形成角反效應。隨著目標姿態(tài)的變化，角反射器旋轉(zhuǎn)，從而產(chǎn)生很強的交叉極化。

(3)表面電磁缺陷散射

金屬目標表面不連續(xù)結(jié)構(gòu)截斷電流從而產(chǎn)生電磁波反射。由于電流矢量方向受到目標介質(zhì)或者結(jié)構(gòu)的影響，產(chǎn)生變極化效應，如圖8(a)所示的不連續(xù)結(jié)構(gòu)。

(4)進氣道和隔道散射

目標進氣道和進氣道與機身之間存在的隔道也會產(chǎn)生較強的后向散射回波。機身與進氣道之間的隔道實際上形成了一個腔體，容易產(chǎn)生諧振引起回波。隨著腔體內(nèi)表面曲率的變化，從而產(chǎn)生復雜的回波極化分量。

圖4　非隱身目標的全極化散射特性Fig.4　Fully polarimetric scattering characteristics of non-stealth aircraft

圖5　較強交叉極化散射區(qū)域雷達散射回波Fig.5　Cross-polarimetric RCS curves

以上的強散射源在雷達波照射下產(chǎn)生很強的主極化散射回波。目標在空中運動時，其相對于雷達的空間位置會不斷改變。類似之前分析的二面角和金屬偶極子產(chǎn)生交叉極化的原理，這些強散射源旋轉(zhuǎn)后，雷達波入射后產(chǎn)生較強的交叉極化回波。

為了進一步分析主極化散射和交叉極化散射關于目標姿態(tài)照射角度的關系，將交叉極化RCS值與主極化RCS值相減。僅選取主極化RCS值大于–20 dBsm的結(jié)果進行作圖。從圖9中可以看出，非隱身飛機的交叉極化較強區(qū)域的分布較廣。

4.2隱身飛機目標全極化散射分析

對于圖3(b)中的隱身飛機目標采用同樣的方法進行分析，得到如圖10所示的全極化散射回波結(jié)果?？梢钥闯觯?jīng)過外形隱身設計，威脅角度范圍內(nèi)的回波RCS得到了較大縮減，即使在米波頻段的主極化分量也在0 dBsm左右。強散射分布在側(cè)向和尾向，如圖10(a)和圖10(d)中的黑框部分角度。對于交叉極化散射分量，較強交叉極化散射點相比非隱身飛機大為減少，僅在(θ, φ)∈[0°, 80°]×[60°, 100°] (圖10(b)和圖10(c)中虛線框部分)等側(cè)向以及斜后方角度范圍內(nèi)分布，主要是由機翼、垂尾與機身形成的二面角在電磁波斜入射時產(chǎn)生。同時強交叉極化散射點分布比較離散，數(shù)目較少并且強度比非隱身飛機弱。如圖11所示，在θ=46°, φ∈[80°, 100°]的角度內(nèi)，交叉極化較高，但是其仍然不超過主極化RCS值。經(jīng)過外形隱身優(yōu)化的目標，其強散射源較少，即使目標位置變化，產(chǎn)生的交叉極化分量比非隱身目標小很多。

圖6　目標翼面Fig.6　Wings

圖7　目標外掛物(導彈)Fig.7　Cruise missle

圖8　目標表面不連續(xù)結(jié)構(gòu)和隔道Fig.8　Discontinuous surface and bump inlet

圖9　交叉極化與主極化RCS差值(主極化大于–20 dBsm)Fig.9　The difference between the cross-polarimetric and the main-polarimetric RCS of non-stealth aircraft (The strength of main-polarimetric RCS is greater than –20 dBsm)

圖10　隱身目標的全極化散射特性Fig.10　Stealth aircraft full-polarimetric scattering characteristics

圖12對隱身飛機的交叉極化與主極化分量差值進行分析?？梢钥闯?，隱身目標在前向和側(cè)向(主要威脅方向)幾乎沒有形成強的交叉極化散射，特別是對于φ極化。較強交叉極化散射多數(shù)集中在目標側(cè)后向區(qū)域。

圖11　隱身目標RCS曲線(θ=46°，φ[80°, 100°])Fig.11　Stealth aircarft RCS curves

圖12　隱身目標交叉極化與主極化RCS差值(主極化大于–20 dBsm)Fig.12　The difference between the cross-polarimetric and the main-polarimetric RCS of steath aircraft (The strength of main-polarimetric RCS is greater than –20 dBsm)

從兩類飛機極化散射仿真結(jié)果可以看出，非隱身目標外形結(jié)構(gòu)復雜，構(gòu)成多種強散射源。在其運動姿態(tài)角度變化時能夠激勵起較強的全極化回波，其中產(chǎn)生強交叉極化的原因可以認為是典型散射結(jié)構(gòu)在電磁波斜入射時等效為原散射矩陣與旋轉(zhuǎn)因子相乘，使得次對角線元素不為零，增強了交叉極化散射。對于隱身目標，氣動外形結(jié)構(gòu)簡單平滑，強散射結(jié)構(gòu)較少。不僅威脅方向RCS極低，交叉極化分量也僅在目標側(cè)向和上、下側(cè)后方向有一定分布。

5　結(jié)論

本文首先分析了典型結(jié)構(gòu)的極化散射產(chǎn)生機理，運用積分方程數(shù)值方法對隱身和非隱身飛機目標在米波頻段進行極化散射仿真計算。可以看出，通常情況下目標散射場的主極化分量大于交叉極化分量。但是對于非隱身目標，電磁波在一定角度照射時，其交叉極化散射回波較強。而經(jīng)過外形優(yōu)化設計的隱身目標，減少了能夠引起強散射的幾何結(jié)構(gòu)，交叉極化散射特性較弱。本文研究內(nèi)容為目標全極化信息利用和特征識別提供重要參考。

[1]陳強.雷達極化中若干理論問題研究[D].[博士論文], 國防科技大學, 2010.Chen Q.Study on theoretical problems in radar polarimetry[D].[Ph.D.dissertation], National University of Defense Technology, 2010.

[2]李金梁, 李永禎, 王雪松.米波極化雷達的反隱身研究[J].雷達科學與技術, 2005, 3(6): 321–326.Li J L, Li Y Z, and Wang X S.Study on anti-stealth with meter-band polarimetric radar[J].Radar Science and Technology, 2005, 3(6): 321–326.

[3]邢世其.人造目標極化雷達三維成像理論與方法研究[D].[博士論文], 國防科技大學, 2013: 3–4.Xing S Q.Study on the 3D imaging of manmade targets based on polarimetric radar[D].[Ph.D.dissertation], National University of Defense Technology, 2013: 3–4.

[4]王雪松.寬帶極化信息處理的研究[D].[博士論文], 國防科技大學, 1999: 4–6.Wang X S.Study on wideband polarization information processing[D].[Ph.D.dissertation], National University of Defense Technology, 1999: 4–6.

[5]崔一.基于SAR圖像的目標檢測研究[D].[博士論文], 清華大學, 2011: 37–57.Cui Y.Study on target detection in synthetic aperture radar images[D].[Ph.D.dissertation], Tsinghua University, 2011: 37–57.

[6]代大海.極化雷達成像及目標特征提取的研究[D].[博士論文],國防科技大學, 2008: 49–71.Dai D H.Study on polarimetric radar imaging and target feature extration[D].[Ph.D.dissertation], National University of Defense Technology, 2008: 49–71.

[7]Jakob Van Zyl.Synthetic Aperture Radar Polarimetry[M].Hoboken, NJ, USA: Wiley, 2011.

[8]安文韜.基于極化SAR的目標極化分解與散射特征提取研究[D].[博士論文], 清華大學, 2010: 11–24.An W T.The polarimetric decomposition and scattering charasteristic extraction of polarimetric SAR[D].[Ph.D.dissertation], Tsinghua University, 2010: 11–24.

[9]劉俊.米波雷達低仰角估計方法研究[D].[博士論文], 西安電子科技大學, 2012: 3–5.Liu J.Study on low-angle estimation in vhf radar[D].[Ph.D.dissertation], Xidian University, 2012: 3–5.

[10]李謙, 林昌祿.雷達極化理論中常見極化散射矩陣特性分析[J].電子科技大學學報, 1993, 22(5): 490–500.Li Q and Lin C L.Analysis of characteristic properties of the common polarization scattering matrices in radar polarization theory[J].Journal of University of Electronic Science and Technology of China, 1993, 22(5): 490–500.

[11]劉志文, 柯有安.雷達反隱身的若干問題與技術途徑[J].現(xiàn)代雷達, 1992, 14(3): 1–8.Liu Z W and Ke Y A.On radar anti-stealth techniques and its problems[J].Modern Radar, 1992, 14(3): 1–8.

[12]林世明, 楊健.目標散射矩陣的特征值理論和雷達天線的最優(yōu)極化[J].電波科學學報, 1995, 10(1/2): 11–15.Lin S M and Yang J.Eigenvalue theory of scattering matrix and optimal polarization of radar antenna[J].Chinese Journal of Radio Science, 1995, 10(1/2): 11–15.

[13]常宇亮, 文玲, 戴幻堯.目標交叉極化散射特性的強度分析[J].雷達科學與技術, 2012, 10(5): 544–548.Chang Y L, Wen L and Dai H Y.Analysis of crosspolarization scattering strength from target[J].Radar Science and Technology, 2012, 10(5): 544–548.

[14]Jong-Sen Lee and Eric Pottier.Polarimetric Radar Imaging: From Basis to Applications[M].Taylor & Francis Group, 2012: 93–95.

[15]聶在平, 胡俊, 姚海英, 等.用于復雜目標三維矢量散射分析的快速多極子方法[J].電子學報, 1999, 27(6): 104–109.Nie Z P, Hu J, Yao H Y, et al..The fast multipole methods for vector analysis of scattering from 3-Dimensional objects with complex structure[J].Acta Electronica Sinica, 1999, 27(6): 104–109.

[16]盛新慶.計算電磁學要論[M].第2版, 合肥: 中國科學技術大學出版社, 2008: 31–39.Sheng X Q.A Brief Treatise on Computational Electronmagnetics[M].The Second Edition, Hefei: Press of University of Science and Technology of China, 2008: 31–39.

[17]桑建華.飛行器隱身技術[M].北京: 航空工業(yè)出版社, 2013: 63–68.Sang J H.Low-Obserbable Technologiese of Aircraft[M].Beijing: Aviation Industry Publishing House, 2013: 63–68.

路東偉(1992–)，男，甘肅蘭州人，于2014年獲得電子科技大學電子信息工程專業(yè)學士學位，并于同年成為電子科技大學微波工程系在讀碩士生。目前從事電磁散射與逆散射、雷達目標特性、極化雷達成像等方面的研究。

E-mail: bookwormldw@qq.com

闕肖峰(1977–)，男，重慶梁平人，于2001年和2008年獲得電子科技大學電磁場與微波技術學士和博士學位，目前為電子科技大學副教授，IEEE會員。主要從事計算電磁學、雷達目標特性、電磁散射與逆散射、天線理論與技術等的研究。

E-mail: quexiaofeng@163.com

綦鑫(1990–)，男，山東青島人，于2014年獲得電子科技大學電子信息工程專業(yè)學士學位，同年保送進入電子科技大學微波工程系攻讀博士學位。目前從事計算電磁學、雷達目標特性、電磁散射與逆散射等的研究。

E-mail: qixin9098@qq.com

聶在平(1946–)，男，陜西西安人，教授，博士生導師，中國電子學會以及IEEE會士，1981年獲電子科技大學碩士學位，1987～1989年在美國伊利諾伊大學電磁研究室從事研究工作。長期致力于電磁散射與輻射、非均勻介質(zhì)中的場與波、計算電磁學的研究，先后主持并且圓滿完成國家重點科技攻關項目、國家863以及自然科學基金重點項目，國防973項目等數(shù)十個研究課題，并且取得了多項重要成果。發(fā)表學術論文超過300篇，同時也是我國《天線工程手冊》副主編。

E-mail: zpnie@uestc.edu.cn

Simulation and Analysis of the Fully Polarimetric Scattering Characteristics of Aircraft in UHF Band

Lu DongweiQue XiaofengQi XinNie Zaiping
(School of Electronic Engineering, University of Electronic Science and Technology of China, Chengdu 610054, China)

This study is based on the application requirements of fully polarimetric radar for aircraft detection and target recognition.We focus on fully polarimetric scattering characteristics of air targets, particularly the cross-polarimetric scattering mechanism and its availability.In the ultra-high frequency band, we conduct numerical modeling and simulation of stealth and non-stealth aircraft targets.The spatial distribution characteristics of polarimetric scattering for different aircraft targets are studied and compared, and the structures that could cause strong cross polarization are analyzed.The results of this study suggest that nonstealth aircraft have more polarimetric characteristics; this will help people use polarimetric radar for detecting aircraft.

Fully polarimetric radar; Electromagnetic scattering; Cross-polarimetric; Stealth aricraft

TN958

A

2095-283X(2016)02-0182-08

10.12000/JR16030

2016-01-31；改回日期：2016-04-07；網(wǎng)絡出版：2016-04-20

路東偉bookwormldw@qq.com

國家自然科學基金(61490695, 61231001)

Foundation Items: The National Natural Science Foundation of China (61490695, 61231001)