基于已知傳播模式數(shù)目的海雜波抑制方法研究

李雪 郭曉彤 王岳松 李吉寧

(中國電波傳播研究所,青島 266107)

?

基于已知傳播模式數(shù)目的海雜波抑制方法研究

李雪 郭曉彤 王岳松 李吉寧

(中國電波傳播研究所,青島 266107)

天波超視距雷達海雜波譜受電離層多模傳播影響大幅展寬,需進行雜波抑制才可實現(xiàn)艦船目標檢測,利用電離層返回散射探測可獲得海雜波譜先驗模式數(shù)目信息. 針對常用的兩種海雜波抑制方法——特征值分解法和循環(huán)對消法,仿真分析了傳播模式數(shù)目識別準確性對海雜波抑制效果的影響. 結(jié)果表明:若傳播模式數(shù)目識別錯誤,則會造成目標漏檢或虛警. 在傳播模式數(shù)目已知的前提下,利用仿真和試驗數(shù)據(jù)研究了目標與一階海雜波多普勒相對位置和多普勒分辨率對上述兩種方法海雜波抑制效果的影響. 目標與一階海雜波越接近,海雜波抑制效果越差,當二者相對多普勒位置小于分辨率時,兩種方法均失效;多普勒分辨率越低,海雜波抑制效果越差;目標與海雜波越接近,對多普勒頻率分辨率要求越高.

天波超視距雷達海雜波譜;傳播模式數(shù)目;奇異值分解法;循環(huán)對消法

DOI 10.13443/j.cjors.2015110902

引 言

天波超視距雷達(Over-The-Horizon Radar, OTHR)利用電離層作為傳輸介質(zhì),可超視距覆蓋幾千千米范圍,實現(xiàn)對彈道導彈、隱身飛機等目標的早期預警,應用前景十分廣闊. 而對艦船等低速目標而言,其多普勒頻譜與海雜波的一階譜分量十分接近,目標檢測難度較大. 又由于電離層的分層特性,信號回波會出現(xiàn)多模效應,導致雜波譜展寬比較嚴重,且模式數(shù)目不定,進一步加大了低速目標的檢測難度[1-2].

為抑制海雜波對目標檢測的影響,目前主要有兩大類方法:第一類是海雜波循環(huán)對消法,將海洋回波的多普勒信號近似認為呈正弦變化,估計一階海雜波的頻率、幅度、相位,在時域上構(gòu)造對應的正弦信號,與回波相減消去一階海雜波,凸顯目標,文獻[3]提出了基于快速傅里葉變換(Fast Fourier Transform,FFT)相位分析的海雜波對消方法,進一步提高了一階海雜波的頻率、相位估計精度,取得了較好效果;第二類是特征值分解法,其中使用頻率最高的是奇異值分解(Singular Value Decomposition,SVD)法,對回波信號進行特征值分解,將反映一階海雜波的奇異值置零,然后重構(gòu)回波信號,實現(xiàn)一階海雜波抑制,凸顯目標信號[4]. 然而,上述兩類方法要求比較準確地估計傳播數(shù)目,即一階海雜波譜峰數(shù)目,否則不能選擇合適的對消次數(shù)或置零個數(shù),嚴重影響海雜波抑制效果,從而導致無法有效實現(xiàn)艦船目標檢測. 針對傳播模式數(shù)目估計問題,文獻[5]提出基于先驗返回散射探測的傳播模式區(qū)提取方法,實現(xiàn)了任意距離單元的傳播模式數(shù)目準確提取.

本文針對常用的兩種海雜波抑制方法——SVD法和基于FFT相位分析雜波循環(huán)對消法,仿真分析了傳播模式數(shù)目識別準確性對海雜波抑制效果的影響. 并在假設傳播模式數(shù)目已知的前提下,利用仿真和試驗數(shù)據(jù)研究了目標與一階海雜波多普勒相對位置和多普勒分辨率對上述兩種方法海雜波抑制效果的影響.

1 海雜波模型

1955年,Crombie提出用諧振散射機理描述高頻無線電波與海浪的一階作用. 根據(jù)該理論,海洋表面波可看成一系列正弦信號的疊加. 假設其波長為L,雷達發(fā)射的無線電磁波波長為λ,入射角為α,當波長滿足Lcosα=λ/2時,各反射點回波的波程差正好為λ的整數(shù)倍,此時,回波信號同相相加,產(chǎn)生一階Bragg散射,且一階Bragg峰對應的多普勒頻率為

(1)

式中: g為重力加速度; α為海浪與雷達照射方向夾角; f0為雷達載頻,MHz; ±號為朝向及背離雷達波束的諧振海浪所產(chǎn)生的正負Bragg峰[6]. 譜峰的相對強度反映了引起海浪的表面風風向. 一階譜峰的強度很大,是海洋回波的主要成分,在高頻海洋回波譜中,表現(xiàn)為一對關(guān)于零頻對稱的譜峰,若存在電離層多模傳播效應,回波譜中還會存在多對一階Bragg譜峰. 本文將利用頻率為fB的兩個正弦信號模擬一階海雜波,同時通過添加不同的多普勒頻移模擬高頻雷達多模傳播下海洋回波譜.

2 海雜波抑制方法

2.1 奇異值分解法

海雜波可以看成是窄帶的時變信號,用時變信號處理技術(shù)能進行頻率跟蹤[7]. 因此,將雷達回波的時域信號構(gòu)造成Hankel矩陣,對其進行奇異值分解,產(chǎn)生的奇異值與一階海雜波的頻率分量一一對應,通過把對應一階海雜波的奇異值置零實現(xiàn)海雜波的抑制.

(2)

矩陣H經(jīng)奇異值分解后表示為

H=USVT.

(3)

式中: S=diag(σ1,σ2,…,σC),σ1≥σ2≥…≥σC≥0是矩陣H的奇異值; U、V為奇異值對應的左、右奇異矢量. 奇異值反映了信號、雜波干擾和噪聲的能量集中情況,前r個較大的奇異值依次主要反映大的信號能量,較小的奇異值對應噪聲,中間部分則主要對應于二階海雜波、慢目標及其他雜波. 由于一階海雜波集中在正負Bragg峰附近,且遠高于目標及其他雜波強度,與較大的奇異值對應,因此,將一階海雜波奇異值置零,重構(gòu)信號后能較好地消除一階海雜波,實現(xiàn)艦船目標檢測.

2.2 基于FFT相位分析的海雜波循環(huán)對消

海雜波循環(huán)對消思想是利用FFT找到一階海雜波的主要分量,通過估計該分量的幅度、頻率和初相重構(gòu)正弦信號,在回波信號中將其減去,循環(huán)迭代此過程直到一階海雜波去除,凸顯出目標信號為止. 因此,這種方法要求對正弦信號的幅度、頻率、相位的估計非常精確. 但是傳統(tǒng)的海雜波對消法頻率分辨率低,實際的一階譜峰對應的頻率可能不在離散頻率點上,導致估計精度下降,剩余雜波產(chǎn)生擴散.

(4)

式中,T為積累時間. 當Δf超出變化范圍[-0.5/T, 0.5/T]時,需要進行

(5)

(6)

φm=φ0+(1-1/N)ΔfTπ.

(7)

3 仿真研究

3.1 傳播模式數(shù)目識別準確性對海雜波抑制影響研究

假設雷達工作頻率為10 MHz,脈沖重復周期為0.2 s,相干積累次數(shù)為512,一階Bragg峰對應的多普勒頻率為0.322 6 Hz. 同時假設當前電離層存在E、F兩個傳播模式,其中E模式引起的相對多普勒頻移為0.05 Hz,其回波雜噪比為35 dB,F模式引起的相對多普勒頻移為0.07 Hz,其回波雜噪比為30 dB. 目標多普勒位于0.35 Hz,信噪比25 dB,該回波信號如圖1黑色線所示. 單從譜圖上只能看到一對一階譜峰,然而該譜峰中包含著E、F兩種傳播模式下的一階海雜波和目標,這就存在傳播模式數(shù)目識別問題.

(a) SVD抑制結(jié)果

(b) 基于FFT相位分析循環(huán)對消抑制結(jié)果圖1 傳播模式數(shù)目判少的抑制結(jié)果

利用基于SVD抑制海雜波和基于FFT相位分析循環(huán)對消兩種方法對上述雷達回波信號進行海雜波抑制. 若按照傳播模式數(shù)目為1對海雜波進行抑制,則只能消除其中一種模式下的一階海雜波,如圖1所示,此時剩余三個較大的譜峰,會均被視為目標,而實際上真實目標只有一個,另外兩個均為模式F的一階譜峰,從而造成虛警. 若按照傳播模式數(shù)目為3對海雜波進行抑制,則會連目標一并消除,如圖2所示,此時會認為沒有目標存在,產(chǎn)生漏警. 因此,當傳播模式數(shù)目識別與真實的傳播模式數(shù)目不一致時,無法實現(xiàn)有效的目標檢測.

若按照傳播模式數(shù)目為2分別對海雜波進行抑制,則結(jié)果如圖3所示,E、F模式下的一階海雜波均被消除,僅凸顯出真實目標位置.

通過以上仿真可看出,只有當傳播模式數(shù)目識別與真實傳播模式數(shù)目一致時,基于SVD海雜波抑制法和基于FFT相位分析循環(huán)對消法才能有效檢測出真實目標. 因此,傳播模式數(shù)目識別的準確性對海雜波抑制效果影響非常大,抑制前必須獲得正確的傳播模式數(shù)目,文獻[5]已提出有效方法,利用基于先驗返回散射探測的傳播模式區(qū)提取傳播模式數(shù)目.

(a) SVD抑制結(jié)果

(b) 基于FFT相位分析循環(huán)對消抑制結(jié)果圖2 傳播模式數(shù)目判多的抑制結(jié)果

(a) SVD抑制結(jié)果

(b) 基于FFT相位分析循環(huán)對消抑制結(jié)果圖3 傳播模式數(shù)目一致的結(jié)果

3.2 目標與海雜波之間的位置關(guān)系對抑制效果影響

假設雷達工作頻率、脈沖重復周期、相干積累次數(shù)、一階Bragg峰對應的多普勒頻率與3.1節(jié)仿真條件一致,并假設當前電離層僅存在一個傳播模式E,其回波雜噪比為35 dB,目標多普勒頻率在[0.29, 0.36] Hz范圍內(nèi)以0.01 Hz步進,信噪比為25 dB.

利用基于SVD抑制海雜波和基于FFT相位分析循環(huán)對消兩種方法對上述雷達回波信號進行海雜波抑制,以目標與一階海雜波多普勒之差為相對位置,統(tǒng)計出目標與海雜波相對位置關(guān)系對抑制效果影響的曲線圖,如圖4所示.

從圖4可以看出,當相對位置為-0.03 Hz時,兩種方法抑制后的信雜比在20 dB左右,隨著目標與一階海雜波相對位置減小,抑制后的信雜比逐漸降低,當目標與海雜波相距小于頻率分辨率時,兩種方法均不能正確檢測出目標,而隨著目標與一階海雜波相對位置增大,抑制后的信雜比再次逐漸升高. 因此,當目標接近一階海雜波時,兩者相對位置越遠,海雜波抑制效果越好,且基于FFT相位分析循環(huán)對消的抑制效果略好.

3.3 多普勒分辨率對海雜波抑制效果影響

假設雷達工作頻率、脈沖重復周期、一階Bragg峰對應的多普勒頻率、傳播模式與3.2節(jié)仿真條件一致,減少采樣點數(shù),使其在[670, 220]范圍內(nèi)以-30遞減.

設目標多普勒為0.36 Hz,與海雜波相對位置為0.03 Hz,信噪比25 dB,采用SVD和基于FFT相位分析循環(huán)對消法分別進行海雜波抑制,結(jié)果如圖5(a)所示. 改變目標位置,使其位于0.35 Hz和0.34 Hz處,與海雜波的相對位置為0.02 Hz和0.01 Hz,抑制后信雜比分別如圖5(b)、(c)所示.

從圖5(a)可以看出,當多普勒頻率間隔為0.007 Hz時,抑制后的信雜比為21 dB,當多普勒頻率間隔增加到0.022 Hz時,SVD抑制海雜波后的信雜比降到15 dB,基于FFT相位分析循環(huán)對消抑制后的信雜比降到16.2 dB. 因此,隨著多普勒分辨率下降,兩種方法的信雜比都隨之降低,相比之下,基于FFT相位分析循環(huán)對消的效果較好,圖5(b)、圖5(c)均能反映出該變化過程.

(b) 目標相對位置0.02 Hz

(c) 目標相對位置0.01 Hz圖5 隨頻率分辨率變化的兩種方法抑制后的信雜比

對比圖5的三幅圖可知:當目標與海雜波相對位置為0.03 Hz時,多普勒頻率間隔大于0.022 Hz時就無法凸顯出真實目標;當目標與海雜波相對位置為0.02 Hz時,頻率間隔大于0.018 Hz時SVD法失效,大于0.022 Hz時基于FFT相位分析循環(huán)對消法失效;當目標與海雜波相對位置為0.01 Hz時,頻率間隔大于0.014 Hz時SVD法失效,大于0.018 Hz時基于FFT相位分析循環(huán)對消法失效. 因此,隨著目標與海雜波相對位置的減小,需要的多普勒分辨率隨之提高.

4 實測數(shù)據(jù)驗證

考慮到OTHR數(shù)據(jù)來源受限,同時高頻窄波束返回散射定頻探測數(shù)據(jù)與OTHR回波數(shù)據(jù)特性相同,均為雙程傳播,且均受電離層多模傳播影響,因此可利用其進行上述算法驗證.

獲取2009年4月19日5時14分某高頻返回散射探測站定頻海雜波回波功率譜,工作頻率20.528 MHz,群距離1 550 km,脈沖重復周期0.05 s,相干積累次數(shù)256,采樣頻率80 kHz. 在多普勒頻率-0.6 Hz和-0.5 Hz處分別添加目標信號,合成的信號譜如圖6、圖7黑色線所示. 結(jié)合先驗返回散射探測的傳播模式區(qū)提取到傳播模式數(shù)目為1,兩種方法的抑制結(jié)果分別如圖6、圖7紅色線所示.

目標位于-0.6 Hz時,SVD抑制后的信雜比為43.249 2 dB,基于FFT相位分析法消除后的信雜比為44.003 0 dB,兩種方法抑制效果相近,基于FFT相位分析循環(huán)法的效果略優(yōu)于SVD. 而當目標位于-0.5 Hz時,兩種方法均無法實現(xiàn)目標檢測. 因此,目標與海雜波越接近,抑制效果越差,與仿真結(jié)果一致.

仿真條件不變,僅將相干積累次數(shù)減少,依次更改為192和128,兩種方法的抑制結(jié)果分別如圖8、圖9所示.

當積累次數(shù)為192時,SVD抑制海雜波后的信雜比為41.070 3 dB,基于FFT相位分析法抑制后的信雜比為42.457 5 dB,較積累次數(shù)為256,兩種方法的抑制效果均下降,但基于FFT相位分析的效

(a) 目標位于-0.6 Hz處

(b) 目標位于-0.5 Hz 處 圖6 SVD抑制海雜波(實測數(shù)據(jù))

果仍然相對較好. 當積累次數(shù)為128時,兩種方法均無法正確檢測目標. 因此,隨著分辨率下降,海雜波抑制效果逐漸變差,同時結(jié)合圖6、圖7可以看出,當目標與海雜波相對位置越近,對多普勒分辨率的要求越高,與仿真結(jié)果一致.

(a) 目標位于-0.6 Hz處

(b) 目標位于-0.5 Hz 處圖7 基于FFT相位分析法消除海雜波(實測數(shù)據(jù))

(a) 積累次數(shù)192

(b) 積累次數(shù)128圖8 不同分辨率下SVD抑制結(jié)果

(a) 積累次數(shù)192

(b) 積累次數(shù)128圖9 不同分辨率下基于FFT相位分析抑制結(jié)果

5 結(jié) 論

本文通過仿真多模效應下奇異值分解法和基于FFT相位分析法抑制海雜波的結(jié)果,反映出已知傳播模式數(shù)目的重要性,并從目標與海雜波相對位置關(guān)系和多普勒頻率分辨率對海雜波抑制結(jié)果的影響兩個方面對比分析兩種方法的目標檢測能力. 經(jīng)試驗驗證,目標與海雜波相對位置越近,抑制效果越不好;多普勒頻率分辨率越低,抑制效果越不好. 目標與海雜波越接近,對多普勒頻率分辨率要求越高. 在同等條件下,基于FFT相位分析法的抑制效果較好.

[1] 周文瑜, 焦培南. 超視距雷達技術(shù)[M]. 北京:電子工業(yè)出版社, 2008: 224.

[2] 凡俊梅, 焦培南, 肖景明. 海洋雜波對高頻雷達檢測海面上低速目標的影響[J]. 電波科學學報, 1997, 12(2): 205-210.

FAN J M, JIAO P N, XIAO J M. The sea clutter effect on the low Doppler targets detection by HF radar[J]. Chinese journal of radio science, 1997, 12(2): 205-210.(in Chinese)

[3] 郭欣, 倪晉麟, 劉國歲. 短相干積累條件下天波超視距雷達的艦船檢測[J]. 電子與信息學報, 2004,26(4):613-618.

GUO X, NI J L, LIU G S. The ship detection of sky wave over-the-horizon radar with short coherent interaction time[J]. Journal of electronics information technology, 2004, 26(4): 613-618. (in Chinese)

[4] 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.

[5] 李雪, 李吉寧, 馮靜, 等. 天波超視距雷達多模傳播抑制[J].電波科學學報, 2015, 30(1):57-62.

LI X, LI J N, FENG J, et al. OTHR multi-mode propagation suppression[J]. Chinese journal of radio science, 2015, 30(1): 57-62. (in Chinese)

[6] ANDERSON S J, MAHONEY A R, TURLEY M D E. Applications of superresolution techniques to HF radar sea echo analysis[C]//Proceedings of the Pacific Ocean Remote Sensing Conference, Melbourne, 1994.

[7] KHAN R H. Ocean-clutter model for high frequency radar[J]. IEEE journal of oceanic engineering, 1991,16(2):181-188.

李雪 (1981-),男,黑龍江人,博士,高級工程師,目前主要研究方向為電離層回波信號處理、雷達信號處理、電波傳播等.

郭曉彤 (1989-),女,河北人,助理工程師,目前主要研究方向為電離層電波傳播、雷達信號處理.

王岳松 (1981-),男,山東人,工程師,目前主要研究方向為計算機算法及應用、大數(shù)據(jù)分析、模式識別算法以及頻率選擇方法.

李吉寧 (1982-),男,山東人,工程師,目前主要研究方向為高頻雷達信號處理和電波傳播.

Sea clutter suppression algorithm based on ionospheric propagation mode number

LI Xue GUO Xiaotong WANG Yuesong LI Jining

(ChinaResearchInstituteofRadioWavePropagation,Qingdao266107,China)

Since the sea clutter spectrum of over-the-horizon radar(OTHR) always has large spreads due to the ionospheric multimode propagation, the clutter suppression needs to be done for detection of ship targets. The information of a priori mode number for the sea clutter spectrum can be acquired via ionospheric backscatter sounding. The effects of the identification for propagation mode number on the suppression of sea clutter are simulated and analyzed based on two common methods of sea clutter suppression, namely, the singular value decomposition (SVD) algorithm and the phase information of Fourier transform algorithm. The results show that target false-alarm or false-dismissal will occur if the identification of propagation mode number is not correct. Furthermore, the influence of the distance between the target and the first-order sea clutter and the Doppler resolution on the sea clutter suppression using the above two methods is studied through both simulated and experimental data with the precondition of known propagation mode number, and consequently two conclusions are obtained. The first conclusion is that the closer the target is to the first-order sea clutter, the worse the suppression of sea clutter operates, and when the distance of the two is within the Doppler resolution range, neither of the two methods works. The second one is that the sea clutter suppression effect is poorer with the lower Doppler resolution, and closer distance between target and clutter will require higher Doppler frequency resolution.

OTHR sea clutter spectrum; propagation mode number; SVD; the phase information of Fourier transform

10.13443/j.cjors.2015110902

2015-11-09

國家自然科學基金青年基金(No.61302006)； 中國電科技術(shù)創(chuàng)新基金(JJ-QN-2013-28)

TN011.2

A

1005-0388(2016)04-0700-07

李雪, 郭曉彤, 王岳松, 等. 基于已知傳播模式數(shù)目的海雜波抑制方法研究[J]. 電波科學學報,2016,31(4):700-705.

LI X, GUO X T, WANG Y S, et al. Sea clutter suppression algorithm based on ionospheric propagation mode number[J]. Chinese journal of radio science,2016,31(4):700-705. (in Chinese). DOI: 10.13443/j.cjors.2015110902

聯(lián)系人： 郭曉彤 E-mail: xtguo1206@sina.cn