胡進(jìn)峰，李萬閣，艾　慧，李會(huì)勇，李　軍，夏　威，徐　威

（1.電子科技大學(xué)電子工程學(xué)院，四川　成都 611731；2.空軍大連通信士官學(xué)校，遼寧　大連 116600）

基于頻域干擾剔除的OTHR射頻干擾抑制算法

胡進(jìn)峰1，李萬閣1，艾慧1，李會(huì)勇1，李軍1，夏威1，徐威2

（1.電子科技大學(xué)電子工程學(xué)院，四川成都 611731；2.空軍大連通信士官學(xué)校，遼寧大連 116600）

目前，天波超視距雷達(dá)（over-the-horizon radar，OTHR）大多是在空域或匹配濾波后的時(shí)頻域抑制射頻干擾，往往會(huì)同時(shí)抑制掉射頻干擾附近的目標(biāo)?？紤]到在匹配濾波后，射頻干擾與目標(biāo)信號(hào)都是近似單頻信號(hào)難以區(qū)分；而在匹配濾波前，目標(biāo)信號(hào)是寬帶的線性調(diào)頻信號(hào)，而射頻干擾是近似單頻的窄帶信號(hào)，射頻干擾與目標(biāo)信號(hào)顯著不同，易于分離?；谠撎攸c(diǎn)，提出了一種基于匹配濾波前的頻域干擾置零的射頻干擾抑制算法。該算法在頻域中檢測(cè)射頻干擾的頻點(diǎn)位置并將該頻點(diǎn)的信號(hào)成分置零，最后變換到時(shí)域?qū)崿F(xiàn)射頻干擾抑制。與傳統(tǒng)方法相比，本文所提算法在抑制射頻干擾時(shí)基本不損害目標(biāo)信號(hào)。實(shí)測(cè)數(shù)據(jù)的處理結(jié)果驗(yàn)證了所提方法的優(yōu)越性。

天波超視距雷達(dá)；射頻干擾；抑制

網(wǎng)址：www.sys-ele.com

0　引　言

射頻干擾抑制是天波雷達(dá)的關(guān)鍵問題之一［1-6］。天波超視距雷達(dá)（over-the-horizon radar，OTHR）工作在3～30 MHz的高頻波段，該波段內(nèi)充滿了各種短波通信和廣播信號(hào)等有源的射頻干擾［7-13］，這些射頻干擾進(jìn)入雷達(dá)接收機(jī)經(jīng)過混頻和雷達(dá)回波信號(hào)混在一起，影響目標(biāo)信號(hào)的檢測(cè)和參數(shù)估計(jì)。

射頻干擾的主要成分為確定性載頻信號(hào)，其通常為單頻信號(hào)且能量較強(qiáng)，在進(jìn)入雷達(dá)接收機(jī)經(jīng)過混頻（去雷達(dá)載頻）后仍然為一單頻信號(hào)。經(jīng)過匹配濾波和多普勒處理后，其在距離多普勒譜中以平行于距離軸的脊背形式存在于某些特定的多普勒單元上，會(huì)遮蔽位于這些多普勒頻率單元上的目標(biāo)，造成目標(biāo)漏檢，影響雷達(dá)系統(tǒng)性能。

目前的射頻干擾抑制方法主要有兩大類：一類是采用自適應(yīng)波束形成方法［14-15］在波束形成時(shí)對(duì)射頻干擾進(jìn)行抑制，此類算法的有效性主要依賴于射頻干擾的空間特性。然而受天波超視距雷達(dá)復(fù)雜工作環(huán)境的影響，其射頻干擾通常在空間上表現(xiàn)出較強(qiáng)的非平穩(wěn)性，此外受天線陣列的限制，該方法不能有效抑制主瓣內(nèi)的射頻干擾和存在范圍過大的射頻干擾；另一類是在匹配濾波之后的距離多普勒域中利用射頻干擾在距離維上的強(qiáng)相關(guān)性，通過特征值分解［16］、子空間投影［17］或平均距離譜法［18］實(shí)現(xiàn)射頻干擾抑制。由于匹配濾波后，目標(biāo)信號(hào)與射頻干擾等都是近似單頻的窄帶信號(hào)，射頻干擾與目標(biāo)回波的可分性較差，因此這類方法在抑制射頻干擾的同時(shí)通常會(huì)對(duì)目標(biāo)造成不可恢復(fù)的影響，尤其是對(duì)于射頻干擾所在多普勒頻率單元上的目標(biāo)，此類方法大多不能保證射頻干擾所在多普勒頻率單元上目標(biāo)的安全性。

針對(duì)這些問題，本文對(duì)天波超視距雷達(dá)中射頻干擾的特性進(jìn)行了研究，發(fā)現(xiàn)在匹配濾波前，射頻干擾和雷達(dá)回波信號(hào)更具有可分性：射頻干擾為窄帶信號(hào)甚至單頻信號(hào)，而天波雷達(dá)回波信號(hào)為具有一定帶寬的線性調(diào)頻信號(hào)?；诖?，本文提出了一種頻域干擾置零的射頻干擾抑制算法，首先將雷達(dá)接收信號(hào)進(jìn)行多普勒處理（進(jìn)行折疊變換到快慢時(shí)間域并對(duì)慢時(shí)間域做快速傅里葉變換（fast Fourier transform，F(xiàn)FT）），在多普勒域中進(jìn)行簡(jiǎn)單的海雜波抑制并做逆多普勒處理變換到時(shí)域，然后對(duì)抑制海雜波后的雷達(dá)接收數(shù)據(jù)做FFT變換到頻域，并設(shè)置門限檢測(cè)射頻干擾位置，最后在原始雷達(dá)接收數(shù)據(jù)的頻域中將所檢測(cè)到的射頻干擾頻點(diǎn)抑制掉，即可實(shí)現(xiàn)射頻干擾抑制。該方法充分考慮射頻干擾在雷達(dá)接收前端為單頻信號(hào)的特點(diǎn)，將天波雷達(dá)中射頻干擾抑制問題轉(zhuǎn)換為一個(gè)寬帶信號(hào)中的窄帶干擾抑制問題，并采用檢測(cè)剔除的思想對(duì)其進(jìn)行抑制，在抑制射頻干擾的同時(shí)最大限度地保護(hù)了目標(biāo)不受影響。實(shí)測(cè)數(shù)據(jù)的處理結(jié)果表明，相對(duì)于傳統(tǒng)的平均距離譜法，本文所提算法可以在盡量不破壞目標(biāo)信號(hào)的前提下很好地抑制射頻干擾。

1　問題分析

本節(jié)將主要對(duì)射頻干擾的來源、特性和雷達(dá)回波的特性進(jìn)行對(duì)比，分析它們之間的差別，為后續(xù)的射頻干擾抑制提供理論基礎(chǔ)。

1.1射頻干擾特性分析

天波超視距雷達(dá)中的射頻干擾主要是由高頻段的短波通訊信號(hào)和廣播電臺(tái)信號(hào)引起的，也被稱為電臺(tái)干擾。在天波雷達(dá)工作的高頻段中充滿了各種廣播電臺(tái)和通信信號(hào)，所以射頻干擾出現(xiàn)得比較頻繁。電臺(tái)信號(hào)的調(diào)制方式有很多，但是大量的實(shí)測(cè)數(shù)據(jù)處理表明，雷達(dá)接收機(jī)中接收到的主要干擾成分是信號(hào)的載波。載波信號(hào)成分是確定的，其在時(shí)域中的相關(guān)性非常強(qiáng)，在信號(hào)處理后期的距離多普勒譜中具有明顯的特征；而電臺(tái)干擾中的隨機(jī)信號(hào)成分在時(shí)域中的相關(guān)性比較弱，對(duì)雷達(dá)回波中的有用信號(hào)影響相對(duì)較小。

為了方便進(jìn)行分析，利用式（1）所示的諧波信號(hào)模型對(duì)射頻干擾進(jìn)行建模

式中，fc為雷達(dá)的工作頻率；ak為第k個(gè)射頻干擾的幅度；fk為第k個(gè)射頻干擾相對(duì)于fc的頻率偏差；φk為第k個(gè)射頻干擾的初始相位；K為射頻干擾的個(gè)數(shù)。

圖1　天波雷達(dá)信號(hào)處理流程

如圖1所示，射頻干擾被雷達(dá)接收機(jī)接收，經(jīng)過混頻處理去掉載頻后的射頻干擾可以表示為

將式（2）用快慢時(shí)間形式進(jìn)行表示，可得

式中，n為掃頻周期序列（慢時(shí)間）；t為某個(gè)脈沖周期內(nèi)的采樣序列（快時(shí)間）；T為雷達(dá)掃頻周期。

式（3）中的射頻干擾經(jīng)過匹配濾波器后，可以得到

式中，mfk（τ）=MF｛exp（j2πfkt）｝；τ為匹配濾波器輸出對(duì)應(yīng)的回波延遲，乘以c/2即可以得到對(duì)應(yīng)的距離；MF｛·｝為匹配濾波算子。

由式（4）可以看出，經(jīng)過匹配濾波器后，射頻干擾在慢時(shí)間域仍然為一單頻信號(hào)。

對(duì)式（4）中的射頻干擾做多普勒處理（對(duì)每個(gè)距離單元中的慢時(shí)間域信號(hào)做FFT）可得

式中，F(xiàn)FT｛·｝指的是快速傅里葉變換算子。由式（5）可以看出，在距離多普勒處理后，射頻干擾存在于某個(gè)特定多普勒頻率的所有距離單元上。

由于射頻干擾的頻偏fk可能超出多普勒頻率觀測(cè)范圍，所以射頻干擾最終在多普勒域中的頻率位置為

式中，fr=1/T為脈沖掃描頻率；T為脈沖重復(fù)周期；mod（·）為求模算子。

由本節(jié)分析可知，射頻干擾在進(jìn)入雷達(dá)接收機(jī)后為一單頻信號(hào)，經(jīng)過匹配濾波處理后其在慢時(shí)間域中仍為一單頻信號(hào)，經(jīng)過多普勒處理后，其在距離多普勒譜中具有明顯的特征：具有很強(qiáng)的距離相關(guān)性，表現(xiàn)為平行于距離維的脊背干擾。

1.2雷達(dá)回波信號(hào)模型

OTHR的發(fā)射信號(hào)通常為線性調(diào)頻信號(hào)，可以表示為

式中，A為發(fā)射信號(hào)的幅度；fc為雷達(dá)工作頻率；T為發(fā)射信號(hào)的重復(fù)周期；N為發(fā)射信號(hào)的脈沖數(shù)；α=B/T為信號(hào)的調(diào)頻率；B為信號(hào)的帶寬。

發(fā)射信號(hào)在遇到目標(biāo)或者海面反射回的雷達(dá)回波，可以看作是發(fā)射信號(hào)在幅度上衰減并增加了時(shí)延和多普勒頻移的信號(hào)。此類信號(hào)是雷達(dá)發(fā)射信號(hào)的反射回波，統(tǒng)稱之為雷達(dá)回波，其可以表示為

式中，sci（t）為第i個(gè)海面的回波信號(hào)；ηci為第i個(gè)海面回波的反射系數(shù)；τci為第i個(gè)海面回波的時(shí)延（其隨時(shí)間t變化，包含了海雜波的多普勒頻移信息）；I為海面回波的總個(gè)數(shù)；stj（t）為第j個(gè)目標(biāo)的回波；ηtj為第j個(gè)目標(biāo)回波的反射系數(shù)；τtj為第j個(gè)目標(biāo)回波的時(shí)延（其隨時(shí)間t變化，包含了目標(biāo)的多普勒頻移信息）；J為目標(biāo)的總個(gè)數(shù)。

此處，將海面回波看作多個(gè)點(diǎn)目標(biāo)的回波綜合。

由式（8）可以看出，雷達(dá)回波是發(fā)射信號(hào)取不同時(shí)延后的加權(quán)疊加，其仍然是具有帶寬B的線性調(diào)頻信號(hào)集。

在匹配濾波后，目標(biāo)信號(hào)和海雜波信號(hào)在慢時(shí)間域中也表現(xiàn)為較單一的單頻信號(hào)（此處特指常規(guī)勻速目標(biāo)，不考慮機(jī)動(dòng)目標(biāo)，機(jī)動(dòng)目標(biāo)的頻譜會(huì)有所展寬）。

由本節(jié)分析可知，在匹配濾波后的多普勒域中，海雜波和目標(biāo)回波表現(xiàn)為一系列單頻信號(hào)的組合，和射頻干擾難以區(qū)分；而在匹配濾波前，雷達(dá)回波信號(hào)中的海雜波和目標(biāo)回波均為具有一定帶寬的線性調(diào)頻信號(hào)，它們?cè)陬l率模式上和射頻干擾具有明顯的差別。下文將基于此，提出相應(yīng)的射頻干擾抑制算法。

2　所提射頻干擾抑制算法

通過第1節(jié)的分析可知，在進(jìn)入雷達(dá)接收機(jī)時(shí)射頻干擾為一近乎單一的單頻信號(hào)，而目標(biāo)回波和海雜波則均為有一定帶寬的線性調(diào)頻信號(hào)。針對(duì)它們之間這種明顯的差異，本節(jié)提出一種基于頻域干擾剔除的射頻干擾抑制算法。該算法首先在多普勒域中抑制海雜波并利用抑制海雜波后的數(shù)據(jù)檢測(cè)射頻干擾的頻點(diǎn)，然后在雷達(dá)原始接收信號(hào)中將檢測(cè)到的頻點(diǎn)置零，實(shí)現(xiàn)射頻干擾的抑制。該方法在匹配濾波前利用射頻干擾和雷達(dá)回波在頻域中的明顯差別抑制射頻干擾，可以盡量避免對(duì)目標(biāo)信號(hào)的破壞。

天波雷達(dá)中進(jìn)入雷達(dá)接收機(jī)的信號(hào)可以表示為

式中，sc（t）為海面回波；st（t）為目標(biāo)回波；j（t）為進(jìn)入接收機(jī)的射頻干擾；w（t）為接收機(jī)外部噪聲。

由于天波雷達(dá)的照射環(huán)境多為廣闊的海面，所以通常情況下雷達(dá)接收信號(hào)中海面回波的能量很強(qiáng)，如圖2（a）所示。此時(shí)，無法直接在頻譜中檢測(cè)出射頻干擾，因此在檢測(cè)射頻干擾前需要對(duì)接收信號(hào)中的海雜波進(jìn)行簡(jiǎn)單地抑制。

圖2　雷達(dá)接收信號(hào)的頻譜分析

由海雜波的產(chǎn)生機(jī)理［19-20］和天波雷達(dá)的工作原理可知，海雜波在多普勒域中存在于零頻附近的低頻區(qū)域中，因此可以通過簡(jiǎn)單的多普勒頻域置零方法對(duì)海雜波進(jìn)行抑制。

如圖3所示，先將雷達(dá)接收信號(hào)變換到多普勒域中對(duì)海雜波進(jìn)行抑制，再通過逆多普勒處理將信號(hào)恢復(fù)到時(shí)域中去，這樣接收信號(hào)中的強(qiáng)海雜波就得到了有效地抑制。

圖3　海雜波抑制流程圖

此時(shí)，如圖2（b）所示，在剩余的接收信號(hào)中射頻干擾相對(duì)于目標(biāo)信號(hào)和噪聲具有一定的能量?jī)?yōu)勢(shì)。將s′（t）利用FFT變換到頻域中，設(shè)置適當(dāng)?shù)臋z測(cè)準(zhǔn)則即可有效地檢測(cè)出射頻干擾所在頻點(diǎn)。

下面對(duì)本節(jié)中射頻干擾的檢測(cè)準(zhǔn)則進(jìn)行介紹，本節(jié)采用“排序CFAR檢測(cè)”［21-22］的思想在頻域中檢測(cè)射頻干擾。如圖4所示，首先設(shè)置參考單元Nref和保護(hù)單元Npro，然后對(duì)參考單元內(nèi)的數(shù)據(jù)從小到大進(jìn)行排序，并取第0.75Nref個(gè)數(shù)據(jù)的2～5倍為檢測(cè)門限，其中·為向上取整操作。

圖4　檢測(cè)門限示意圖

對(duì)海雜波抑制后的雷達(dá)接收信號(hào)做FFT轉(zhuǎn)換到頻域中，設(shè)置好檢測(cè)門限，若頻譜幅度超過檢測(cè)門限則認(rèn)為其是射頻干擾。若檢測(cè)到射頻干擾，則在原始雷達(dá)接收信號(hào)的頻譜中將其置零，即可實(shí)現(xiàn)射頻干擾抑制。

該算法的具體流程如下：

步驟1將雷達(dá)接收信號(hào)轉(zhuǎn)換為快慢時(shí)間域信號(hào)，并對(duì)慢時(shí)間域做傅里葉變換實(shí)現(xiàn)多普勒處理；

步驟2在雷達(dá)接收信號(hào)的多普勒域中，對(duì)處于零頻區(qū)域附近的海雜波頻譜進(jìn)行置零操作；

步驟3對(duì)多普勒域海雜波置零后的信號(hào)，做逆多普勒處理，即對(duì)慢時(shí)間域做逆傅里葉變換；

步驟4對(duì)海雜波抑制后的接收數(shù)據(jù)，做FFT，在頻域中設(shè)置門限檢測(cè)射頻干擾的位置；

（4）公犬的去勢(shì)。治療前列腺膿腫的最為常用的方法。具體操作過程為陰囊部剪毛、清洗、消毒。左手沿著陰囊頸部握住犬的睪丸，將其輕輕壓向陰囊底部，使兩個(gè)睪丸正好位于陰囊縫際的兩側(cè)，固定住睪丸；切口分別位于陰囊縫際的兩側(cè)0.5 cm處，后手持刀平行于陰囊縫切開陰囊的皮膚和總鞘膜，且勿傷及睪丸實(shí)質(zhì)，切口長(zhǎng)3～4 cm，將睪丸輕輕的擠出陰囊。在睪丸上方的4 cm左右處，貫穿結(jié)扎精索，結(jié)扎要確實(shí)，以防術(shù)后出血。在結(jié)扎線下方1～2 cm處切斷精索，除去睪丸，精索斷端用碘酊消毒，傷口不用縫合，消毒即可。

步驟5對(duì)原始雷達(dá)接收信號(hào)做FFT變換到頻域，并利用上一步中得到的射頻干擾位置，在頻域?qū)⑵渲昧悖?/p>

步驟6將射頻干擾置零后的頻域信號(hào)通過逆快速傅里葉變換（inverse FFT，IFFT）變換到時(shí)域中。

該方案的具體實(shí)施框圖如圖5所示。

圖5　所提射頻干擾抑制算法流程

3　仿真分析

本節(jié)將用傳統(tǒng)的平均距離譜法［18］和本文所提算法分別對(duì)OTHR實(shí)測(cè)數(shù)據(jù)進(jìn)行處理，分析處理結(jié)果，驗(yàn)證該算法的有效性。

實(shí)測(cè)數(shù)據(jù)說明：該實(shí)測(cè)數(shù)據(jù)由我國(guó)某實(shí)驗(yàn)型OTHR所錄取，為雷達(dá)接收信號(hào)經(jīng)過去載頻后的基帶數(shù)據(jù)，主要針對(duì)空中飛行目標(biāo)檢測(cè)，發(fā)射信號(hào)的脈沖重復(fù)周期為12 ms，脈沖重復(fù)周期個(gè)數(shù)為512，距離分辨率為15 km，觀測(cè)多普勒范圍為-41.67～41.67 Hz。

圖6為某幀含有射頻干擾的實(shí)測(cè)數(shù)據(jù)的距離多普勒譜。已知在該幀數(shù)據(jù)中含有2個(gè)目標(biāo)，它們的距離多普勒位置如表1所示。

圖6　射頻干擾抑制前的距離多普勒譜

表1　距離多普勒位置

射頻干擾所在的多普勒單元為-5.86 Hz。

由于目標(biāo)1和射頻干擾在同一個(gè)多普勒單元內(nèi)，所以目標(biāo)1被掩蓋，無法直接發(fā)現(xiàn)。下面首先采用本文所提算法對(duì)其進(jìn)行處理。

如圖7所示，對(duì)海雜波進(jìn)行抑制，在頻譜上設(shè)置門限可以有效地檢測(cè)出射頻干擾所在頻點(diǎn)位置。

圖7　海雜波抑制后數(shù)據(jù)的頻譜和檢測(cè)門限

然后在原始數(shù)據(jù)中將檢測(cè)出的射頻干擾頻點(diǎn)置零，再進(jìn)行后續(xù)的距離多普勒處理，可以得到如圖8所示的距離多普勒譜。在圖8中可以看出，經(jīng)過所提算法處理后，射頻干擾得到了有效地抑制，目標(biāo)1顯露出來，清晰可見。

圖9為目標(biāo)1所在距離單元數(shù)據(jù)經(jīng)過所提算法處理前的多普勒譜。從中可以看出，射頻干擾和目標(biāo)1同處于一個(gè)多普勒單元中，無法直接分辨。

圖10為目標(biāo)1所在距離單元采用所提算法處理后的結(jié)果。從中可以看出，經(jīng)過所提算法處理后，在多普勒單元-5.86 Hz中射頻干擾被抑制掉，只剩下目標(biāo)信號(hào)。

圖8　所提算法處理后的距離多普勒譜

圖9　目標(biāo)1所在距離單元處理前的多普勒譜

圖10　目標(biāo)1所在距離單元處理后的多普勒譜

下面給出傳統(tǒng)平均距離譜法［18］的處理結(jié)果。圖11為傳統(tǒng)距離譜法處理后的距離多普勒譜。從中可以看出，由于目標(biāo)1回波信號(hào)正好被射頻干擾信號(hào)所覆蓋，經(jīng)過平均距離譜法處理后，目標(biāo)1信號(hào)也被抑制掉了。

圖11　傳統(tǒng)平均距離譜法處理后的距離多普勒譜

圖12為目標(biāo)1所在距離單元經(jīng)過傳統(tǒng)平均距離譜法處理后的多普勒譜。從圖12可以看出，目標(biāo)1所在多普勒單元-5.86 Hz處沒有明顯尖峰，表明目標(biāo)1也被抑制掉了。

圖12　目標(biāo)1所在距離單元經(jīng)過傳統(tǒng)平均距離譜法處理后的多普勒譜

經(jīng)過對(duì)上述處理結(jié)果的分析可得，傳統(tǒng)的平均距離譜法在對(duì)射頻干擾進(jìn)行抑制時(shí)，會(huì)將射頻干擾所在多普勒單元中的目標(biāo)同時(shí)抑制掉；而本文所提算法在匹配濾波前，在信號(hào)頻譜中抑制單頻的射頻干擾，可以在盡量不破壞目標(biāo)的前提下有效地抑制射頻干擾，即使目標(biāo)與射頻干擾在相同的多普勒頻率單元上也不會(huì)被抑制掉。

4　結(jié)　論

本文主要針對(duì)當(dāng)前天波雷達(dá)中最普遍、影響最大的短波通信和廣播電臺(tái)的射頻干擾信號(hào)的抑制進(jìn)行研究。通過研究分析發(fā)現(xiàn)，相對(duì)于匹配濾波后的數(shù)據(jù)，在匹配濾波前的雷達(dá)接收數(shù)據(jù)中，射頻干擾和雷達(dá)回波信號(hào)更具有可分性。基于此提出了一種基于頻域干擾置零的天波OTHR射頻干擾抑制算法，該算法在海雜波抑制后的雷達(dá)接收數(shù)據(jù)中檢測(cè)射頻干擾并直接將該干擾頻率成分置零，所述方法可以在不破目標(biāo)信號(hào)的前提下有效地抑制射頻干擾。實(shí)測(cè)數(shù)據(jù)的處理結(jié)果表明，即使目標(biāo)與射頻干擾在相同的多普勒頻率單元上也不會(huì)被抑制掉。

［1］Li H Y，Xu D W，Hu J F，et al.FRFT based algorithm for maneuvering target detection with OTH radar［J］.Systems Engineering and Electronics，2014，36（9）：1725-1730.（李會(huì)勇，許丁文，胡進(jìn)峰，等.基于FRFT的天波雷達(dá)多機(jī)動(dòng)目標(biāo)檢測(cè)［J］.系統(tǒng)工程與電子技術(shù)，2014，36（9）：1725-1730.）

［2］Bo C，Gu H，Su WM.OTHR sea clutter suppression algorithm based on higher order singular value decomposition［J］.Systems Engineering and Electronics，2014，36（5）：872-878.（薄超，顧紅，蘇衛(wèi)民.基于高階奇異值分解的OTHR海雜波抑制算法［J］.系統(tǒng)工程與電子技術(shù)，2014，36（5）：872-878.）

［3］Wu X，Chen J W，Bao Z，et al.Study on key factor and technology for detecting ballistic target in sky wave over-the-horizon radar［J］.Systems Engineering and Electronics，2013，35（11）：2297-2302.（吳瑕，陳建文，鮑拯，等.天波超視距雷達(dá)探測(cè)彈道目標(biāo)關(guān)鍵因素與技術(shù)研究［J］.系統(tǒng)工程與電子技術(shù)，2013，35（11）：2297-2302.）

［4］Lu K，Liu X.Enhanced visibility of maneuvering targets for high-frequency over-the-horizon radar［J］.IEEE Trans.on Antennas and Propagation，2005，53（1）：404-411

［5］Guo X，SunH，Yeo T S.Transient interference excision in over-the-horizon radar using adaptive time-frequency analysis［J］. IEEE Trans.on Geoscience and Remote Sensing，2005，43（4）：722-735

［6］Abramovich Y I，F(xiàn)razer G J，Johnson B A.Principles of modeselective MIMO OTHR［J］.IEEE Trans.on Aerospace and Electronic Systems，2013，49（3）：1839-1868.

［7］Yi J，Wan X，Cheng F，et al.Computationally efficient RF interference suppression method with closed-Form maximum likelihood estimator for HF surface wave over-the-horizon radars［J］. IEEE Trans.on Geoscience and Remote Sensing，2013，51（4）：2361-2372.

［8］Chen G，Zhao Z，Zhu G，et al.HF radio-frequency interference mitigation［J］.IEEE Trans.on Geoscience and Remote Sensing Letters，2010，7（3）：479-482.

［9］Wang W，Wyatt L R.Radio frequency interference cancellation for sea-state remote sensing by high-frequency radar［J］.IET Radar，Sonar&Naυigation，2011，5（4）：405-415.

［10］Meyer F J，Nicoll J B，Doulgeris A P.Correction and characterization of radio frequency interference signatures in l-band synthetic aperture radar data［J］.IEEE Trans.on Geoscience and Remote Sensing，2013，51（10）：4961-4972.

［11］Guo X，Sun H B，Yeo T S.Interference cancellation for highfrequency surface wave radar［J］.IEEE Trans.on Geoscience and Remote Sensing，2008，46（7）：1879-1891.

［12］Zhou H，Wen B Y.Radio frequency interference suppression in small-aperture high-frequency radars［J］.IEEE Trans.on Geoscience and Remote Sensing，2012，9（4）：788-792.

［13］Zhou H，Wen B Y，Wu S C.Dense radio frequency interference suppression in HF radars［J］.IEEE Signal Processing Letters，2005，12（5）：361-364.

［14］Fabrizio GA，Gershman A B，Turley MD.Robust adaptive beamforming for HF surface wave over-the-horizon radar［J］. IEEE Trans.on Aerospace and Electronic Systems，2004，40（2）：510-525.

［15］Su H，Liu H，Shui P，et al.Adaptive beamforming for nonstationary HF interference cancellation in skywave over-the-horizon radar［J］.IEEE Trans.on Aerospace and Electronic Systems，2013，49（1）：312-324.

［16］Deng W，Wen B Y，Zhou H.Suppression of radio-frequency interference in HFGW radar based on eigen decomposition［J］. Chinese Journal of Radio Science，2004，19（5）：565-559.（鄧巍，文必洋，周浩.基于特征分解的高頻地波雷達(dá)抗射頻干擾研究［J］.電波科學(xué)學(xué)報(bào)，2004，19（5）：565-568.）

［17］Niu J，Li M，Zhang L，et al.Radio frequence interference suppression of high frequence groundwave radar based on range-Doppler spectrum［J］.Periodical of Ocean Uniυersity of China，2014，44（12）：111-114.（牛炯，黎明，張玲，等.基于距離多普勒譜的高頻地波雷達(dá)的射頻干擾抑制［J］.中國(guó)海洋大學(xué)學(xué)報(bào)，2014，44（12）：111-114.）

［18］Luo Y Y，Wen B Y，Zhou H.Suppressing of radio-frequency interferences in chirp radar based on averaging range spectrum［J］.Wuhan Uniυersity Journal，2003，19（5）：649-652.（羅郁楊，文必洋，周浩.利用平均距離譜法抑制線性調(diào)頻雷達(dá)射頻干擾［J］.武漢大學(xué)學(xué)報(bào)：理學(xué)版，2003，49（5）：649-652.）

［19］Rosenberg L.Characterization of high grazing angle x-band seaclutter doppler spectra［J］.IEEE Trans.on Aerospace and Electronic Systems，2014，50（1）：406-417.

［20］Li Y J，Wei Y S，Xu R Q，et al.Simulation analysis and experimentation study on sea clutter spectrum for high-frequency hybrid sk ysurface wave propagation mode［J］.IET Radar，Sonar&Naυigation，2014，8（8）：917-930.

［21］Deng X B，Shi C H，Gao C.Subspace signal detection in compound-Gaussian clutter［J］.Systems Engineering and Electronics，2013，35（9）：1836-1840.（鄧曉波，施長(zhǎng)海，高超，等.復(fù)合高斯雜波中子空間信號(hào)檢測(cè)［J］.系統(tǒng)工程與電子技術(shù)，2013，35（9）：1836-1840.）

［22］Zhang Y R，Gao M G，Li Y J.Performance of cooperative suppress jamming on LFM pulsed radar's ML-CFAR detection［J］. Systems Engineering and Electronics，2014，36（11）：2156 2163.（張養(yǎng)瑞，高梅國(guó)，李云杰.LFM脈沖雷達(dá)恒虛警檢測(cè)的協(xié)同壓制干擾［J］.系統(tǒng)工程與電子技術(shù)，2014，36（11）：2156 2163.）

胡進(jìn)峰（1976-），通訊作者，男，副教授，博士，主要研究方向?yàn)殡S機(jī)信號(hào)雷達(dá)、天波雷達(dá)與MIMO雷達(dá)信號(hào)處理。

E-mail：hujf@uestc.edu.cn

李萬閣（1990），男，碩士研究生，主要研究方向?yàn)樘觳ɡ走_(dá)信號(hào)處理。

E-mail：liwange@126.com

艾慧（1990-），女，碩士研究生，主要研究方向?yàn)樘觳ɡ走_(dá)信號(hào)處理。

E-mail：1285483757@qq.com

李會(huì)勇（1975-），男，教授，博士，主要研究方向?yàn)殛嚵行盘?hào)與自適應(yīng)信號(hào)處理。

E-mail：hyli@uestc.edu.cn

李軍（1977-），男，副教授，博士，主要研究方向?yàn)槔走_(dá)系統(tǒng)與信號(hào)處理、自適應(yīng)及陣列信號(hào)處理。

E-mail：lijunsc@uestc.edu.cn

夏威（1980-），男，副教授，博士，主要研究方向?yàn)楦咚賹?shí)時(shí)信號(hào)處理和雷達(dá)系統(tǒng)與信號(hào)處理。

E-mail：wx@uestc.edu.cn

徐威（1974-），男，副教授，碩士，主要研究方向?yàn)楣饫w通信。

E-mail：xw_ld@sina.com

Radio frequency interference suppression algorithm based on interference elimination of frequency domain in OTHR

Hu Jin-feng1，Li Wan-ge1，Ai Hui1，Li Hui-yong1，Li Jun1，Xia Wei1，Xu Wei2
（1.School of Electronic Engineering，Uniυersity of Electronic Science and Technology of China，Chengdu 611731，China；2.Air Force Communications Officer School，Dalian 116600，China）

In over-the-horizon radar（OTHR），the traditional methods suppress the radio frequency interference（RFI）in space or time domain after matched filtering，and they tend to suppress the target around the RFI at the same time.We notice that RFI and the target signal are the similar narrowband signal after matched filtering，so it is difficult to distinguish them in the time domain or frequency domain；while before matched filtering，the target signal is a wideband signal，and the RFI is a narrowband signal.Based on these features，a RFI suppression algorithm based on interference elimination in the frequency domain is proposed.The algorithm firstly suppresses the sea clutter in the Doppler domain and transforms the signal into the frequency domain，detects the frequency position of RFI in the frequency domain，and then sets the frequency position to zero in the frequency domain of the original signal，finally transforms the signal to the time domain.Compared with the traditional methods，the proposed algorithm can effectively suppress the RFI and protect the target signal.Experimental results verify the superiority of the proposed algorithm.

over-the-horizon radar（OTHR）；radio frequency interference（RFI）；suppression

TN 958.93；TN 957

A

10.3969/j.issn.1001-506X.2016.05.12

1001-506X（2016）05-1046-06

2015-05-14；

2015-09-11；網(wǎng)絡(luò)優(yōu)先出版日期：2016-01-12。

網(wǎng)絡(luò)優(yōu)先出版地址：http://www.cnki.net/kcms/detail/11.2422.TN.20160112.1434.004.html

國(guó)家自然科學(xué)基金項(xiàng)目（61101172，61101173，61201280，61371184）；中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金（ZYGX2014J021）資助課題