曹懷剛，任群言，郭圣明，馬力

(1.中國(guó)科學(xué)院 水聲環(huán)境特性重點(diǎn)實(shí)驗(yàn)室，北京 100190；2.中國(guó)科學(xué)院大學(xué)，北京 100049)

矢量水聽器由1個(gè)聲壓通道和3個(gè)具有偶極子指向性的加速度通道組成，單矢量水聽器即可實(shí)現(xiàn)聲源的方位估計(jì)，具有指向性跟聲源頻率無(wú)關(guān)、不存在左右舷模糊的優(yōu)點(diǎn)，同時(shí)還可以有效抑制各向同性噪聲的優(yōu)點(diǎn)。近年來(lái)基于單矢量水聽器方位估計(jì)的研究和應(yīng)用有了很大的進(jìn)展。Nehorai等[1]提出了2種基于單矢量水聽器的方位估計(jì)方法(基于聲強(qiáng)的算法和基于振速協(xié)方差矩陣的算法)并對(duì)其方位估計(jì)的性能進(jìn)行了分析。Levin等[2-3]提出了一種基于最大似然估計(jì)的單矢量水聽器方位估計(jì)方法，并研究了混響環(huán)境下的方位估計(jì)性能。梁國(guó)龍等[4]結(jié)合矢量水聽器自身陣列流型的特點(diǎn)，將陣列信號(hào)處理中的MUSIC算法應(yīng)用到單矢量水聽器上，用單個(gè)矢量水聽器實(shí)現(xiàn)了窄帶信號(hào)和寬帶信號(hào)的高分辨方位估計(jì)。張維等[5]采用量子粒子群求解聲壓和質(zhì)點(diǎn)振速組成的非線性相關(guān)方程組，實(shí)現(xiàn)多目標(biāo)方位估計(jì)。笪良龍等[6]基于海上實(shí)驗(yàn)數(shù)據(jù)分析了平均聲強(qiáng)器和復(fù)聲強(qiáng)器2種方法的方位估計(jì)性能。胡承彥等[7]提出了一種利用單矢量水聽器聲壓和加速度互譜的方位估計(jì)方法。由于單矢量水聽器的指向性指數(shù)最高為6 dB[8]，傳統(tǒng)的基于單矢量水聽器的方位估計(jì)方法對(duì)信噪比都有較大的依賴。因此，一種能在相干噪聲干擾下實(shí)現(xiàn)較高精度的方位估計(jì)的方法是目前所急需的。機(jī)器學(xué)習(xí)是近年來(lái)的研究熱點(diǎn)，其在聲源定位方面表現(xiàn)出了優(yōu)于傳統(tǒng)方法的良好性能。NIU 等[9-10]研究了基于垂直陣的前饋神經(jīng)網(wǎng)絡(luò)、支持向量機(jī)和隨機(jī)森林3種機(jī)器學(xué)習(xí)方法在聲源測(cè)距上的性能，并用實(shí)驗(yàn)數(shù)據(jù)進(jìn)行了驗(yàn)證。HUANG等[11]將卷積神經(jīng)網(wǎng)絡(luò)的輸出作為前饋神經(jīng)網(wǎng)絡(luò)的輸入，研究了深度神經(jīng)網(wǎng)絡(luò)在垂直陣聲源測(cè)距中的性能?；跈C(jī)器學(xué)習(xí)的垂直陣水面水下目標(biāo)分辨、水平陣方位估計(jì)等也都取得了一定的進(jìn)展[12-14]?；跈C(jī)器學(xué)習(xí)的方法具有比傳統(tǒng)的定位方法更高的定位精度和更好的環(huán)境適應(yīng)能力。本文利用深度學(xué)習(xí)的方法，不依賴實(shí)驗(yàn)數(shù)據(jù)，僅利用KRAKEN仿真的矢量聲場(chǎng)數(shù)據(jù)來(lái)訓(xùn)練卷積神經(jīng)網(wǎng)絡(luò)，利用訓(xùn)練好的神經(jīng)網(wǎng)絡(luò)模型來(lái)估計(jì)目標(biāo)方位，以期獲得比傳統(tǒng)方法更高的估計(jì)精度。

1 矢量信號(hào)模型和數(shù)據(jù)的預(yù)處理

1.1 矢量信號(hào)模型

矢量水聽器能同時(shí)共點(diǎn)地測(cè)量聲場(chǎng)的聲壓和加速度的3個(gè)正交分量，考慮聲場(chǎng)中的1個(gè)微分單元，根據(jù)牛頓第二定律可得：

(1)

(2)

(3)

(4)

式中：ax、ay、az分別為加速度a的3個(gè)正交分量;α為極角;θ為方位角，是本文方位估計(jì)值。

1.2 數(shù)據(jù)的預(yù)處理

在聲壓和加速度輸入到神經(jīng)網(wǎng)絡(luò)之前需要對(duì)其進(jìn)行預(yù)處理，以便神經(jīng)網(wǎng)絡(luò)可以更有效的提取特征，首先對(duì)聲壓和加速度進(jìn)行歸一化處理：

(5)

圖1 加速度a及其3個(gè)正交分量ax、ay、az 的投影圖Fig.1 Acceleration a and projection of three orthogonal components ax, ay and az

其次，將歸一化后的聲壓和加速度進(jìn)行互譜運(yùn)算后再進(jìn)行快拍平均：

(6)

由于本文的深度學(xué)習(xí)方法屬于監(jiān)督學(xué)習(xí)，所以在訓(xùn)練卷積神經(jīng)網(wǎng)絡(luò)的時(shí)候需要用到標(biāo)簽(真實(shí)方位角)，標(biāo)簽數(shù)據(jù)的預(yù)處理如下：

(7)

式中：tn為第n個(gè)采樣點(diǎn)的標(biāo)簽，d=[1,2,3,…,360]；dnr為第n個(gè)采樣點(diǎn)的真實(shí)角度；σ表示方位角的模糊度范圍，在本文中σ為15。

1.3 估計(jì)結(jié)果的度量標(biāo)準(zhǔn)

本文中設(shè)估計(jì)值θe和真實(shí)值θr之差的絕對(duì)值小于10°(|θe-θr|<10°)的估計(jì)值為正確的估計(jì)結(jié)果，用估計(jì)準(zhǔn)確率作為方法性能的度量標(biāo)準(zhǔn)之一，估計(jì)準(zhǔn)確率為：

(8)

式中：NC表示估計(jì)正確的估計(jì)結(jié)果的個(gè)數(shù)；N為總的估計(jì)結(jié)果數(shù)。

均方根誤差ERMSE為方法性能的度量標(biāo)準(zhǔn)，均方根誤差表示為：

(9)

2 卷積神經(jīng)網(wǎng)絡(luò)的設(shè)計(jì)和訓(xùn)練

2.1 神經(jīng)網(wǎng)絡(luò)的設(shè)計(jì)

本文采用的深度神經(jīng)網(wǎng)絡(luò)為卷積神經(jīng)網(wǎng)絡(luò)。卷積神經(jīng)網(wǎng)絡(luò)整體上由輸入層、隱藏層和輸出層組成，其中隱藏層包括卷積層、激活函數(shù)、池化層和全連接層。卷積層是卷積神經(jīng)網(wǎng)絡(luò)的核心，卷積的工作原理可以理解為：在三維輸入數(shù)據(jù)上滑動(dòng)n×n的窗口，在每個(gè)可能的位置停止并提取周圍特征的三維數(shù)據(jù)塊。每個(gè)三維數(shù)據(jù)塊與學(xué)到的同一個(gè)權(quán)重矩陣(卷積核)做向量積，轉(zhuǎn)化為一維的向量。對(duì)所有這些向量進(jìn)行空間重組，使其轉(zhuǎn)化為三維輸出特征數(shù)據(jù)。輸出特征數(shù)據(jù)中的每個(gè)空間位置都對(duì)應(yīng)于輸入特征數(shù)據(jù)中相同的位置。該卷積層的輸出即為下一個(gè)卷積層的輸入。卷積的計(jì)算過程為：

(10)

式中：K表示卷積核；I表示滑動(dòng)窗口提取的數(shù)據(jù)；i、j表示數(shù)據(jù)的位置。一般來(lái)說(shuō)，計(jì)算機(jī)在處理數(shù)據(jù)時(shí)，數(shù)據(jù)會(huì)被離散化，所以卷積運(yùn)算由積分變?yōu)殡x散求和。設(shè)最后一個(gè)卷積層的輸出為x=[x1,x2,…,xk]，卷積層的輸出即為全連接層的輸入，對(duì)于全連接層，其計(jì)算的過程為：

(11)

式中：s為神經(jīng)元的序數(shù)；wsk、ws0為權(quán)重和偏置；as經(jīng)過激活函數(shù)的轉(zhuǎn)化后即可得到該全連接層的輸出bs：

bs=f(as)

(12)

該全連接層的輸出即為下一個(gè)全連接層的輸入，并重復(fù)式(11)、(12)的運(yùn)算，其中激活函數(shù)f(·)采用線性整流函數(shù)(ReLU)，其表達(dá)式為：

f(x)=max(0,x)

(13)

本文所采用的卷積神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)如圖2所示，卷積層中“3×3”表示卷積核的大小，第3個(gè)數(shù)字表示卷積核的個(gè)數(shù)(如第1卷積層“3×3×128”中“128”表示該卷積層共有128個(gè)卷積核)。由于輸入數(shù)據(jù)的維度較小(3×202)，在設(shè)計(jì)的神經(jīng)網(wǎng)絡(luò)中并沒有添加池化層。為降低訓(xùn)練過程中神經(jīng)網(wǎng)絡(luò)過擬合造成的影響，在卷積層和全連接層之間，全連接層和輸出層之間加入隨機(jī)失活層。圖2所示為本文所用卷積神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)全連接層一共有4個(gè)，每個(gè)全連接層包含2 048個(gè)神經(jīng)元。

2.2 神經(jīng)網(wǎng)絡(luò)的訓(xùn)練

神經(jīng)網(wǎng)絡(luò)的訓(xùn)練和學(xué)習(xí)過程主要圍繞以下4個(gè)方面：1)由多個(gè)層組成網(wǎng)絡(luò)層；2)輸入數(shù)據(jù)和相應(yīng)的目標(biāo)值；3)損失函數(shù)，即用于學(xué)習(xí)的反饋信號(hào)；4)決定學(xué)習(xí)過程如何進(jìn)行的優(yōu)化器。關(guān)系如圖3所示。

圖2 卷積神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)Fig.2 The structure of the CNN

圖3 卷積神經(jīng)網(wǎng)絡(luò)的訓(xùn)練流程Fig.3 The processing diagram for CNN training

神經(jīng)網(wǎng)絡(luò)的訓(xùn)練和學(xué)習(xí)過程為通過損失函數(shù)和優(yōu)化器不斷更新權(quán)重使神經(jīng)網(wǎng)絡(luò)的預(yù)測(cè)值不斷接近真實(shí)目標(biāo)值的過程。本文采用的訓(xùn)練數(shù)據(jù)為KRAKEN模型仿真的數(shù)據(jù)，海洋環(huán)境為典型的淺海波導(dǎo)，海底為分層介質(zhì)包括沉積層和基底。仿真的聲源、接收和環(huán)境參數(shù)如表1所示。共仿真輸入數(shù)據(jù)360×25×2×4×6=432 000組，將利用表1中參數(shù)仿真的數(shù)據(jù)訓(xùn)練的神經(jīng)網(wǎng)絡(luò)命名為CNN-3。

3 實(shí)驗(yàn)仿真及數(shù)據(jù)處理

本文通過仿真數(shù)據(jù)研究CNN-3在不同信噪比下的方位估計(jì)性能。同時(shí)引入加權(quán)直方圖法[6,15]來(lái)跟CNN-3進(jìn)行對(duì)比。

仿真的環(huán)境參數(shù)如下：在1°～360°等間隔的取100個(gè)采樣點(diǎn)，聲源深度為4 m，聲源距離2 km，信號(hào)頻段為800～900 Hz，采樣間隔為1 Hz，海深61 m，沉積層厚度為6 m，沉積層聲速為1 600 m/s，沉積層密度為1.6 g/cm3，沉積層衰減為0.3 dB/λ，基底聲速為1 750 m/s，基底密度為1.8 g/cm3，基底衰減為0.3 dB/λ。利用KRAKEN模型仿真聲壓和質(zhì)點(diǎn)加速度，并加入色噪聲，設(shè)置-5、0、5、10和15 dB 5個(gè)信噪比，分別利用加權(quán)直方圖法和CNN-3對(duì)仿真的帶噪聲數(shù)據(jù)進(jìn)行方位估計(jì)，得到的2種方法在不同信噪比下的估計(jì)準(zhǔn)確率和均方根誤差如表2所示。

表1 仿真環(huán)境的參數(shù)Table 1 The parameters of simulated environment

由表2的結(jié)果可知，2種方法的方位估計(jì)性能隨信噪比的增加而增加，當(dāng)SNR≥10 dB時(shí)，2種方法的估計(jì)準(zhǔn)確率均達(dá)到100%，均方根誤差均小于2°。而當(dāng)SNR<10 dB時(shí)，CNN-3的方位估計(jì)性能要明顯的優(yōu)于加權(quán)直方圖的方法，特別是在SNR=-5 dB時(shí)，CNN-3的估計(jì)準(zhǔn)確率為87%，均方根誤差為6.9°；而加權(quán)直方圖法的估計(jì)準(zhǔn)確率只有66%，均方根誤差為28.7°。隨著信噪比的增加，兩者的性能差距在逐漸縮小?？傮w來(lái)說(shuō)，卷積神經(jīng)網(wǎng)絡(luò)(CNN-3)具有良好的方位估計(jì)性能，特別是當(dāng)信噪比較低時(shí)，其估計(jì)結(jié)果依然十分可靠。

表2 加權(quán)直方圖法和CNN-3在不同信噪比下方位估計(jì)的準(zhǔn)確率和均方根誤差

2018年在中國(guó)南海進(jìn)行了單矢量水聽器的被動(dòng)方位估計(jì)實(shí)驗(yàn)。實(shí)驗(yàn)過程中將船的輻射噪聲作為聲源，聲源船船長(zhǎng)60 m，船寬26 m，圍繞矢量潛標(biāo)逆時(shí)針作半徑約為1 km的圓周運(yùn)動(dòng)，船速保持在10 kN左右，其運(yùn)動(dòng)的航跡如圖4(a)所示，聲源船與矢量潛標(biāo)的距離隨時(shí)間的變化如圖4(b)所示。實(shí)驗(yàn)海區(qū)海深61 m，矢量水聽器以潛標(biāo)的形式布放在水下44 m的深度。矢量水聽器的采樣率為12 kHz，有效工作頻段為20～3 000 Hz。實(shí)驗(yàn)過程中，矢量水聽器的航向角變化如圖4(c)所示，航向角的變化代表矢量水聽器自身水平方向的轉(zhuǎn)動(dòng)，從圖中來(lái)看其隨時(shí)間有較快且幅度較大的抖動(dòng)，本文中的方位估計(jì)結(jié)果都根據(jù)航向角進(jìn)行了補(bǔ)償。矢量水聽器的各通道時(shí)頻圖如圖5所示，分別是聲壓和加速度通道在100～1 000 Hz頻帶內(nèi)的時(shí)頻圖。

圖4 實(shí)驗(yàn)過程中聲源和矢量水聽器的主要參數(shù)變化Fig.4 The main parameter changes of source ship and vector sensor during the experimen

實(shí)驗(yàn)海區(qū)距離三亞港口較近，有大量的航船來(lái)往，由于遠(yuǎn)處航船噪聲的能量主要集中在低頻段(800 Hz以下)，但由于其強(qiáng)度較大，對(duì)估計(jì)近處的目標(biāo)艦船的方位角有較大的干擾，因此在進(jìn)行實(shí)驗(yàn)數(shù)據(jù)處理時(shí)選擇800～900 Hz的頻段。

將矢量水聽器測(cè)得的聲壓和質(zhì)點(diǎn)加速度按式(5)、(6)進(jìn)行預(yù)處理后輸入訓(xùn)練好的卷積神經(jīng)網(wǎng)絡(luò)CNN-3，得到的方位估計(jì)結(jié)果如圖6(b)所示，圖6(a)為加權(quán)直方圖法得到的方位估計(jì)結(jié)果，其結(jié)果作為深度學(xué)習(xí)法的對(duì)照。

將圖6的估計(jì)結(jié)果分3段來(lái)分析：

1)在0～200 s的時(shí)間段內(nèi)，圖6(a) 基于加權(quán)直方圖的估計(jì)結(jié)果的準(zhǔn)確率為29%，均方根誤差為13.3°；圖6(b) 基于CNN-3的估計(jì)結(jié)果的準(zhǔn)確率為85%，均方根誤差為7.6°。在這個(gè)時(shí)間段內(nèi)卷積神經(jīng)網(wǎng)絡(luò)的性能要明顯的優(yōu)于加權(quán)直方圖法，但2種方法均存在較大誤差，誤差出現(xiàn)的原因可能是由于在這個(gè)時(shí)間段內(nèi)有較強(qiáng)的干擾出現(xiàn)。

圖5 矢量水聽器接收航船噪聲信號(hào)的時(shí)頻Fig.5 The time frequency diagram of ship noise received by a vector sensor

圖6 實(shí)驗(yàn)數(shù)據(jù)的方位估計(jì)結(jié)果Fig.6 The azimuth of experimental data estimated

2)在200～400 s的時(shí)間段內(nèi)，圖6(a)基于加權(quán)直方圖法的估計(jì)結(jié)果的準(zhǔn)確率為24%，均方根誤差為52.1°；圖6(b) 基于CNN-3的估計(jì)結(jié)果的準(zhǔn)確率為24%，均方根誤差為22.9°。在這個(gè)時(shí)間段內(nèi)，2種方法的估計(jì)準(zhǔn)確率都非常低，但結(jié)合圖6和其均方根誤差來(lái)看，基于加權(quán)直方圖的估計(jì)結(jié)果已經(jīng)嚴(yán)重的偏離真實(shí)值，而基于CNN-3的估計(jì)結(jié)果雖然也有較大的誤差，但沒有嚴(yán)重偏離真實(shí)值，仍然有一定的參考意義。這一段較大誤差出現(xiàn)的原因?yàn)椋河蓤D4(b)可知，在200～400 s的時(shí)間段內(nèi)聲源距離增加，聲源距離的增加導(dǎo)致信噪比的下降從而造成了誤差的增大。

3)在400～1 331 s的時(shí)間段內(nèi)，圖6(a) 基于加權(quán)直方圖的估計(jì)結(jié)果準(zhǔn)確率為93%，均方根誤差為5.7°；圖6(b) 基于CNN-3的估計(jì)結(jié)果的準(zhǔn)確率為97%，均方根誤差為4.8°。在這個(gè)時(shí)間段內(nèi)2種方法都有較好的表現(xiàn)，CNN-3的估計(jì)結(jié)果略好于加權(quán)直方圖法。

綜上，在各個(gè)時(shí)間段內(nèi)，基于卷積神經(jīng)網(wǎng)絡(luò)方法的性能要優(yōu)于加權(quán)直方圖法，尤其是在0～400 s干擾較強(qiáng)、信噪比較低時(shí)，卷積神經(jīng)網(wǎng)絡(luò)的優(yōu)勢(shì)更加明顯，進(jìn)一步驗(yàn)證了本文方法的可靠性。

對(duì)于實(shí)驗(yàn)中出現(xiàn)的誤差，除了上面分析的幾點(diǎn)外還存在以下3點(diǎn)：1)作為聲源的船只并不是一個(gè)點(diǎn)聲源，船長(zhǎng)為60 m，而記錄GPS的天線位于船艉，在1 km的聲源距離上由此造成的誤差最大能到3.4°；2)估計(jì)結(jié)果根據(jù)羅經(jīng)記錄到的航向角進(jìn)行了補(bǔ)償，但羅經(jīng)記錄的航向角本身可能存在一定的誤差，進(jìn)行補(bǔ)償時(shí)就把誤差引入了結(jié)果之中；3)矢量水聽器x、y通道的指向性并不是嚴(yán)格的正交，不滿足加權(quán)直方圖法和仿真的訓(xùn)練數(shù)據(jù)正交的前提假設(shè)，從而引入了誤差。

4 結(jié)論

1)與很多機(jī)器學(xué)習(xí)定位方法相比，本文方法不需要實(shí)測(cè)數(shù)據(jù)作為訓(xùn)練數(shù)據(jù)，僅利用KRAKEN仿真數(shù)據(jù)對(duì)卷積神經(jīng)網(wǎng)絡(luò)進(jìn)行訓(xùn)練，降低了訓(xùn)練成本。

2)仿真結(jié)果表明利用卷積神經(jīng)網(wǎng)絡(luò)的方法可以獲得精確的方位估計(jì)結(jié)果，特別是在低信噪比的情況下(-5 dB)其依然有較高的估計(jì)準(zhǔn)確率EACCU=87%，ERMSE=6.9°，高于加權(quán)直方圖法EACCU=66%，ERMSE=28.7°。實(shí)驗(yàn)數(shù)據(jù)處理結(jié)果說(shuō)明在各個(gè)時(shí)間段基于卷積神經(jīng)網(wǎng)絡(luò)的方法都要優(yōu)于傳統(tǒng)的加權(quán)直方圖法，進(jìn)一步證明了卷積神經(jīng)網(wǎng)絡(luò)法的可靠性。

3)由實(shí)驗(yàn)數(shù)據(jù)分析結(jié)果可知，其他船只的干擾、矢量水聽器各通道的指向性不嚴(yán)格正交、聲源距離增加導(dǎo)致信噪比的下降、聲源船不滿足點(diǎn)聲源條件以及羅經(jīng)存在誤差等都是造成方位估計(jì)結(jié)果不準(zhǔn)確的原因，在以后的實(shí)驗(yàn)中應(yīng)盡量避免這些因素的干擾。