雙基地高頻地波雷達(dá)海雜波抑制

仇永斌，張 寧，張樹春

(1.哈爾濱工業(yè)大學(xué)電子與信息工程學(xué)院，150001哈爾濱，chouyongbin@126.com;2.空軍第一飛行學(xué)院飛行仿真技術(shù)研究所，150001哈爾濱)

無論單基地還是雙基地高頻地波雷達(dá)，海雜波在多普勒譜上都表現(xiàn)為兩個明顯幅度占優(yōu)的譜峰，即眾所周知的Bragg峰以及Bragg峰之間的連續(xù)區(qū)，被分別定義為一階海雜波和二階海雜波.其產(chǎn)生機(jī)理和特性在諸多文獻(xiàn)[1－4]中都有詳細(xì)分析.強(qiáng)功率的海雜波會淹沒附近的目標(biāo)信息，在多普勒譜上形成檢測盲區(qū).因而針對海雜波，尤其是一階海雜波，開發(fā)出很多抑制算法.經(jīng)過實(shí)測數(shù)據(jù)檢驗(yàn)有效的算法如奇異值分解法[5－6]、線性預(yù)測法[7]、相鄰雜波單元雜波對消[8]、雜波預(yù)白化法[9]、Root循環(huán)對消法[10]等，甚至 Music 子空間算法[11]也能在一定程度上抑制海雜波.這些算法是針對單基地海雜波提出的.單基地海雜波只與雷達(dá)工作頻率相關(guān)，對于平穩(wěn)的完全發(fā)展的海，一階Bragg峰多普勒頻率是緩時變的，在洋流作用下，與理論值會發(fā)生偏移，但這種時變特性和偏移特性都不影響把海雜波當(dāng)作單頻信號來模型化.而對于雙基地海雜波，不但與雷達(dá)工作頻率有關(guān)，還同雜波單元的雙基地角有關(guān)，在距離和趨近于基線長度時，海雜波頻率具有向零頻收斂的特性，同時在雙基地角和洋流的作用下，具有明顯的展寬效應(yīng)，在多普勒(RD)譜上，一階海雜波譜表現(xiàn)為兩條收斂的且逐漸展寬的曲線，因而不同于單基地的直線譜形態(tài).這種特性，直接排除了基于線性預(yù)測和相鄰雜波單元對消等雜波抑制方法應(yīng)用的可行性.同時，對雙基地而言，雙基地角很難精確估計(jì)，這使得雙基地海雜波抑制變得難以實(shí)現(xiàn).Music子空間算法把海雜波當(dāng)作無方向性的類噪聲信號進(jìn)行處理，可以實(shí)現(xiàn)海雜波抑制，也可以用于提取海面洋流信息[12－13].奇異值分解法和Root循環(huán)對消算法，單雙基地都可以應(yīng)用，但需要事先檢測出海雜波，才能有針對性地抑制雜波.但對于批量處理數(shù)據(jù)來說，先檢測后抑制顯然意義不大.對于雙基地布局，洋流與地貌狀況直接影響雜波譜形態(tài)，導(dǎo)致譜峰分裂以及理論值與預(yù)測值發(fā)生偏差，使得海雜波背景下目標(biāo)檢測變得異常困難，一個直接的影響是雜波會被當(dāng)作目標(biāo)被檢測和跟蹤，形成大量的虛假目標(biāo)和虛假航跡.針對這一問題，本文根據(jù)海雜波譜的先驗(yàn)知識，結(jié)合海面的平穩(wěn)特性，在RD譜上事先確定海雜波頻率范圍，然后利用Clean算法，在波束形成之前，利用天線陣元的時域數(shù)據(jù)，把各通道的海雜波抑制掉.這種方法盡管損失了一些出現(xiàn)在雜波頻率區(qū)域的目標(biāo)信息，但能夠?qū)崟r處理數(shù)據(jù)，為后續(xù)正常的目標(biāo)自動檢測提供便利.

1 雙基地海雜波先驗(yàn)知識

1.1 雙基地海雜波譜的影響因素

雙基地一階Bragg峰多普勒頻移的理論值為

式中:g為重力加速度;λ為雷達(dá)工作波長;φ為雙基地角;fMB為同頻率單基地海雜波多普勒頻移.雙基地角在雜波單元上是一個變量，因而，在距離和趨近于基線長度的過程中，一階海雜波譜具有漸進(jìn)展寬的現(xiàn)象;洋流的作用，不但會加劇這種展寬效應(yīng)，還會造成雜波頻率與理論值發(fā)生偏移.圖1是一次雙基地測量實(shí)驗(yàn)的布局，接收天線為8陣元線性陣列，工作頻率4.92 MHz.在接收天線陣列法向+35°方向有密集排列的島嶼暗礁.

在島嶼暗礁的方向上進(jìn)行波束形成，測量的海面回波譜如圖2，從圖2上可以看到，海雜波展寬現(xiàn)象非常明顯，同時，從第30個距離單元開始，譜峰呈現(xiàn)嚴(yán)重的分裂和急遽展寬形態(tài).圖3是距離門12上的功率譜，從圖3上可以看到，Bragg峰分裂成3個獨(dú)立譜峰，圖3中虛線為Bragg線理論值.這種形態(tài)分布在島嶼方向附近，方向波束的RD譜上持續(xù)距離長，各距離單元上形態(tài)一致，只是譜峰獨(dú)立程度有所差別，該形態(tài)在測量時段內(nèi)一直是持續(xù)的，因而，排除了這種形態(tài)有目標(biāo)和其他雜波的可能.從RD譜上看，具有雙基地海雜波特有的曲率特征.因而可以斷定，該形態(tài)譜峰是由Bragg峰分裂形成.這種形態(tài)對雜波檢測會遭成嚴(yán)重影響，基于理論預(yù)測值以及譜峰形態(tài)特征的檢測方法無法分辨具有獨(dú)立譜峰分裂形態(tài)且與理論預(yù)測值差別較大的海雜波，導(dǎo)致Bragg峰經(jīng)常被當(dāng)作目標(biāo)檢測出來.

圖1 測量實(shí)驗(yàn)雙基地布局

圖2 展寬和譜峰分裂的雙基地海雜波譜

圖3 距離門12上的功率譜

圖4是另一次實(shí)驗(yàn)的結(jié)果，波束方向?yàn)樘炀€陣列法向+48°(實(shí)驗(yàn)布局和工作條件如圖5所示)，譜展寬效應(yīng)表現(xiàn)不明顯，但Bragg線的理論值與實(shí)測值在近距離門端存在較大偏差，分析表明，該偏差和曲率不一致性是由地貌遮擋以及洋流作用共同產(chǎn)生的.地貌遮擋造成雙基地角估計(jì)不準(zhǔn)確，而真實(shí)的雙基地角卻難以估計(jì);洋流引起的頻率偏移也無法準(zhǔn)確估計(jì).同時，一個非常明顯的特征是:正負(fù)Bragg峰向?qū)τ诶碚擃A(yù)測值的偏移量不對稱，由于風(fēng)向的影響，正負(fù)Bragg峰譜寬度也不對稱.

圖4分析表明，雙基地海雜波譜峰位置和形態(tài)受多重因素影響，利用理論預(yù)測值來檢測雜波變得非常困難.實(shí)際處理結(jié)果也表明，海雜波經(jīng)常被當(dāng)作目標(biāo)來檢測和跟蹤，造成大量的虛假目標(biāo)，后續(xù)的航跡關(guān)聯(lián)變得異常復(fù)雜，形成大量虛假航跡.因而，海雜波對正常的目標(biāo)檢測和跟蹤造成嚴(yán)重影響，尋求能夠?qū)崟r處理的雜波抑制方法顯得非常必要.

圖4 理論值與實(shí)測值發(fā)生較大偏差的海雜波譜

1.2 雙基地海雜波譜分析

對于完全發(fā)展的海，海面可以看做平穩(wěn)的隨機(jī)過程，海雜波多普勒頻移是緩時變的.海態(tài)平穩(wěn)條件下，海雜波譜線位置通常比較固定.海雜波這一特性，可以用來事先在RD譜上確定海雜波的頻率范圍.該頻率范圍內(nèi)，可以全部當(dāng)做雜波來處理，事實(shí)上，對于復(fù)雜的雙基地雜波譜，即便目標(biāo)出現(xiàn)在此區(qū)域，也很難做出準(zhǔn)確的檢測.以第2次測量實(shí)驗(yàn)為例，說明雜波頻率范圍的確定方法.

首先分析該實(shí)驗(yàn)系統(tǒng)的布局.雙基地實(shí)驗(yàn)系統(tǒng)的布局示意圖如圖5所示.

圖5 雙基地布局和地貌示意圖

該系統(tǒng)距離和分辨率為6 km;基線長度為47.76 km;接收天線陣列由8個陣元組成，間距15 m;法向與基線的夾角為51.35°，雷達(dá)工作頻率3.776 MHz.圖5中，機(jī)動接收站所處地形復(fù)雜，導(dǎo)致理論上是全向的接收天線，由于地物遮擋變得覆蓋范圍有限，而且波束邊界也難以精確定位.

考察圖5中接收天線陣列陣元的波束圖形，有兩條頻率線可以確定.第1條頻率線是陣元法向所確定的Bragg線，因?yàn)殛囋ň€方向無地物遮擋，因而一定有一階海雜波譜線與之對應(yīng).第2條是各距離單元上的最大雙基地角所確定的該單元上的最小一階海雜波多普勒頻移.

考察第i個距離門，該距離門的最大距離和為L+i×ΔR，最小距離和為L+(i－1)×ΔR.設(shè)P為距離門內(nèi)任一點(diǎn)，由圖5中三角形TPR可知

式中:α為PR與基線的夾角;φ為P點(diǎn)對應(yīng)的雙基地角;r為該點(diǎn)對應(yīng)的橢圓距離和，er=r/L.P點(diǎn)對應(yīng)的雙基地角可寫為

式(2)即為雙基地角表達(dá)式.當(dāng)P點(diǎn)位于組成距離單元的最小距離和橢圓與基線垂直平分線的交點(diǎn)時，該點(diǎn)對應(yīng)的雙基地角為該距離門內(nèi)最大雙基地角.無論是利用橢圓的幾何圖形性質(zhì)，還是利用雙基地角表達(dá)式對α求偏導(dǎo)數(shù)，并令=0，可求得唯一極值點(diǎn)，都可以得到φ取得最大雙基地角時，有

圖5中接收天線陣列右側(cè)的地物遮擋，只影響第1個距離門上的最大雙基地角，當(dāng)距離門＞1時，最大雙基地角不受地物遮擋影響，是能夠唯一確定的，因而能夠在一階雜波譜中找到對應(yīng)的多普勒頻率分量.由式(2)、式(3)可求出對應(yīng)的雙基地角，由式(1)求出對應(yīng)的Bragg線.

圖6給出雙基地接收天線陣列陣元1的通道RD譜，數(shù)據(jù)參數(shù)為:相干積累時間為152.99 s;頻率分辨率為Δf=0.006 5 Hz;數(shù)據(jù)長度1 024點(diǎn).圖6中一階海雜波多普勒頻率的內(nèi)側(cè)邊界由最小多普勒頻率線平移逼近獲得，外側(cè)邊界由法線確定的多普勒頻率線平移逼近獲得.這兩條頻率線曲率并不一致，利用這一特性，平移這兩條譜線，使其在RD譜上逼近一階海雜波譜的邊界.經(jīng)過平移和目視匹配，如圖6中虛線所示，盡量將一階Bragg峰的頻率范圍包含到兩條虛線之中.由于左右譜峰主瓣有可能不對稱，因而平移的尺度可能不同，視具體的測量結(jié)果而定.

由于海面是平穩(wěn)隨機(jī)過程，在同一次測量過程中，8個陣元的一階海雜波譜峰都能涵蓋在這兩條經(jīng)過平移的頻率線之間.圖7給出了3 h后的陣元8的RD譜，可以看到，同樣的兩條頻率線，依然涵蓋了一階海雜波譜峰的頻率變化范圍.因而，這兩條頻率線所確定的頻率范圍，可以作為雜波抑制的基準(zhǔn)，用于批量處理通道數(shù)據(jù).

圖6 雙基地第1陣元RD譜

圖7 3 h后雙基地第8陣元RD譜

對于完全發(fā)展的海，實(shí)測譜數(shù)據(jù)研究表明，雙基地Bragg線頻率隨時間是緩變的，圍繞均值的波動，因而是一個平穩(wěn)的隨機(jī)變量，可以用緩時變頻率的正弦信號來模型化.同時，穩(wěn)定洋流造成Bragg峰位置偏移的范圍是可預(yù)知的，這一特征，反映到雜波譜隨時間尺度的變化特征就是:盡管隨時間的推移，Bragg峰頻率隨時間是變化的，但一階雜波譜總是在固定的頻率范圍內(nèi)出現(xiàn).同時，在任意時刻，Bragg峰在RD譜上的范圍也是固定的.即Bragg峰在該時刻的RD譜上隨距離尺度的變化是緩變的，同時范圍也是固定的.這兩個特性決定了這樣一種特征:如果用兩條曲線在任意時刻的RD譜上逼近Bragg峰頻率的兩條邊界，并適度留有余量，那么在后續(xù)的任何時刻，只要海態(tài)保持平穩(wěn)，雷達(dá)工作狀態(tài)穩(wěn)定，Bragg峰一定出現(xiàn)在這兩根曲線包含的頻率范圍內(nèi).這兩根直線確定的頻率范圍即Bragg線頻帶.雙基地Bragg線頻帶是恒定的、時不變的.Bragg線恒出現(xiàn)在Bragg線頻帶內(nèi).定義這樣一種特征，即Bragg線頻帶的時不變的特征，為雙基地一階海雜波的平穩(wěn)性特征.

Bragg線的平穩(wěn)性特征表明，既然是平穩(wěn)的隨機(jī)過程，時間平均可以取代集合平均，任何時段內(nèi)的統(tǒng)計(jì)特征都是集合平均的取樣.因而，在測量時段內(nèi)，只要海態(tài)是平穩(wěn)的，在任意時刻對波束或通道RD譜上對Bragg線頻帶進(jìn)行取樣，并適度留有余量，都可以做為整個測量時段的Bragg線頻帶的估計(jì).同時，距離單元內(nèi)的Bragg峰具有局部幅度占優(yōu)特征，在功率譜上一般以局部極大值或最大值點(diǎn)形式存在.由于具有譜峰分裂形態(tài)，極大值點(diǎn)可能不是唯一的.因而，在譜峰分裂等形態(tài)下，只要落入Bragg頻帶內(nèi)的極大值點(diǎn)，即可認(rèn)為是Bragg峰而進(jìn)行抑制.這些特征即為雜波抑制的基礎(chǔ)先驗(yàn)知識.

2 海雜波信號模型和抑制算法

2.1 海雜波信號模型

在實(shí)際處理中，單基地一階海雜波一直被當(dāng)作復(fù)正弦信號來處理.盡管雙基地一階海雜波譜因雙基地角影響而形成具有曲率的曲線，并且因距離單元上雙基地角的變化而形成譜擴(kuò)展，但距離門內(nèi)的譜形態(tài)依然表現(xiàn)為兩個明顯的譜峰，因而一樣可以用復(fù)正弦信號來模型化，至多因譜峰分裂造成頻率分量增加.由上述處理過程知道，事先已經(jīng)在RD譜上劃分出海雜波譜區(qū)域，該區(qū)域內(nèi)的雜波分量和疑似雜波的目標(biāo)分量，都將作為雜波分量被抑制.而出現(xiàn)在雜波譜區(qū)域內(nèi)的非機(jī)動目標(biāo)信號，也一樣可以用復(fù)正弦信號模型化，因此，假設(shè)雜波譜區(qū)域內(nèi)有K個復(fù)頻率分量(其中包括一階海雜波分量和疑似雜波信號的目標(biāo)信號、噪聲信號或者電離層雜波信號等)，雜波譜區(qū)內(nèi)的海雜波信號模型可以寫為

式中:Ai、fi、φi分別為對應(yīng)信號分量的幅度、頻率和初相;Ts為信號采樣周期;N為時域信號序列的長度.其中，fmin≤fi≤ fmax，[fmin，fmax]為確定的雜波頻率范圍.

2.2 海雜波抑制算法

所謂的海雜波抑制，是把位于雜波區(qū)域內(nèi)的能量從整體信號能量中分離出來，而又不改變其他信號的原貌．目的是使雜波能量不被當(dāng)做目標(biāo)檢測出來，從而實(shí)現(xiàn)正常的目標(biāo)檢測，而不是把雜波能量淹沒的信號提取出來．

據(jù)信號模型，結(jié)合海雜波信號具有頻率緩時變特性可知，由時域信號序列構(gòu)造的Hankel矩陣，奇異值分解之后的奇異值能夠跟蹤窄帶緩時變頻率的正弦信號分量，因此，基于奇異值分解的Hankel矩陣降秩法能夠?qū)崿F(xiàn)雜波抑制．但該算法計(jì)算速度慢，無法實(shí)現(xiàn)實(shí)時數(shù)據(jù)處理．

另一種雜波抑制思路，則是基于正弦信號參數(shù)估計(jì)．如果把信號模型中的正弦信號參數(shù)估計(jì)出來，在時域信號中將該分量減去，則消除了該信號分量的能量，從而實(shí)現(xiàn)雜波抑制．這種正弦信號參數(shù)估計(jì)的算法很多，比如Clean算法;改進(jìn)的Clean算法如Relax算法和ROOT循環(huán)對消算法;APES(正弦信號幅度相位估計(jì))算法;Capon算法等．后兩種算法是超分辨率譜估計(jì)算法，對正弦信號參數(shù)估計(jì)的更精確．Clean算法的優(yōu)點(diǎn)是速度快，便于實(shí)時大批量處理數(shù)據(jù)．

2.3 Clean算法

Clean算法［14］是一種基于FFT的快速譜估計(jì)算法．該算法的要求是:組成信號時域序列的各分量信號的持續(xù)時間長度都與進(jìn)行傅立葉變換處理的長度相同，這一點(diǎn)雜波信號模型無疑是能夠滿足的．針對復(fù)正弦信號分量組成的雜波模型，提取雜波信號能量的計(jì)算方法為:

1)對信號的時域序列S(n)作 FFT，得到S(ω)，求取S(ω)落在雜波頻率范圍內(nèi)的最大極大值點(diǎn)．該極值點(diǎn)所對應(yīng)的幅度，相位和頻率分別作為第1個分量信號的參數(shù)估計(jì)，即

第1個被提取的雜波分量可以寫為

將該分量從原時域序列中減去，提取該分量信號能量，即

對剩余信號Sr(n)重復(fù)上述過程，直到滿足收斂條件．這里進(jìn)行雜波抑制的目的是使雜波區(qū)域內(nèi)的能量不被當(dāng)做目標(biāo)檢測出來，因而只要雜波區(qū)域內(nèi)的最大極大值點(diǎn)的功率小于平均噪聲功率，即可認(rèn)為雜波抑制完成．而平均噪聲功率很容易事先求?。虼耍赏ㄟ^計(jì)算雜波區(qū)域內(nèi)最大極大值點(diǎn)的功率來設(shè)定收斂條件．

2)如果把相位的估計(jì)利用剩余能量最小原則來確定，則Clean算法就變成了ROOT循環(huán)對消算法．即在［0，2π］上搜索 φi，使得

取得最小值［10］．

Clean等譜估計(jì)方法用于信號能量提取，一般是從最大值點(diǎn)開始，反復(fù)提取剩余能量的最大值信號．但海面回波譜，尤其是通道數(shù)據(jù)，海雜波一般不是能量最大的信號．實(shí)驗(yàn)數(shù)據(jù)表明，在近距離門端，能量最大信號為地物回波和固定目標(biāo)回波，造成零頻能量最高．其次是航線上的大型艦船目標(biāo)，海雜波沒有表現(xiàn)為能量最強(qiáng)的信號．實(shí)際上，風(fēng)向?qū)ｋs波雙邊譜幅度的影響，也會造成海雜波至少有一側(cè)能量較低，因而，海雜波能量最強(qiáng)的狀況并不是普遍現(xiàn)象．這里把Clean算法應(yīng)用到固定頻率范圍內(nèi)，一方面是由于海雜波在功率譜上表現(xiàn)為極值點(diǎn)尖峰，盡管由于譜峰分裂或者其他信號疊加會造成信號分量增加，但既然無從辨別海雜波頻率范圍內(nèi)的信號屬何種信號類型，不如當(dāng)作海雜波信號抑制掉;另一方面，這樣處理并不違背原算法的信號參數(shù)估計(jì)準(zhǔn)則;再有，雜波頻率范圍內(nèi)頻率點(diǎn)有限，其中極大值點(diǎn)更少，迭代對消無需幾次，即可把雜波能量提取干凈，計(jì)算速度快，能夠?qū)崿F(xiàn)實(shí)時處理．實(shí)際處理結(jié)果也表明，這樣做是可行的．其他譜估計(jì)算法，即便是相對較快的ROOT算法，計(jì)算速度也遠(yuǎn)達(dá)不到實(shí)時處理．

3 雜波抑制處理結(jié)果

雜波抑制的處理過程為:首先利用給出的雜波譜特征先驗(yàn)知識，即Bragg線頻帶在通道RD譜或波束RD譜上的平穩(wěn)性特征，在RD譜上預(yù)先確定Bragg線頻帶;Bragg峰以極大值點(diǎn)的形式只存在于Bragg線頻帶內(nèi)，由于雙基地海雜譜中，正負(fù)Bragg峰在洋流的作用下發(fā)生的譜峰偏移不一定具有對稱性，因而只要落入Bragg頻帶內(nèi)的極大值點(diǎn)，都將被作為雜波抑制掉．對于頻帶內(nèi)的極大值點(diǎn)，利用復(fù)正弦信號作為模型，利用Clean算法對模型參數(shù)進(jìn)行估計(jì)，之后在信號的時域序列中將該信號模型的估值減去，從而抑制掉頻帶內(nèi)所有極大值點(diǎn)對應(yīng)的頻率分量，即提取該分量的信號能量．剩余能量的信號即為雜波抑制后的時域序列．

給出Clean算法得到的雜波抑制結(jié)果，這里采用的是波束形成之前的通道數(shù)據(jù)，在時域完成雜波抑制之后，再進(jìn)行數(shù)字波束形成，得到雜波抑制后波束方向上的RD譜．圖8，圖9給出了雜波抑制前后16°方向上的RD譜.圖9上可以看出雜波能量明顯被抑制掉，甚至遠(yuǎn)低于背景噪聲功率，因而不會被當(dāng)作目標(biāo)檢測出來.圖10，圖11給出了－32°方向上雜波抑制前后的RD譜.從圖形上看，雜波抑制處理效果相當(dāng)明顯，但Bragg頻帶內(nèi)的目標(biāo)信號也同時被抑制掉了，因而損失了這部分信息.

圖8 未抑制雜波16°方向RD譜

圖9 Clean算法處理后16°方向上的RD譜

圖10 未抑制雜波－32°方向RD譜

圖11 Clean算法處理后－32°方向上的RD譜

作為對比，利用奇異值分解算法進(jìn)行同樣的雜波抑制處理，結(jié)果顯示，兩種方法在處理效果上沒有明顯差別，但奇異值分解算法有明顯的能量殘留，因而Clean算法略優(yōu)于奇異值分解算法.但計(jì)算速度上Clean算法遠(yuǎn)比奇異值分解算法快得多.以處理一個通道120個距離單元的速度對比來說，Clean算法在普通個人PC上運(yùn)行，只需要不到0.5 s，而奇異值分解算法則需要6 h以上.因而Clean算法可以實(shí)現(xiàn)實(shí)時處理.

同樣的處理方法也可以應(yīng)用到雜波影響嚴(yán)重的波束方向上，不過需要重新快速的擬合出兩條邊界曲線.本文所謂的實(shí)時，是指在雷達(dá)工作開始的一小段時間，可以通過先驗(yàn)知識確定雜波頻率范圍，之后的處理過程，只要海面是平穩(wěn)的，可以一直延續(xù)下去.

對于Bragg線頻帶內(nèi)損失的信號，由于系統(tǒng)是多頻同時工作的，多頻可以部分解決信號損失的問題.如果大量信號出現(xiàn)在該頻帶內(nèi)，則針對該頻帶內(nèi)的雜波形態(tài)特征進(jìn)行雜波檢測，研究表明，如果雜波分裂，則在相當(dāng)?shù)慕嵌确秶途嚯x范圍內(nèi)具有同樣的形態(tài)特征.

4 結(jié)論

1)由于雙基地海雜波受多重因素影響，比如洋流和島嶼等因素引起譜峰嚴(yán)重展寬和分裂;洋流和地貌等造成雜波頻率理論預(yù)測值與實(shí)測值發(fā)生較大偏差.

2)在距離和較小時，雜波譜展寬等，導(dǎo)致雜波背景下的目標(biāo)檢測變得異常困難，因而，從通道數(shù)據(jù)入手，采用一種類似于強(qiáng)制性“挖除”的方法來抑制海雜波.

3)提出的算法主要實(shí)現(xiàn)的是犧牲少量淹沒在雜波區(qū)內(nèi)的目標(biāo)信號，達(dá)到正常目標(biāo)的檢測.對雙基地來講，海雜波形成速度譜上的檢測盲區(qū)，淹沒于其中的目標(biāo)信號本身確實(shí)難以檢測和分離，無法準(zhǔn)確判定是目標(biāo)信號還是雜波信號.

[1]BARRICK D.First-order theory and analysis of MF/HF/VHF scatter from the sea[J].IEEE Transactions on Antennas and Propagation，1972，20(1):2－10.

[2]TRIZNA D B.A model for the First-order doppler Spectrum for bistatic HF radar surface wave sea scatter[C]//IGARRS 2001.Sydney:Australia，2001:7 －11.

[3]TRIZNA D B.A bistatic HF radar for current mapping and robust ship tracking[C]//OCEANS 2008.Wash-ington，DC:IEEE，2008:1－6.

[4]GILL E，HUANG W，WALSH J.The effect of the bistatic scattering angle on the High-Frequency radar cross sections of the ocean surface[J].Geoscience and Remote Sensing Letters，IEEE，2008，5(2):143 －146.

[5]POON M W Y，KHAN R H，LE-NGOC S.A singular value decomposition(SVD)based method for suppressing ocean clutter in high frequency radar[J].IEEE Transactions on Signal Processing，1993，41(3):1421－1425.

[6]HICKEY K，KHAN RH，WALSH J.Parametric estimation of ocean surface currents with HF radar[J].IEEE Journal of Oceanic Engineering，1995，20(2):139 －144.

[7]KHAN R H.Ocean-clutter model for high-frequency radar[J].IEEE Journal of Oceanic Engineering，1991，16(2):181－188.

[8]Defence R＆D Canada，Raytheon Canada Ltd.Experiments with high frequency surface wave radar and development of a practical ground screen method for HFSWR at cape bonavista[R].Report Number:DRDC-OTTAWA-CR-1999-018;RCL-DND-309;ContractReport(CAN)，1998－11－17.

[9]XIE J，YUAN Y，LIU Y.Super-resolution processing for HF surface wave radar based on Pre-Whitened MUSIC[J].IEEE Journal of Oceanic Engineering，1998，23(4):313－321.

[10]ROOT B.High-frequency over-the-horizon radar ship detection through clutter cancellation:an alternative to high-resolution spectral estimation[C]//Proceedings of Advanced Signal Processing Algorithms，Architectures，and Implementations VII. San Diego， CA:SPIE，1998:488－500.

[11]TONY P，REZA D.System and method for spectral generation in radar:US Patent，6822606[P].2004 －09－16.

[12]Barrick D E，ISAACSON J，LILLEBOE P M，et al.O-cean surface current mapping with bistatic HF radarz:US Patent，6774837[P].2004 －08 －10.

[13]楊紹麟，柯亨玉，侯杰昌，等.OSMAR2000基于MUSIC的超分辨率海洋表面流算法[J].武漢大學(xué)學(xué)報，2001，47(5):601－608.

[14]鄒虹，保錚.多分量線性調(diào)頻信號的時頻分析[D].西安:西安電子科技大學(xué)，2001.