一種天波雷達多徑擴展多普勒雜波抑制方法

羅 歡 陳建文 鮑 拯

?

一種天波雷達多徑擴展多普勒雜波抑制方法

羅 歡*①陳建文②鮑 拯②

①(空軍預警學院研究生管理大隊 武漢 430019)②(空軍預警學院三系 武漢 430019)

該文分析了天波超視距雷達(Over The Horizon Radar, OTHR)多徑擴展多普勒雜波(Spread Doppler Clutter, SDC)的產(chǎn)生機理。由于陣列存在幅相誤差且期望信號的功率大于SDC功率，自適應(yīng)數(shù)字波束形成(Adaptive Digital Beam Forming, ADBF)將降低SDC抑制能力，同時還會導致信號對消，嚴重降低信噪比。針對以上問題，該文提出一種自適應(yīng)抑制SDC的方法。該方法首先采用改進噪聲子空間擬合自校正法消除陣列幅相誤差，得到期望信號和SDC準確的到達仰角，然后采用正交投影權(quán)矢量進行ADBF處理，避免了強期望信號條件下ADBF權(quán)矢量估計不準的問題。理論分析和仿真實驗表明該方法能夠較徹底地抑制多徑SDC。

天波超視距雷達；多徑擴展多普勒雜波；幅相誤差校正；自適應(yīng)數(shù)字波束形成

1 引言

天波超視距雷達(Over The Horizon Radar, OTHR)工作頻段為3~30 MHz，利用電離層對高頻電波的反射可實現(xiàn)對800~3500 km的海面和空中目標進行探測[1,2]。在OTHR某些工作頻率，電磁波可能經(jīng)過不同層高的電離層折射，使得來自于不同距離的回波信號經(jīng)過不同層高的電離層反射后，在同一時間到達接收陣列，由于群路徑相同，接收信號多普勒譜表現(xiàn)為幾組不同反射信號多普勒譜的疊加[2,3]。在實際工程中，由于電磁波傳播時受到非主選電離層(如F層)非線性擾動的影響，經(jīng)其反射回來的多普勒譜被嚴重展寬的幾率很高，展寬的范圍也較大，甚至可覆蓋-40~40 Hz之間的多普勒區(qū)域[4]，將這種經(jīng)過非主選電離層傳播而被接收的重影回波稱為多徑擴展多普勒雜波(Spread Doppler Clutter, SDC)。如果穩(wěn)定電離層(一般為E層，主選電離層)回波與擾動電離層(如F層)的SDC疊加在一起被OTHR陣列同時接收，則目標將被SDC覆蓋。

本文主要研究由接收陣列旁瓣進入的非主選電離層多徑SDC的抑制問題。OTHR是大型陣列，不可避免的陣列幅相誤差將降低仰角DOA的估計精度和ADBF的性能；同時由于期望信號(經(jīng)主選層返回的海雜波)的功率大于SDC功率，嚴重影響協(xié)方差矩陣的性能，在進行ADBF時將導致信號對消，嚴重降低信雜噪比[13]。本文首先提出了改進的陣列幅相誤差自校正方法消除陣列誤差，然后利用正交投影權(quán)矢量進行ADBF處理抑制多徑SDC，可以較好解決以上問題。

2 問題描述

2.1 多徑SDC產(chǎn)生機理分析

多徑SDC的產(chǎn)生原理可通過圖1來說明。理想條件下，位于處的OTHR以穩(wěn)定的電離層E層為信道(主選信道)，經(jīng)由----探測位于處的艦船目標和其上空的飛機目標，記為E-E模式。不失一般性，以擾動的F層為例，經(jīng)過F層傳播產(chǎn)生的多徑SDC有3種，一是發(fā)射信號經(jīng)過擾動F層折射后產(chǎn)生多普勒擴展，在海面反射后經(jīng)過E層返回被接收，記為F-E模式，其傳輸路徑如圖1中----所示；二是發(fā)射信號經(jīng)E層折射到達海面，在海面反射后經(jīng)F層傳播產(chǎn)生多普勒擴展后被接收，記為E-F模式，其傳輸路徑如圖1中----所示；三是發(fā)射和接收兩次都經(jīng)過F層傳播，經(jīng)過兩次展寬影響，記為F-F模式，其傳輸路徑如圖1中----所示。F-F模式的回波受擾動影響大，但是強度較弱，而E-F和F-E模式的回波較強，對E-E模式穩(wěn)定的回波影響較大。

圖1 多徑SDC的產(chǎn)生機理示意圖

由于主選和非主選電離層的高度差異使得不同多徑信號的DOD或者DOA差異明顯。具體來講，E-E與E-F發(fā)射路徑經(jīng)過相同的電離層(E層)，對應(yīng)的仰角DOD近乎相同；接收路徑經(jīng)過不同電離層，對應(yīng)的仰角DOA差異較大。E-E與F-E發(fā)射路徑經(jīng)過不同電離層，對應(yīng)的仰角DOD差異較大；接收路徑經(jīng)過相同電離層，對應(yīng)的仰角DOA幾乎相同。F-F模式與E-E模式的DOD和DOA都有明顯差異。DOD差異明顯的SDC可采用發(fā)射自適應(yīng)空域處理進行抑制；DOA差異明顯的SDC可采用接收自適應(yīng)空域處理進行抑制。本文研究的由接收陣列旁瓣進入的非主選層多徑SDC屬于DOA差異明顯的SDC抑制問題，即抑制E-F傳播模式下的SDC。

2.2 陣列模型

2.3 存在陣列幅相誤差時的多徑SDC抑制問題

若陣列不存在幅相誤差且快拍數(shù)據(jù)只含SDC和噪聲時，采用常規(guī)ADBF方法即可抑制SDC，例如最小方差無失真響應(yīng)(MVDR)自適應(yīng)波束形成法[15]。但是，當陣列存在幅相誤差且快拍數(shù)據(jù)中還含有期望信號時，SDC對目標遮蓋明顯，而且期望信號(主瓣海雜波)的功率一般大于SDC(旁瓣雜波)的功率，直接使用含有期望信號的數(shù)據(jù)構(gòu)造的協(xié)方差矩陣會導致信號相消。另外，陣列幅相誤差不僅會影響DOA估計精度，還會導致ADBF性能下降。下面對此進行說明。

當陣列存在幅相誤差時，式(1)可寫為

3 多徑SDC抑制方法

從第2節(jié)的分析可見，OTHR陣列幅相誤差以及期望信號功率較大的現(xiàn)實情況均會造成SDC抑制效果嚴重降低，本節(jié)在消除陣列幅相誤差的基礎(chǔ)上，采用基于到達仰角DOA估計的正交權(quán)矢量ADBF算法，能夠較徹底地抑制SDC，且不會導致期望信號對消。

3.1 OTHR陣列幅相誤差自校正

陣列幅相誤差校正常用的方法主要有有源校正法和自校正法。有源校正法通過在空間設(shè)置方位精確已知的輔助源對陣列誤差進行離線估計，而自校正方法通常根據(jù)優(yōu)化函數(shù)對陣列誤差進行在線估計。有源校正法需要輔助信源的方位精確已知，如果輔助源的方位信息出現(xiàn)偏差，會降低陣列校正效果，對于天波雷達還必須部署相應(yīng)的專用有源校正目標，因此成本較高且難以實現(xiàn)隨時校正[16]；而自校正方法不需要輔助源，最經(jīng)典的屬文獻[17]提出的噪聲子空間擬合(Noise Subspace Fitting, NSF)自校正法，但該方法對誤差矩陣迭代初值的選取有較高要求，若選擇的迭代初值與實際陣列幅相誤差差異較大，NSF法將失效。針對NSF自校正法的實際問題，本節(jié)對其進行改進，以達到穩(wěn)健校正OTHR陣列幅相誤差的目的。

(4)進行MUSIC空間譜估計

3.2 基于正交投影權(quán)矢量的自適應(yīng)波束形成方法

綜上所述，本文提出的抑制由OTHR接收旁瓣進入的多徑SDC的方法流程為：

(2)采用改進NSF自校正方法消除陣列幅相誤差，并得到期望信號和SDC的DOA估計值；

(3)采用式(17)的正交投影權(quán)矢量進行自適應(yīng)波束形成抑制SDC。

4 仿真實驗與分析

4.1 參數(shù)設(shè)置

4.2 仿真分析

對回波數(shù)據(jù)進行脈沖壓縮和常規(guī)DBF后得到的距離-多普勒譜如圖2所示，圖3是2250 km處的多普勒頻譜圖。從圖2和圖3可以看出位于2235~2265 km處的SDC能量較強，完全遮蓋了位于2250 km處的多普勒頻率為0.8 Hz的艦船目標，另外，由于陣列幅相誤差和期望信號的影響，噪聲基底抬高明顯，期望信號中的海雜波被削弱，這驗證了2.3節(jié)的分析結(jié)論。

圖2 常規(guī)處理后的距離-多普勒譜

圖4是分別采用NSF法與改進NSF法對陣列幅相誤差進行校正的結(jié)果圖，圖4(a)是采用常規(guī)MUSIC法和經(jīng)過Toeplitz預處理得到新的協(xié)方差矩陣后再進行MUSIC估計得到的DOA初值估計。由圖4(a)可見，陣列幅相誤差導致MUSIC譜分辨率嚴重下降且出現(xiàn)多個偽峰，使得估計的DOA初值與真實值差異較大，而且估計出的DOA個數(shù)多于真實DOA個數(shù)；采用Toeplitz預處理后的協(xié)方差矩陣估計的MUSIC譜性能得到很大提升，沒有出現(xiàn)偽峰，DOA初值估計也比較準確，分別為31.6°和44.0°。圖4(b)是分別采用NSF法與改進NSF法消除幅相誤差后估計的期望信號與SDC的DOA值，可以看出NSF法校正失效，采用改進NSF法能夠比較準確地估計出DOA值(經(jīng)過4次迭代后停止)，分別為29.7°和44.8°。

圖3 2250 km處的多普勒譜

圖4 NSF法與改進NSF法的校正結(jié)果比較

圖5 ADBF后的距離-多普勒譜()

圖6 ADBF后的距離-多普勒譜()

圖7 ADBF后的距離-多普勒譜()

圖8 ,和對應(yīng)的方向圖

表1不同仰角估計誤差時的信雜噪比改善量(dB)

仰角仰角估計誤差(°) 0.20.40.60.81.01.21.41.61.8 目標到達仰角18.017.617.216.916.616.416.116.015.8 SDC到達仰角17.616.415.514.813.912.912.111.811.6 目標、SDC到達仰角17.116.014.714.113.312.311.210.610.2

5 結(jié)論

電離層是OTHR進行目標探測時必須依靠的媒質(zhì)，電離層多層結(jié)構(gòu)和時變特性會產(chǎn)生多徑SDC，影響OTHR對慢速目標的檢測，如何有效抑制多徑SDC是OTHR面臨的難題之一。本文綜合考慮陣列存在幅相誤差和回波數(shù)據(jù)中含有強期望信號的實際情況，提出了一種抑制從接收陣列旁瓣進入的多徑SDC的方法。該方法采用改進NSF自校正法消除陣列幅相誤差，然后采用正交投影ADBF方法能夠在完整保留期望信號的同時較好地抑制多徑SDC，易于工程實現(xiàn)。理論分析和實驗仿真結(jié)果表明本文方法可有效抑制多徑SDC，改善慢速艦船目標檢測性能。

[1] Quan Y H, Xing M D, Zhang L,.. Transient interference excision and spectrum reconstruction for OTHR[J]., 2012, 48(1): 42-44.

[2] Anderson S J and Abramovich Y I. A unified approach to detection, classification, and correction of ionospheric distortion in HF sky wave radar systems[J]., 1998, 33(4): 1055-1067.

[3] Su H T, Liu H W, Shui P L,.. Adaptive beamforming for nonstationary HF interference cancellation in skywave Over-the-Horizon radar[J]., 2013, 49(1): 312-324.

[4] Br?ms M, Lundborg B, Jodalen V,.. Doppler effects on high latitude HF paths during an ionospheric disturbance[C]. Seventh International Conference on HF Radio Systems and Techniques, Sweden, 1997: 75-79.

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

Zhou Wen-yu and Jiao Pei-nan. Technology in Over The Horizon Radars[M]. Beijing: Publishing House of Electronics Industry, 2008: 276-277.

[6] Abramovich Y I, Frazer G J, and Johnson B A. Noncausal adaptive spatial clutter mitigation in monostatic MIMO radar: fundamental limitations[J]., 2010, 4(1): 40-54.

[7] Frazer G J. Simplified processing for time-staggered MIMO waveforms with interference reduction properties[C]. IEEE Radar Conference, Atlanta, GA, 2012: 801-805.

[8] 趙志國, 陳建文, 鮑拯. 分時波束空間MIMO天波雷達發(fā)射ADBF方法[J]. 華中科技大學學報(自然科學版), 2013, 41(2): 64-69.

Zhao Zhi-guo, Chen Jian-wen, and Bao Zheng. Transmit ADBF method for time diverse beam space MIMO-OTHR[J].(), 2013, 41(2): 64-69.

[9] 趙志國, 陳建文, 鮑拯. 一種高能量發(fā)射效率MIMO天波雷達方案[J]. 電子學報, 2012, 40(12): 2519-2524.

Zhao Zhi-guo, Chen Jian-wen, and Bao Zheng. A MIMO skywave radar scheme with high transmitting power efficiency[J]., 2012, 40(12): 2519-2524.

[10] Kazanci O, Bilik I, and Krolik J. Wavefront adaptive raymode processing for over-the-horizon HF radar clutter mitigation[C]. Signals, Systems and Computers Conference, Pacific, CA, 2007: 2191-2194.

[11] Bilik I, Kazanci O, and Krolik J. Radar clutter mitigation via space-time wavefront adaptive sensing[C]. IEEE Acoustics, Speech and Signal Processing, Las Vegas, USA, 2008: 2589-2592.

[12] Krolik J, Mecca V, Kazanci O,.. Multipath spread-Doppler clutter mitigation for over-the-horizon radar[C]. IEEE Radar Conference, Rome, 2008: 1-5.

[13] 廖桂生, 保錚, 張林讓. 基于特征結(jié)構(gòu)的自適應(yīng)波束形成新算法[J]. 電子學報, 1998, 26(3): 23-26.

Liao Gui-sheng, Bao Zheng, and Zhang Lin-rang. A new eigenstructure-based algorithm for adaptive beamforming [J]., 1998, 26(3): 23-26.

[14] Frazer G J, Abramovich Y I, and Johnson B A. Multiple-input multiple-output over-the-horizon radar: experimental results[J].&, 2009, 3(4): 290-303.

[15] 王永良, 丁前軍, 李榮鋒. 自適應(yīng)陣列處理[M]. 北京: 清華大學出版社, 2009: 66-67.

Wang Yong-liang, Ding Qian-jun, and Li Rong-feng. Adaptive Array Processing[M]. Beijing: Publishing House of Tsinghua University, 2009: 66-67.

[16] 田文龍, 李高鵬, 許榮慶. 利用自動識別系統(tǒng)信息進行高頻地波雷達天線陣校正[J]. 電子與信息學報, 2012, 34(5): 1065-1069.

Tian Wen-long, Li Gao-peng, and Xu Rong-qing. High frequency surface wave radar antenna array calibration using automatic identification system data[J].&, 2012, 34(5): 1065-1069.

[17] Friedlander B and Weiss A J. Direction finding in the presence of mutual coupling[J]., 1991, 39(3): 273-284.

[18] Chen Y M, Lee J H, Yeh C C,.. Bearing estimation without calibration for randomly perturbed arrays[J]., 1991, 39(1): 194-197.

羅 歡： 男，1986年生，博士生，研究方向為天波超視距雷達信號處理.

陳建文： 男，1964年生，教授，博士生導師，從事天波超視距雷達信號處理、陣列信號處理等的研究.

鮑 拯： 男，1977年生，博士，講師，從事天波超視距雷達信號處理、陣列信號處理等的研究.

A Multipath Spread Doppler Clutter Suppression Method for Skywave Over the Horizon Radar

Luo Huan①Chen Jian-wen②Bao Zheng②

①(,,430019)②(,,430019)

This paper analyzes the production mechanism of the multipath Spread Doppler Clutter (SDC) in skywave Over The Horizon Radar (OTHR). Because the amplitude and phase errors always exist in the array and the power of desired signal (sea clutter) is greater than the power of SDC, the SDC suppression ability declines through conventional Adaptive Digital Beam Forming (ADBF) process, at the same time, the desired signal is weakened and Signal-to-Noise Ratio (SNR) is seriously decreased. To resolve the above problems, an adaptive SDC suppression method is proposed. In this method, an improved Noise Subspace Fitting (NSF) method is utilized to eliminate amplitude and phase errors of OTHR array, and the direction of arrival angle of desired signal and SDC can be accurately obtained. Then, orthogonal projection weight vector is used in ADBF to solve the problem that conventional ADBF weight vector can not be accurately estimated because of the strong desired signal. Theoretical analysis and simulation results show that the scheme can completely suppress multipath SDC.

Skywave Over The Horizon Radar (OTHR); Multipath spread Doppler clutter; Amplitude and phase error correction; Adaptive Digital Beam Forming (ADBF)

TN958

A

1009-5896(2014)11-2691-07

10.3724/SP.J.1146.2013.01910

羅歡 luohuan5566@sina.com

2013-12-06收到，2014-03-31改回

國家自然科學基金(61072132)資助課題