(電子信息系統(tǒng)復(fù)雜電磁環(huán)境效應(yīng)國家重點(diǎn)實(shí)驗(yàn)室,河南洛陽471003)

0 引言

箔條干擾[1]是使用最為廣泛的無源干擾之一,飛機(jī)和艦船等軍事目標(biāo)都經(jīng)常投放箔條彈以對抗來襲的導(dǎo)彈[2-3]。從箔條投入之初就一直有人對其進(jìn)行雷達(dá)散射特性[4-8]和對抗方法[9-18]的研究。箔條的大量使用對雷達(dá)的性能提出了嚴(yán)重挑戰(zhàn),如何對抗箔條干擾是雷達(dá)面臨的重要問題之一。

極化信息是電磁波除時域、頻域和空域信息以外的又一可資利用的重要信息,利用極化信息對抗箔條干擾很早就有人開展相關(guān)研究[13-18],其中當(dāng)箔條干擾和目標(biāo)不可分辨時,通常采用極化濾波的方法[13],當(dāng)箔條干擾和目標(biāo)可分辨時,通常采用極化識別的方法[14-18]。文獻(xiàn)[14-17]給出的極化識別方法得到了廣泛應(yīng)用,這些極化識別方法針對的是未配重的箔條干擾,采用垂直極化發(fā)射,用垂直和水平極化雙極化接收,基于兩極化通道幅度比來識別箔條干擾和艦船目標(biāo)。箔條生產(chǎn)加工技術(shù)的進(jìn)步,使得通過配重技術(shù)以控制箔條干擾的取向成為可能,原有的極化識別方法在對抗新型配重箔條干擾時的識別性能無法令人滿意。

本文正是基于此背景,對利用極化信息的箔條干擾識別方法進(jìn)行了改進(jìn)。首先對配重箔條干擾的取向統(tǒng)計模型進(jìn)行了介紹,基于此對配重箔條干擾的極化統(tǒng)計特性進(jìn)行了研究,指出了原有方法識別箔條干擾時存在的問題?；谂渲嘏c未配重箔條干擾的全極化特性,對原有的極化識別特征量進(jìn)行了改進(jìn),提出了一種改進(jìn)的箔條干擾極化識別方法。通過計算機(jī)仿真對各情況下改進(jìn)方法和原有方法的性能進(jìn)行了對比,證明了對配重后的垂直取向箔條干擾,原有方法基本上無法正確識別,而改進(jìn)方法較原有方法具有更強(qiáng)的適用性,對各種類型的箔條干擾都能以較高的概率正確識別。

1 配重箔條的極化特性分析

1.1 配重箔條干擾的取向統(tǒng)計模型

假設(shè)箔條的形狀為細(xì)長圓柱形,則箔條的取向d可以用(θd,?d)表示,其中θd∈[0,π]表示箔條姿態(tài)的俯仰角,?d∈[0,2π]表示箔條姿態(tài)的方位角。沒有配重的箔條散開后,可以認(rèn)為姿態(tài)服從均勻分布,均勻分布時箔條姿態(tài)(θd,?d)的概率分布密度函數(shù)為

箔條生產(chǎn)加工技術(shù)的進(jìn)步,使得通過配重技術(shù)可以控制箔條干擾的取向。經(jīng)過配重后,箔條的取向表現(xiàn)為俯仰向上集中于某一特定的角度,可以假設(shè)其取向的概率分布為:方位角(也即?d)在0～2π范圍內(nèi)服從均勻分布,θd的概率分布為相對集中于某一傾角θc以及一定的分散程度(也即方差D2)。對箔條而言,θd的概率分布總是關(guān)于θd=π/2水平面上下對稱的,于是可以用“θd在[0,π]范圍內(nèi)以傾角θc和π-θc為中心服從方差為D2的正態(tài)分布”來近似描述箔條取向仰角的概率分布。因此,箔條空間取向的概率密度函數(shù)[6]可以表示為

1.2 配重箔條干擾的極化特性分析

假定雷達(dá)視線方向(LOS)位于水平面內(nèi),建立如圖1如示的雷達(dá)散射坐標(biāo)系,箔條軸向在雷達(dá)散射坐標(biāo)系下的姿態(tài)用(θ,?)表示,則可以認(rèn)為θ=π/2-θd,?=?d。箔條軸向與雷達(dá)視線的夾角為Θ。將箔條等效為半波偶極子,則箔條的極化散射矩陣[8]為

式中,L為半波偶極子的等效長度,η=120πΩ為自由空間的特性阻抗,z0為偶極子的散射阻抗,cosΘ=cosθcos?,d·eh=cosθsin?,d·ev=sinθ。

圖1 雷達(dá)散射坐標(biāo)系

箔條云由大量姿態(tài)不同的箔條偶極子構(gòu)成,根據(jù)大數(shù)定理,可知箔條云單個極化通道的幅度服從瑞利分布,RCS服從指數(shù)分布,且箔條云極化散射矢量k=[SHHSVVSVH]T服從零均值的復(fù)高斯分布,其協(xié)方差矩陣為

式中,〈·〉為對箔條的所有可能的姿態(tài)求平均,N為箔條的數(shù)目。利用數(shù)值積分可以得到箔條各極化通道的平均RCS。

由圖2可以看出,當(dāng)D很大時,箔條的分布趨于均勻分布,此時θc的大小對各極化通道平均RCS的影響較小,當(dāng)D較小時,箔條的平均RCS隨θc起伏較大。當(dāng)θc趨于90°時,箔條的取向分布接近水平,此時H H通道箔條的平均RCS達(dá)到最大,當(dāng)θc和D趨于0°時,箔條的取向接近垂直,此時垂直極化通道箔條的RCS達(dá)到最大,σVV越來越大,σHV越來越小逐漸增大。

原有的極化識別方法沒有考慮箔條的配重,采取雙極化體制,發(fā)射極化為垂直接收,用垂直極化和水平極化來接收,原有識別方法的特征量為

艦船的極化散射以平板和二面角為主,電磁計算和實(shí)測數(shù)據(jù)皆表明:艦船目標(biāo)同極化通道的回波幅度遠(yuǎn)大于交叉極化通道,即

也就是說,艦船目標(biāo)的特征量ρ1一般較大,如果箔條干擾未配重的話,其特征量ρ1一般較小,但是箔條配重后,當(dāng)D較小且θc趨于0°時,ρ1較小的結(jié)論將不再成立,這意味著對于配重箔條,如果能夠?qū)⒉瓧l的取向集中于垂直,則先前的極化識別方法將會識別錯誤,誤將該類型的箔條識別為艦船。下面將提出一種改進(jìn)的極化識別方法,以保證箔條極化識別方法的穩(wěn)健性和普適性。

2 改進(jìn)的箔條干擾極化識別方法

2.1 改進(jìn)的極化識別方法

圖2 各極化通道平均RCS隨θc和D變化的三維圖

艦船目標(biāo)和箔條干擾的結(jié)構(gòu)組成存在很大差異:艦船的結(jié)構(gòu)以鐵板組成的平面散射和兩塊鐵板組成的兩面角散射為主,箔條干擾是由大量數(shù)目的箔條絲組成,箔條的取向分布對箔條云團(tuán)的散射特性具有重要的決定作用。由于單根箔條的散射波是線極化波,而線極化的方向僅與箔條的取向有關(guān),所以無論箔條是否進(jìn)行配重,取向?yàn)楹晤愋?也無法使得兩個同極化通道和交叉極化通道的回波RCS比值(即)同時很大,但對艦船目標(biāo),兩者同時較大?；谠摬町?提出改進(jìn)的箔條識別方法。

圖3 ρ2隨θc和D變化的三維圖

從圖3可以看出,ρ2的取值不超過3,這說明無論箔條的取向?yàn)楹晤愋?也無法使得同時很大,但對艦船目標(biāo)兩者同時較大。因此可以將雷達(dá)導(dǎo)引頭體制由雙極化體制改為全極化體制,在此基礎(chǔ)上對原有極化識別方法進(jìn)行改進(jìn),將識別特征量修正為

2.2 改進(jìn)方法和原有方法的性能對比

假設(shè)艦船的極化回波矢量服從零均值的復(fù)高斯分布,即箔條和艦船的分布形式相同,只是極化協(xié)方差矩陣中元素的取值不同。假設(shè)艦船的極化協(xié)方差滿足如下:

由于同極化通道之間的相關(guān)性,導(dǎo)致難以對改進(jìn)方法的識別性能進(jìn)行理論分析,本文采取計算機(jī)仿真對改進(jìn)方法和原有方法的識別性能進(jìn)行比較分析。選取幾種典型的箔條取向分布參數(shù),進(jìn)行仿真分析。根據(jù)上一節(jié)的數(shù)值積分結(jié)果可以得到箔條云極化協(xié)方差矩陣的元素取值,如表1所示,同時表1中也給出仿真時艦船目標(biāo)極化特性的參數(shù)設(shè)置。基于表1中的參數(shù)進(jìn)行仿真分析。由于式(6)和式(8)給出的識別特征量的取值范圍較大,不便于分析,所以對識別特征量進(jìn)行非線性變換(可以證明,該變換為單調(diào)遞增函數(shù),不影響識別概率),令

表1 箔條和艦船極化特性仿真參數(shù)設(shè)置

用H0假設(shè)表示雷達(dá)目標(biāo)為箔條,用H1假設(shè)表示雷達(dá)目標(biāo)為艦船,根據(jù)統(tǒng)計檢測理論可知,假設(shè)識別特征量為T,識別性能主要取決于p(T|H0)和p(T|H1)。p(T|H0)和p(T|H1)差別越大,識別性能越好。利用連續(xù)5個脈沖的回波,基于計算機(jī)仿真得到典型取向分布箔條干擾和艦船目標(biāo)的識別特征量的統(tǒng)計分布p(T1|H0)、p(T1|H1)、p(T2|H0)、p(T2|H1),如圖4所示,圖中實(shí)線表示艦船,虛線為箔條干擾,分別對應(yīng)假設(shè)H0和H1,每幅子圖的左邊為原有方法,識別特征量為T1,右邊為改進(jìn)方法,識別特征量為T2,圖4(a)、(b)、(c)分別對應(yīng)為箔條取向服從均勻分布、水平分布和垂直分布的情況。

圖4 典型取向分布箔條干擾和艦船目標(biāo)的識別特征量的統(tǒng)計分布

通過圖4可以看出,對原有方法,水平分布時p(T1|H0)、p(T1|H1)差別明顯大于其他三種分布,這說明原有方法對水平分布的識別性能較好,但垂直分布時,原有方法得到的p(T1|H0)、p(T1|H1)差別很小,即對垂直分布箔條云,原有方法得到的特征量,箔條云和艦船目標(biāo)的可分性較差。通過圖4(c)可以看出,改進(jìn)方法得到的特征量艦船和箔條的可分性明顯強(qiáng)于原有方法。對其他的幾種分布形式,原有方法和改進(jìn)方法得到的特征量,艦船和箔條的可分性大致相當(dāng),在圖4中的區(qū)別不大。根據(jù)5個脈沖的回波,選擇相應(yīng)的識別門限,使得對艦船的識別率為99%,對箔條干擾進(jìn)行識別,識別結(jié)果如表2所示。

_表2 原有方法和改進(jìn)方法的識別性能對比

由表2可以看出,對垂直分布的箔條干擾,原有方法的識別率只有0.8%,但改進(jìn)方法的識別率達(dá)到99.8%,對均勻分布和水平分布,改進(jìn)方法的識別率和原有方法相當(dāng)。這說明了本文提出的改進(jìn)方法能夠?qū)Ω鞣N類型的箔條進(jìn)行識別,較原有方法具有更強(qiáng)的適應(yīng)性。

3 結(jié)束語

本文基于配重箔條的極化特性分析,提出了一種箔條極化識別的改進(jìn)方法,通過計算機(jī)仿真證明了改進(jìn)方法較原有方法具有更強(qiáng)的適應(yīng)性。當(dāng)然,識別性能的提高需要以雷達(dá)極化信息的獲取為基礎(chǔ),必須由雙極化體制改為全極化體制。另外,本文的分析主要是基于計算機(jī)仿真,對艦船統(tǒng)計特性的假設(shè)也比較簡單,這些都有待于在下一步工作中改進(jìn)。

[1]陳靜.雷達(dá)箔條干擾原理[M].北京:國防工業(yè)出版社,2007:1-40.

[2]焦淑瑜,蓋添.火控雷達(dá)抗無源干擾能力試驗(yàn)與評估[J].雷達(dá)科學(xué)與技術(shù),2013,11(2):141-144,149.

[3]蔣志彪,駱魯秦,李鑫,等.機(jī)載箔條彈投放系統(tǒng)作戰(zhàn)效能評估[J].雷達(dá)科學(xué)與技術(shù),2013,11(3):246-250.JIANG Zhi-biao,LUO Lu-qin,LI Xin,et al.Operational Effectiveness Evaluation of Airborne Chaff Cartridge Delivery Subsystem[J].Radar Science and Technology,2013,11(3):246-250.(in Chinese)

[4]PANDEY A K.Modeling and Simulation of Chaff Cloud with Random Orientation and Distribution[C]∥IEEE MTT-S International Microwave and RF Conference,New Delhi,India:IEEE,2013:1-4.

[5]SEO D W,NAM H J,KWON O J,et al.Dynamic RCS Estimation of Chaff Clouds[J].IEEE Trans on Aerospace and Electronic Systems,2012,48(3):2114-2127.

[6]MARCUS S W.Dynamics and Radar Cross Section Density of Chaff Clouds[J].IEEE Trans on Aerospace and Eletronic Systems,2004,40(1):93-102.

[7]TANG B,SHENG X.The Simulation of Statistical Characteristics of Dense Chaff Clouds[C]∥IEEE Global High Tech Congress on Electronics,Shenzheng:IEEE,2013:191-194.

[8]李金梁,曾勇虎,申緒澗,等.均勻取向箔條云的RCS極值研究[J].雷達(dá)科學(xué)與技術(shù),2012,10(3):316-319.

[9]LI X,LIN L,SHAO X.A Target Polarization Recognition Method for Radar Echoes[C]∥International Conference on Microwave and Millimeter Wave Technology,Chengdu:IEEE,2010:1644-1647.

[10]YANG Yong,FENG De-jun,ZHANG Wen-ming,et al.Detection of Chaff Centroid Jamming Aided by GPS/INS[J].IET Radar,Sonar and Navigation,2013,7(2):130-142.

[11]ZHAN G,LU J,HU S,et al.Research on K-L Transform Based on Wavelet Transform for Anti-Chaff-Jamming[C]∥International Conference on Transportation,Mechanical and Electrical Engineering,Changchun:IEEE,2011:1500-1503.

[12]LAI Q,DAI H,ZHAO J,et al.A Novel Approach of Countering Centroid Jamming by Using INS Information in Terminal Guidance[C]∥IEEE 10th International Conference on Signal Processing,Beijing:IEEE,2010:2121-2124.

[13]楊勇,肖順平,馮德軍,等.雷達(dá)導(dǎo)引頭斜投影抗質(zhì)心干擾性能分析[J].電子學(xué)報,2014,42(3):439-444.

[14]沈允春,謝俊好,劉慶普.識別箔條云新方案[J].系統(tǒng)工程與電子技術(shù),1995,17(4):60-63.

[15]劉慶普,沈允春.箔條云極化識別方案性能分析[J].系統(tǒng)工程與電子技術(shù),1996,18(11):1-7.

[16]SHAO X,DU H,XUE J.Theoretical Analysis of Polarization Recognition Between Chaff Cloud and Ship[C]∥IEEE International Workshop on Anti-Counterfeiting,Security and Identification,Xiamen:IEEE,2007:125-129.

[17]SHAO X,DU H,XUE J.A New Method of Ship and Chaff Polarization Recognition Under Rain and Snow Cluster[C]∥IEEE International Workshop on Anti-Counterfeiting,Security and Identification(ASID),Xiamen:IEEE,2007:142-147.

[18]TANG B,LI H M,SHENG X Q.Jamming Recognition Method Based on the Full Polarisation Scattering Matrix of Chaff Clouds[J].IET Microwaves Antennas and Propagation,2012,6(13):1451-1460.