周 彤，宋立臣

(大連測控技術(shù)研究所，遼寧 大連 116013)

0 引 言

在艦船水下輻射噪聲測量的過程中，安全、可靠地引導(dǎo)被測目標(biāo)通過測量區(qū)域，并保證被測目標(biāo)與測量系統(tǒng)之間的有效距離以滿足測試需求是極其重要的，也是一直關(guān)心的問題。現(xiàn)有的被測目標(biāo)引導(dǎo)方法，往往依靠單一導(dǎo)航信標(biāo)配合測距系統(tǒng)實(shí)現(xiàn)，通常僅提供被測目標(biāo)與導(dǎo)航信標(biāo)之間的單點(diǎn)距離信息，缺乏方位信息，使得目標(biāo)實(shí)際航行軌跡與試驗(yàn)要求軌跡有一定的偏差，造成測量單程無效，從而降低了測試效率，延長了試驗(yàn)時間[1]。

目前方位估計算法的通常假設(shè)目標(biāo)在觀測時間內(nèi)為靜止?fàn)顟B(tài)，即靜態(tài)方位算法，同時采用大計算量和多快拍數(shù)的批處理形式進(jìn)行方位估計[2]。然而在實(shí)際情況下，目標(biāo)通常為運(yùn)動狀態(tài)，批處理方式的靜態(tài)方位估計方法不適用于實(shí)時連續(xù)跟蹤。同時，靜態(tài)方位估計方法是互不相關(guān)地估計2個連續(xù)時刻的方位信息，并沒有考慮到2個連續(xù)時刻狀態(tài)的關(guān)聯(lián)性，進(jìn)而造成動態(tài)方位估計偏差較大。隨著對運(yùn)動目標(biāo)的測向精確度和實(shí)時處理的要求越來越高，靜態(tài)方位估計算法已不能滿足日益增長的需求，因此研究適用于動態(tài)目標(biāo)的方位跟蹤算法成了一個新趨勢。本文將靜態(tài)方位估計中的MUSIC算法和粒子濾波跟蹤算法相結(jié)合，既提高了方位估計的精度更減小了計算量以適用于對目標(biāo)的實(shí)時方位跟蹤。

1 AVS 和 MUSIC 方位估計算法

1.1 AVS 接收信號模型

通常在實(shí)際應(yīng)用中矢量水聽器接收的水下目標(biāo)輻射聲場為平面波聲場，因此假設(shè)聲源為平面波。則矢量水聽器接收的聲壓和振速的關(guān)系為

實(shí)際應(yīng)用中處理的信號為離散信號，一般采用快拍數(shù)對每時刻數(shù)據(jù)進(jìn)行批處理。假設(shè)快拍數(shù)為N在k時刻的聲源信號為

接收信號包含方位角和俯仰角信息，進(jìn)而可以對數(shù)據(jù)進(jìn)行方位估計。

1.2 MUSIC 方位估計

MUSIC算法作為子空間特征分解相關(guān)算法的代表，實(shí)現(xiàn)了方位估計技術(shù)向高分辨的發(fā)展。MUSIC算法的核心思想是利用數(shù)據(jù)的自相關(guān)矩陣的內(nèi)在結(jié)構(gòu)特性，將數(shù)據(jù)的協(xié)方差矩陣分解為和陣列流形一致的信號子空間，以及與信號子空間正交的噪聲子空間，并利用2個子空間的正交性，通過對譜峰進(jìn)行搜索來估計信源的方位。

工程實(shí)際中，為提高信號處理的精度，采用多快拍數(shù)的觀測數(shù)據(jù)。設(shè)快拍數(shù)為L，則接收數(shù)據(jù)的自相關(guān)矩陣可表示為

其中，G為與噪聲相關(guān)的特征值矩陣對應(yīng)的噪聲子空間。

因此，通過譜估計式（11）的峰值搜索估計和峰值相對應(yīng)的角度值即為入射的方向信息。

在信噪比（SNR）相對較高的噪聲環(huán)境中，MUSIC方法能夠通過如圖1所示的尖銳峰值來呈現(xiàn)聲源DOA。聲源真實(shí)方位為（20.5°，40.5°），然而，當(dāng)SNR低時，峰值可能失真并且估計的DOA可能遠(yuǎn)離真實(shí)位置，如圖2所示。

圖1 MUSIC 方法 DOA 估計性能 SNR=10 dBFig.1 DOA estimated performance based on the MUSIC algorithm SNR=10 dB

圖2 MUSIC 方法 DOA 估計性能 SNR=-10 dBFig.2 DOA estimated performance based on the MUSIC algorithm SNR=-10 dB

2 基于粒子濾波算法的方位跟蹤方法

為了給出DOA跟蹤問題的總體框架，首先定義狀態(tài)空間模型。假設(shè)目標(biāo)聲源當(dāng)前方位為，并以（rad/s）的速度移動。因此，目標(biāo)聲源的狀態(tài)可以由當(dāng)前方位和運(yùn)動速度構(gòu)成，即。采用勻速模型CV（Constant Velocity）模擬單信源的方位角度變化情況[5]，DOA的運(yùn)動速度為常數(shù)。則狀態(tài)空間模型為

矢量水聽器的輸出矩陣是天然的觀測模型，即為式（7）。采用基于遞歸方式的貝葉斯重要性采樣，并以預(yù)測和更新的方式實(shí)現(xiàn)狀態(tài)估計[6]。

因此粒子權(quán)重在算法的濾波性能上具有重要的作用。當(dāng)粒子接近真實(shí)位置時，粒子的似然函數(shù)具有較大的權(quán)重，則為了讓粒子在更新時使權(quán)重大的粒子能夠在重采樣時取代權(quán)重小的粒子，由式（20）可知似然函數(shù)的確定是算法性能的關(guān)鍵。假設(shè)測量的過程噪聲為高斯噪聲，則似然函數(shù)的最大似然估計為[9]

似然函數(shù)的主瓣通常比較平緩，尤其在低信噪比的條件下，因此不能有效的通過權(quán)重衡量真實(shí)狀態(tài)周圍的粒子并且可能出現(xiàn)虛假峰值。同時當(dāng)快拍數(shù)很大時，將會導(dǎo)致粒子的權(quán)重接近于0，進(jìn)而導(dǎo)致對粒子的權(quán)重估計失效。

因此引入MUSIC方法作為似然函數(shù)，即

則似然函數(shù)被重塑，以增強(qiáng)在高似然區(qū)域采樣粒子的權(quán)重。這一步非常重要，因?yàn)樗軌驇椭罄m(xù)重采樣算法更有效地選擇和復(fù)制粒子。

給出最終PF-MUSIC方法的流程為

步驟1 初始化：采用隨機(jī)方法估計初始角度，從均勻分布中采樣，粒子數(shù)為N，。初始化各粒子權(quán)重為。

步驟7 令k=k+1，重復(fù)步驟2～步驟7。直到k＞T時，終止計算。

3 方位跟蹤方法性能驗(yàn)證分析

3.1 仿真分析方位跟蹤方法性能

圖3 SNR=10 條件下的方位跟蹤Fig.3 DOA tracking under the condition of SNR=10 dB

圖3 給出了MUSIC方法和PF-MUSIC跟蹤方法的DOA跟蹤性能。MUSIC方法的性能在多個時間步長處，DOA估計結(jié)果遠(yuǎn)離真實(shí)軌跡。然而，PF-MUSIC方法能夠結(jié)合來自聲源的動力學(xué)模型的時間信息以及來自當(dāng)前測量時刻的空間信息，從而能夠一致地跟蹤聲源方位。它能始終鎖定聲源，并呈現(xiàn)令人滿意的DOA跟蹤。另外，雖然初始DOA對于PF-MUSIC方法是未知的并且被假定為均勻分布，但PF-MUSIC跟蹤方法能夠快速地收斂到真實(shí)軌跡。

為了充分評估跟蹤性能，采用50次蒙特卡羅平均模擬均方根誤差（RMSE）。如果絕對誤差小于1°，則被認(rèn)為是正確的估計。圖4給出不同SNR下對RMSE的影響。采用從-10 dB～0 dB 的不同 SNR，以 2 dB的增量產(chǎn)生噪聲環(huán)境。

由于結(jié)合了時間信息，所提出的PF-MUSIC跟蹤算法在估計DOA方面比MUISC方位方法執(zhí)行得更好。PF-MUSIC跟蹤算法即使在非常低的SNR環(huán)境下（例如，SNR=-10 dB），依然能夠保持對聲源DOA的鎖定。同時，PF-MUSIC方法不需要對角度進(jìn)行三維搜索，在計算量上遠(yuǎn)低于MUSIC方法，更適用于對目標(biāo)的實(shí)時方位跟蹤的需求。

圖4 不同快拍數(shù)和 SNR 下的 RMSEFig.4 RMSE under the conditions of different snapshots and SNR

3.2 海上試驗(yàn)驗(yàn)證

2017年8月，在大連海域水深60 m區(qū)域，試驗(yàn)驗(yàn)證實(shí)時方位跟蹤方法的性能。信標(biāo)的導(dǎo)航信號為4 kHz的單頻CW脈沖信號，脈沖周期為T=2 s。水下目標(biāo)和導(dǎo)航信標(biāo)正橫距離大于100 m，則聲場可視為平面波聲場。

為驗(yàn)證方位跟蹤方法的實(shí)際性能，這里采用短基線方位方法作為對比方案。分別在目標(biāo)首部和尾部各安裝2個接收器，用以接收導(dǎo)航信標(biāo)信號?；€長度為 53 m。

試驗(yàn)所采用的矢量水聽器為二維球型，典型的矢量水聽器接收的聲壓及振速的功率譜如圖5所示。

短基線和PF-MUSIC方法的跟蹤效果如圖6所示。

可知PF-MUSIC方法的方位跟蹤效果良好，并且與短基線方位相比曲線更平滑，適用于水下目標(biāo)的實(shí)時方位跟蹤。同時，僅需在艦船上加裝單矢量水聽器即可滿足測試要求，更適用于工程實(shí)際。

4 結(jié) 語

將粒子濾波算法和MUSIC方法相結(jié)合，應(yīng)用于矢量水聽器接收的水下運(yùn)動目標(biāo)數(shù)據(jù)的處理，有效地對運(yùn)動目標(biāo)的方位進(jìn)行跟蹤。采用CV模型對聲源運(yùn)動狀態(tài)進(jìn)行建模，將MUSIC方法的估計函數(shù)作為粒子濾波算法中的似然函數(shù)，有效改進(jìn)了在低信噪比條件下的跟蹤性能。通過結(jié)合時間和空間信息，PF-MUSIC跟蹤方法在二維DOA跟蹤中優(yōu)于傳統(tǒng)的MUSIC方法，并且不需要空間搜索降低了計算量，適用于對目標(biāo)的實(shí)時跟蹤，同時在較低信噪比的條件下也能夠獲得較好的精度，對水下目標(biāo)的實(shí)時方位跟蹤上有較好的應(yīng)用前景。

圖5 矢量水聽器接收聲壓和振速信號的頻譜Fig.5 Spectrums of Sound pressure and vibration speed received by the vector hydrophone

圖6 方位跟蹤效果Fig.6 The effect of DOA tracking