兩種干擾抵消方法在拖曳陣應(yīng)用中的對比*

孫春艷　申和平　劉乙樂

(北京神州普惠科技股份有限公司　北京　100085)

?

兩種干擾抵消方法在拖曳陣應(yīng)用中的對比*

孫春艷申和平劉乙樂

(北京神州普惠科技股份有限公司北京100085)

強干擾抵消是大孔徑光纖拖曳線列陣探測系統(tǒng)中的重要研究內(nèi)容。常見的強干擾抵消方法是基于自適應(yīng)濾波的干擾抵消，需要將干擾信號作為參考輸入，屬于時域處理方法?；贐artlett權(quán)的波束零陷形成方法是頻域處理方法，通過在干擾的方向形成零陷對干擾進行抑制。論文將這兩種方法用于抵消實際拖曳數(shù)據(jù)中的近距離強目標，對比兩種方法的抵消效果以及應(yīng)用復雜度，為工程實踐提供參考。實驗結(jié)果表明，基于Bartlett權(quán)的波束零陷形成方法實現(xiàn)簡單，可同時抵消多個干擾，且抵消效果更好，更適合工程應(yīng)用。

強干擾抵消;自適應(yīng)濾波;Bartlett波束零陷

Class NumberTB556

1　引言

大孔徑光纖拖曳線列陣關(guān)注的是低頻遠距離目標，但由于減震降噪技術(shù)的發(fā)展，加上遠距離的傳播，陣列接收到的信號具有極低的信噪比。要成功地實現(xiàn)對弱目標的檢測，光纖拖曳探測系統(tǒng)僅通過加大其陣列的孔徑是不夠的，本艦強干擾和近距離的強目標都會影響對弱目標的檢測。對于拖線陣，本艦干擾可以通過加長拖纜的長度，使接收陣列遠離本艦來減弱本艦干擾的影響。而對于近距離的強目標，只能通過改進信號處理算法來解決。

強干擾抵消方法有很多，比較常用的兩種方法是基于自適應(yīng)濾波的干擾抵消[1]和自適應(yīng)零陷波束形成[2～3]。前者是時域處理方法，需要將強干擾信號作為參考輸入，進行自適應(yīng)濾波后獲得干擾抵消后的輸出。關(guān)鍵問題是參考信號的選取，被動的情況下參考信號未知，需要在干擾方向先進行預成波束，得到對干擾信號的估計。后者是頻域處理方法，采用零陷權(quán)的方式來抑制干擾。本文研究的基于Bartlett權(quán)的零陷波束形成方法，僅需給定要抵消的干擾方位，在干擾方向上形成零陷，達到對干擾抵消的目的。

本文旨在尋找一種簡單實用的干擾抵消方法。利用上述兩種方法分別處理相同的拖曳陣探測數(shù)據(jù)，對比二者的實現(xiàn)復雜度和抵消效果。

2　頻域塊自適應(yīng)濾波方法

基于自適應(yīng)濾波的干擾抵消算法核心部分是自適應(yīng)濾波，其將干擾信號作為參考輸入進行濾波，得到抵消干擾后的輸出。自適應(yīng)濾波器[4]是以最小均方誤差為準則的最優(yōu)濾波器。傳統(tǒng)的時域LMS是按點進行濾波的，在數(shù)據(jù)量大時運算效率低。頻域塊LMS自適應(yīng)濾波方法使濾波過程一塊一塊的進行，改善了計算的復雜度，提高了運算的速率。這種實現(xiàn)方式最初是由Clark[5～6]和Ferrara[7～8]各自獨立提出來的。FLMS算法的主要思想是，將數(shù)據(jù)分成若干塊，每塊含有N個點。每處理完一塊后修改一次權(quán)值，通過傅立葉變換在頻率域上完成濾波器系數(shù)的自適應(yīng)。同時，算法中還用到了重疊存儲的方法為快速卷積提供了強有力的工具。在文獻[9]中對FLMS算法進行了詳細的論述，在此僅給出算法實現(xiàn)的主要步驟。

在數(shù)據(jù)重疊50%存儲的情況下FLMS算法的步驟如下：

1)初始化權(quán)向量為2M×1維零向量，估計功率為2M×1維全1向量。

2)輸入信號為2M個點，每次更新M個點。對輸入信號濾波：

3)誤差估計：

e(k)=d(k)-y(k)

(3)

(4)

4)估計功率更新：

P(k)=λP(k-1)+(1-λ)|X(k)|2

(5)

η(k)=1/(1+P(k))

(6)

5)權(quán)值更新：

φ(k)=IFFT[η(k)XH(k)E(k)]的最前M個元素

(7)

(8)

算法中有兩個參數(shù)需要設(shè)置，步長μ和遺忘因子λ，一般取值在0～1之間，對于不同的數(shù)據(jù)需要調(diào)節(jié)。

3　基于Bartlett權(quán)的波束零陷方法

當有多個相干聲源同時存在于測量區(qū)域范圍內(nèi)時，希望可以找到一種權(quán)，使感興趣的目標方位上波束形成輸出功率最大，而強干擾或不感興趣的目標波束輸出為零。基于Bartlett權(quán)的零陷波束[10～11]形成方法就是基于這一準則。

上述準則利用公式進行描述，可以表述為

max(P1)=max(E[|y(t)|2])=max(WHE[X(t)XH(t)]W)

(9)

式中，max表示取最大值，W為待求解的最優(yōu)權(quán)。X(t)為接收數(shù)據(jù)，是發(fā)射信號與目標方向的導向矢量的乘積。

式(9)中的權(quán)受下述條件約束：

(10)

式中θ1,θ2,…,θN為干擾的方位。a(θ1),a(θ2),…,a(θN)為干擾的導向矢量，經(jīng)過求解可得波束零陷的最佳權(quán)向量為

(11)

式中，B=[I-H(HHH)-1HH]，在B中H=[a(θ1),a(θ2),…,a(θN)]。

得到加權(quán)向量后，利用加權(quán)向量對接收數(shù)據(jù)進行加權(quán)，即可得到零陷后的波束輸出。

Popt(θ)=WHRW

(12)

4　拖曳實驗數(shù)據(jù)處理

對2014年在某海域進行的海上拖曳陣試驗數(shù)據(jù)進行處理。主要研究兩種方法對近距離強目標的抵消效果，以及算法實現(xiàn)的復雜程度。拖線陣包含128個基元，間距2.4m。

對某時段的實驗數(shù)據(jù)進行波束形成處理，得到的方位歷程如圖1所示。

圖1　遠場波束形成歷程圖

從圖1中可以看到，該時段在120°附近存在一個強目標。

1)自適應(yīng)干擾抵消方法

針對上述方位歷程中120°方向強目標在某一時刻的數(shù)據(jù)進行自適應(yīng)抵消，不同參數(shù)下的抵消結(jié)果如圖2、3所示。

圖2　抵消113°方向-步長0.8，遺忘因子0.5

圖3　抵消113°方向-步長0.8，遺忘因子0.8

從圖2和圖3中可以看出，在步長為0.8，遺忘因子為0.5時，113°方向的信號被抵消7.24dB左右，而在步長為0.8，遺忘因子為0.8時，113°方向的信號僅被抵消2.73dB。說明不同的參數(shù)下抵消效果是不同的，實際中要想獲得好的抵消效果需要對參數(shù)不斷進行調(diào)節(jié)。

相同的參數(shù)對49°方向的目標進行抵消，抵消結(jié)果如圖4、5所示。

圖4　抵消49°方向-步長0.8，遺忘因子0.5

圖5　抵消49°方向-步長0.8，遺忘因子0.8

從圖4和圖5中可以看出，對于49°方向的目標，兩對參數(shù)下均只能抵消1dB左右，要想獲得較好的抵消效果，還需不斷調(diào)整參數(shù)進行嘗試。且對于多個目標不能同時進行抵消，需循環(huán)處理。

2)Bartlett波束零陷

對(1)節(jié)的數(shù)據(jù)進行Bartlett波束零陷抵消，分別抵消113°和49°方向的目標，抵消結(jié)果如圖6、7所示。

圖6　抵消113°方向

圖7　抵消49°方向

從圖6和圖7中可以看出，抵消前后113°方向目標的主瓣下降很多，在15.5dB左右，49°方向的目標抵消后也下降9dB左右。該方法沒有需要調(diào)節(jié)的參數(shù)，僅需輸入要抵消的方位即可。應(yīng)用簡單，效果較好。

要想同時抵消49°和113°的目標，只需在要抵消的方位處輸入兩個角度，就可以實現(xiàn)同時抵消。結(jié)果如圖8所示。

圖8　同時抵消兩個方位

從圖8中可以看出，同時抵消與分別抵消效果相同。

5　結(jié)語

本文將兩種抵消方法應(yīng)用于實際的海試拖曳數(shù)據(jù)處理中，對比在實際應(yīng)用中兩種方法的可行性。基于自適應(yīng)濾波的干擾抵消方法，需要先在要抵消的方位預成波束獲得參考輸入，再進行抵消。不同的參數(shù)效果差別較大，且同一組參數(shù)并不適用于所有的干擾，調(diào)節(jié)起來較復雜，效果不佳。而基于Bartlett權(quán)的波束零陷形成算法不需要進行參數(shù)調(diào)節(jié)，且針對多個干擾可以同時抵消，抵消效果較好，簡單易行。綜上，基于Bartlett權(quán)的零陷波束形成方法更適合在抑制強干擾檢測弱目標的實際中進行應(yīng)用。

[1]梅繼丹.水聲聲圖測量技術(shù)研究[D].哈爾濱:哈爾濱工程大學,2010:70-75.

[2]G.W.Elko，A.N.Pong.Simple adaptive first-order differential microphone[C]//NJ,USA:IEEE Piscataway,1995:169-172.

[3]F.Luo，J.Yang，C.Pavlovic，A.Nehorai.Adaptive null-forming scheme in digital hearing aids[J].IEEE Trans On Signal Processing,2002,50(7):1583-1590.

[4]Simon Haykin.Adaptive Filter Theory[M].NJ,USA:Prentice Hall.Inc，1991:2-4.

[5]G.A.Clark，S.K.Mitra，S.R.Parker.Block implementation of adaptive digital filter[J].IEEE Trans Circuits Syst,1981,28(6):584-592.

[6]G.A.Clark，S.R.Parker，S.K.Mitra.A unified approach to timeand frequencydomain realization of FIR adaptive digital filter[J].IEEE Trans Acoust Speech Signal Process,1983,31(5):1073-1083.

[7]E.R.Ferrara.Fast implementation of LMS adaptive filter[J].IEEE Trans Acoust Speech Signal Process,1980,28(4):474-475.

[8]J.J.Shynk.Frequency-domain and multirate adaptive filtering[J].IEEE Signal Process,1992,9(1):14-37.

[9]田超.頻域塊LMS自適應(yīng)濾波算法的研究[D].西安:西安理工大學,2010:15-23.

[10]梅繼丹,惠俊英,王逸林等.Bartlett波束形成的波束零陷權(quán)設(shè)計[J].哈爾濱工程大學學報,2008,29(12):1315-1318.

[11]梅繼丹,王新勇,惠俊英等.近場聲聚焦波束形成與波束零陷研究[J].大連海事學報,2009,35(3):21-24.

Contrast of Two Interference Cancellations in the Application of Towed Array

SUN ChunyanSHEN HepingLIU Yile

(Appsoft Technology Co.Ltd,Beijing100085)

Strong interference cancellation is an important research aspect in the detection system of large-scale optical towed linear array.As a processing method in time domain,interference cancellation based on adaptive filtering,which is used commonly,requires interference signals as reference inputs.And as a processing method in frequency domain,null-forming based on Bartlett weight can restrain interference by null-forming in the direction of interference.In this article,both of two methods mentioned above are applied to cancel nearby strong target in practical data of towed array,the effects and application complexities of both methods are contrasted,and the results provide a reference for engineering practice.The experiment results show that null-forming based on Bartlett weight,which is able to cancel multiple targets simultaneously and more effectively,is more suitable for engineering application.

strong interference cancellation,adaptive filtering,Bartlett null-forming

2016年3月12日,

2016年4月28日

國家重大科學儀器設(shè)備開發(fā)專項資助項目(編號：2013YQ140431)資助。

孫春艷，女，工程師，研究方向：水聲信號處理。申和平，男，高級工程師，研究方向：水聲信號處理。劉乙樂，男，工程師，研究方向：水聲信號處理。

TB556DOI：10.3969/j.issn.1672-9730.2016.09.029