崔向陽(yáng) 劉 穎 付學(xué)斌 田 歡 陳 亮 王愛(ài)榮
(1.西安電子工程研究所 西安 710100;2.西安現(xiàn)代控制技術(shù)研究所 西安 710065)
隨著電子戰(zhàn)技術(shù)的發(fā)展,雷達(dá)工作的電磁環(huán)境日益復(fù)雜,在檢測(cè)目標(biāo)的同時(shí),需要高效的對(duì)抗空間干擾。對(duì)有源干擾的感知和測(cè)量是雷達(dá)抗干擾的前提,干擾源定向技術(shù)通過(guò)對(duì)空間信號(hào)進(jìn)行采樣處理,估計(jì)出輻射源的方位、俯仰信息,為抗干擾提供先驗(yàn)信息。
對(duì)干擾源定向的常用方法主要有:
1)比幅測(cè)向法,即通過(guò)對(duì)偵查空間進(jìn)行波束的全掃描,當(dāng)波束掃描到干擾源時(shí),利用比幅方法進(jìn)行測(cè)角,這種方法對(duì)波束掃描及占用的系統(tǒng)資源都有要求;
2)相位干涉儀測(cè)向,利用鑒相方法來(lái)估計(jì)干擾源方向,這種方法需要基線長(zhǎng)度有一定要求;
3)空間譜估計(jì)及超分辨的方法,方法以MUISC,ESPRIT等算法為代表,具有多角度分辨和高精度的特點(diǎn),但是實(shí)現(xiàn)方式較為復(fù)雜。
與此同時(shí),大規(guī)模數(shù)字處理器件的快速發(fā)展,更大規(guī)模的數(shù)字陣列雷達(dá),越來(lái)越多地進(jìn)入到了工程研制階段,為干擾測(cè)向提供了更多的可能性。干擾入射方向的估計(jì),屬于波達(dá)方向估計(jì)的范疇。傳統(tǒng)的雷達(dá),由于系統(tǒng)的自由度不夠,能夠采用的波達(dá)方向估計(jì)技術(shù)并不多。在數(shù)字陣列中對(duì)于波達(dá)方向的估計(jì),如MUISC、ESPRIT等算法,這些算法雖然測(cè)角精度高,但是考慮到雷達(dá)工作的波束駐留時(shí)間很短,工程應(yīng)用中受限于運(yùn)算時(shí)間的要求難以實(shí)現(xiàn),而對(duì)于兩維數(shù)字陣列,其復(fù)雜度進(jìn)一步提高。
對(duì)于一均勻線陣,其陣元個(gè)數(shù)為,其陣元間距為,為載波信號(hào)的波長(zhǎng)。對(duì)于遠(yuǎn)場(chǎng)個(gè)窄帶入射信號(hào)的時(shí)間平穩(wěn)信號(hào),其如入射角為,其第個(gè)陣元的輸出為
(1)
其中,=1,2,…,-1;()為入射信號(hào)。
對(duì)比FFT的表達(dá)式為
(2)
其中,F(xiàn)FT后第個(gè)譜峰對(duì)應(yīng)的空間采樣點(diǎn)為;_FFT為FFT的點(diǎn)數(shù),通過(guò)對(duì)應(yīng)與FFT的關(guān)系,獲得空間采樣頻率,第譜峰對(duì)應(yīng)的入射角度為
(3)
對(duì)于FFT估計(jì)的角度分辨率,其瑞利限(均勻線陣)由3 dB波束寬度決定為
(4)
兩維數(shù)字陣列其,分別為行、列陣元間距,行列陣元數(shù)分別為,,其全陣面接收的第個(gè)采樣快拍的數(shù)據(jù)矩陣為
(5)
本文中所提到的干擾測(cè)向方法,首先需要進(jìn)行接收數(shù)據(jù)的選取。需要對(duì)兩維子陣平面中,選取一個(gè)“十字架”,見(jiàn)圖1所示,即中間的一行及一列,進(jìn)行干擾測(cè)向。對(duì)于單個(gè)陣元,陣元的輻射的特性可近似認(rèn)為為全向天線。這樣選取某一列的所有行,見(jiàn)式(6)所示。
(6)
則在方位維可近似認(rèn)為全向波束,俯仰維為一維線陣,可以進(jìn)行俯仰維干擾角度測(cè)向。同樣的,選取某一列的所有行,即
()=(0,()1,()…-1,())
(7)
則在俯仰維可近似認(rèn)為全向波束,方位維為一維線陣,可以進(jìn)行方位維干擾角度測(cè)向。
圖1 陣面陣元選取平臺(tái)示意圖
其對(duì)應(yīng)的方位角和俯仰角分別為式(8)和式(9)所述,此角度為陣面坐標(biāo)系下的錐掃描角,n_FFT為FFT點(diǎn)數(shù)。
其中方位角為
(8)
其中俯仰角為
(9)
為了提高估計(jì)的可信度,干擾估計(jì)需要對(duì)多次快拍進(jìn)行處理,假設(shè)共采樣個(gè)快拍,(),=0,1,…,,為天線陣面“列”采樣點(diǎn);(),=0,1…,,為天線陣面“行”采樣點(diǎn)。通常情況下,由于快拍采樣過(guò)程極短,入射干擾角度幾乎沒(méi)有變化,那么多次快拍之間干擾入射角度應(yīng)該一致,即個(gè)干擾的入射角度為式(10)所示。
=[(0)(1) …()]
(10)
此外則測(cè)得的FFT后的峰值位置也應(yīng)一致。個(gè)快拍個(gè)干擾入射對(duì)應(yīng)的譜峰位置的信號(hào)如式(11)所示。
(11)
計(jì)算獲得方位維個(gè)干擾方向的次采樣快拍的記錄值為式(12)所示。
=[(0)(1) …()]
(12)
對(duì)應(yīng)的譜峰位置的信號(hào)為式(13)所示。
(13)
對(duì)俯仰維譜峰信號(hào)的矩陣×,及方位維譜峰信號(hào)的矩陣×,進(jìn)行自相關(guān)運(yùn)算。獲得相關(guān)矩陣為式(14)所示。
(14)
構(gòu)造與相關(guān)矩陣對(duì)應(yīng)的角度匹配矩陣,該矩陣并不計(jì)算具體數(shù)值,而是通過(guò)對(duì)矩陣1×,1×進(jìn)行與×矩陣相同的矩陣運(yùn)算過(guò)程,生成角度查找表。如式(15)所示。
(15)
即(,)該點(diǎn)的方位角為(),俯仰角為()。
相關(guān)矩陣×幅度越大則表明相關(guān)性越大,即就是俯仰與方位角度配對(duì)程度越高。針對(duì)個(gè)干擾通過(guò)對(duì)相關(guān)矩陣×幅度排序后,獲得個(gè)極值。通過(guò)個(gè)極值在×矩陣中的獲得該極值數(shù)據(jù)的下標(biāo)(,),對(duì)比相關(guān)矩陣×與查找表×的對(duì)應(yīng)關(guān)系,如圖2所示,從矩陣×中,根據(jù)坐標(biāo)(,),并找到個(gè)配對(duì)的干擾角度(,),完成干擾角度配對(duì)。
圖2 角度配對(duì)示意圖
首先,設(shè)計(jì)雷達(dá)系統(tǒng)工作時(shí)序,增加靜默偵查區(qū),如圖3所示。干擾測(cè)向是為干擾對(duì)消提供角度先驗(yàn)知識(shí),所以在干擾測(cè)量時(shí),如果有己方雷達(dá)的輻射信號(hào)進(jìn)入干擾測(cè)量的通道,在對(duì)消時(shí)就會(huì)消去己方有用信號(hào)。所以需要在雷達(dá)發(fā)射周期中增加靜默區(qū),此時(shí)不輻射電磁信號(hào),僅接收信號(hào)(干擾)。此時(shí)雷達(dá)靜默工作,雷達(dá)陣面處于不發(fā)射僅接收狀態(tài)。
圖3 雷達(dá)工作時(shí)序示意圖
其次,在確定周圍環(huán)境無(wú)干擾時(shí)(如天線暗室環(huán)境),對(duì)行(或列)采樣的快拍進(jìn)行FFT,記錄雷達(dá)正常工作時(shí),統(tǒng)計(jì)雷達(dá)靜默區(qū)的接收機(jī)噪聲電平,并根據(jù)系統(tǒng)需求,當(dāng)正常工作時(shí)靜默區(qū)接收機(jī)噪聲電平大于某一值時(shí)認(rèn)為系統(tǒng)受到干擾,其中值可由上位機(jī)參數(shù)控制或設(shè)置經(jīng)驗(yàn)值。正常工作時(shí),如雷達(dá)靜默區(qū)偵聽(tīng)到的信號(hào)高于門(mén)限一定數(shù)量時(shí),判定為系統(tǒng)受到干擾。假設(shè)干擾的個(gè)數(shù)為個(gè),通常設(shè)置為系統(tǒng)能夠處理的最大干擾個(gè)數(shù),這取決于系統(tǒng)能夠?qū)垢蓴_的自由度和運(yùn)算能力。
再次,按照第一,二章所述的方法。利用FFT方法進(jìn)行干擾測(cè)向,充分利用數(shù)字陣列靈活的優(yōu)勢(shì),單獨(dú)處理“行”,“列”單次采樣快拍進(jìn)行FFT運(yùn)算,可以獲得方位、俯仰維指向干擾角度的FFT峰值點(diǎn),記錄其幅度是否大于門(mén)限值。
最后,采用第三章節(jié)所述方法,對(duì)測(cè)出的兩維角度進(jìn)行角度匹配,獲得最終的處理結(jié)果。
其處理流程見(jiàn)圖4所示,采用本文所述方法,在某研制雷達(dá)中取得了較好的測(cè)量結(jié)果。
圖4 處理流程示意圖
這里通過(guò)仿真設(shè)置兩個(gè)干擾源,其入射方向分別為:干擾1入射方位俯仰角為[-5,2],干擾2入射方位俯仰角為[15,-15]。首先通過(guò)干擾估計(jì)可以獲得方位俯仰的分別指向見(jiàn)圖5所示,之后經(jīng)過(guò)角度匹配的處理之后,獲得其真實(shí)指向,分別為[-4.99,2.08],[15.02,-15.06]。
圖5 多次采樣快排進(jìn)行FFT角度估計(jì)
實(shí)驗(yàn)中,設(shè)置兩個(gè)干擾源,且分別使用點(diǎn)頻連續(xù)波和噪聲干擾兩種形式的干擾。由于實(shí)現(xiàn)條件限制,干擾源位置無(wú)法靈活變動(dòng),這里通過(guò)轉(zhuǎn)動(dòng)天線伺服,測(cè)量獲得與天線法線偏離的角度后,通過(guò)讀取伺服轉(zhuǎn)動(dòng)角度值并加到天線法線偏移值上,來(lái)獲得測(cè)量真值。這里俯仰測(cè)試0°和5°,方位測(cè)試-5°和5°,即天線法線與伺服碼盤(pán)0°對(duì)齊時(shí),兩個(gè)干擾角度為[-5,0],[5,5]。之后伺服方位每5°轉(zhuǎn)動(dòng)測(cè)試一次,具體測(cè)試結(jié)果見(jiàn)表1所示。從表1能夠看到,對(duì)噪聲的干擾角度估計(jì)誤差大于對(duì)點(diǎn)頻的估計(jì)。噪聲類型干擾角度的估計(jì)誤差大約在0.3°,對(duì)點(diǎn)頻類型干擾角度的估計(jì)誤差小于0.1°。
表1 實(shí)際測(cè)試結(jié)果
本文通過(guò)增加雷達(dá)發(fā)射靜默周期,偵聽(tīng)環(huán)境干擾,同時(shí)通過(guò)對(duì)接收的數(shù)據(jù)進(jìn)行分析,根據(jù)設(shè)置的門(mén)限閾值,判定是否有干擾進(jìn)入,同時(shí)通過(guò)FFT的方法快速形成方向圖,進(jìn)行空域?yàn)V波搜索,獲得方位、俯仰維度的干擾指向角度,之后通過(guò)對(duì)干擾入射方向的數(shù)據(jù)進(jìn)行相關(guān)估計(jì),完成兩個(gè)維度的干擾角度配對(duì)處理,獲得配對(duì)后的干擾入射角的方位、俯仰兩維角度。具有相應(yīng)速度快的特點(diǎn),在對(duì)測(cè)向精度要求不高的場(chǎng)景,可以快速獲得干擾的兩維指向。