矢量圓陣時域解析 MVDR 算法研究

王易川，李海濤，高 鑫，陳 鵬

(1. 海軍潛艇學院，山東 青島 266042；2. 海軍青島雷達聲吶修理廠，山東 青島 266000)

矢量圓陣時域解析 MVDR 算法研究

王易川1，李海濤1，高 鑫1，陳 鵬2

(1. 海軍潛艇學院，山東 青島 266042；2. 海軍青島雷達聲吶修理廠，山東 青島 266000)

為提高矢量圓陣的方位估計性能，將 VTAMVDR（MVDR based on Time-domain Aanalysis signals of Vector Sensor）應用在矢量圓陣中，提出矢量圓陣時域解析 MVDR 算法，研究了矢量圓陣時域解析 MVDR 算法的原理及實現流程，然后通過仿真分析對比矢量圓陣時域解析 MVDR 算法與矢量圓陣常規波束形成的波束寬度、方位估計性能和分辨率。仿真結果表明，矢量圓陣時域解析 MVDR 算法方位估計性能和多目標分辨能力優于常規波束形成算法。

矢量圓陣；時域解析；波束形成

0 引 言

矢量水聽器可以同時測量聲場的聲壓和質點振速信息的傳感器，通過聯合處理聲壓和振速得到不同的組合指向性，并為最簡單的線列陣帶來了分辨左/右舷模糊的能力[1–3]。目前對于矢量線陣的空間譜估計已經有較多的文獻發表，而對于矢量圓陣的研究較少[4–6]。圓陣和直線陣相比有諸多優點，可以進行 360°無模糊測向。在矢量線陣空間譜估計研究方面，文獻[7]提出了 VTAMVDR（MVDR based on Time-domain Aanalysis signals of Vector Sensor）空間譜估計方法，VTAMVDR的基本思想是應用 TAMVDR（MVDR based on Timedomain Aanalysis signals）的理念將接收到的聲壓振速時域實數信號通過希爾伯特變換構造出復數形式的解析信號，再利用線性約束條件使陣列輸出功率最小的準則確定復數形式的最優權向量，從而獲得來波方向的波束輸出功率。在構造時域解析信號的基礎上，對每一路時域解析信號引入復權，利用期望方向的輸出相應不變的約束使陣列輸出功率最小而獲得權向量的最優解。利用 TAMVDR，可在單次快拍下獲得穩定優化解，信號檢測和方位分辨性能較優。相對于其他算法需要很多快拍數據才能收斂而言，該算法相對降低了對條件的依賴，同時運算量相對較小。本文研究了VTAMVDR 在矢量圓陣中的應用。首先從矢量圓陣的陣列模型出發，提出了矢量圓陣時域解析 MVDR 算法，研究了矢量圓陣時域解析 MVDR 算法原理及實現流程，然后通過仿真分析對比了矢量圓陣時域解析MVDR 算法與矢量圓陣常規波束形成的波束寬度和分辨率。

1 矢量圓陣時域解析 MVDR 算法原理

圖 1 給出了一個在圓周上均勻分布的平面離散圓陣圖，設離散圓陣的半徑為 r，陣元為 N 個二維矢量水聽器，把時間的參考點選在圓心 O 上，設入射信號來自 0°方向，相鄰 2 個陣元的夾角為 α =2π/N，假設圓陣每次波束掃描時，有 2M 個基元參與工作，記 2M個矢量水聽器 p，vx，vy接收信號經時延以及濾波之后的輸入信號分別為：

式中：Sp，Sx，Sy為期望目標信號矢量；Vp，Vp，Vy為陣列加性各向同性噪聲矢量與干擾之和。

利用 Hilbert 變換 H (·) 構造時域復數信號為：

陣列輸出功率為：

式中

對權向量的求解可以表述為有約束的優化問題：

式中： A=[1,1,···,1]T為 2M 個 1 組成的列矩陣；θ 為掃描方位。

先構造1個代價函數：

上式對 w 求導，并令導數為0，得到：

其輸出功率為：

圖 2 給出了矢量圓陣時域解析自適應波束形成流程圖。

從圖 2 可看出，矢量圓陣時域解析 MVDR 完全基于時域處理，在數據長度滿足條件時，一般只需要1次快拍，不需要進行子帶分解，減小了計算量。當陣元個數較多時，能大大減小運算量，方便使用 DSP 實現矢量水聽器陣列的工程化應用。

2 數值仿真分析

均勻圓陣，半徑為 0.8 m，由 32 個矢量水聽器陣元組成，采樣率為 40 kHz，處理頻帶范圍 [2 kHz，4 kHz]，1次快拍點數為 1 024 點，目標在觀測時間內方位不變，信噪比為 5 dB，加性噪聲為高斯白噪聲，采用一次快拍，目標在觀測時間內方位不變，加性噪聲為高斯白噪聲。考察矢量圓陣常規波束形成和矢量圓陣時域解析自適應波束形成多次方位估計的均方根誤差（root mean square error, RMSE）[8]。

定義 N 次目標方位估計的 RMSE 為：

其中 θij和ij分別為第 i 個目標第 j 次試驗的真值和估計值。

圖 3 給出了單目標條件下矢量圓陣常規波束形成和矢量圓陣時域解析MVDR算法的波束輸出對比。

從圖 3 可看出，矢量圓陣常規波束形成下降 3 dB，波束寬度為 10.87°，本文提出的矢量圓陣時域解析MVDR 算法為 1.11°，明顯波束寬度較窄，本文算法的波束寬度為常規波束形成的 10.21%。

圖 4 給出了 2 種方法 RMSE 隨信噪比變化示意圖，可以看出隨著信噪比的變化，矢量圓陣常規波束形成方法的性能弱于矢量圓陣時域解析 MVDR 算法。

仿真條件不變，有 2 個寬帶信號入射，分別位于90°和 98°方位，目標在觀測時間內方位不變，信噪比為 5 dB，加性噪聲為高斯白噪聲，假設估計的兩目標方位為和，如果滿足下式：

則認為是一次成功雙目標分辨。

圖 5 給出了雙目標條件下矢量圓陣常規波束形成和矢量圓陣時域解析自適應波束形成的波束輸出，可以看出矢量圓陣時域解析自適應波束形成可以明顯的分辨 2 個目標而常規波束形成不能分辨。

圖 6 給出了不同角度間隔下 2 種方法的分辨能力對比圖，可以看出本文方法與常規波束形成方法相比分辨性能較優。

3 結 語

本文提出了矢量圓陣時域解析 MVDR 算法，理論和仿真分析了本文方法的性能，通過和常規波束形成進行對比，可以看出本文方法在方位估計性能和目標分辨能力上明顯優于常規波束形成，具有較高的估計精度，較低的雙目標分辨門限。

[1]惠俊英, 惠娟. 矢量聲信號處理基礎[M]. 北京:國防工業出版社, 2009.

[2]孫貴青, 張春華, 黃海寧, 等. 聲矢量傳感器線陣的左右舷分辨[J]. 哈爾濱工程大學學報, 2010, 31(7):848–855. SUN Gui-qing, ZHANG Chun-hua, HUANG Hai-ning, et al. Left-right resolution of acoustic vector senor line arrays[J]. Journal of Harbin Engineering University. 2010, 31(7):848–855.

[3]楊秀庭. 矢量水聽器及其在聲吶系統中的應用研究[D]. 北京：中國科學院, 2009.

[4]楊德森,朱中銳,時勝國. 矢量圓陣測向方法[J]. 哈爾濱工程大學學報, 2012, 33(10):1259–1264. YANG De-sen, ZHU Zhong-rui, SHI Sheng-guo. Directon-ofarrival estimation for a circular vector-sensor array[J]. Journal of Harbin Engineering University, 2012, 33(10):1259–1264.

[5]康曉濤, 王志洋, 康博宇, 等. 四元數在均勻圓形矢量傳感器陣列信號參數估計中的應用[J].吉林大學學報 (工學版). 2013, 42(S1):154–159. KANG Xiao-tao, WANG Zhi-yang, KANG Bo-yu, et al. Application of quaternion vector sensor in the uniform circular array signal parameter estimation[J]. Journal of Jiin University:Enginerring Science Edition. 2013, 42(S1):154–159.

[6]李曉青, 趙亞衛, 李新波, 等. 均勻圓型聲矢量陣的四元數二維波達方向估計[J].吉林大學學報. 2014, 32(5):451–457. LI Xiao-qing, ZHAO Ya-wei, LI Xin-bo, et al. Twodimensional DOA estimation using quaternion based on UCA of acoustic vector sensor[J]. Journal of Jiin University. 2014, 32(5):451–457.

[7]李海濤, 李智忠, 宋志杰, 等. 一種新的矢量線列陣波束形成方法研究[J] 傳感器與微系統, 2015, 34(10):144–150. LI Hai-tao, LI Zhi-zhong, SONG Zhi-jie, et al. A new algorithm for vector sensor line array beamforming[J]. Transducer and Microsystem Technologies, 2015, 34(10):144–150.

[8]梁國龍, 馬巍, 范展, 等. 矢量聲吶高速運動目標穩健高分辨方位估計[J]. 物理學報, 2013, 62(14):282-290. LIANG Guo-long, MA Wei, FAN Zhan, et al. A high resolution rbust localization approach of high speed target based on vector sonar[J].Acta Phys.Sin. 2013, 62(14):282-290.

Studying of circular vector sensor array time-domain analysis mvdr algorithm

WANG Yi-chuan1, LI Hai-tao1, GAO Xin1, CHEN Peng2
(1. Navy Submarine Academy, Qingdao 266042, China; 2. Navy Qingdao Radar and Sonar Repair Factory, Qingdao 266000, China)

To improve the resolution of direction estimation of circular vector sensor array, a new algorithm based on VTAMVDR is used in circular vector sensor array signal processing. Circular vector sensor array time-domain analysis MVDR algorithm is proposed. Studied the principle and process flow of circular vector sensor array time-domain analysis MVDR algorithm. Compared the beam width and resolution capability of circular vector sensor array time-domain analysis MVDR with circular vector sensor conventional beamforming algorithm. Theoretical analysis and computer simulations show that the proposed algorithm is more effective in improving resolution of multiple targets than the traditional algorithm

circular vector sensor array；time-domain analysis；beamforming

TN911.7

A

1672–7619(2017)03–0131–04

10.3404/j.issn.1672–7619.2017.03.027

2016–07–18；

2016–10–08

國家自然科學基金資助項目（60901054）

王易川(1976–)，男，碩士，副教授，研究方向為水聲目標識別。