羅光成 杜馬千里 茆琳 李智忠

（1.91001部隊(duì)，北京，100036；2.92196部隊(duì)，青島，262100）

（3.海軍潛艇學(xué)院，青島，262199）

大間隔非等間距陣 DOA估計(jì)技術(shù)可應(yīng)用在以下場(chǎng)景中：多個(gè)聲吶浮標(biāo)節(jié)點(diǎn)組合成陣探測(cè)，水下固定式聲吶在同樣范圍內(nèi)使用盡可能少的陣元以降低成本或部分陣元損壞等。當(dāng)陣元間距與所要探測(cè)信號(hào)的波長(zhǎng)不滿足空間奈奎斯特采樣定律，即陣元間距大于二分之一波長(zhǎng)時(shí)（定義為大間隔），經(jīng)典的基于二次型的波束搜索算法，如CBF、MVDR等易出現(xiàn)柵瓣虛警現(xiàn)象，因而大間隔非等間距陣 DOA估計(jì)技術(shù)需具備柵瓣抑制能力。作為新的采樣理論，壓縮感知（或者壓縮采樣）可通過(guò)開發(fā)信號(hào)的稀疏特性，在遠(yuǎn)小于Nyquist采樣率的條件下，用隨機(jī)采樣獲取信號(hào)的離散樣本，然后通過(guò)非線性重建算法完美地重建信號(hào)[1]?；谶@一思想，本文擬利用稀疏空域信號(hào)重構(gòu)方法開展大間隔不等間距陣的 DOA估計(jì)技術(shù)研究。

稀疏信號(hào)重構(gòu)技術(shù)近些年來(lái)得到快速發(fā)展，并被廣泛應(yīng)用到各個(gè)領(lǐng)域，包括圖像重建與恢復(fù)、小波降噪、雷達(dá)成像、聲吶目標(biāo)定位等[2-5]，在頻譜估計(jì)和陣列處理的背景下也出現(xiàn)了一些新的算法，包括 l1范數(shù)最小化[6]、匹配追蹤[7]、凸優(yōu)化方法[8]等。比較經(jīng)典的是Gorodnitsky等人利用迭代最小范數(shù)加權(quán)稱為FOCUSS進(jìn)行 DOA估計(jì)[9]，后來(lái)Chen提出了基于追蹤準(zhǔn)則的稀疏信號(hào)估計(jì)方法 (Basis Pursuit Denoising ，BPDN)[10]，通過(guò)最小化空間信號(hào)的l1范數(shù)進(jìn)行稀疏性約束，結(jié)合噪聲功率限制的陣列數(shù)據(jù)的稀疏重構(gòu)擬合約束，利用凸優(yōu)化工具獲得稀疏的 DOA信號(hào)估計(jì)結(jié)果。

但是上述這些方法存在諸如不適用于多快拍相干處理、超參數(shù)的選擇沒(méi)有規(guī)則、計(jì)算量大無(wú)法實(shí)際應(yīng)用等一系列問(wèn)題。本文基于文獻(xiàn)[11]利用SBL，假設(shè)聲源幅度、接收陣元信號(hào)均滿足復(fù)高斯分布，基于證據(jù)最大推導(dǎo)求解分布參數(shù)，提出了一種使用由多個(gè)測(cè)量快拍組成協(xié)方差矩陣，進(jìn)行大間距陣 DOA估計(jì)。該方法可以在高分辨 DOA估計(jì)時(shí)帶來(lái)多項(xiàng)益處，包括：（1）能對(duì)相干信號(hào)進(jìn)行 DOA估計(jì)；（2）相對(duì)于基于特征分解的波束形成方法可以使用任意數(shù)量的快拍數(shù)；（3）可實(shí)現(xiàn)在線迭代算法；（4）具有比MVDR、MUSIC等經(jīng)典算法更高的方位分辨率。

1 算法原理

假設(shè)傳感器個(gè)數(shù)為N，采樣快拍數(shù)為L(zhǎng)，暫時(shí)只考慮窄帶信號(hào)，則得到陣列接收頻域信號(hào)

2 仿真驗(yàn)證

仿真參數(shù)：陣元位置坐標(biāo)如圖1所示，5個(gè)陣元非等間距近似布置在一條直線上，總長(zhǎng)255 m，平均陣間距對(duì)應(yīng)的半波長(zhǎng)頻率約為11.8 Hz，假設(shè)聲源頻率64 Hz，不滿足奈奎斯特采樣定理，因而存在柵瓣模糊現(xiàn)象。假設(shè)目標(biāo)方位 90°，輸入信噪比 6 dB，估計(jì)協(xié)方差矩陣所用快拍數(shù)為 40個(gè)，進(jìn)行了100次蒙特卡洛仿真。結(jié)果如圖2、3所示。

圖1 陣元位置示意圖

從圖2(a)、(b)的處理結(jié)果可以看出，CBF和MVDR都在60°和120°方向上出現(xiàn)明顯的柵瓣，相對(duì)CBF，MVDR具備一定的柵瓣抑制能力，而圖2(c)SBL結(jié)果中無(wú)明顯的柵瓣軌跡，說(shuō)明稀疏處理能夠較好的抑制柵瓣模糊。

圖2 三種算法仿真結(jié)果

3 海試數(shù)據(jù)驗(yàn)證

數(shù)據(jù)說(shuō)明：陣元位置同圖1，存在兩個(gè)目標(biāo)。處理頻率64 Hz，2 s產(chǎn)生一個(gè)頻域快拍，50個(gè)快拍組成采樣協(xié)方差矩陣，20 min數(shù)據(jù)長(zhǎng)度，不同算法的處理結(jié)果如圖3所示。

對(duì)比不同算法的海試數(shù)據(jù)處理結(jié)果，可以看到：圖3(a)、(b)中CBF和MVDR雖然有較高的目標(biāo)1和目標(biāo)2的空間能量譜輸出，但是如果沒(méi)有聲源數(shù)目的先驗(yàn)知識(shí)，很難排除柵瓣的真實(shí)性，從圖3(c)可以看到，SBL算法具有柵瓣抑制能力，能夠得到更清晰軌跡。

圖3 三種算法仿真結(jié)果

4 總結(jié)

本文將空域稀疏處理應(yīng)用于大間距陣 DOA估計(jì)，提出了基于SBL算法的大間隔非等間距陣DOA 估計(jì)技術(shù)，利用仿真和實(shí)際數(shù)據(jù)對(duì)多種算法進(jìn)行了對(duì)比分析，最后得出本文所提算法具有最佳的柵瓣抑制效果。實(shí)際工程應(yīng)用需考慮陣元位置誤差、快拍間相干性以及處理帶寬等其它更多因素的影響，需要開展更加全面深入的研究。