李子高，李淑秋，聞疏琳

(中國(guó)科學(xué)院 聲學(xué)研究所，北京 100190)

一種基于獨(dú)立成分分析的水下寬帶信號(hào)波達(dá)方位估計(jì)方法

李子高，李淑秋，聞疏琳

(中國(guó)科學(xué)院 聲學(xué)研究所，北京 100190)

在水下目標(biāo)特征信息缺乏、聲納陣列孔徑有限的條件下，對(duì)水下寬帶信號(hào)開展高分辨率、高精度的波達(dá)方位估計(jì)研究，具有明顯的研究意義和應(yīng)用價(jià)值。利用水下信源統(tǒng)計(jì)獨(dú)立的特點(diǎn)，采用基于獨(dú)立成分分析的狀態(tài)一致性變換算法進(jìn)行方位估計(jì)；考慮到水下加性噪聲的高斯分布特性，利用小波降噪的預(yù)處理技術(shù)，消減混合信號(hào)中的高斯噪聲成分；針對(duì)狀態(tài)一致性變換在低頻段的性能退化，提出一種空間重采樣修正方法。使用上述基于空間重采樣修正的狀態(tài)一致性變換算法對(duì)仿真及海試數(shù)據(jù)進(jìn)行處理，目標(biāo)方位估計(jì)結(jié)果表明：該方法具有角度分辨率高、需要陣元數(shù)少、對(duì)先驗(yàn)知識(shí)依賴小、低頻段性能穩(wěn)定的優(yōu)點(diǎn)，相比常規(guī)波束形成、子空間法、傳統(tǒng)狀態(tài)一致性變換3類方法更具性能優(yōu)越性與實(shí)用性。

信息處理技術(shù)；寬帶信號(hào)；小波降噪；空間重采樣；頻域獨(dú)立成分分析；波達(dá)方位估計(jì)

0 引言

波達(dá)方位(DOA)估計(jì)是水下陣列信號(hào)處理中一個(gè)重要的研究方向，廣泛運(yùn)用于目標(biāo)探測(cè)、水聲通信等領(lǐng)域。目前常用的DOA方法包括基于常規(guī)波束形成(CBF)的多波束方法，以及基于信號(hào)子空間估計(jì)的高精度DOA方法。前者性能穩(wěn)健但是方位分辨率受到瑞利限約束，后者具有很高精度，但是在實(shí)際應(yīng)用中對(duì)子空間的估計(jì)誤差會(huì)導(dǎo)致算法性能退化[1-3]。

獨(dú)立成分分析(ICA)是一種常用的盲源分離算法，其原理是利用非高斯源信號(hào)之間的統(tǒng)計(jì)獨(dú)立特性，在缺乏信號(hào)特征、信道模型等先驗(yàn)信息的情況下，完成對(duì)源信號(hào)的估計(jì)以及對(duì)混合模型的提取[4]。

由于水下被動(dòng)聲信號(hào)往往基于卷積混合模型[5]，本文采用復(fù)數(shù)快速獨(dú)立成分分析(FastICA)方法[6]，對(duì)陣列采樣信號(hào)在頻域進(jìn)行分解，直接得到頻域信號(hào)混合模型，進(jìn)而獲得入射聲源的DOA估計(jì)，這種方法稱為狀態(tài)一致性變換(SCT)[7]。本文采用的FastICA方法，具有收斂速度快、步長(zhǎng)對(duì)分離結(jié)果影響小、需要陣元數(shù)少等優(yōu)點(diǎn)[8]。由于ICA算法模型限制，混合源信號(hào)中最多只能有一項(xiàng)具有高斯分布特性。因此，針對(duì)水下加性噪聲的高斯分布特性，本文采用小波降噪技術(shù)消減陣列采樣信號(hào)中的高斯噪聲成分[9-10]。此外，留意到在低頻段的分解性能較差，本文提出基于空間重采樣(SR)[11]的修正方法，改善SCT在低頻段的算法表現(xiàn)，本文稱這種方法為SR_SCT. 仿真以及試驗(yàn)數(shù)據(jù)證實(shí)了上述方法在處理寬帶信源以及克服本艦噪聲源干擾中的優(yōu)越性。

1 信號(hào)模型

以均勻線形陣列(ULA)為例，如圖1所示。實(shí)際水下陣列被動(dòng)采樣信號(hào)往往基于卷積混合模型，因此，信號(hào)在頻點(diǎn)fk上的窄帶混合模型可以寫為

Xk=AkSk+Zk，

(1)

式中：k=1,2,…，K為頻點(diǎn)坐標(biāo)；聲源信號(hào)Sk=[s1,k,s2,k,…,sN,k]T，sn,k表示第n號(hào)聲源在頻點(diǎn)fk上的分量；陣列流形矢量Ak=[a1,k,a2,k,…,aN,k]，其中，an,k=[an,k(1),an,k(2),…,an,k(M)]T，an,k(m)=exp (-j2πfk(m-1)dsin (θn)/c)，d為陣元間距，c為聲速；Zk為加性高斯噪聲在頻點(diǎn)fk上的分量。

圖1 均勻線列陣信號(hào)接收模型Fig.1 Received signal model of ULA

2 算法介紹

2.1 算法框圖

本文算法流程圖如圖2所示。算法分為4步：1)基于小波變換對(duì)原始數(shù)據(jù)進(jìn)行降噪處理，減少混合信號(hào)中高斯噪聲成分；2)對(duì)陣列信號(hào)時(shí)頻分析，對(duì)各子頻帶信號(hào)進(jìn)行空間重采樣處理；3)對(duì)各頻帶分別進(jìn)行窄帶ICA處理，得到相應(yīng)的信號(hào)混合模型估計(jì)；4)從各頻點(diǎn)混合矩陣中提取DOA信息，綜合各頻點(diǎn)結(jié)果，得到寬帶入射信號(hào)的DOA估計(jì)。

2.2 小波降噪處理

根據(jù)ICA的前提，信號(hào)源中最多只能有一個(gè)服從高斯分布[4]?？紤]到接收陣列的加性噪聲以及水下環(huán)境噪聲往往呈現(xiàn)為高斯分布，因此在對(duì)采樣信號(hào)進(jìn)行ICA處理之前，加入一個(gè)預(yù)處理流程來(lái)降低混合信號(hào)中的高斯白噪聲成分。

多尺度小波變換常用來(lái)進(jìn)行高斯噪聲抵消[9-10]，本文對(duì)采樣數(shù)據(jù)的噪聲縮減處理分為以下3步：

1) 采樣數(shù)據(jù)經(jīng)離散小波變換(DWT)獲得小波參數(shù)集合；

2) 估算各層小波系數(shù)的高斯噪聲閾值，基于縮減規(guī)則處理各層系數(shù)；

3) 經(jīng)離散小波逆變換，得到降噪后的采樣數(shù)據(jù)。

采樣信號(hào)中高斯噪聲閾值估計(jì)由(2)式給出：

(2)

式中：σN表征采樣數(shù)據(jù)的噪聲能量，本文以絕對(duì)中位差進(jìn)行估計(jì)，

σN=1/0.676 7median(|X-medina(X)|)，

(3)

閾值縮減規(guī)則由(4)式給出：

(4)

式中：CX(i)表示信號(hào)X的第i層小波系數(shù)；α為一正常量，用以調(diào)節(jié)噪聲縮減的幅度，α值越小，降噪幅度越小，反之α值越大，高斯噪聲成分存留越少，但同時(shí)會(huì)導(dǎo)致有效信號(hào)的破壞程度更高。

2.3 基于ICA的窄帶信號(hào)DOA估計(jì)方法

ICA的算法模型為y=As+na，y為混合信號(hào)，A為隨機(jī)混合矩陣，s為源信號(hào)，na為加性噪聲。ICA算法的目標(biāo)是在僅知道y的前提下，獲取混合模型，解混矩陣，以及對(duì)源信號(hào)s的估計(jì)=Hy.

對(duì)比(1)式給出的窄帶陣列信號(hào)混合模型，可以看出與ICA算法模型具有相同形式。在水聲領(lǐng)域,考慮到被動(dòng)聲吶接收的不同艦船輻射噪聲，以及不同聲源發(fā)出的主動(dòng)聲納信號(hào)之間相互獨(dú)立，并與環(huán)境噪聲獨(dú)立，聲源信號(hào)滿足ICA模型的統(tǒng)計(jì)獨(dú)立的假設(shè)[8]，因此ICA方法可以用于處理水聲信號(hào)。在不具備源信號(hào)數(shù)目等先驗(yàn)知識(shí)的條件下，本文采用基于縮進(jìn)模式復(fù)數(shù)FastICA算法[12]，對(duì)獨(dú)立源信號(hào)逐一分解，步驟如下：

1) 白化處理，消除采樣信號(hào)的二階相關(guān)性。

2)選取列向量相互正交的隨機(jī)矩陣作為解混矩陣迭代初始值。

3)對(duì)列向量wn進(jìn)行更新，迭代公式如下：

(5)

式中：wn為解混矩陣W的第n列向量；目標(biāo)函數(shù)g(y)=1/(ε+y)；ε為一常量。

4)去除wn與wp,p=1,2,…,n-1方向上的投影值，從而保證各列向量相互正交。

(6)

5)判斷wn是否收斂，如不收斂，返回第3步繼續(xù)迭代。若收斂，進(jìn)行下一步。

6)根據(jù)當(dāng)前解混向量wn，得到混合矩陣中對(duì)應(yīng)元素n的估計(jì)值。由于不一定是方陣，不能直接以W-1作為A的估計(jì)值，我們以Moore-Penrose偽逆矩陣替代[13]，如(7)式所示：

n=Q-1wn.

(7)

由于信號(hào)混合模型A中包含了入射信號(hào)到達(dá)不同陣元的時(shí)差信息(TDOA)，分析上文得到的信道混合模型估計(jì)中各陣元之間的TDOA，可以得到來(lái)波方位的估計(jì)，這種參數(shù)估計(jì)方法稱為SCT[7].

信源n在陣元a、b之間的TDOA信息包含在中對(duì)應(yīng)元素的相位差中：

(8)

式中：γ為中對(duì)應(yīng)元素幅值比；f為窄帶頻點(diǎn)。

對(duì)比ULA陣列流形矢量Ak，得到入射角θ的估計(jì)值為

(9)

(10)

2.4 常規(guī)寬帶SCT算法

2.3節(jié)給出了針對(duì)單一頻點(diǎn)的SCT參數(shù)估計(jì)方法。對(duì)于寬帶入射信號(hào)，常規(guī)的SCT算法流程如下：

先對(duì)陣列采樣信號(hào)進(jìn)行短時(shí)傅里葉變換(STFT)變換至頻率- 時(shí)間域，從而將信號(hào)頻帶分解為若干窄頻帶(子頻帶)，窄帶信號(hào)模型為

X(fk)=A(fk,Θ)S(fk)+Z(fk),k=1,2,…,K,

(11)

式中：Θ為聲源方位角，Θ=[θ1,θ2,…,θN].

2.5 基于空間重采樣的改進(jìn)SCT算法

實(shí)際使用中發(fā)現(xiàn)，SCT算法在某些頻段性能較差[13]：在高頻段，空域混疊現(xiàn)象導(dǎo)致錯(cuò)誤方位譜峰的出現(xiàn)；在較低頻率處，陣元間距所引起的信號(hào)相差較小，由于算法性能對(duì)于ICA分解誤差非常敏感，導(dǎo)致低頻段處DOA估計(jì)誤差較大。

SR_SCT算法步驟如下：

1)對(duì)陣列時(shí)域采樣信號(hào)進(jìn)行短時(shí)傅里葉變換，分解為多個(gè)子頻帶信號(hào)。

2)選取適當(dāng)?shù)膮⒖碱l率：fmin

3)獲取各個(gè)頻點(diǎn)處的重采樣矩陣：

(12)

式中：u=1,2,…,M′，v=1,2,…,M，M′=?Mfk/f0」，ψk=min (π,πfk/f0).

4)依次計(jì)算各頻點(diǎn)對(duì)應(yīng)的重采樣信號(hào)：XR(fk)=T(fk)X(fk).

5)對(duì)各頻點(diǎn)的重采樣信號(hào)進(jìn)行窄帶DOA估計(jì)，將所得SCT譜進(jìn)行功率疊加，得到寬帶信號(hào)的空間譜估計(jì)。

3 數(shù)據(jù)處理結(jié)果

我們采用以下4個(gè)實(shí)驗(yàn)對(duì)本文提出的方法進(jìn)行驗(yàn)證。

3.1 SCT算法與SR_SCT算法對(duì)比實(shí)驗(yàn)

仿真實(shí)驗(yàn)環(huán)境參數(shù)如表1所示。其中ULA陣間距為0.5 m，因此不出現(xiàn)空域混疊的頻率上限為1 500/(2×0.5)=1 500 Hz. 本文選取的目標(biāo)頻段為150～1 500 Hz，短時(shí)傅里葉變換窗長(zhǎng)度為250點(diǎn)。

實(shí)驗(yàn)對(duì)比了常規(guī)SCT算法與SR_SCT的算法表現(xiàn)。源信號(hào)選取3組不相干的寬帶線性調(diào)頻信號(hào)，入射方位角為-45°、-10°、30°. 選取連續(xù)3個(gè)陣元的采樣數(shù)據(jù)，分別使用兩種方法進(jìn)行處理，圖3(a)、圖3(b)給出了在各個(gè)頻點(diǎn)處的DOA估計(jì)結(jié)果的對(duì)比?？梢钥闯鰞煞N方法在低頻段算法性能

表1 仿真實(shí)驗(yàn)參數(shù)表Tab.1 Simulation arguments

有明顯區(qū)別：常規(guī)SCT算法在低頻段存在大量無(wú)效估計(jì)值,估計(jì)精度無(wú)法保證，SR_SCT算法在全頻段均表現(xiàn)出良好的估計(jì)性能。圖3(c)給出兩種方法的寬帶SCT空間譜對(duì)比，可以看出SR_SCT算法具有更陡的峰值以及更低的旁瓣能量，其中常規(guī)SCT算法中較高的旁瓣能量主要來(lái)源就是低頻段的錯(cuò)誤估計(jì)值。

圖3 SCT算法與SR_SCT算法估計(jì)結(jié)果對(duì)比 Fig.3 Comparison of estimated results of SCT and SR_SCT

此外，為評(píng)估兩種算法的計(jì)算復(fù)雜度，統(tǒng)計(jì)每次獨(dú)立實(shí)驗(yàn)中兩種算法的運(yùn)行時(shí)間，結(jié)果如圖5所示。觀察發(fā)現(xiàn)，結(jié)果與預(yù)期相符合，由于多出空間重采樣步驟，SR_SCT算法相比常規(guī)SCT算法運(yùn)行時(shí)間更長(zhǎng)，計(jì)算復(fù)雜度更高。在此代價(jià)下SR_SCT算法獲得更高DOA估計(jì)精度以及穩(wěn)定性。在對(duì)目標(biāo)分辨精度要求更高的場(chǎng)合，SR_SCT算法更具競(jìng)爭(zhēng)力。

同時(shí)需要注意到，隨著信噪比提升，兩種算法的耗時(shí)差距逐步縮小，在30 dB處趨近于相同，由于空間重采樣的計(jì)算復(fù)雜度與信噪比不相干，為常量值，因此可以看出在低信噪比環(huán)境下，SR_SCT算法相比常規(guī)SCT算法收斂時(shí)間更長(zhǎng)。但是考察收斂結(jié)果可以發(fā)現(xiàn)，在低于0 dB的情況下，兩者DOA估計(jì)結(jié)果均存在大量錯(cuò)誤值，收斂峰值存在偏差；當(dāng)信噪比環(huán)境大于0 dB時(shí)，SR_SCT算法運(yùn)行時(shí)間變化趨于穩(wěn)定，RMSE值也迅速下降至穩(wěn)定值附近，具備較好的收斂精度，相比之下，雖然常規(guī)SCT算法具有更快的收斂速度，但是其估計(jì)結(jié)果存在大量誤差，實(shí)際情況是其更快的收斂于錯(cuò)誤的峰值。

圖4 均方根誤差分布Fig.4 RMSEs of SCT and SR_SCT

圖5 運(yùn)行時(shí)間Fig.5 Running times of SCT and SR_SCT

3.2 CBF算法、MUSIC算與SR_SCT算法對(duì)比實(shí)驗(yàn)

3.1節(jié)實(shí)驗(yàn)證實(shí)了在陣元數(shù)與源信號(hào)數(shù)相等的情況下，兩種基于ICA的DOA估計(jì)方法的適用性，以及二者的算法性能比較。本節(jié)實(shí)驗(yàn)將給出SR_SCT算法與常規(guī)波束形成(CBF)以及基于子空間思想的MUSIC算法的性能對(duì)比。當(dāng)聲場(chǎng)中存在3個(gè)目標(biāo)時(shí)，我們將參與DOA估計(jì)的水聽器信號(hào)數(shù)目設(shè)置為4個(gè)，從而構(gòu)成滿足MUSIC算法的工作條件。為同時(shí)考察信源入射方位角鄰近和分散的情況，我們?cè)O(shè)定信源DOA為-30°、-20°、30°，其中，前兩組信源波達(dá)角接近，第3組分布較遠(yuǎn)。圖6(a)～圖6(c)分別給出3種算法在各個(gè)頻點(diǎn)的窄帶DOA估計(jì)結(jié)果，圖6(d)給出將各頻點(diǎn)窄帶空間譜疊加得到的寬帶信號(hào)空間譜估計(jì)?？梢钥闯?，CBF由于陣元數(shù)較少，分辨率受瑞利限約束，難以分辨方位接近的兩組信源；MUSIC方法在高頻段具有很高的角度分辨率，但在低頻段表現(xiàn)較差，在空間譜中存在很明顯的錯(cuò)誤峰值，并且存在一定程度的角度估計(jì)偏差；本文提出的SR_SCT方法在全頻段均有良好的角度分辨率，空間譜目標(biāo)對(duì)應(yīng)峰值相比MUSIC算法具有更陡峭分布。綜上，在信源方位分布接近時(shí)，CBF與MUSIC算法性能均呈現(xiàn)不同程度的退化，SR_SCT算法相比前二者具有更高的穩(wěn)健型和準(zhǔn)確性。

圖6 3種DOA估計(jì)算法性能對(duì)比Fig.6 Performance comparison of CBF,MUSIC and SR_SCT

3.3 MUSIC算法與SR_SCT算法對(duì)比實(shí)驗(yàn)(信源數(shù)欠估計(jì))

信源數(shù)目是信源先驗(yàn)知識(shí)中很重要的內(nèi)容，對(duì)信源數(shù)的估計(jì)誤差往往會(huì)極大影響DOA估計(jì)結(jié)果。為驗(yàn)證在缺乏信源數(shù)信息的條件下，本文所提出方法的有效性，本節(jié)實(shí)驗(yàn)中設(shè)計(jì)的算法工作參數(shù)為：4組傳感器對(duì)3組信源進(jìn)行采樣，估計(jì)的信源數(shù)為2個(gè)，信源DOA與實(shí)驗(yàn)2相同，圖7(a)和圖7(b)給出MUSIC與SR_SCT的處理結(jié)果。觀察發(fā)現(xiàn)在信源數(shù)欠估計(jì)的情況下，MUSIC算法性能明顯惡化，缺失目標(biāo)峰值且DOA估計(jì)出現(xiàn)偏移，相比之下SR_SCT仍然給出準(zhǔn)確的方位估計(jì)，盡管出現(xiàn)一定程度的峰值衰減。出現(xiàn)這種現(xiàn)象的原因是ICA分解過程中各信源的估計(jì)順序存在隨機(jī)性，雖然在單一頻點(diǎn)上由于信源數(shù)目欠估計(jì)而缺失部分目標(biāo)峰值，但在不同頻點(diǎn)上缺失的目標(biāo)并不相同，通過全頻帶累加仍然得到了全部目標(biāo)對(duì)應(yīng)的峰值。信源數(shù)欠估計(jì)僅導(dǎo)致目標(biāo)峰值的降低，而不會(huì)丟失目標(biāo)，ICA收斂順序的隨機(jī)性增強(qiáng)了在先驗(yàn)知識(shí)缺乏的情況下算法的魯棒性。

圖7 信源數(shù)欠估計(jì)情形下MUSIC與SR_SCT的性能對(duì)比Fig.7 DOA estimates of MUSIC and SR_SCT with underestimating of sources number

3.4 海試數(shù)據(jù)處理

選取一段海試數(shù)據(jù)進(jìn)行算法驗(yàn)證。試驗(yàn)中采用一段32元等間距分布線列水聽器陣進(jìn)行信號(hào)采集，陣元間距1 m，通過拖曳形式保持陣型正直。水聽器列陣所處環(huán)境中包含3個(gè)主要聲源：主動(dòng)聲納信號(hào)、附近水域路過船只輻射噪聲以及線列陣拖船輻射噪聲。其中，主動(dòng)聲源距離較遠(yuǎn)可視為直達(dá)波，入射方位在-14°左右；拖船以及附近船只的輻射噪聲由于在淺海環(huán)境中且距離接收陣列較近，入射信號(hào)呈現(xiàn)多途相干特征，其中拖船噪聲方位在60°～80°之間分布，路過船只方位大致在25°. 圖8給出了這段數(shù)據(jù)的時(shí)頻譜圖，從圖中可以看出聲場(chǎng)中入射信號(hào)的頻率和時(shí)間分布特征。從圖8中可以看出主動(dòng)聲納信號(hào)為線性調(diào)頻信號(hào)，帶寬為200～1 000 Hz，此外，船艦輻射噪聲主要分布于600 Hz以下，在400 Hz處具有較強(qiáng)的線譜。

圖8 時(shí)間頻率譜Fig.8 LoFAR spectrum

分別采用常規(guī)SCT和SR_SCT對(duì)上述數(shù)據(jù)進(jìn)行處理，獲得如下處理結(jié)果，如圖9所示。圖9(a)和圖9(b)給出兩種算法在各個(gè)頻點(diǎn)處的窄帶DOA估計(jì)結(jié)果，橫線代表實(shí)際目標(biāo)方位。可以看出，兩種方法在500～1 000 Hz頻段性能接近，但是SR_SCT在低頻段性能更優(yōu)。由于另外兩組船體輻射信號(hào)主要分布在600 Hz以下，在低頻段常規(guī)SCT無(wú)法有效分解出兩組目標(biāo)的方位，相比之下SR_SCT具有更好的分辨能力，更多的估計(jì)值聚集在目標(biāo)方位附近。此結(jié)論在圖9(c)的方位譜中有更明確體現(xiàn)。圖9(c)給出兩種方法得到的寬帶空間譜估計(jì)結(jié)果，可以看出，對(duì)于主動(dòng)聲納信號(hào)，兩種方法均得出正確的目標(biāo)峰值，但是SR_SCT得到更明顯的另兩組局部峰值，對(duì)應(yīng)于兩組分布于低頻段的船體輻射噪聲。試驗(yàn)數(shù)據(jù)處理結(jié)果與仿真結(jié)果具有一致性，符合對(duì)算法性能的預(yù)期。

圖9 SCT與SR_SCT的DOA估計(jì)結(jié)果對(duì)比Fig.9 Comparison of estimated results of SCT and SR_SCT

本文所研究的陣列采樣信號(hào)，主要由艦船輻射噪聲、主動(dòng)聲納信號(hào)以及加性環(huán)境噪聲組成。通過對(duì)比圖9中船艦輻射噪聲與主動(dòng)聲納信號(hào)的DOA估計(jì)結(jié)果可以發(fā)現(xiàn)，后者具有更準(zhǔn)確的估計(jì)結(jié)果，這是由于兩者的組成成分具有不同的統(tǒng)計(jì)特性。艦船輻射噪聲多為寬帶譜疊加窄帶線譜的形式[14]。據(jù)學(xué)者們對(duì)SCBX數(shù)據(jù)庫(kù)中艦船輻射噪聲的統(tǒng)計(jì)規(guī)律分析，窄帶線譜成分呈亞高斯分布，寬帶譜成分呈高斯分布[15]。對(duì)于主動(dòng)聲納信號(hào)，除寬帶白噪聲信號(hào)以外，大多數(shù)都服從非高斯分布。海洋環(huán)境噪聲以及傳感器加性噪聲均呈現(xiàn)高斯分布。因此對(duì)于本文研究的SR_SCT方法能夠有效作用于艦船輻射噪聲中的線譜成分，以及絕大多數(shù)主動(dòng)聲納信號(hào)成分。相比之下，由于主動(dòng)聲納信號(hào)具備更充分的非高斯特性，相比艦船輻射噪聲，估計(jì)結(jié)果具有更好的精確性和穩(wěn)定性。

船體輻射噪聲的DOA估計(jì)值與真實(shí)值存在一定的偏移，這是由于船體輻射噪聲在淺海聲場(chǎng)環(huán)境中體現(xiàn)為多途相干特性，從而導(dǎo)致ICA分解過程中對(duì)混合模型的估計(jì)產(chǎn)生誤差。對(duì)于相干信源的估計(jì)誤差，是基于ICA的DOA估計(jì)方法的固有問題。

4 結(jié)論

針對(duì)水下寬帶信號(hào)的卷積混合模型，本文提出基于頻域獨(dú)立分量分析的SR_SCT算法。針對(duì)常規(guī)SCT算法在低頻段算法性能的缺陷，本文提出將空間重采樣技術(shù)與SCT相結(jié)合的思路，仿真以及試驗(yàn)數(shù)據(jù)的處理結(jié)果證實(shí)了SR_SCT相比常規(guī)SCT，在低頻段具有更高的估計(jì)精度以及更穩(wěn)定的工作性能。此外，通過與常規(guī)波束形成以及MUSIC的對(duì)比，SR_SCT算法體現(xiàn)出分辨率高，陣元數(shù)需求少，以及對(duì)源信號(hào)先驗(yàn)知識(shí)依賴性低的優(yōu)點(diǎn)。上述優(yōu)點(diǎn)證實(shí)了SR_SCT在水下寬帶目標(biāo)檢測(cè)中的應(yīng)用價(jià)值。但是如何克服淺海環(huán)境中常見的多途現(xiàn)象對(duì)DOA估計(jì)結(jié)果的干擾，有待進(jìn)一步研究。

References)

[1] 李啟虎. 聲吶信號(hào)處理引論[M]. 北京：科學(xué)出版社, 2012. LI Qi-hu. Introduction to sonar signal processing [M]. Beijing: Science Press,2012. (in Chinese)

[2] Ahmad Z, Ali I. Three decades of development in DOA estimation technology[J]. Indonesian Journal of Electrical Engineering and Computer Science, 2014, 12(8): 6297-6312.

[3] Suryaprakash R T, Nadakuditi R R. The performance of MUSIC in white noise with limited samples and missing data[C]∥2014 IEEE Radar Conference. Philadelphia, PA,US: IEEE, 2014: 0940-0944.

[4] Hyv?rinen A. Independent component analysis: recent advances[J]. Phillosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 2013, 371(1984): 20110534.

[5] 章新華, 范文濤, 康春玉, 等. 波束形成與獨(dú)立分量分析融合的寬帶高分辨方位估計(jì)方法[J]. 聲學(xué)學(xué)報(bào), 2009, 34(4): 303-310. ZHANG Xin-hua, FAN Wen-tao, KANG Chun-yu, et al. Broad-band high-resolution bearing estimation method combined with independent component analysis and beamforming [J]. Chinese Journal of Acoustics, 2009, 34(4): 303-310. (in Chinese)

[6] Ruan Z L, Li L P, Qian G B. On extending the complex FastICA algorithms to noisy data[J]. Neural Networks, 2014, 60(C): 194-202.

[7] Nesta F, Omologo M. Generalized state coherence transform for multidimensional TDOA estimation of multiple sources[J]. IEEE Transactions on Audio, Speech, and Language Processing, 2012, 20(1): 246-260.

[8] Du K L, Swamy M N S. Independent component analysis[M]∥Neural Networks and Statistical Learning. London: Springer, 2014: 419-450.

[9] Jansen M. Noise reduction by wavelet thresholding[M]. Berlin: Springer Verlag, 2012.

[10] 林杰, 付夢(mèng)印, 李道平. 自適應(yīng)小波閾值去噪算法及在圖像處理中的應(yīng)用[J]. 兵工學(xué)報(bào), 2011, 32(7): 896-900. LIN Jie, FU Meng-yin, LI Dao-ping. Self-adaptive wavelet threshold de-noising method and its application in image processing [J]. Acta Armamentarii, 2011, 32(7): 896-900. (in Chinese)

[11] Goodwin M M, Elko G W. Constant beamwidth beamforming[C]∥IEEE International Conference on Acoustics, Speech, and Signal Processing.. Minneapolis, MN, US: IEEE, 1993, 1: 169-172.

[12] Wei T. A study of the fixed points and spurious solutions of the deflation-based FastICA algorithm[J]. Neural Computing and Applications, 2015: 1-12.

[13] Sawada H, Mukai R, Araki S, et al. A robust and precise method for solving the permutation problem of frequency-domain blind source separation[J]. IEEE Transactions on Speech and Audio Processing, 2004, 12(5): 530-538.

[14] Yang X, Mu L, Liu Y. Study on line-spectrum stability of ship underwater radiated noise[C]∥ IEEE/OES China Ocean Acoustics Symposium. Harbin, China: IEEE, 2016.

[15] 倪晉平. 水聲信號(hào)盲分離技術(shù)研究 [D]. 西安：西北工業(yè)大學(xué), 2002. NI Jin-ping. A study on blind separation of underwater acoustic signals [D]. Xi’an: Northwestern Polytechnical University, 2002. (in Chinese)

A Novel ICA-based DOA Estimation Method for Underwater Wideband Source

LI Zi-gao, LI Shu-qiu, WEN Shu-lin

(Institute of Acoustics, Chinese Academy of Sciences, Beijing 100190, China)

The estimation of high-resolution and high-accurate direction-of-arrive (DOA) of underwater wideband acoustic sources in the cases of limited array aperture and lack of target feature information is studied. Considering the independent statistical characteristics of underwater bandwidth acoustic sources, a state coherence transform (SCT) method based on independent component analysis (ICA) is utilized to obtain DOA estimation. For the Gaussian distribution characteristics of additive noise, the wavelet denoising pre-processing is used to reduce the Gaussian components in the mixed signals. A spatial resampling method is proposed to fix the performance degradation of SCT in low-frequency range. In the computer simulation and sea trial, the proposed algorithm has the advantages of higher resolution, fewer sensors, less dependence on the priori knowledge and better low-frequency performance compared with other traditional algorithms. The signal processing results confirm the feasibility and performance superiority of the proposed algorithm.

information processing technology; wideband signal; wavelet denoising; spatial resampling; cFastICA; DOA estimation

2016-03-17

國(guó)家“863”計(jì)劃項(xiàng)目(Y212051211)

李子高(1988—)，男，博士研究生。E-mail：328828794@qq.com

李淑秋(1963—)，女，研究員，博士生導(dǎo)師。E-mail：lsq@mail.ioa.ac.cn

TN929.3

A

1000-1093(2017)01-0114-09

10.3969/j.issn.1000-1093.2017.01.015