復(fù)雜電磁環(huán)境下微弱目標探測技術(shù)研究?

(中國電子科技集團公司第三十八研究所,安徽合肥230088)

0 引言

外輻射源雷達本身并不發(fā)射信號,主要利用非合作的廣播電臺信號[1-3]、電視信號[4]、衛(wèi)星信號及其他輻射源作為目標照射源,不容易被敵方所偵察和干擾,提高了雷達的隱蔽性。同時雷達照射源大多是低頻段的信號,增加了探測隱身飛機的能力。外輻射源雷達因具有反干擾強、抗反輻射導(dǎo)彈、反隱身及成本低等特點,逐漸受到重視,國內(nèi)外高校及研究機構(gòu)也進行了大量的研究和試驗[1-5]。

雷達系統(tǒng)利用民用電臺發(fā)射的直達波信號與目標反射的回波信號進行多普勒相關(guān)處理來進行目標的檢測和定位。由于系統(tǒng)必須接收來自電臺的直達波信號,系統(tǒng)與電臺間不能有遮擋,那么天線在接收目標回波的同時,必然會受到直達波和經(jīng)高山及建筑物等反射的多徑信號所干擾。雖然天線的法線能避開電臺方向,但是通過反射和繞射所接收到的直達波信號仍遠遠大于目標所反射的回波信號;同時由于雷達站周圍的環(huán)境可能比較復(fù)雜,來自近距離不同固定物所反射的多徑信號同樣也強于回波信號;一般情況下,國家對電臺的分布是有一定限制的,不會出現(xiàn)電臺相互干擾的問題,但是高靈敏度的接收機經(jīng)常能接收到較遠地方發(fā)射的電臺信號,在海邊和國境附近也會接收到對方的電臺干擾信號。由于電臺干擾信號遠遠強于目標所反射的回波信號,若不采取有效措施則限制了雷達系統(tǒng)的探測性能。

1 復(fù)雜電磁環(huán)境下的干擾信號分析

外輻射源雷達是一種特殊體制下的雙/多基地非合作雷達系統(tǒng),其發(fā)射站為位置明確的民用廣播、電視信號和自主廣播等發(fā)射臺。發(fā)射塔全向輻射的電磁波信號連續(xù)不間斷地照射廣闊的空域,經(jīng)電磁波信號照射的空中目標產(chǎn)生相應(yīng)的反射回波信號。雷達系統(tǒng)的接收站采用兩個通道,其中回波接收通道采用陣列天線形式,用于接收目標所反射的回波信號;參考接收通道采用單天線形式,指向發(fā)射臺方向,用于接收直達波參考信號。其定位原理圖如圖1所示。

圖1 外輻射源雷達的定位圖

其中發(fā)射站的位置為T,雷達接收站的位置為R,目標的位置為G,收發(fā)站之間的距離為L,目標和接收站之間的距離為Rr,目標和發(fā)射站之間的距離為Rt,在目標探測中通過回波接收通道采集的信號與參考接收通道采集的信號進行長時相干積累處理,得到兩者之間的時差估計τ,即

式中,S=Rr+Rt為距離和,c為光速。根據(jù)幾何關(guān)系可得到Rr為

式中,φr為通過多波束比幅技術(shù)所測得的方位角。

從上面的外輻射源定位系統(tǒng)中可以看出,雷達所處的電磁環(huán)境較為復(fù)雜,經(jīng)目標反射的回波信號完全淹沒在很強的雜波干擾中,特別是弱目標和處于探測遠區(qū)的目標回波。圖2顯示了在復(fù)雜電磁環(huán)境下回波的頻譜特性,在頻譜圖中不僅含有直達波和多徑信號,而且同頻和鄰頻電臺信號的干擾較為嚴重。在這種環(huán)境下進行目標檢測,根本不能發(fā)現(xiàn)微弱目標,嚴重時系統(tǒng)都無法檢測到目標。因此必須利用數(shù)字陣列雷達的優(yōu)勢和信號濾波處理算法,從空域[6]和時域[7-10]等方面采取綜合抗干擾處理方法,實現(xiàn)干擾雜波的抑制,同時采用長時相干積累處理提高信號的增益,特別是微弱目標的檢測,實現(xiàn)在復(fù)雜電磁環(huán)境下的連續(xù)跟蹤探測。

圖2 復(fù)雜電磁環(huán)境下回波的頻譜特性

2 空時自適應(yīng)處理的干擾抑制技術(shù)

2.1 基于空域的自適應(yīng)干擾置零處理

其原理是根據(jù)不同的最優(yōu)化準則,建立相應(yīng)的數(shù)學(xué)模型,通過調(diào)整各陣元的加權(quán)系數(shù)在雜波干擾處形成零點,得到干擾抑制的目的。針對無源雷達采用最大信干噪比(MSINR)準則,使陣列輸出信號與干擾加噪聲之比最大。其陣列接收信號為

式中,a(θ0)為信號來波方向的導(dǎo)向矢量,s0(t)為信號的復(fù)包絡(luò),Xi+n(t)為干擾加噪聲向量,與信號不相關(guān)。信號協(xié)方差矩陣Rs和干擾噪聲協(xié)方差矩陣Ri+n分別為

最大信干噪比準則是使輸出信干噪比最大,即

式中,W為權(quán)值,用拉格朗日乘子法,其目標函數(shù)為

對式(7)求導(dǎo)可得

式中,λ=SINR,最優(yōu)權(quán)矢量Wopt是與最大特征值λ對應(yīng)的特征向量,將式(4)代入式(8),解得

根據(jù)空域自適應(yīng)置零算法的分析,利用雷達在復(fù)雜電磁環(huán)境下采集的實際回波數(shù)據(jù),從圖2的回波頻譜中可以看出,回波中既有直達波和多徑信號的干擾,同時存在相鄰電臺的干擾。為此采用空域自適應(yīng)干擾抑制處理,其結(jié)果如圖3所示。

圖3 空域干擾抑制前后的處理結(jié)果

2.2 基于時域的改進型RLS自適應(yīng)干擾抑制處理

外輻射源雷達主要是利用回波信號與參考信號的長時相干積累處理來探測目標,由于回波信號很弱,而回波通道接收到的信號含有強直達波、多徑信號以及電臺干擾信號,嚴重影響了目標信號的檢測,抑制干擾信號尤為重要,直接影響著雷達系統(tǒng)的性能。為此采用自適應(yīng)信號處理的方法,通過實時更新自適應(yīng)濾波的系數(shù)得到抑制干擾信號的目的,而自適應(yīng)處理算法有很多種,比如最小均方(LMS)誤差算法[9-10]、遞推最小二乘(RLS)算法[11-12]、卡爾曼濾波(Kalman)算法、盲自適應(yīng)算法等,其中RLS算法的處理效果及工程實現(xiàn)得到了很好的應(yīng)用。

自適應(yīng)RLS算法在第k個采樣點對應(yīng)的一組輸入信號為xk(i),i=0,1,…,N,那么該算法對應(yīng)的輸出為

式中,N為濾波器的階數(shù),ek為誤差項,w(i)為濾波器的第i項系數(shù)。

式中,Pk為計算的遞歸項,Gk=。

在自適應(yīng)RLS算法中,遺忘因子λ同LMS誤差算法中步長μ的作用相似,在常規(guī)處理算法中根據(jù)經(jīng)驗值來選擇λ,不能實時調(diào)整來滿足系統(tǒng)要求,因此通過引入遺忘因子函數(shù)來解決這一問題。根據(jù)λ與誤差的關(guān)系來實時調(diào)整λ值,當(dāng)輸出誤差e(n)較大時,應(yīng)調(diào)整λ值使其較小而加快收斂速度;當(dāng)輸出誤差e(n)較小時,應(yīng)調(diào)整λ值使其較大而獲得較高的收斂精度;但是若只利用誤差信號e(n)來調(diào)整遺忘因子λ(n),可能會出現(xiàn)誤差信號使濾波系數(shù)偏離理想值。因此還需要調(diào)整λ,在系統(tǒng)初始狀態(tài)下使λ值較小而加快算法的收斂速度;當(dāng)算法接近收斂時使λ值較大而獲得較好的收斂精度。因此改變這一設(shè)計函數(shù),得到遺忘因子函數(shù)為

根據(jù)式(13)對遺忘因子函數(shù)的修改,對常規(guī)的RLS濾波算法進行修正,得到的改進型RLS算法計算流程如下所示:

利用在實驗系統(tǒng)中采集的實際數(shù)據(jù)進行仿真分析,其中圖4(a)表示按常規(guī)LMS處理算法對干擾抑制的處理結(jié)果,圖4(b)表示按改進型RLS處理算法對干擾抑制的處理結(jié)果。從兩個圖的比較可以看出,改進型RLS處理算法相對于LMS處理算法在干擾抑制上提高了10 dB以上,同時在收斂速度和穩(wěn)態(tài)精度上也都得到一定的提高。改進型RLS處理算法采用閉環(huán)自適應(yīng)處理算法,根據(jù)干擾的強度來快速調(diào)整濾波器的權(quán)系數(shù),抑制回波中的干擾信號。對于不同環(huán)境下的直達波和多徑信號等干擾信號都具有較好的抑制作用,具有較好的穩(wěn)健性。在實際數(shù)據(jù)的仿真中發(fā)現(xiàn)只有干擾強度發(fā)生變化時干擾抑制比變化較大,而與環(huán)境和頻率的變化影響不大。

圖4 LMS和改進RLS時域干擾抑制的處理結(jié)果比較

3 微弱目標的長時相干積累技術(shù)

外輻射源雷達在對低信噪比的微弱目標進行搜索檢測時,接收到目標的回波信噪比一般都很低,雖然采取空域和時域等多種干擾抑制措施,對這些低信噪比的目標也很難檢測到,因此對目標回波的能量進行長時有效積累,增加有用信號的能量,提高檢測目標的信噪比,從而保證雷達的探測性能。長時相干積累處理主要是利用陣列天線接收到的回波信號與輔助天線接收到的參考信號進行相關(guān)運算。由于該算法運算量大,工程上不易實現(xiàn),因此采用二維處理方法獲得目標的距離和多普勒信息,其處理方法如圖5所示。

圖5 距離和多普勒的二維處理技術(shù)

雷達接收到的參考信號相當(dāng)于外輻射源所發(fā)射的信號,若發(fā)射信號(即參考信號)為

根據(jù)二維處理技術(shù),參考信號分成了n部分,則每部分的參考信號為

同時相應(yīng)的回波信號為

因此根據(jù)參考信號與目標回波信號的積累處理,將參考信號延遲Δt后與接收到的目標回波信號共軛相乘,則輸出信號為

對式(21)中n部分的輸出結(jié)果,在實現(xiàn)頻域的FFT濾波處理[14],可搜索到在不同距離信息Δt和多普勒信息fd時與目標回波信號完全達到匹配,目標回波信號的能量得到了有效積累,實現(xiàn)了目標的檢測。

4 實驗仿真和結(jié)果分析

在該雷達實驗系統(tǒng)中,以電臺發(fā)射的廣播信號為輻射源,采取陣列天線接收目標所反射的回波信號,輔助天線指向電臺發(fā)射信號,其中廣播信號的頻率為97.6 MHz,周邊環(huán)境干擾比較強,如圖2所示。通過該實驗采集的數(shù)據(jù)進行仿真分析,其中常規(guī)的LMS處理和積累技術(shù)得到的結(jié)果如圖6所示,其中圖6(a)為目標的距離維顯示結(jié)果,圖6(b)為目標的速度維顯示結(jié)果。從圖中可知,只能發(fā)現(xiàn)兩個較大的目標信息,其中目標1位于41.3 km處,徑向速度為-154.7 m/s;目標2位于63 km處,徑向速度為-200.6 m/s。而采取上面分析的針對復(fù)雜電磁環(huán)境下采取的處理方法得到結(jié)果如圖7所示,其中圖7(a)為目標的距離維顯示結(jié)果,圖7(b)為目標的速度維顯示結(jié)果。不僅發(fā)現(xiàn)了兩個近距離的大目標信息,而且還發(fā)現(xiàn)了遠距離處的另外兩個微弱目標信息,其中目標3位于123.8 km處,徑向速度為193.8 m/s;目標4位于225.8 km處,徑向速度為-20 m/s,同時大目標的信雜比也相應(yīng)提高了近5 dB左右。因此在復(fù)雜電磁環(huán)境下,特別是多徑信號和強電臺信號的干擾,針對微弱目標的探測方法,大大提高目標的檢測概率,實現(xiàn)了目標的連續(xù)跟蹤和定位。

圖6 復(fù)雜電磁環(huán)境下常規(guī)處理檢測結(jié)果

圖7 復(fù)雜電磁環(huán)境下微弱目標檢測結(jié)果

5 結(jié)束語

基于外輻射源的無源雷達系統(tǒng)中,分析了在復(fù)雜電磁環(huán)境下的直達波、多徑信號和其他電臺信號的干擾,嚴重影響了無源雷達系統(tǒng)的探測性能,特別是微弱目標的檢測和定位。為此,提出了從空域自適應(yīng)處理技術(shù)和時域自適應(yīng)處理技術(shù)來抑制強干擾信號,利用實測數(shù)據(jù)進行處理和對比分析,驗證了算法的有效性,同時分析了低信噪比下目標回波的長時相干積累技術(shù),降低了實時處理的運算量。最后在實驗中進行了實際數(shù)據(jù)的仿真和應(yīng)用,針對在復(fù)雜電磁環(huán)境下的弱目標檢測技術(shù)提高了大目標的信噪比,而且探測到遠距離的微弱目標,在實際系統(tǒng)的應(yīng)用中取得了很好的效果。

[1]OLSEN K E,BAKER C J.FM-Based Passive Bistatic Radar as a Function of Available Bandwidth[C]∥2008 IEEE Radar Conference,Rome:IEEE,2008:1266-1273.

[2]王俊,保錚,張守宏.無源探測與跟蹤雷達系統(tǒng)技術(shù)及其發(fā)展[J].雷達科學(xué)與技術(shù),2004,2(3):129-135.WANG Jun,BAO Zheng,ZHANG Shou-hong.Passive Detecting/Tracking Radar System Technologies and Its Development[J].Radar Science and Technology,2004,2(3):129-135.(in Chinese)

[3]ZHENG Heng,ZHAO Hong-li,LI Fei.Accidental Radio Jamming Suppression in Passive Radar[C]∥2008 IEEE Radar Conference,Rome:IEEE,2008:294-298.

[4]萬顯榮.基于低頻段數(shù)字廣播電視信號的外輻射源雷達發(fā)展現(xiàn)狀與趨勢[J].雷達學(xué)報,2012,1(2):109-123.

[5]汪玲,卜曉楠,龐鑠,等.多徑環(huán)境下的無源分布式孔徑雷達成像方法[J].系統(tǒng)工程與電子技術(shù),2013,35(1):61-67.

[6]PARIHAR M S,BESU A,KOUL S K.Efficient Spurious Rejection and Null Steering Using Slot Antennas[J].IEEE Antennas and Wireless Propagation Letters,2011,10(1):207-210.

[7]李飛,關(guān)堂新,趙洪立,等.基于DTV廣播的外輻射源雷達對消算法研究[J].雷達科學(xué)與技術(shù),2011,9(5):397-400.LI Fei,GUAN Tang-xin,ZHAO Hong-li,et al.Research on Clutter Canceling Algorithm in Passive Radar Based on DTV Broadcast[J].Radar Science and Technology,2011,9(5):397-400.(in Chinese)

[8]KIM H W,PARK H S,LEE S K,et al.Modified-Filtered-u LMS Algorithm for Active Noise Control and Its Application to a Short Acoustic Duct[J].Mechanical Systems and Signal Processing,2011,25(1):475-484.

[9]許德剛,朱子平,洪一.自適應(yīng)濾波在無源探測中對雜波抑制的應(yīng)用[J].系統(tǒng)工程與電子技術(shù),2006,28(2):202-204.

[10]田福慶,羅榮,李克玉,等.基于改進的雙曲正切函數(shù)變步長LMS算法[J].系統(tǒng)工程與電子技術(shù),2012,34(9):1758-1763.

[11]ZAKHAROV Y V,WHITE G P,LIU J.Low-Complexity RLS Algorithms Using Dichotomous Coordinate Descent Iterations[J].IEEE Trans on Signal Processing,2008,56(7):3150-3161.

[12]劉梅,金龍.基于自適應(yīng)RLS算法的系統(tǒng)辨識應(yīng)用[J].中山大學(xué)研究生學(xué)刊(自然科學(xué)、醫(yī)學(xué)版),2013,34(1):80-89.