矢量圆阵时域解析 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–)，男，硕士，副教授，研究方向为水声目标识别。