路翠華，李國(guó)林，熊 波

（海軍航空工程學(xué)院 a.七系；b.兵器科學(xué)與技術(shù)系，山東 煙臺(tái) 264001）

0 引言

調(diào)頻引信根據(jù)回波信號(hào)與發(fā)射信號(hào)的頻差來(lái)提取目標(biāo)信息，但是壓制性干擾的存在可以阻礙引信對(duì)目標(biāo)信息的正常檢測(cè)。壓制性干擾的最佳干擾波形就是隨機(jī)性最強(qiáng)的波形，在平均功率限定的條件下，正態(tài)分布噪聲的隨機(jī)性最大，為最佳壓制性干擾波形。自適應(yīng)濾波技術(shù)發(fā)展至今，已成為干擾抑制[1-6]、回波消除[7]、泄露抵消[8]等方面重要的技術(shù)手段。本文將自適應(yīng)濾波技術(shù)應(yīng)用到線性調(diào)頻引信中，用來(lái)對(duì)壓制性高斯白噪聲進(jìn)行抑制。

自適應(yīng)濾波算法的選擇影響濾波器的收斂速度和精度。LMS算法是一種用瞬時(shí)值估計(jì)梯度矢量的自適應(yīng)算法，具有算法簡(jiǎn)單、計(jì)算量少的優(yōu)點(diǎn)，因而得到了廣泛的應(yīng)用，但其存在收斂速度、跟蹤速度和收斂精度的矛盾。為了解決這一矛盾，人們提出了變步長(zhǎng)自適應(yīng)濾波算法。文獻(xiàn)[9]提出其步長(zhǎng)因子可隨迭代次數(shù)n的增加而逐漸減少的自適應(yīng)濾波算法；文獻(xiàn)[10]提出一種時(shí)間平均估值梯度的自適應(yīng)濾波算法；文獻(xiàn)[11]提出一種步長(zhǎng)因子μ 正比于誤差信號(hào)e (n)大小的自適應(yīng)濾波算法；文獻(xiàn)[12]提出一種步長(zhǎng)因子μ 與e(n)和 x (n)的互相關(guān)函數(shù)的估值成正比的自適應(yīng)濾波算法；文獻(xiàn)[13]給出了Sigmoid函數(shù)變步長(zhǎng)LMS算法；文獻(xiàn)[14]也提出了一種變步長(zhǎng)自適應(yīng)濾波算法。本文采用文獻(xiàn)[14]提出的變步長(zhǎng)自適應(yīng)濾波算法，該算法滿足變步長(zhǎng)自適應(yīng)濾波算法的步長(zhǎng)調(diào)整原則[15]，并克服了Sigmoid 函數(shù)變步長(zhǎng)LMS算法在自適應(yīng)穩(wěn)態(tài)階段仍有較大步長(zhǎng)變化的不足。

1 線性調(diào)頻引信的工作原理

線性調(diào)頻引信為鋸齒波調(diào)制，發(fā)射波的頻率按鋸齒波規(guī)律變化。其發(fā)射、反射和差頻信號(hào)變化規(guī)律如圖1所示。

圖1 鋸齒波調(diào)制時(shí)間—頻率曲線

圖1a）為發(fā)射與回波信號(hào)的時(shí)間—頻率曲線。實(shí)線為發(fā)射信號(hào)頻率 ft，虛線為回波信號(hào)頻率 fr，f0為載波頻率，?F為最大頻偏，Tm為調(diào)制信號(hào)周期。圖1b）為混頻器輸出端差頻信號(hào)的時(shí)間—頻率曲線，其中，F(xiàn)b表示差頻頻率，為nT ～nT +τ 期間的差頻。圖1b）中虛線為考慮多普勒效應(yīng)影響時(shí)的情況。圖1c）為混頻器輸出端差頻信號(hào)電壓Ui與時(shí)間的曲線。由圖1可求得發(fā)射信號(hào)頻率 ft與回波信號(hào)頻率 fr的表達(dá)式為：

式中：c為電磁波的傳播速度。

差頻bF為：

即

由式(3)、(4)看出，當(dāng)調(diào)制參數(shù)mT 和 ?F 一定時(shí)，差頻bF 與距離R 成正比。

若存在壓制性干擾，干擾信號(hào)為正態(tài)分布的噪聲。通過線性調(diào)頻引信通帶的壓制性噪聲太強(qiáng)時(shí)，引信將無(wú)法正確測(cè)得差頻信號(hào)頻率bF，從而無(wú)法正確得到導(dǎo)彈與目標(biāo)的距離R，此時(shí)需要對(duì)引信通帶內(nèi)的噪聲進(jìn)行抑制。

2 噪聲抑制原理

線性調(diào)頻信號(hào)中含有高斯白噪聲時(shí)，調(diào)頻信號(hào)是可預(yù)測(cè)的，但高斯白噪聲的可預(yù)測(cè)性比較差，利用兩者可預(yù)測(cè)性的差異進(jìn)行噪聲抑制。噪聲抑制原理如圖2所示。

圖2 噪聲抑制原理

設(shè)含有噪聲的信號(hào)為x，x=s+v，其中s為線性調(diào)頻信號(hào)，v為高斯白噪聲。利用前向預(yù)測(cè)濾波器可得到信號(hào)s的前向預(yù)測(cè)信號(hào)。誤差e=s +v?，根據(jù)誤差e 調(diào)整前向預(yù)測(cè)濾波器的權(quán)值，使→s，從而使e → v，達(dá)到噪聲抑制的目的。

因?yàn)橐殴ぷ鲿r(shí)間比較短，所以要選計(jì)算量少、收斂速度快的自適應(yīng)濾波算法。本文采用文獻(xiàn)[14]提出的變步長(zhǎng)自適應(yīng)濾波算法。該算法初始收斂階段較大，對(duì)應(yīng)的μ (n)較大，算法收斂速度較快；當(dāng)算法進(jìn)入穩(wěn)態(tài)時(shí)，達(dá)到最小，此時(shí) μ (n)也達(dá)到最小，由此得到最佳Wiener 解。算法如下：

式中：d (n)為期望響應(yīng)；X (n)為輸入信號(hào)矢量；W (n)為權(quán)值；μ (n)為步長(zhǎng)因子；參數(shù)α＞0控制函數(shù)的形狀；參數(shù)β＞0控制函數(shù)的取值范圍。

3 噪聲抑制性能仿真與分析

對(duì)式(8)兩邊取數(shù)學(xué)期望，則

為深入貫徹中央全面從嚴(yán)治黨要求和鄉(xiāng)村振興戰(zhàn)略，培養(yǎng)優(yōu)秀基層黨建和社會(huì)管理人才，提高基層經(jīng)濟(jì)發(fā)展管理質(zhì)量，11月10日-16日，十二師二二二團(tuán)首批28名由連隊(duì)、社區(qū)“兩委”和致富帶頭人、合作社社長(zhǎng)組成的考察團(tuán)赴對(duì)口援疆省市，山西長(zhǎng)治市開展了為期一周的觀摩學(xué)習(xí)。

實(shí)際上很難達(dá)到自適應(yīng)濾波的理想情況，下面通過仿真來(lái)分析線性調(diào)頻引信采用變步長(zhǎng)LMS算法進(jìn)行噪聲抑制的性能。取變步長(zhǎng)LMS算法中的參數(shù)α、β為：α=1 000，β=0.001 54（濾波參數(shù)α、β 選取依據(jù)經(jīng)驗(yàn)值）。當(dāng)SNB=0 dB時(shí)，線性調(diào)頻引信采用變步長(zhǎng)LMS算法進(jìn)行噪聲抑制的效果如圖3所示。

圖3 SNB=0 dB時(shí)噪聲抑制效果

當(dāng)SNB=?5 dB時(shí)，線性調(diào)頻引信采用變步長(zhǎng)LMS算法進(jìn)行噪聲抑制的效果如圖4所示。

圖4 SNB=?5 dB時(shí)噪聲抑制效果

圖3a）和圖4a）為原始信號(hào)，圖3b）和圖4b）為含有噪聲的信號(hào)，圖3c）和圖4c）為噪聲抑制后的輸出信號(hào)。由圖3c）可以看出，SNB=0 dB時(shí)，線性調(diào)頻引信采用變步長(zhǎng)LFM算法很好地抑制了線性調(diào)頻信號(hào)中的噪聲。由圖4c）可以看出，SNB=?5 dB時(shí)，噪聲抑制后的信號(hào)雖然仍舊含有部分噪聲，但信噪比得到了明顯的改善。

當(dāng)SNB=?10 dB時(shí)，線性調(diào)頻引信采用變步長(zhǎng)LMS算法進(jìn)行噪聲抑制的效果如圖5所示。

圖5 SNB=?10 dB時(shí)噪聲抑制效果

由圖5c）可以看出，當(dāng)SNR=?10 dB時(shí)，變步長(zhǎng)LMS算法的參數(shù)為α=1 000，β=0.001 54時(shí)，算法不收斂，達(dá)不到很好的噪聲抑制效果。此時(shí)，需要調(diào)整變步長(zhǎng)LFM算法的參數(shù)α、β。

表1 噪聲抑制前后信噪比及改善量

由表1可以看出，線性調(diào)頻引信利用變步長(zhǎng)LMS算法進(jìn)行噪聲抑制明顯改善了信噪比，并且在算法收斂的條件下，信噪比越小，信噪比的改善量越大。

4 結(jié)論

將自適應(yīng)濾波技術(shù)應(yīng)用到線性調(diào)頻引信中，用來(lái)對(duì)壓制性噪聲進(jìn)行抑制。根據(jù)線性調(diào)頻信號(hào)和噪聲可預(yù)測(cè)性的差異，利用前向預(yù)測(cè)濾波器對(duì)噪聲進(jìn)行抑制，分析了噪聲抑制原理，自適應(yīng)濾波算法采用文獻(xiàn)[14]提出的變步長(zhǎng)LMS 自適應(yīng)濾波算法。取變步長(zhǎng)LMS算法中參數(shù) α=1 000，β=0.001 54，通過仿真分析可以發(fā)現(xiàn)，在SNB=?5dB時(shí)仍然可以達(dá)到很好的噪聲抑制效果，在算法收斂的條件下，信噪比越小，信噪比的改善量越大。

[1]謝春健.基于自適應(yīng)濾波的擴(kuò)頻通信系統(tǒng)抗干擾研究[D].西安:西北工業(yè)大學(xué),2005.

[2]胡軍鋒.寬帶通信中的窄帶干擾抑制技術(shù)研究[D].西安:西安電子科技大學(xué),2004.

[3]周紅彬.直接序列擴(kuò)頻通信中窄帶干擾的識(shí)別和抵消[D].西安:西安電子科技大學(xué),2003.

[4]李紅霞.直接序列擴(kuò)頻信號(hào)中多窄帶干擾的自適應(yīng)抑制技術(shù)[D].成都:電子科技大學(xué),2001.

[5]SAMI AROMAA,PERTTI HENTTU,MARKKU JUNTTI.Transform-selective interference suppression algorithm for spread-spectrum communications[J].IEEE SIGNAL PROCESSING LETTERS,2005,12(1):49-51.

[6]HAOSHEN,ANTONIA PAPANDREOU-SUPPAPPOLA.Wideband time-varying interference suppression using matched signal transforms[J].IEEE Trans.Signal Process.,2005,53(7):2607-2612.

[7]林耀榮.自適應(yīng)濾波理論及其在回波消除中的應(yīng)用研究[D].廣州:華南理工大學(xué),1999.

[8]韓俊杰,黨懷鎖,焦光龍.FTLMS算法在引信泄漏信號(hào)對(duì)消器設(shè)計(jì)中的應(yīng)用[J].彈箭與制導(dǎo)學(xué)報(bào),2002,22(4):66-68.

[9]GITLIN R D,WEINSTEIN S D.On the design of gradient algorithms for digitally implemented adaptie filters[J].IEEE Trans.on CT,1973,20(2):125-136.

[10]GITLIN R D,WEINSTEIN S D.The effects of large interference on the tracking capability of digitally implemented echo cancellers[J].IEEE Trans.on COM,1978,26(6):833-839.

[11]YASUKAWA H,SHIMADA S,FURUKRAWA I.Acoustic echo canceller with high speech quality[C]//ICASSP,1987:2125-2128.

[12]葉華,吳伯修.變步長(zhǎng)自適應(yīng)濾波算法的研究[J].電子學(xué)報(bào),1990,18(4):63-69.

[13]吳光弼,祝琳瑜.一種變步長(zhǎng)LMS 自適應(yīng)濾波算法[J].電子學(xué)報(bào),1994,22(1):55-60.

[14]高鷹,謝勝利.一種變步長(zhǎng)LMS 自適應(yīng)濾波算法及分析[J].電子學(xué)報(bào),2001,29(8):1094-1097.

[15]譚景繁,歐陽(yáng)景正.一種新的變步長(zhǎng)自適應(yīng)濾波算法[J].數(shù)據(jù)采集與處理,1997,12(3):171-174.