基于極化特性分析的箔條干擾識別方法

焦丹丹， 劉 崢， 呂方方， 張艷艷

(西安電子科技大學雷達信號處理國家重點實驗室， 陜西西安 710071)

0 引言

箔條具有廉價、制造簡單、使用方便、強干擾性等優點，這使得箔條成為最早使用和使用最廣泛的無源干擾物[1]。箔條在戰場上，尤其是海戰場上的使用為艦船目標免遭反艦導彈的襲擊迎來了曙光，但與此同時大大降低了導彈的命中率。因此，箔條的使用推動了抗箔條干擾技術的發展。

目前，抗箔條干擾主要是在分析箔條干擾和艦船回波在不同維度特性差異的基礎上，提取出不同的特征量進行識別[2]。目標的極化特性在識別真假目標方面可以提供更加豐富的信息，如何充分利用目標的極化特性，根據箔條干擾和艦船回波在極化域的特性差異，選取出性能優良的特征量是研究者迫切需要解決的問題。對于現有的極化識別方法，選取的特征量大多針對某種特定類型的箔條干擾或者某些特定條件下的箔條干擾時才有效，特征量的穩健性和魯棒性較差，這是由于現有特征量會隨雷達參數、環境參數、干擾自身參數的變化而改變，從而在一定程度上影響雷達識別目標的性能。因此，系統分析所選特征量的極化特性，定量評估特征量在區分真假目標時的可分性性能就顯得尤為重要。

針對上述問題，本文做了如下工作：首先，對現有極化角均值特征量進行了特性分析，指出了現有極化角均值特征量存在的缺點與不足，在對現有極化角均值特征量修正的基礎上，得到最小極化角均值特征量；其次，對最小極化角均值特征量進行了特性分析，系統討論了該特征量隨敏感參數變化的一般性規律；最后，用可分性測度函數定量分析了最小極化角均值特征量在區分艦船回波和箔條干擾時的可分性性能。研究結果表明：最小極化角均值特征量比現有極化角均值特征量性能良好，從而有望為新型雷達系統設計與抗干擾技術的研究提供理論參考與技術支持。

1 特征量的選取

1.1 研究對象

常用箔條干擾類型有圓柱形箔條、鋁箔條和V形箔條[1]，如圖1所示。本文研究高頻波段下的雷達面臨的箔條干擾，由于高頻波段雷達使用頻率多在X波段及以上，而V形箔條適用于S波段以下，因此，本文對該類型的箔條不予以研究。

圖1 箔條種類

表1為本文的箔條研究對象。實際應用中使用最多的是圓柱形箔條，特別是鍍鋁玻璃絲箔條，下文若無特殊說明，均以鍍鋁玻璃絲箔條為研究案例展開相關研究。

表1 箔條研究對象

1.2 研究模型

箔條回波特性受箔條空間取向統計分布的影響。常用箔條姿態服從球面均勻分布和雙峰正態分布兩種模型[3-4]。假設箔條空間取向矢量為d=(θ，φ)，其中θ為箔條在空中的姿態俯仰角，φ為箔條在空中的姿態方位角。

當d服從球面均勻分布時，其分布函數為

(1)

實際使用中，由于箔條參數精心設計及箔條投擲方式各異等原因，箔條姿態傾角分布更貼近于雙峰正態分布，即箔條取向矢量d中的姿態方位角φ服從[0，2π]上的均勻分布，另一個姿態傾斜角θ用“θ在[0，π]范圍內以傾角θc和π-θc為中心，服從標準差為D的正態分布”來刻畫[5]。

當d服從雙峰正態分布時，其分布函數為

(2)

式中，ζ為

(3)

實際上，箔條服從上述兩種分布是對立統一的，這是由于當式(2)中的方差D趨于無窮大時，式(2)刻畫的雙峰正態分布趨向于式(1)刻畫的球面均勻分布。

由于某一種箔條空間取向的統計模型不能準確代表所有箔條在空中分布的情況，因此，本文選取了4種典型的箔條姿態進行研究，如表2所示。表2是箔條服從4種典型空間取向模型時的具體參數。

表2 箔條取向研究模型

1.3 極化角均值特征量特性分析

共垂直極化通道的極化角γ1[6]和共水平極化通道的極化角γ2定義為

(4)

(5)

式中，π表示圓周率，arctan表示反正切函數，AVV，AHV，AHH依次表示共垂直極化通道、交叉極化通道和共水平極化通道的回波振幅。

共垂直極化通道的極化角均值ξ1[6]和共水平極化通道的極化角均值ξ2為

(6)

(7)

顯然，要研究箔條極化角均值特征量，首要任務就是研究極化雷達體制下的箔條平均RCS特性。以箔條中心為原點，正南方向為x軸、正東方向為y軸、豎直方向為z軸，建立地理坐標系下箔條散射示意圖，如圖2(a)所示。

設箔條取向矢量為d=(θ,φ)，姿態傾角θ是箔條與z軸方向的夾角，姿態方位角φ是箔條在xoy平面內的投影量與x軸的夾角。類似地，設散射場nR在地理坐標系中的方向矢量為nR=(θ1,φ1)，入射電場nT方向為nT=-nR，θr為雷達入射余角，則有關系式[5]：

(8)

式中，eh，ev依次代表單位水平極化方向和單位垂直極化方向。

以箔條中心為原點，以單位水平極化方向eh為x軸，以單位垂直極化方向ev為y軸，以入射電場方向nT為z軸，建立極化坐標系下箔條散射示意圖，如圖2(b)所示。

圖2 箔條散射示意圖

箔條在圖2(b)中的方向矢量為(Θ,α)，其中，Θ是箔條與z軸方向的夾角，α是箔條在xoy平面內的投影量與x軸的夾角。則有關系式[5]：

(9)

通常用Stokes參數刻畫箔條的散射極化特性[5]。由于成熟期的箔條RCS很大，具備保護真實目標的能力，因此本文研究對象為成熟期箔條，該時期箔條云散射回波Stokes矢量為單根箔條Stokes矢量之和，箔條云團平均RCS是單根箔條平均RCS之和，即

(10)

式中，N為箔條云團中箔條的數目，σvv，σhv，σhh依次表示單根箔條共垂直極化通道、交叉極化通道和共水平極化通道的散射回波Stokes矢量平均值[5]。根據箔條雷達回波的互易性，σHV=σVH。

若雷達的收/發極化方式為垂直或水平時，則有[5]

(11)

式中，Mij為單根箔條平均Mueller矩陣中第i行第j列的元素，計算公式為[5]

Mij=?M′ij(θ,φ)f(θ,φ)dΩ

(12)

式中，f(θ,φ)為箔條空間取向為(θ,φ)時的概率密度函數，表達式為式(1)和式(2)，dΩ=sinΘdθdφ，M′ij為單根箔條極化散射Mueller矩陣中第i行第j列的元素，表達式為[5]

(13)

式中，η為自由空間特性阻抗，η=120π Ω，L為半波偶極子的等效長度，L=λ/π，λ為雷達入射波長，z0為半波偶極子的輻射阻抗，z0=731 Ω。

根據上述分析過程，可得極化雷達體制下，共垂直極化通道極化角均值ξ1和共水平極化通道極化角均值ξ2的理論值的計算步驟：

1) 根據式(8)和式(9)解出Θ(θ,φ)，α(θ,φ)的數值表達式。

2) 將Θ(θ,φ)，α(θ,φ)代入式(13)計算單根箔條極化散射Mueller矩陣中第i行第j列的元素M′ij。

3) 將箔條空中取向分布模型的概率密度函數f(θ,φ)以及M′ij代入式(12)中，求解單根箔條平均Mueller矩陣中第i行第j列的元素Mij。

4) 將Mij代入式(11)求解單根箔條共垂直極化通道、交叉極化通道和共水平極化通道的散射回波Stokes矢量平均值，即σvv，σhv和σhh。

5) 將σvv，σhv，σhh代入式(10)可得整個箔條云團共垂直極化通道、交叉極化通道和共水平極化通道的平均RCS，即σVV，σHV，σHH。

6) 將σVV，σHV，σHH代入式(6)和式(7)中，得共垂直極化通道極化角均值ξ1和共水平極化通道極化角均值ξ2的理論值。

1.4 最小極化角均值特征量的選取

根據上述1.3節的理論分析過程可知，影響極化角均值的因素有：雷達入射余角θr、雷達入射波長λ和箔條空中取向分布函數f(θ,φ)(簡稱箔條姿態)。下面根據單一變量原則，結合極化角均值ξ1,ξ2的理論計算步驟，討論共垂直極化通道極化角均值ξ1和共水平極化通道極化角均值ξ2隨上述參數變化的理論變化規律。

圖3是箔條服從模型3時，極化角均值ξ1和ξ2分別隨雷達入射余角、雷達入射波長、箔條姿態的變化圖。其中，虛線代表極化角均值ξ1，實線代表極化角均值ξ2，箔條根數為1萬根，圖3(a)中雷達載頻為17 GHz，圖3(b)中雷達入射余角為45°，圖3(c)中雷達載頻為17 GHz，雷達入射余角為45°。

圖3 箔條服從模型3時極化角均值隨參數變化規律

由圖3(a)可知，對于垂直極化箔條(即箔條服從模型3)，在低入射余角時，現有極化角均值ξ1取值很大，但此時極化角均值ξ2取值很小；由圖3(c)可知，當影響箔條姿態分布的參數θc和D都趨向于0，即箔條趨向垂直正態分布時，現有極化角均值ξ1取值高達88°以上。若基于箔條極化角均值小于艦船極化角均值的判斷準則，用現有極化角均值ξ1來區分真假目標，則會出現判斷失誤的情況。因此我們需要對現有極化角均值特征量ξ1進行修正。

實測數據及文獻資料皆表明[7]，艦船極化角均值主要在60°以上，而通過上述討論可知，無論雷達入射余角、雷達入射波長、箔條姿態如何變化，極化角均值ξ1和ξ2無法同時很大。基于這種差異，可用最小極化角均值特征量ξ來區分箔條和艦船。ξ定義為

ξ=min{ξ1,ξ2}

(14)

2 最小極化角均值特征量特性分析

2.1 最小極化角均值特征量的特性分析

下面采用“理論分析+仿真驗證”的方法，討論箔條的最小極化角均值特征量ξ(簡稱ξ)隨雷達入射余角、雷達入射波長、箔條姿態的變化規律。其中，理論分析方法與1.3節中的分析方法一致，仿真驗證中箔條回波數據采用計算機仿真方法獲取[4]。

1) 雷達入射余角

圖4 (a)～(d)是箔條服從模型1～模型4時ξ隨雷達入射余角θr的變化圖。其中，雷達載頻為17 GHz。從圖4可知，隨θr的增大，ξ不變或是緩慢增大，ξ取值范圍為[18,58)。

圖4 ξ隨雷達入射余角變化圖

2) 雷達入射波長

圖5是箔條服從模型1～模型4時ξ隨雷達入射波長λ的變化圖,其中，雷達入射余角θr為45°。從圖5可知，ξ不隨λ的變化而變化，ξ取值范圍為(24,58]。

圖5 ξ隨入射波長變化圖

3) 箔條姿態

圖6是ξ隨箔條姿態的變化圖。其中，雷達載頻為17 GHz，入射余角為45°。從圖6可知，ξ隨箔條姿態的改變而變化，ξ取值范圍為[5,58)。

圖6 ξ隨箔條姿態(θc和D)變化圖

在圖6中，選取了4種典型的箔條姿態(模型1～模型4)進行仿真驗證，理論分析結果與仿真驗證結果基本吻合，如表3所示。

表3 箔條最小極化角均值實驗理論對比

另外，圖7為箔片1、箔片2服從模型4時，最小極化角均值特征量ξ隨雷達入射余角的變化情況。

圖7 箔片ξ隨雷達入射余角變化圖

從圖7可知，無論是鍍鋁玻璃絲箔條、箔片1或是箔片2，最小極化角均值特征量ξ隨特定參數的變化規律基本一致，與箔條型號幾乎無關。

2.2 最小極化角均值特征量特性分析總結

由于2.1節研究過程僅分析了最小極化角均值特征量ξ隨單一參數的變化規律，這種變化規律是否具有普適性，需要進一步深入研究。因此本節基于大量的理論分析以及仿真驗證歸納總結了最小極化角均值特征量ξ隨參數變化的一般性規律，如表4所示。

表4 箔條最小極化角均值特征量特性分析總結結果

從表4可知，無論雷達入射余角為多少、雷達載頻為多大、箔條空中取向如何，箔條回波的最小極化角均值特征量ξ的取值范圍主要在60°以下，與理論分析結果吻合。相比之下，艦船回波的最小極化角均值特征量ξ的取值范圍主要在60°以上。

上述差異是由二者的結構組成造成的。對于艦船目標，其極化散射分量主要由平板和二面角組成[7]，艦船回波同極化與交叉極化的比值遠大于1，即式(6)和式(7)中的ρ1?1，ρ2?1，通過式(6)和式(7)計算可知，艦船回波ξ1和ξ2的取值主要在60°以上，進一步根據式(14)可知，箔條最小極化角均值ξ的取值主要在60°以上，與現有電磁計算和實測數據結果吻合。同理，箔條云團中單根箔條的散射是線極化波[7]，從而箔條回波同極化與交叉極化的比值無法同時很大，即箔條回波的ρ1和ρ2取值無法同時很大，計算可知，箔條干擾ξ1和ξ2的取值無法都在60°以上，進一步根據式(14)可知，箔條最小極化角均值ξ的取值主要在60°以下，與現有箔條統計特性和實測數據結果吻合。

綜上所述，可利用箔條干擾和艦船回波在最小極化角均值特征量ξ上的取值差異，區分箔條和艦船。

3 最小極化角均值特征量可分性分析

為了觀察最小極化角均值特征量ξ在區分艦船和箔條時的區分效果，本文以鍍鋁玻璃絲箔條為研究案例，給出了雷達入射余角為10°，雷達載頻為17 GHz，艦船回波和箔條回波用最小極化角均值特征量ξ處理后的對比結果，如圖8、圖9所示。其中，箔條回波數據采用計算機仿真方法獲取[4]，艦船回波數據用FEKO電磁計算軟件獲得[8]。

圖8 箔條與艦船極化角均值處理之后的對比圖

對比圖8和圖9可知，若基于箔條極化角均值小于艦船極化角均值的判斷準則，現有極化角均值特征量確實在對抗垂直極化箔條干擾時判斷失誤，而最小極化角均值特征量ξ可區分任何取向的箔條干擾和艦船回波。

圖9 箔條與艦船本文特征量處理之后的對比圖

特征可分性測度是定量描述不同特征量優劣程度的數學工具，為衡量最小極化角均值特征量在區分真假目標中的貢獻程度，本文用可分性測度函數J2作為度量工具，定量分析本文特征量在區分艦船回波和箔條干擾時的可分性性能。J2定義為

(15)

式中，tr稱為矩陣的跡，Sb代表類間離差，Sw代表類內離差。

類內離差Sw是各類模式協方差矩陣的先驗概率加權值平均，其表達式為[9]

(16)

類間距離Sb是每一類模式均值向量與模式總體均值向量之間平方距離的先驗概率加權和[9]，表達式為

(17)

式中，i=1,2,…,c0。tr(Sb)是每一類均值與全局均值之間平均距離的一種測度。

用J2作為可分性測度函數，對圖8(b)和圖9中最小極化角均值特征量ξ的可分性性能進行定量分析，可分性結果如表5所示。原則上，選取的特征量應該使類內分散度盡量小，類間分散度盡可能大，即特征量的可分性測度值越大，特征量的可分性性能越好。從表5可知,最小極化角均值特征量的可分性測度值在箔條服從不同模型時取值不同，但可分性測度值最小也在20以上。結合圖8、圖9以及表5可知，最小極化角均值特征量的可分性較好。

表5 箔條最小極化角均值特征量的可分性測度值

6 結束語

本文根據箔條空中取向分布模型以及箔條極化信息，給出了共垂直極化通道和共水平極化通道箔條極化角均值的理論計算步驟，通過對現有極化角均值特征量和最小極化角均值特征量特性分析的基礎上，基于大量的仿真回波數據驗證了最小極化角均值特征量特性分析結果的正確性和可靠性，同時也證明了最小極化角均值特征量在區分真假目標時的穩健性和魯棒性更好。關于采用可分性測度函數定量分析特征量在區分箔條干擾和艦船回波時的貢獻程度時，可分性測度值多大，特征量的可分性性能好，這屬于抗箔條干擾識別算法的研究范疇，也是后續工作的一個重要研究方向。