基于迭代對(duì)消的外輻射源雷達(dá)目標(biāo)檢測方法

周子鉑, 王彬彬, 張朝偉, 劉建衛(wèi), 徐穎鑫, 王志會(huì)

(1． 空軍預(yù)警學(xué)院， 湖北武漢 430019; 2． 海軍工程大學(xué)， 湖北武漢 430033)

0 引言

外輻射源雷達(dá)(Passive Bistatic Radar，PBR)沒有發(fā)射機(jī)，也就不輻射電磁波，其通過接收空間已經(jīng)存在的電磁波實(shí)現(xiàn)目標(biāo)檢測。因?yàn)橥廨椛湓蠢走_(dá)工作在靜默狀態(tài)，所以敵方很難實(shí)施精確的定向電磁干擾與物理攻擊。鑒于其卓越的抗干擾、反隱身以及較低的制造成本等特點(diǎn)，近些年吸引了越來越多的研究興趣。外輻射源雷達(dá)的目標(biāo)檢測主要依靠相干積累來實(shí)現(xiàn)。相干積累基于雙通道實(shí)現(xiàn)，這意味著外輻射源雷達(dá)至少包含有兩個(gè)接收通道，參考通道用于接收直達(dá)波信號(hào)；觀測通道用于接收目標(biāo)回波，當(dāng)然也會(huì)摻雜著雜波干擾。由于輻射信號(hào)功率較低、目標(biāo)波動(dòng)等因素的影響，觀測通道內(nèi)的目標(biāo)回波能量遠(yuǎn)低于直達(dá)波和多徑雜波。因此，雜波消除是外輻射源雷達(dá)目標(biāo)檢測的關(guān)鍵步驟。此外，即使在雜波消除之后，目標(biāo)回波仍然遠(yuǎn)遠(yuǎn)比噪聲更弱，所以需要對(duì)觀測信號(hào)進(jìn)行相干積累以提升信號(hào)功率，最終實(shí)現(xiàn)目標(biāo)檢測。

目前廣泛使用的調(diào)頻廣播(Frequency Modulation, FM)、電視(Television，TV)、第四代移動(dòng)通信系統(tǒng)(Long Time Evolution, LTE)、全球定位導(dǎo)航系統(tǒng)(Global Navigation Satellite System, GNSS)以及各類商業(yè)衛(wèi)星使用的信號(hào)都可以作為外輻射源雷達(dá)的信號(hào)源。鑒于此類信號(hào)的輻射源安裝位置海拔比較高且信號(hào)覆蓋范圍廣等特點(diǎn)，外輻射源雷達(dá)便可以利用這些信號(hào)實(shí)現(xiàn)隱身低空目標(biāo)的探測，彌補(bǔ)主動(dòng)雷達(dá)存在探測盲區(qū)的不足。在眾多的機(jī)會(huì)輻射源中，GNSS具備建設(shè)完整、工作穩(wěn)定、全球覆蓋等特點(diǎn)，是一種非常適合用于目標(biāo)監(jiān)視的信號(hào)源。GNSS主要包括有美國的全球定位系統(tǒng)(Global Position System, GPS)、歐盟的伽利略系統(tǒng)(Galileo System)、俄羅斯的格洛納斯系統(tǒng)(GLONASS)以及中國的北斗導(dǎo)航系統(tǒng)(Beidou)，各方信號(hào)都可以被外輻射源雷達(dá)利用。因此基于GNSS信號(hào)的外輻射源雷達(dá)目標(biāo)檢測具備相當(dāng)?shù)目尚行浴Ｔ缭?993年Tsui等人就提出基于GNSS外輻射源的鏡面反射目標(biāo)模型，利用定向天線分別接收直達(dá)波和目標(biāo)回波，然后經(jīng)過信號(hào)處理檢測到目標(biāo)，并估計(jì)目標(biāo)的雙基地距離和多普勒頻率。

目前已有的外輻射源目標(biāo)檢測方法主要從時(shí)域、頻域和空域著手解決雜波對(duì)目標(biāo)檢測的干擾。Colone等在文獻(xiàn)[22]提出的擴(kuò)展相消算法(Extensive Cancellation Algorithm, ECA)利用參考天線接收直達(dá)波信號(hào)并構(gòu)造雜波子空間，基于目標(biāo)回波空間與雜波空間正交的特征，將觀測信號(hào)投影于雜波子空間的正交補(bǔ)完成雜波抑制，最后計(jì)算投影信號(hào)與直達(dá)波信號(hào)的互模糊函數(shù)實(shí)現(xiàn)相干積累獲得目標(biāo)檢測結(jié)果。鑒于GNSS信號(hào)互模糊函數(shù)多普勒頻帶較寬的特點(diǎn)，零多普勒頻率的固定雜波可能會(huì)掩蓋低多普勒頻率的微弱目標(biāo)，甚至強(qiáng)目標(biāo)也會(huì)掩蓋弱目標(biāo)，影響該算法的目標(biāo)檢測效果。并且ECA算法需要矩陣求逆，在數(shù)據(jù)比較長時(shí)會(huì)帶來很大的計(jì)算量，難以滿足實(shí)時(shí)信號(hào)處理的要求，尤其在較低信噪比(Signal-to-Noise-Ratio, SNR)時(shí)，該算法的性能會(huì)受到嚴(yán)重影響。Zhao等提出利用觀測信號(hào)的頻域互相關(guān)先檢測出功率比較強(qiáng)的雜波信號(hào)，然后利用估計(jì)得到的雜波參數(shù)重構(gòu)雜波信號(hào)并將其從觀測信號(hào)中減除，相比于時(shí)域處理方法而言，該算法可有效降低計(jì)算量，但是該算法需要主觀設(shè)置閾值以判定雜波，閾值選擇的準(zhǔn)確性嚴(yán)重影響到算法的魯棒性。Zhao等首先分析回波多普勒頻率響應(yīng)，根據(jù)多普勒頻率將信號(hào)分離，并在相應(yīng)多普勒頻率響應(yīng)的估計(jì)過程中消除零多普勒頻率信號(hào)實(shí)現(xiàn)雜波抑制，該算法對(duì)數(shù)字調(diào)頻廣播信號(hào)效果比較明顯，但是對(duì)其他信號(hào)體制不具備普適性。基于空域的雜波抑制方法主要通過降低天線副瓣的方法壓制雜波干擾，雖然這樣可以在信號(hào)處理之前從根本上減少雜波的干擾，但這顯然增加了外輻射源雷達(dá)的系統(tǒng)復(fù)雜度。

針對(duì)目標(biāo)反射GNSS回波十分微弱的特點(diǎn)，本文首先根據(jù)回波模型分析觀測信號(hào)的組成成分，然后根據(jù)雜波強(qiáng)烈而目標(biāo)回波微弱的特點(diǎn)，將信號(hào)的雜波成分從觀測信號(hào)中減去。最后利用脈沖壓縮實(shí)現(xiàn)觀測矩陣的相干積累，這樣可以在一定程度上避免長時(shí)間觀測過程中目標(biāo)運(yùn)動(dòng)狀態(tài)對(duì)檢測結(jié)果的影響。

1 信號(hào)模型

GNSS信號(hào)是將導(dǎo)航信息調(diào)制在周期產(chǎn)生的擴(kuò)頻序列碼上產(chǎn)生的，然后再混頻到系統(tǒng)所使用的頻段之上。其中使用的擴(kuò)頻序列碼信號(hào)是以1 ms為周期循環(huán)產(chǎn)生的，此處為簡化數(shù)學(xué)推導(dǎo)過程忽略調(diào)制在其上的導(dǎo)航信息，將GNSS信號(hào)建模為簡單的以1 ms為周期的擴(kuò)頻序列碼信號(hào)。

基于GNSS信號(hào)的外輻射源雷達(dá)目標(biāo)檢測模型如圖1所示，在地球軌道飛行的衛(wèi)星被視為輻射源，雷達(dá)接收機(jī)放置于參考坐標(biāo)系的原點(diǎn)處。衛(wèi)星和目標(biāo)的坐標(biāo)分別為(,,)和(,,)，而且目標(biāo)速度為(,,)。表示衛(wèi)星與雷達(dá)之間的直線距離,和則表示目標(biāo)回波的發(fā)射距離和接收距離。為雙基地角,表示目標(biāo)速度與雙基地角的平分線之間的夾角。雷達(dá)不僅能接收到衛(wèi)星廣播的直達(dá)波信號(hào)，還會(huì)接收到來自目標(biāo)反射的回波，以及地面反射的回波。此處設(shè)置兩個(gè)接收通道，參考通道負(fù)責(zé)接收直達(dá)波信號(hào)，觀測通道的波束則指向目標(biāo)接收觀測信號(hào)。假設(shè)衛(wèi)星運(yùn)動(dòng)已補(bǔ)償完成，并將信號(hào)建模為不含導(dǎo)航信息的擴(kuò)頻序列碼，用函數(shù)()表示，混頻之后的基帶觀測信號(hào)可建模為

(1)

式中，為載波波長，為光速，為地雜波的散射系數(shù)，為目標(biāo)的散射系數(shù)，為時(shí)間，和分別為雜波路徑數(shù)和目標(biāo)數(shù)，,分別為地雜波的發(fā)射距離和接收距離，()表示隨時(shí)間變化的高斯白噪聲。

圖1 GNSS外輻射源目標(biāo)檢測模型

本文為提高信號(hào)處理的效率，將雷達(dá)接收的連續(xù)波信號(hào)按周期=1 ms進(jìn)行分段截取，并將截取的小段回波信號(hào)作為整體進(jìn)行處理。在忽略導(dǎo)航信號(hào)的情況下，可以近似每一段截取的信號(hào)在衛(wèi)星發(fā)射的時(shí)候都是相同，接收回波的差異主要由目標(biāo)散射引起，因此可以根據(jù)接收回波進(jìn)行目標(biāo)特征提取，進(jìn)而檢測并識(shí)別目標(biāo)。

這樣可以就將連續(xù)波信號(hào)重新建模為脈寬為，脈沖重復(fù)頻率為1 kHz的脈沖信號(hào)。鑒于目標(biāo)相對(duì)較低的飛行速度和極高的衛(wèi)星高度(兩萬余公里)，可近似目標(biāo)在周期時(shí)間內(nèi)是靜止的，而在各周期之間目標(biāo)是運(yùn)動(dòng)的，式(1)所示連續(xù)波信號(hào)可建模為下列脈沖信號(hào)：

(2)

式中，,分別表示快時(shí)間和慢時(shí)間，且0<<，即目標(biāo)在快時(shí)間域被認(rèn)為是靜止的，而目標(biāo)的運(yùn)動(dòng)狀態(tài)僅與慢時(shí)間相關(guān)。經(jīng)上述處理，連續(xù)波信號(hào)就被分段截取成以為脈寬的觀測信號(hào)矩陣。那么，數(shù)字采樣之后的觀測信號(hào)模型可重寫如下：

(3)

同樣，數(shù)字采樣之后的參考通道信號(hào)如下式所示：

[(1)(2) …(-1)()]

(4)

式中，為參考信號(hào)幅度，為周期時(shí)間內(nèi)直達(dá)波信號(hào)的長度，與采樣率相關(guān)。

2 基于迭代對(duì)消的外輻射源雷達(dá)目標(biāo)檢測

2.1 ECA雜波抑制算法

作為一種經(jīng)典算法，ECA利用接收回波在雜波正交補(bǔ)空間的投影實(shí)現(xiàn)雜波抑制。經(jīng)上述分析可知，接收回波主要由目標(biāo)回波和雜波組成，并且目標(biāo)回波子空間與雜波子空間相互正交。通過將雷達(dá)接收回波投影于直達(dá)波及其時(shí)延信號(hào)張成的雜波子空間的正交補(bǔ)，可以實(shí)現(xiàn)接收回波的雜波去除。雜波子空間如下式所示：

(5)

式中，表示距離抑制的階次，[·]表示矩陣轉(zhuǎn)置。矩陣中的第一行到第列分別表示帶有0到-1次平移的多路徑樣本。

根據(jù)目標(biāo)回波子空間與雜波子空間的正交性，ECA算法雜波子空間系數(shù)的求解可以轉(zhuǎn)化為下列最優(yōu)化問題：

(6)

式中，為子空間系數(shù)。

詳細(xì)考察式(6)可知，該式是一個(gè)標(biāo)準(zhǔn)的二階凹函數(shù)。求解該準(zhǔn)則的共軛梯度，并令其等于零：

(7)

式中，[·]表示矩陣的共軛轉(zhuǎn)置。此時(shí)，可求得

=()

(8)

因此，經(jīng)過投影之后的剩余目標(biāo)回波信號(hào)如下式所示：

=-

(9)

2.2 基于迭代對(duì)消的雜波抑制

因?yàn)樾l(wèi)星較高的飛行高度，圖1所示的目標(biāo)檢測模型表明衛(wèi)星發(fā)射的信號(hào)需要經(jīng)歷很遠(yuǎn)的傳輸路徑才能照射到目標(biāo)，然后經(jīng)目標(biāo)散射被雷達(dá)接收到，一般認(rèn)為低于噪聲40～50 dB。一方面，由于地面復(fù)雜的地理環(huán)境，雷達(dá)接收到的地雜波不僅強(qiáng)烈，甚至可能包含很多傳輸路徑，即目標(biāo)散射系數(shù)遠(yuǎn)小于地雜波散射系數(shù)?。另一方面，觀測通道也會(huì)接收到來自衛(wèi)星發(fā)射的直達(dá)波，此部分為觀測信號(hào)中能量最強(qiáng)的分量。

經(jīng)上述分析，觀測信號(hào)能量主要被雜波占據(jù)，而目標(biāo)回波很微弱，即信號(hào)的主要成分為雜波(非目標(biāo)回波均為雜波)，目標(biāo)回波可近似歸于噪聲進(jìn)行考慮。那么，雷達(dá)回波可歸納如下：

(10)

式中:

(11)

重構(gòu)之后的觀測信號(hào)矩陣可用表示：

(12)

觀測矩陣的奇異值能有效反映觀測回波中各分量的權(quán)值，較大的奇異值反映回波中的主要分量，而其余較小的奇異值則反映噪聲的功率。如式(10)所示，觀測信號(hào)由雜波和噪聲組成，其中雜波功率強(qiáng)度遠(yuǎn)大于噪聲，可以近似回波的主要分量為雜波。因此，可對(duì)觀測信號(hào)矩陣進(jìn)行矩陣分解以獲得觀測矩陣的奇異值及相關(guān)組成分量，如式(13)所示：

(13)

若矩陣的秩為，那么，式中=diag(,,,…,)，而(=1,2,…,)為觀測矩陣的全部非零奇異值，且滿足≥≥…≥≥+1≈…≈。根據(jù)奇異值分解理論可知，酉矩陣由的特征向量(,=1,2,…,)構(gòu)成，而酉矩陣則由的特征向量(,=1,2,…,)構(gòu)成。那么，觀測矩陣可進(jìn)一步分解為

(14)

如前所述，觀測信號(hào)主要能量由雜波占據(jù)，且目標(biāo)回波功率遠(yuǎn)低于噪聲，因此觀測矩陣的前個(gè)大奇異值(,,,…,)表示觀測矩陣中的雜波分量，將回波矩陣大奇異值對(duì)應(yīng)的主成分消除掉就可以實(shí)現(xiàn)雜波消除。雖然已知目標(biāo)回波功率遠(yuǎn)低于雜波功率，但是多徑雜波數(shù)量和目標(biāo)特性等先驗(yàn)信息是未知的，從而無法在起始階段即確定上述大特征值個(gè)數(shù)。因此，為有效消除觀測矩陣中的雜波成分，在執(zhí)行雜波消除的過程中，采用式(15)中逐步迭代的方式。

(15)

(16)

圖2 觀測矩陣奇異值的比值

如前所述，基于外輻射源的目標(biāo)回波功率遠(yuǎn)低于噪聲功率，并且可能受剩余雜波的影響。因此，雜波抑制之后仍需對(duì)回波進(jìn)行相干積累以提升目標(biāo)回波功率，然后才能實(shí)現(xiàn)比較穩(wěn)定的目標(biāo)檢測。將觀測矩陣消除雜波之后，分別在距離維和多普勒維進(jìn)行二維脈沖壓縮實(shí)現(xiàn)觀測矩陣的相干積累，首先用匹配濾波進(jìn)行距離維的脈沖壓縮獲取目標(biāo)的距離信息：

((,),))

(17)

式中，(,)表示直達(dá)波信號(hào)。在獲得目標(biāo)距離信息之后，再進(jìn)行方位向壓縮獲取目標(biāo)多普勒維信息，如下所示：

(,)=((,),)

(18)

為獲得最佳雜波消除效果，此處利用信息熵來表征相干積累結(jié)果的無序或有序狀態(tài)，相干積累結(jié)果雜亂無章而無明顯目標(biāo)顯示的時(shí)候?yàn)闊o序狀態(tài)，此時(shí)積累結(jié)果的熵較大；若相干積累結(jié)果為有顯著目標(biāo)顯示的狀態(tài)即為有序狀態(tài)，且目標(biāo)峰值高于底噪越多，熵就越小。

()=-(,)ln(,)

(19)

Δ()=()-(-1)

(20)

其中(0)表示未進(jìn)行迭代，而直接對(duì)觀測矩陣進(jìn)行相干積累后獲得結(jié)果的熵。

綜上所述，算法執(zhí)行框圖如圖3所示，具體包括的步驟有：

1) 采集觀測矩陣和直達(dá)信號(hào)向量；

2) 計(jì)算觀測矩陣的奇異值、特征向量酉矩陣和，并逐步對(duì)觀測矩陣進(jìn)行迭代雜波對(duì)消；

3) 對(duì)雜波對(duì)消之后的觀測信號(hào)進(jìn)行相干積累獲得距離-多普勒(Range Doppler, RD)圖譜；

4) 計(jì)算迭代過程中相干積累結(jié)果熵的變化，熵值下降最劇烈的相應(yīng)迭代次數(shù)為最佳迭代次數(shù)，并且利用當(dāng)次迭代之后的相干積累結(jié)果進(jìn)行目標(biāo)檢測。

圖3 算法執(zhí)行框架

2.3 算法復(fù)雜度分析

表1 乘法計(jì)算量比較

(a) 不同脈沖采樣點(diǎn)數(shù)條件下的計(jì)算量

3 實(shí)驗(yàn)結(jié)果

此部分利用GPS L1波段仿真信號(hào)進(jìn)行實(shí)驗(yàn)，采樣頻率為2.046 MHz，信號(hào)截取長度=1 ms，即1ms的采樣點(diǎn)數(shù)為2 046，脈沖重復(fù)頻率為1 kHz，觀測時(shí)間為200 ms。仿真實(shí)驗(yàn)中，觀測信號(hào)包含一路直達(dá)波，五路多徑雜波以及三路微弱目標(biāo)回波，具體仿真參數(shù)如表2所示。

表2 相關(guān)仿真參數(shù)

3.1 仿真實(shí)驗(yàn)效果分析

為有效說明本文所述方法的有效性，實(shí)驗(yàn)中不僅對(duì)所提方法的目標(biāo)檢測效果進(jìn)行討論，并且與經(jīng)典的ECA方法效果進(jìn)行比較。

本實(shí)驗(yàn)中最佳迭代次數(shù)由迭代過程中相干積累結(jié)果熵值變化的最小值確定，圖5呈現(xiàn)了迭代過程中熵的變化。圖5(a)顯示迭代過程中熵值的變化，由圖可見，隨著迭代次數(shù)的增加，熵呈現(xiàn)先減后增的趨勢，在1～4次迭代之間，相干積累結(jié)果的熵逐漸下降，且在第3次和第4次迭代之間熵值減小最為顯著，那么，當(dāng)次迭代可獲得最佳目標(biāo)檢測結(jié)果；此外，觀察熵變的整個(gè)過程可見，由于此實(shí)驗(yàn)中包含3個(gè)目標(biāo)，在4～6次迭代之間，迭代過程中目標(biāo)回波也被逐漸對(duì)消，到第6次迭代之后，只剩下一個(gè)目標(biāo)回波，所以此時(shí)相干積累結(jié)果的熵最小；當(dāng)?shù)螖?shù)大于6時(shí)，3個(gè)目標(biāo)都被完全消除，相干積累無明顯峰值出現(xiàn)，所以熵值會(huì)逐漸增大到趨于穩(wěn)定，如圖所示。

因?yàn)榇颂幍鷮?duì)消的目的是消除雜波并實(shí)現(xiàn)目標(biāo)檢測，所以熵減的趨勢是關(guān)注的重點(diǎn)，且在這個(gè)過程中熵的變化率小于零。圖5(b)則反映了整個(gè)迭代對(duì)消過程中熵的變化率，可見，第4次迭代造成相干積累結(jié)果的熵減最顯著，與圖5(a)吻合。因此可以判定在本實(shí)驗(yàn)中迭代次數(shù)為4次即可獲得最佳目標(biāo)檢測結(jié)果，雜波抑制之后的相干積累結(jié)果如圖6所示。

(a) 相干積累結(jié)果的熵變

(a) 原始觀測矩陣的檢測結(jié)果

圖6給出了本文所提方法的目標(biāo)檢測效果，并與經(jīng)典的ECA方法進(jìn)行了比較。圖6(a)表示未進(jìn)行雜波抑制，直接利用觀測矩陣進(jìn)行相干積累的結(jié)果，從圖中可看到，由于強(qiáng)雜波的存在，相干積累結(jié)果的峰值出現(xiàn)在接近零多普勒頻率處，即雜波被積累得到很高的功率。因?yàn)槟繕?biāo)回波本身比較微弱，在相干積累增益相同的情況下，目標(biāo)自然會(huì)被強(qiáng)雜波掩蓋，從而無法實(shí)現(xiàn)目標(biāo)檢測。但是在利用本文所提方法消除雜波之后就可以有效實(shí)現(xiàn)目標(biāo)檢測，如圖6(b)所示。在雜波被消除之后，相干積累增強(qiáng)了目標(biāo)回波功率，目標(biāo)峰值即出現(xiàn)在相應(yīng)距離單元和多普勒單元處，而無明顯的雜波剩余峰值。從圖中讀取3個(gè)目標(biāo)的雙基地距離和多普勒頻率如表3所示。與表2中的參數(shù)相比，參數(shù)估計(jì)精度比較高，可見本文所提方法可以有效實(shí)現(xiàn)目標(biāo)檢測，并估計(jì)獲得比較精確的目標(biāo)雙站距離和多普勒頻率參數(shù)。

表3 不同方法的目標(biāo)檢測結(jié)果

為有效說明所提方法的效果，圖6還給出了經(jīng)典的ECA算法對(duì)雜波的抑制性能。圖6(c)為分別對(duì)200次脈沖進(jìn)行ECA雜波抑制，然后進(jìn)行相干積累的結(jié)果。圖中顯示了4個(gè)峰值，其中3個(gè)多普勒頻率遠(yuǎn)遠(yuǎn)偏離于零的峰值為目標(biāo)，在接近零頻處的峰值為強(qiáng)烈的雜波剩余，具體4個(gè)峰值的參數(shù)如表3所示。相比于表2，表3中不僅顯示了3個(gè)目標(biāo)的參數(shù)估計(jì)結(jié)果，還給出了ECA算法雜波剩余的相關(guān)參數(shù)。此處由于數(shù)據(jù)長度有限，導(dǎo)致相干積累結(jié)果的多普勒頻譜發(fā)生頻譜泄露，原本應(yīng)當(dāng)位于零頻處的雜波偏移到接近零頻的0.625 8 Hz處。由圖6(c)可見，剩余雜波的峰值強(qiáng)度顯然高于目標(biāo)峰值19 dB以上，因此該雜波剩余將會(huì)對(duì)后續(xù)的目標(biāo)檢測構(gòu)成嚴(yán)重影響。

圖6(d)為利用ECA算法對(duì)一次脈沖的觀測信號(hào)進(jìn)行雜波抑制之后的互模糊函數(shù)，可見模糊函數(shù)在零頻處存在一個(gè)很深的凹口，表明ECA算法可以實(shí)現(xiàn)對(duì)零頻雜波的抑制。然而GNSS信號(hào)互模糊函數(shù)的多普勒頻帶較寬，可見在零多普勒頻率附近仍然存在較嚴(yán)重的雜波剩余，甚至?xí)蜎]目標(biāo)回波。因此，圖6(d)中并無明顯目標(biāo)指示，即無法實(shí)現(xiàn)有效的目標(biāo)檢測。

根據(jù)表2可知，雜波平均比目標(biāo)回波強(qiáng)60 dB以上，再綜合圖6(c)和圖6(d)目標(biāo)檢測效果，ECA算法能起到一定的雜波衰減作用，但是其殘留的雜波剩余仍然比較強(qiáng)烈，在200 ms的相干積累之后，雜波強(qiáng)度依然遠(yuǎn)遠(yuǎn)強(qiáng)過目標(biāo)回波，導(dǎo)致目標(biāo)檢測效果受到嚴(yán)重影響。

為充分證明本方法對(duì)復(fù)雜環(huán)境的適應(yīng)性，在不同的信噪比和信雜比(Signal-Clutter-Ratio, SCR)條件下分別做50次蒙特卡洛仿真實(shí)驗(yàn)，計(jì)算不同的雜波噪聲環(huán)境下獲得最佳檢測結(jié)果需要迭代的次數(shù)的平均值，如圖7所示。圖7(a)所示為不同信噪比條件下需要的迭代次數(shù)，圖中橫軸為目標(biāo)2的信噪比，圖中標(biāo)記點(diǎn)處標(biāo)出的為表1所示實(shí)驗(yàn)條件所需的迭代次數(shù)。根據(jù)結(jié)果可知，在該設(shè)定的SNR范圍內(nèi)，平均迭代4～5次就可以獲得最佳檢測結(jié)果，因此可以認(rèn)為本方法在低SNR條件下仍然能獲得比較穩(wěn)健的目標(biāo)檢測效果。圖7(b)為不同信雜比(目標(biāo)2與直達(dá)波的功率比值，SCR)條件下的迭代次數(shù)，同樣，圖中標(biāo)記點(diǎn)處標(biāo)出的為表1所示實(shí)驗(yàn)條件下的所需迭代次數(shù)，與圖7(a)類似，在該寬泛的SCR范圍內(nèi)，4～5次即可獲得目標(biāo)的檢測結(jié)果。綜上，即使在較強(qiáng)的噪聲和雜波干擾條件下，本方法具備較強(qiáng)的魯棒性。

(a) 不同SNR條件所需迭代次數(shù)

3.2 算法計(jì)算量分析

表4所示為經(jīng)典ECA算法和本文提出的基于迭代對(duì)消的雜波抑制方法對(duì)200 ms觀測信號(hào)處理時(shí)間的比較，仿真實(shí)驗(yàn)執(zhí)行處理器為Intel Core(TM)i5-7200U CPU@2.5GHz，仿真軟件為MATLAB 2017a。

表4 算法執(zhí)行時(shí)間比較

綜上所示，本方法執(zhí)行時(shí)間遠(yuǎn)小于經(jīng)典的ECA算法，而且本方法對(duì)于雜波的衰減量遠(yuǎn)大于ECA算法，因此本文所提方法更適用于實(shí)時(shí)外輻射源目標(biāo)檢測。

4 結(jié)束語

針對(duì)GNSS信號(hào)用于機(jī)會(huì)輻射源的應(yīng)用，本文提出一種基于迭代對(duì)消的外輻射源目標(biāo)檢測算法，該算法針對(duì)目標(biāo)反射回波十分微弱的特征，對(duì)構(gòu)建的觀測矩陣進(jìn)行奇異值分解，將較大特征值對(duì)應(yīng)的雜波成分從中消去進(jìn)而抑制回波信號(hào)中的強(qiáng)烈雜波，然后再對(duì)消除雜波后的觀測矩陣進(jìn)行脈沖壓縮實(shí)現(xiàn)相干積累完成目標(biāo)檢測。在雜波消除過程中，通過搜索每一次迭代之后相干積累結(jié)果熵變的最小值判定雜波消除完成，以獲得最佳檢測結(jié)果。最后利用GPS系統(tǒng)L1波段的仿真信號(hào)進(jìn)行實(shí)驗(yàn)，證明了本方法可在較寬泛的信噪比和信雜比范圍內(nèi)獲得穩(wěn)定的目標(biāo)檢測效果。