鄭 博, 杜 科

(上海無線電設(shè)備研究所,上海200090)

0 引言

戰(zhàn)場電磁環(huán)境日益復(fù)雜,干擾技術(shù)迅猛發(fā)展,不斷向數(shù)字化、多樣化、智能化方向前進(jìn)。種類繁多的干擾已經(jīng)嚴(yán)重威脅PD體制雷達(dá)制導(dǎo)系統(tǒng)的精確制導(dǎo)。為使我國PD雷達(dá)導(dǎo)引頭在復(fù)雜電子戰(zhàn)環(huán)境中充分發(fā)揮作用,開展雷達(dá)導(dǎo)引頭抗干擾技術(shù)的研究是十分緊迫和必要的[1]。

自衛(wèi)式欺騙干擾是精確模擬真實目標(biāo)回波的一種干擾信號。敵機通過欺騙導(dǎo)彈的量測距離、量測速度等信息影響制導(dǎo)精度,從而達(dá)到保護(hù)自身目的?，F(xiàn)階段對于自衛(wèi)式欺騙干擾的研究僅限于單重速度、距離欺騙干擾[2～4],而多重速度、距離欺騙干擾的研究較少。

本文提出一種PD雷達(dá)導(dǎo)引頭抗欺騙干擾信息處理新技術(shù),利用回波提取的多維信息(距離、速度、角度、幅度)以及多幀信號之間的相關(guān)性,通過數(shù)據(jù)關(guān)聯(lián)、多點跡建立、點跡特征提取、目標(biāo)/干擾識別等過程有效抵抗距離、速度欺騙干擾(單重/多重)。

1 自衛(wèi)式速度、距離欺騙干擾特征

設(shè)PD雷達(dá)導(dǎo)引頭的發(fā)射波形為

式中:f0為載頻頻率;μ為線性調(diào)頻斜率(μ=0則為普通脈沖)。

多重距離拖引干擾,其m個拖引干擾回波信號波形為

式中:τi(t)為距離欺騙干擾的延時量,在檢測過程中體現(xiàn)在時域上出現(xiàn)多個距離不同的目標(biāo)。圖1為多拖速距離拖引干擾。

圖1 多拖速距離拖引干擾

而多重速度拖引干擾,其n個拖引干擾回波信號波形為

式中:f di為速度欺騙干擾的頻偏量,在檢測過程中體現(xiàn)在多普勒頻域上出現(xiàn)多個速度不同的目標(biāo)。圖2為多拖速速度拖引干擾。

現(xiàn)有導(dǎo)引頭抗干擾效果評估中,對于距離、速度欺騙干擾,只要導(dǎo)引頭可以實施角跟蹤就認(rèn)為是有效抗干擾。但是錯誤的距離、速度信息仍然會影響制導(dǎo)精度。也有學(xué)者利用干擾信號能量較大的特點,干擾出現(xiàn)時,AGC電壓突變,將大信號刪除來抗干擾。然而,當(dāng)多個干擾出現(xiàn)時,現(xiàn)有抗干擾方法將束手無策。

圖2 多拖速速度拖引干擾

2 信息關(guān)聯(lián)抗干擾算法

多重欺騙干擾在其距離或速度曲線中可以明顯區(qū)分(如圖1、圖2所示)。故將雷達(dá)航跡理論融入到導(dǎo)引頭數(shù)據(jù)處理中,利用PD雷達(dá)導(dǎo)引頭各個時刻測量到的導(dǎo)彈-目標(biāo)相對距離、導(dǎo)彈-目標(biāo)相對速度、目標(biāo)相對導(dǎo)引頭天線的角誤差等信息進(jìn)行多幀信息關(guān)聯(lián),形成目標(biāo)和導(dǎo)彈的相對點跡。

在PD雷達(dá)導(dǎo)引頭跟蹤目標(biāo)過程中,將目標(biāo)回波信息進(jìn)行關(guān)聯(lián),形成目標(biāo)和導(dǎo)彈的相對點跡。當(dāng)欺騙干擾出現(xiàn)時,干擾的信號攜帶的導(dǎo)彈-目標(biāo)相對距離、導(dǎo)彈-目標(biāo)相對速度、目標(biāo)相對導(dǎo)引頭天線的角誤差等信息與目標(biāo)點跡不同,從而得到不同點跡。從點跡上可以很容易確定干擾的出現(xiàn)。

由于點跡攜帶了豐富的目標(biāo)信息,利用目標(biāo)與干擾之間特征差異,從而識別目標(biāo),剔除干擾。

2.1 信息關(guān)聯(lián)抗干擾流程

在脈沖多普勒雷達(dá)導(dǎo)引頭中利用信息關(guān)聯(lián)思想抗干擾,加強了數(shù)據(jù)處理部分功能,而信號檢測部分保持不變,主要包括:回波信號采樣,距離門重排,按距離門FFT,時頻二維恒虛警等模塊[5]。而數(shù)據(jù)處理部分則見下圖虛線部分,主要包括:信號關(guān)聯(lián)模塊、跟蹤起始/終止模塊、濾波預(yù)測模塊、跟蹤波門計算等部分。圖3所示為信息關(guān)聯(lián)抗干擾流程圖。

圖3 信息關(guān)聯(lián)抗干擾流程圖

2.2 信息關(guān)聯(lián)抗干擾功能模塊

2.2.1 信息關(guān)聯(lián)模塊

(1)信號關(guān)聯(lián)

信號關(guān)聯(lián)主要任務(wù)是將每幀檢測到的量測與已有點跡關(guān)聯(lián);在沒有點跡的情況下,則將同一目標(biāo)的多幀回波信號進(jìn)行關(guān)聯(lián),進(jìn)行點跡建立。

目前常用的信號關(guān)聯(lián)算法主要有[6]:聯(lián)合概率數(shù)據(jù)關(guān)聯(lián)算法(JPDA),多假設(shè)方法(MH T)以及密集雜波環(huán)境下的數(shù)據(jù)關(guān)聯(lián)快速算法(FAFDA)。不同算法因其假設(shè)前提不同,存在不同的適應(yīng)條件,以及各自的利弊,在實際應(yīng)用中則應(yīng)該按具體情況合理選擇。表1所示為關(guān)聯(lián)算法特點總結(jié)。

表1 關(guān)聯(lián)算法比較

(2)跟蹤波門

本文采用的相關(guān)波門形成方法是在跟蹤空間內(nèi)定義一個矩形區(qū)域,即矩形波門。

假定量測維數(shù)為nz,新息向量v(k)、量測向量z(k)及量測的預(yù)測向量的第i個分量分別用表示,新息協(xié)方差S(k)的第i行第j列元素用Sij表示,當(dāng)量測z(k)的所有分量均滿足關(guān)系時,稱該量測為候選回波。其中:K G為波門常數(shù),在實際應(yīng)用中往往取較大的值(K G≥3.5)。

新息向量vi(k)的范數(shù)為

在一定的假設(shè)條件下,可以證明,g(k)是服

假定高斯誤差模型成立,在nz維量測和新息誤差模型相互獨立的條件下,正確量測落入矩形跟蹤門概率為式中:V G為跟蹤門的體積

正確量測位于跟蹤門外的概率為

當(dāng)所選PˉG較小時,有近似關(guān)系PˉG≈1-nz PˉG。

nz維矩形波門的面(體)積為

利用信息協(xié)方差矩陣的標(biāo)準(zhǔn)差進(jìn)行歸一化后得其面(體)積為

若不同分量對應(yīng)的波門常數(shù)K G互不相同,則式(9)和式(10)分別為

本文仿真時選用矩形跟蹤門+速度門輔助約束的方法對目標(biāo)量測進(jìn)行初選,以進(jìn)一步減少相關(guān)波門內(nèi)無關(guān)點跡的個數(shù)。

對一個被跟蹤的目標(biāo)來說,前后兩幀之間目標(biāo)的速度變化是連續(xù)的,不會產(chǎn)生大的突變。因此,可用一個速度門限將其限定在某個范圍之內(nèi)。

式中:a max表示目標(biāo)機動加速度可能的最大值,可根據(jù)實際情況選擇不同的值。仿真中a max=500 m/s2,矩形跟蹤門常數(shù)K G=2.6,相應(yīng)的門概率為P G=0.990 7。

(3)點跡起始/終止

本文在仿真時,采用基于順序處理技術(shù)的一種點跡起始算法,可稱之為直觀滑窗方法。

假設(shè)(ri,vi),i=1,2,…,N為N次連續(xù)掃描獲得的彈目距離—速度觀測值,如果這N次掃描中有某M對觀測值滿足以下約束規(guī)則,那么該方法就認(rèn)定應(yīng)起始一條點跡。否則,滑窗右移一幀繼續(xù)進(jìn)行判定。

由于彈目距離觀測值量化誤差Er的存在,相鄰兩幀的彈目距離觀測值之差為0或E r,即有

此時,目標(biāo)速度的最大值為vmax=Er/T,T為幀間隔時間。

由于目標(biāo)速度變化是連續(xù)的。因此,相鄰兩幀之間的目標(biāo)速度不會產(chǎn)生突變,即有

式中:amax表示目標(biāo)機動加速度可能的最大值,可根據(jù)實際情況選擇不同的值。

直觀滑窗方法計算簡單,在沒有真假目標(biāo)先驗信息的情況下,是較好的選擇。

本文在仿真時采用點跡終結(jié)判斷準(zhǔn)則是連續(xù)K次掃描丟失目標(biāo)決定點跡終結(jié)[7]。

2.2.2 干擾判決模塊

干擾判決模塊主要實現(xiàn)干擾信號的邏輯判斷,若干擾出現(xiàn)時伴隨著能量的急劇變化,導(dǎo)致能量較小的目標(biāo)被淹沒;若干擾出現(xiàn)時,能量沒有出現(xiàn)明顯變化,則同一時刻能夠檢測到目標(biāo)和干擾兩個信號。因此,可以根據(jù)AGC電平突變、及點跡個數(shù)與測量值個數(shù)的關(guān)系變化為判定依據(jù),來確定濾波判斷開始時刻。

從該時刻開始,根據(jù)每條點跡的濾波速度和測量速度是否匹配[8],作為點跡是目標(biāo)或干擾的標(biāo)志(速度-距離同步拖引干擾的判決方法將在后續(xù)工作中給出)。

3 實驗仿真與結(jié)果分析

3.1 多重距離拖引干擾

設(shè)置目標(biāo)初始距離25 km,速度900 m/s;每幀數(shù)據(jù)間隔15 m s,目標(biāo)在每幀之間是連續(xù)的,干擾在每幀之間存在缺失。干擾1拖速為600 m/s,即實際速度為300 m/s,拖引時間為2.4 s;干擾2拖速為450 m/s,即實際速度為450 m/s,拖引時間為2.4 s;干擾3拖速為300m/s,即實際速度為600m/s,拖引時間為2.4 s。圖4為三重距離拖引干擾點跡濾波。

圖4 三重距離拖引干擾點跡濾波圖

由圖可見,在干擾出現(xiàn)前,目標(biāo)建立出穩(wěn)定點跡。而在第68幀時刻同時出現(xiàn)三個干擾,在第88、104、116幀時刻,三個干擾先后被識別并跟蹤。經(jīng)過一段時間的點跡起始,三個干擾也形成穩(wěn)定點跡。圖5為目標(biāo)和干擾的速度濾波結(jié)果。

圖5 目標(biāo)和干擾速度濾波結(jié)果

從圖中可以得知,目標(biāo)的濾波速度與點跡量測速度基本相同,而三個干擾的濾波速度與點跡量測速速度都存在不同程度的分離。所以區(qū)別干擾,只需通過濾波速度和量測速度的差值對比來識別,進(jìn)而得到真實目標(biāo)信號角度、速度、距離的信息并進(jìn)行跟蹤。

3.2 多重速度拖引干擾

設(shè)置目標(biāo)初始距離25 km,速度600 m/s;同時施放三個速度拖引干擾,拖速分別為300 m/ss,150m/ss,225m/ss,間隔 5 s,拖5 s;每幀數(shù)據(jù)間隔10ms。圖6為三重速度拖引干擾示意圖。

圖6 三重速度拖引干擾示意圖

在第501幀時刻,同時出現(xiàn)三個干擾,由于初始階段三個干擾速度和目標(biāo)速度相互之間比較接近,關(guān)聯(lián)濾波器無法區(qū)分干擾和目標(biāo),只能通過對每幀跟蹤波門內(nèi)的距離測量值進(jìn)行概率加權(quán)獲得一個等效的距離測量值,并利用該等效的距離測量值對原目標(biāo)點跡進(jìn)行更新。隨著目標(biāo)速度和干擾速度之間的差值逐漸增大,濾波器在503幀時刻開始區(qū)分出目標(biāo)和干擾,并開始出現(xiàn)多條點跡。圖7為跟目標(biāo)和干擾點跡的濾波速度和點跡量測速速度之間的關(guān)系。

圖7 目標(biāo)和干擾速度濾波結(jié)果

雖然目標(biāo)濾波曲線有一定波動,但是仍然可以清楚看到目標(biāo)速度濾波曲線與速度量測曲線的差值在零附近波動;而速度拖引干擾的點跡濾波曲線與量測曲線存在較大差異。從上圖可以清楚的區(qū)分目標(biāo)與干擾,實現(xiàn)干擾剔除,進(jìn)而得到真實目標(biāo)信號角度、速度、距離的信息并進(jìn)行跟蹤。

4 結(jié)束語

本文通過信息關(guān)聯(lián)數(shù)據(jù)處理方法抗自衛(wèi)式欺騙干擾。將雷達(dá)航跡理論引入導(dǎo)引頭中,利用目標(biāo)回波信號多幀之間的關(guān)聯(lián)特性,通過點跡建立分離干擾與目標(biāo);再通過目標(biāo)與干擾在點跡上的一些特性差異加以判斷,從而識別目標(biāo)與干擾。

[1] 楊海林,郭愛芳,侯民勝.PD雷達(dá)抗速度欺騙干擾的計算機仿真[J].電子工程師,2005.

[2] 趙雪飛,趙川,劉崢.PD雷達(dá)抗距離-速度同步拖引干擾的頻譜識別法[J].雷達(dá)對抗,2005(3).

[3] 侯民勝.雷達(dá)抗距離欺騙干擾的計算機仿真研究[J].電子工程師,2003,(9).

[4] 曾茂生.脈沖多普勒雷達(dá)技術(shù)特征及干擾方法的探討[J].艦船電子對抗.2003,26(1):16-18.

[5] 高烽.雷達(dá)導(dǎo)引頭概論[M].電子工業(yè)出版社,2010.

[6] Samuel.S.Blackman.M ultip le Target Tracking w ith Radar App lications[M].USA:A rtech House,1986.

[7] 陳瓊,邱峰,王俊.一種實用的多目標(biāo)航跡跟蹤算法[J].火控雷達(dá)技術(shù).2008(1).

[8] Faruk Kural.A M ethod for Detecting RGPO/VGPO Jamm ing[J].Signal Processing and Communiaction App lications Conference,IEEE,2004.