賈 愚，賀 青，羅來源

（盲信號處理重點實驗室，成都　610041）

小孔徑天波超視距雷達多徑模糊污染消除及海雜波抑制*

賈 愚**，賀 青，羅來源

（盲信號處理重點實驗室，成都610041）

由于電離層的結(jié)構(gòu)分層性及厚度的不均勻，在超視距探測中會出現(xiàn)多徑模糊的現(xiàn)象，使同一個目標(biāo)出現(xiàn)在多個距離單元上，造成目標(biāo)識別的不準(zhǔn)確。針對傳統(tǒng)天波超視距雷達（OTHR）接收陣列不具備俯仰分辨能力的不足，利用小孔徑OTHR在俯仰方向的分辨能力，將幾何估計結(jié)果與測向結(jié)果相結(jié)合，在俯仰方向自適應(yīng)波束形成，消除了多徑模糊帶來的影響。同時改進了基于奇異值分解（SVD）的海雜波抑制算法抑制海雜波，提高了目標(biāo)信雜比以便于后續(xù)檢測。仿真實驗結(jié)果顯示所提方法能夠很好地消除多徑模糊帶來的影響，同時抑制了海雜波，信雜比提高了9 dB。

天波超視距雷達；多徑模糊；小孔徑；海雜波抑制；二維陣列；自適應(yīng)波束形成

1　引　言

天波超視距雷達（Over-The-Horizon Radar， OTHR）是一種利用電離層對高頻電磁波的反射作用進行超視距探測的大范圍警戒雷達。由于利用電離層反射探測，因此在某些工作頻率，電磁波的反射點位于不同電離層，造成相同目標(biāo)擁有多個不同路徑的回波，在距離維上形成多個偽峰。傳統(tǒng)天波超視距雷達通常利用返回散射探測設(shè)備、垂直探測設(shè)備、斜向探測設(shè)備對需要工作的區(qū)域進行聯(lián)合評估［1］選擇最佳的工作頻率，以避免多模和多徑模糊，但是當(dāng)電離層變化劇烈，無法選擇單一傳播模式時，多徑模糊難以避免。

近年來，國內(nèi)外學(xué)者對OTHR二維陣列做了許多研究∶文獻［2-3］介紹和分析了法國的二維超視距雷達系統(tǒng)NOSTRADAMUS；文獻［4］對未來天波超視距雷達的發(fā)展進行了展望；文獻［5］介紹了二維OTHR的模式選擇方法；文獻［6］對二維陣列的優(yōu)化設(shè)計進行了研究討論；文獻［7］研究了利用二維陣列的新的選頻方法。小孔徑OTHR是一種新體制的二維OTHR，它利用小孔徑圓陣作為其接收陣［8］。針對多徑模糊問題，可以利用文獻［7］的方法進行最佳工作頻率的自適應(yīng)選擇，另外還可以通過在信號域?qū)夭〝?shù)據(jù)進行信號處理，減輕多徑模糊帶來的影響。本文正是從這一角度進行研究，通過粗細(xì)兩次俯仰角估計，并在此基礎(chǔ)上做自適應(yīng)波束形成，以消除多徑模糊對目標(biāo)探測的影響。為了進一步提高信雜比，以便后續(xù)檢測，本文改進了基于奇異值分解（Singular Value Decomposition，SVD）的海雜波抑制的方法，在幅度和峰值進行二維搜索，自適應(yīng)地選擇奇異值，更好地消除了海雜波對信號的影響。

2　問題描述

2.1多徑模糊原理

電離層按電子濃度的高度變化可分為D、E、F層，其中天波超視距雷達主要利用E、F層的反射作用進行探測。E層高度一般在90～140 km之間，F(xiàn)層在140 km以上［1］。如圖1所示，發(fā)射信號由線陣發(fā)出探測到目標(biāo)后發(fā)生反射，其中一部分回波經(jīng)過電離層反射被小孔徑圓陣接收陣所獲取。但是由于電離層各區(qū)域濃度的差異，因此有一部分回波經(jīng)過路徑A-B-C沿在E層反射被接收，另一部分回波經(jīng)過路徑A-D-C沿F層反射被接收，由于回波時長的不同，因此同一目標(biāo)會被判定為不同距離上的兩個目標(biāo)，使目標(biāo)檢測造成困難。

圖1　多徑模糊原理示意圖Fig.1 Illustration of multipath fuzzy

由圖1所示，通過不同電離層路徑傳播的回波在到達接收陣時，其回波在俯仰上到達角（Direction of Arrival，DOA）是明顯不同的，本文正是從俯仰DOA的差異出發(fā)，在俯仰方向上采用接收波束自適應(yīng)空域處理來解決回波的多徑模糊問題。

對于傳統(tǒng)OTHR，由于采用線陣作為接收陣，無法獲取俯仰方向上的信息，因此，面對多徑模糊的問題時，處理手段比較被動，通常是利用探測設(shè)備獲取電離層的狀態(tài)參數(shù)，根據(jù)電離層參數(shù)來選擇最佳的工作頻率，使雷達工作于單一傳播模式，從而避免多徑模糊帶來的干擾。小孔徑OTHR采用小孔徑圓陣作為接收陣，不僅可以減少資源消耗，降低系統(tǒng)復(fù)雜度，還能獲取回波俯仰方向上的信息，用于主動地消除多徑模糊。

2.2信號模型

小孔徑OTHR采用圓陣作為接收陣［8］。設(shè)圓陣半徑為R，陣元數(shù)為N，OTHR工作波長為λ，回波信號俯仰角為θ，方位角為φ。

在時刻t，陣列接收信號可以表示為［9］

式中∶p表示不同路徑的回波；γp是不同路徑衰減因子；a（φp，θp）是圓陣導(dǎo)向矢量，

sp（t）=［sp1（t），sp2（t）…spN（t）］T；n（t）為噪聲矢量。假設(shè)（φ0，θ0）為期望信號的來波方向，則其余方向回波是通過不同路徑來波方向的回波。sp（t）中spN為第N個陣元接收到的p路徑的回波復(fù)包絡(luò)∶

式中∶fd為目標(biāo)多普勒頻率；fb為一階海雜波Bragg頻率，它與工作頻率有關(guān)。在OTHR多徑模糊問題中，若考慮只有E、F兩層傳播，發(fā)射模式為單模，那么回波路徑只有A-B-C和A-D-C兩種，兩者回波在方位上角度相同，因此區(qū)別只存在于俯仰角θk和回波復(fù)包絡(luò)sp（t）。

3　多徑模糊消除及海雜波抑制算法

3.1多徑模糊消除算法原理

由2.2節(jié)討論可知，由E、F兩層造成的多徑模糊可以簡化為A-B-C和A-D-C量路徑的回波，兩者回波方向在方位方向相同、俯仰上差異明顯，因此接收信號可以表示為

如果在沒有進行模式識別的前提下，將獲得的數(shù)據(jù)假定為E層傳播，進行常規(guī)處理，由于有F層的目標(biāo)回波，因此，相同目標(biāo)會在距離維上形成兩個回波峰。此時結(jié)合垂直探測設(shè)備和斜向探測設(shè)備，可以獲得當(dāng)前電離層的等效高度及可能傳播的模式。由等效高度、傳播模式，結(jié)合目標(biāo)回波峰的位置可以得到估計俯仰角度θf、θe，其幾何關(guān)系如圖2所示。

圖2　粗估計仰角幾何關(guān)系示意圖Fig.2 Geometrical relationship of pitch angle in rough estimation

式中∶OC、PC由電離層垂測設(shè)備獲得；PA、OA可由傳播時延進行估算。將所得數(shù)據(jù)做二維DOA估計，可以獲得比較精確的DOA估計值（φ，θ1）、（φ，θ2）。由于E、F兩層傳播的仰角差異明顯，將DOA估計值與粗估計值θf進行對比，如果θ1≈θf，即可判斷θ1為F層傳播過來的多徑污染。將多徑傳播的回波當(dāng)作干擾信號，再將式（6）的接收數(shù)據(jù)做自適應(yīng)波束形成，即可消除多徑污染。

二維DOA估計采用二維MUSIC方法，它是一維MUSIC的擴展。根據(jù)式（6），寫為矩陣形式有式中∶H表示共軛轉(zhuǎn)置，對R進行特征分解。假設(shè)有m個目標(biāo)，則其中較小的N-m個特征值張成噪聲空間EN，m個較大的特征值張成目標(biāo)子空間Es，由于噪聲和目標(biāo)子空間是正交的，在θ、φ上進行二維搜索，則有譜函數(shù)∶

進行譜峰搜索得到估計值（φ，θ1）、（φ，θ2），與粗估計值θf進行對比，如若θ1≈θf，則將（φ，θ1）作為干擾，（φ，θ2）當(dāng)作期望信號做MVDR波束形成［10］，即

利用Lagrange乘子法求解這一優(yōu)化問題，可得自適應(yīng)權(quán)矢量為

多徑模糊消除后的陣列輸出為

總結(jié)消除多徑模糊的算法流程如圖3所示。

圖3　多徑模糊消除算法流程Fig.3 The algorithm flow of multipath fuzzy contamination elimination

3.2改進型的海雜波抑制算法

經(jīng)過多徑模糊消除處理后，由于海雜波的一階Bragg峰能量較大，會給目標(biāo)檢測帶來困難，經(jīng)過電離層相位污染的海雜波譜甚至?xí)箤捯灾卵蜎]速度較慢的艦船目標(biāo)，因此有必要對海雜波進行抑制。

一階Bragg峰可以看成是兩個單頻的正弦信號，經(jīng)過快時間匹配濾波后的目標(biāo)信號也可視為單頻的正弦信號，因此可以利用SVD分解，將信號的不同頻率區(qū)分開。

經(jīng)過快時間匹配濾波的信號表示為y1（m），m= 1，2，…，M，將其構(gòu)造Hankel矩陣H（y）∶

式中∶M為數(shù)據(jù)長度；c由信號中正弦信號數(shù)目r決定，通常取c=3r+1。

對H進行奇異值分解∶

式中∶S為H的奇異值，S=diag［σ1，σ2，…，σc］，且σ1≥σ2≥…≥σc≥0；U、V為對應(yīng)奇異值的左右奇異矢量組成的矩陣。較大的r個奇異值對應(yīng)r個正弦分量，剩下c-r個較小的奇異值對應(yīng)小的雜波和噪聲分量。因此，傳統(tǒng)方法通過將前兩個對應(yīng)Bragg分量的奇異值置零，再構(gòu)造回波信號進行雜波抑制。構(gòu)造新的Hankel矩陣如下∶

式中∶S1=diag［0，0，σ3，…，σc］。由H1重構(gòu)時域信號∶

式中∶m=i+j-1；n為符合要求的元素；y＇（m）為雜波抑制后的時域信號。

有時由于受到電離層污染的影響，一部分海雜波能量擴展到信號能量空間，另外為了更好地消除噪聲的影響，從而進一步提高信雜噪比，本文改進了SVD雜波抑制的方法。由于信號的頻率信息保留在U中，對U的列向量做快速傅里葉變換（Fast Fourier Transform，F(xiàn)FT）可以得到對應(yīng)奇異值的信號分量。本文在將奇異值置零時，改進了判決準(zhǔn)則，做進一步篩選。具體步驟如下∶

（1）對U的每一列做FFT得到矩陣UFFT；

（2）目標(biāo)信號頻率和Bragg頻率已經(jīng)是已知的先驗信息，因此可以快速在fd處進行譜峰搜索，如果有峰值則記搜得的峰值大小為v（d）；

（3）求頻譜的平均幅度記為ˉv；

（4）求v（d）/ˉv與門限δ做比較，如果大于門限則可判斷為信號分量，小于門限則可能是噪聲或毛刺；

（5）將判斷為信號分量的對應(yīng)奇異值保留，其余分量全部置零，得到S2重構(gòu)數(shù)據(jù)矢量。

4　仿真實驗與分析

參考國內(nèi)某次小孔徑OTHR實驗參數(shù)，OTHR系統(tǒng)參數(shù)為∶接收陣元N=80，半徑R=190 m，信號體制采用線性調(diào)頻信號，工作頻率f=15 MHz，信號采用線性調(diào)頻信號，信號帶寬為B=10 kHz，脈寬τ =20 ms，每個周期采400個點，積累256個周期。一階Bragg峰的頻率fb=0.395 Hz，目標(biāo)多普勒頻率fd=0.8 Hz，假設(shè)E層傳播層高為120 km，F(xiàn)層傳播層高為254 km，E層傳播目標(biāo)在第94個距離單元，F(xiàn)層傳播目標(biāo)在第107個距離單元，即分別位于705 km處和802.5 km處；通過估算，E層回波仰角為60°，F(xiàn)層回波仰角為50°，方位角均為30°。

在未進行多徑模糊消除時，常規(guī)處理后得到距離-多普勒圖如圖4所示?？梢钥吹接捎陔婋x層分層的影響，通過E、F層傳播的回波在距離軸上形成兩個目標(biāo)峰。將回波數(shù)據(jù)做二維MUSIC，得到結(jié)果如圖5所示。通過譜峰搜索，可以得到比較精確的回波仰角，實驗測得仰角為49.5°和60.5°。

圖4　多徑下距離多普勒圖Fig.4 Range-Doppler diagram of multipath

圖5　二維MUSIC測向圖Fig.5 Direction finding diagram of 2-D MUSIC

利用仰角信息得到期望信號俯仰方向，做MVDR自適應(yīng)波束形成，如圖6所示，可以看到在49.5°左右形成了零陷，而在期望的60°方向形成波束，整體旁瓣比常規(guī)波束形成的稍低。

圖6　MVDR方向圖Fig.6 Directional diagram of MVDR

利用式（15）得到自適應(yīng)波束形成后效果如圖7所示，可以看到在802.5 km處的偽目標(biāo)已經(jīng)被抑制了，此時只剩下705 km處的真實目標(biāo)，但是過高的海雜波仍然影響檢測，為此，根據(jù)3.2節(jié)的海雜波抑制算法進行雜波抑制。

圖7　多徑污染消除后的距離多普勒圖Fig.7 Range-Doppler diagram after multipath fuzzy contamination elimination

如圖8所示，可以看到改進型的SVD算法由于準(zhǔn)確地篩選了對應(yīng)的奇異值，把其他雜波、噪聲分量都抑制了，因此和傳統(tǒng)的SVD算法相比，雜波抑制效果更好，噪底更低，噪聲部分更加平滑。

圖8　改進后SVD算法和傳統(tǒng)SVD算法對比Fig.8 Contrast of improved SVD algorithm withconventional SVD algorithm

圖9是經(jīng)過雜波抑制后的距離-多普勒圖，雖然仍有部分一階Bragg峰的余量，但是目標(biāo)已經(jīng)可以很清楚地分辨出來。

圖9　雜波抑制后的距離-多普勒圖Fig.9 Range-Doppler diagram after sea clutter suppression

圖10是針對多目標(biāo)單干擾情況的仿真。在多普勒頻率fd1=0.8 Hz處由于回波角度剛好達到E、F層傳播的條件，產(chǎn)生了一個多徑干擾，位置在802.5 km附近。在與偽峰相同距離處有一多普勒fd2=2.5 Hz的快速目標(biāo)，距離多普勒如圖10（a）所示。仿真設(shè)定fd2=2.5 Hz的快速目標(biāo)也是通過E層傳播，但是并沒有產(chǎn)生多徑干擾，在抑制F層慢目標(biāo)干擾時，快速目標(biāo)會被保留下來。從圖10（b）可以看到，E層快速目標(biāo)經(jīng)過本文算法仍被很好地抑制掉，這與理論推導(dǎo)是一致的。圖10（c）是抑制海雜波以后的距離-多普勒圖，由結(jié)果可以看到本文算法對于多目標(biāo)情況下仍然適用。

圖10　多目標(biāo)單干擾距離多普勒Fig.10 Range-Doppler diagram of multiple target with one interference

圖11是針對多目標(biāo)多干擾情況的仿真，目標(biāo)都位于705 km處，一個快速目標(biāo)，一個低速目標(biāo)，兩目標(biāo)都在F層產(chǎn)生了多徑干擾，在經(jīng)過本文方法處理后也獲得了比較理想的效果。

圖11　多目標(biāo)多干擾距離-多普勒Fig.11 Range-Doppler diagram of multiple target with multiple interference

5　結(jié)束語

本文利用小孔徑OTHR接收陣在俯仰角上的分辨能力，在先驗粗估計信息的基礎(chǔ)上，在俯仰維上進行自適應(yīng)波束形成，很好地抑制了多徑模糊干擾。與傳統(tǒng)天波超視距雷達相比，提供了一種利用波束形成直接消除多徑模糊的方法，經(jīng)過仿真驗證有較好的效果，證明了算法的有效性。另外，為了進一步提高信雜比，本文通過合理篩選與期望信號對應(yīng)的奇異值，改進了SVD海雜波抑制法，理論分析與實驗結(jié)果表明改進的SVD海雜波抑制法能夠較好地抑制海雜波，并將噪底降低，信雜比提高了9 dB，改善了目標(biāo)檢測性能。

然而本文算法較為依賴估計得到的俯仰角度信息，如果粗估計角度信息有誤，在此基礎(chǔ)上所做的自適應(yīng)波束形成必然會產(chǎn)生誤差。因此后續(xù)工作中會將高分辨率二維DOA估計、穩(wěn)健波束形成技術(shù)納入到研究中來，以提高算法的適應(yīng)能力。

［1］ 周文瑜，焦培南.超視距雷達技術(shù)［M］.北京∶電子工業(yè)出版社，2008.

ZHOU Wenyu，JIAO Peinan.Over the horizon radar technology［M］.Beijing∶Publishing House of Electronics Industry，2008.（in Chinese）

［2］ SAILLANT S，AUFFRAY G，DOREY P.Exploitation of elevation angle control for a 2-D hf skywave radar［C］// Proceedings of 2003 International Radar Conference.Adelaide，Australia∶IEEE，2003∶662-666.

［3］ BAZIN V，MOLINIE J P，MUNOZJ，et al.A general presentation about the OTH-Radar NOSTRADAMUS［C］// Proceedings of 2006 IEEE International Radar Conference.Machester，UK∶IEEE，2006∶634-642.

［4］ 吳瑕，陳建文，鮑拯，等.新體制天波超視距雷達的發(fā)展與研究［J］.宇航學(xué)報，2013，34（5）∶671-678.

WU Xia，CHEN Jianwen，BAO Zheng，et al.Development and research on the new system sky-wave over-the-horizon radar［J］.Journal of Astronautics，2013，34（5）∶671 -678.（in Chinese）

［5］ ABRAMOVICH Y I，F(xiàn)RAZERM G J，JOHNSON B A. Principles of mode-selective MIMO OTHR［J］.IEEE Transactions on Aerospace and Electronic Systems，2013，49（3）∶1839-1868.

［6］ 嚴(yán)韜，陳建文，鮑拯.基于改進遺傳算法的天波超視距雷達二維陣列稀疏優(yōu)化設(shè)計［J］.電子與信息學(xué)報，2014，36（12）∶3014-3020.

YAN tao，CHEN Jianwen，BAO Zheng，et al.Optimization design of sparse 2-D array for over-the-horizon radar（OTHR）based on improved genetic algorithm［J］.Journal of Electronics&Information Technology，2014，36（12）∶3014-3020.（in Chinese）

［7］ 羅忠濤，何子述，盧琨，等.二維陣列超視距雷達的選頻方法［J］.電子與信息學(xué)報，2015，37（4）∶975-981.

LUO Zhongtao，HE Zishu，LU Kun，et al.Operating frequency selection for sky-wave over-the-horizon radar with 2-D array［J］.Journal of Electronics&Information Technology，2015，37（4）∶975-981.（in Chinese）

［8］ HE Q，LUO L Y.Analysis of over-the-horizon target detection based on the small aperture array［C］//Proceedings of the 2014 7th International Congress on Image and Signal Processing.Dalian∶IEEE，2014∶1100-1104.

［9］ DOU D X，LI M，HE Z S.Multi-mode clutter suppression of multiple-input-multiple-output over-the-horizon radar based on blind source separation［J］.IET Radar Sonar&Navigation，2015，9（8）∶956-966.

［10］ 張賢達.現(xiàn)代信號處理［M］.2版.北京∶清華大學(xué)出版社，2002∶126-134.

ZHANG Xianda.Modern Signal Processing［M］.2nd ed.Beijing∶Tsinghua University Press，2002∶126-134.（in Chinese）

賈 愚（1990—），男，四川西昌人，2013年于電子科技大學(xué)獲工學(xué)學(xué)士學(xué)位，現(xiàn)為盲信號處理重點實驗室碩士研究生，主要研究方向為雷達信號處理和陣列信號處理；

JIA Yu was born in Xichang，Sichuan Province，in 1990.He received the B.S.degree from University of Electronic Science and Technology of China in 2013.He is now a graduate student.His research concerns radar signal processing and array signal processing.

Email∶jiayu_9012@163.com

賀 青（1984—），男，四川成都人，博士研究生，主要研究方向為雷達信號處理和陣列信號處理；

HE Qing was born in Chengdu，Sichuan Province，in 1984. He is currently working toward the Ph.D.degree.His research concerns radar signal processing and array signal processing.

羅來源（1956—），男，湖南長沙人，博士，高級工程師，主要研究方向為系統(tǒng)總體技術(shù)、短波信號處理。

LUO Laiyuan was born in Changsha，Hunan Province，in 1956.He is now a senior engineer with the Ph.D.degree.His research concerns system design and HF signal processing.

Multipath Fuzzy Contamination Elimination and Sea Clutter Suppression for Small Aperture Sky-wave Over-the-horizon Radars

JIA Yu，HE Qing，LUO Laiyuan
（National Key Laboratory of Science and Technology on Blind Signal Processing，Chengdu 610041，China）

∶Multipath propagation，caused by the layered structure of ionosphere，makes one target appear in different range cells，which leads to an inaccurate result of target identification.The receiving array of conventional sky-wave over-the-horizon radar（OTHR）is uniform linear array which doesn＇t have the resolution capability of elevation direction while small aperture OTHR system has it，because of its uniform circular receiving array.To solve the problem of multipath propagation，this paper proposes to make adaptive beamforming in elevation direction with the result of geometric estimation and spectrum estimation under small aperture system.Moreover，to develop signal-to-clutter ratio（SCR）of echo signal，it improves the ocean clutter suppression algorithm that is based on singular value decomposition（SVD）.Theoretical analysis and simulation result prove the proposed algorithm can eliminate multipath effect efficiently and avoid false targets.Besides，the improved SVD suppression algorithm increases the quality of ocean clutter suppression.

∶sky-wave over-the-horizon（OTHR）radar；multipath fuzzy；small aperture；sea clutter suppression；2-D array；adaptive beamforming

TN958.93

A

1001-893X（2016）05-0568-07

10.3969/j.issn.1001-893x.2016.05.016

賈愚，賀青，羅來源.小孔徑天波超視距雷達多徑模糊污染消除及海雜波抑制［J］.電訊技術(shù)，2016，56（5）∶568-574.［JIA Yu，HE Qing，LUO Laiyuan.Multipath fuzzy contamination elimination and sea clutter suppression for small aperture sky-wave over-the-horizon radars［J］.Telecommunication Engineering，2016，56（5）∶568-574.］

2015-10-08；

2016-01-12Received date：2015-10-08；Revised date：2016-01-12

**通信作者：jiayu_9012@163.comCorresponding author：jiayu_9012@163.com