• CN:11-2187/TH
  • ISSN:0577-6686

机械工程学报 ›› 2020, Vol. 56 ›› Issue (22): 46-55.doi: 10.3901/JME.2020.22.046

• 仪器科学与技术 • 上一篇    下一篇

扫码分享

二维动态网格压缩波束形成声源识别方法

樊小鹏1, 余立超2,3, 褚志刚2,3, 杨洋2,3, 李丽1   

  1. 1. 广东电网有限责任公司电力科学研究院 广州 510082;
    2. 重庆大学机械传动国家重点实验室 重庆 400044;
    3. 重庆大学汽车工程学院 重庆 400044
  • 收稿日期:2020-02-22 修回日期:2020-09-08 出版日期:2020-11-20 发布日期:2020-12-31
  • 通讯作者: 褚志刚(通信作者),男,1978年出生,博士,教授,博士研究生导师。主要研究方向为振动噪声测量分析技术,数字信号处理,噪声源识别技术。E-mail:zgchu@cqu.edu.cn
  • 作者简介:樊小鹏,男,1986年出生,博士。主要研究方向为电力系统噪声环境保护、电气设备故障诊断。E-mail:283128493@qq.com;余立超,男,1994年出生。主要研究方向为阵列声源识别技术及其应用。E-mail:810177064@qq.com;杨洋,女,1988年出生,博士,讲师。主要研究方向为噪声源识别技术理论及其应用,工程信号处理。E-mail:yangyang911127@cqu.edu.cn;李丽,女,1971年出生,博士,教授级高级工程师。主要研究方向为电力系统噪声环境保护、电气设备故障诊断。E-mail:liligz@163.com
  • 基金资助:
    国家自然科学基金(11874096,11704040)和中国南方电网科技(GDKJXM20180152)资助项目。

Two-dimensional Dynamic Grid Compressive Beamforming for Acoustic Source Identification

FAN Xiaopeng1, YU Lichao2,3, CHU Zhigang2,3, YANG Yang2,3, LI Li1   

  1. 1. Electric Power Research Institute of Guangdong Power Grid Limited Liability Corporation, Guangzhou 510082;
    2. State Key Laboratory of Mechanical Transmissions, Chongqing University, Chongqing 400044;
    3. College of Automotive Engineering, Chongqing University, Chongqing 400044
  • Received:2020-02-22 Revised:2020-09-08 Online:2020-11-20 Published:2020-12-31

摘要: 基于平面传声器阵列测量和固定网格的传统二维压缩波束形成通过建立阵列传声器测量声压信号和假想声源网格点未知源强之间的欠定线性方程组,基于主声源通常为稀疏分布这一事实,利用稀疏促进算法求解上述方程组从而获得声源波达方向和源强的准确估计,进而准确识别声源。然而,当声源波达方向与网格点不一致、即发生基不匹配时,性能会劣化。为克服该问题,提出二维动态网格压缩波束形成声源识别方法。定义网格坐标和源强分布矢量为变量并采用对数求和罚函数构建目标函数以促进解的稀疏性;基于优化最小化框架在目标函数的基础上构造合适的替代函数以降低优化复杂度;通过梯度下降法对替代函数进行迭代优化求解,从而使网格坐标和源强分布矢量逐渐收敛至真实值附近。仿真和试验均表明,相较于传统固定网格的二维压缩波束形成,该方法能够克服基不匹配问题、获得更高的定位精度和量化精度;该方法能够适用于传声器随机分布的平面阵列,且无需先验的信噪比(噪声干扰)及声源稀疏度等信息,即使在传声器数量较少的情况下也能得到低污染高分辨率的声源成像,保证了高精度的二维波达方向估计和源强量化。

关键词: 压缩波束形成, 动态网格, 平面传声器阵列, 声源识别, 波达方向

Abstract: Conventional two-dimensional compressive beamforming, which is based on planar microphone array measurements and fixed discrete grids, establishes an underdetermined linear system of equations between the sound pressure signals measured by microphones and the unknown source strengths of the assuming acoustic sources corresponding to grids. Based on the fact that the main sources are generally sparse, the above equations can be solved by the sparsity promotion algorithms, then accurate estimation of the source direction-of-arrivals(DOAs) and quantification of the source strengths are achieved, thus identifying sources accurately. However, its performance deteriorates when sources do not coincide with the grids, namely the basis mismatch occurs. To solve this issue, a two-dimensional dynamic grid compressive beamforming for acoustic source identification is developed. First, the grid coordinates and the source strength distribution vector are defined as variables, and the objective function is constructed using a log-sum penalty function to promote the sparsity of the solution. Subsequently, a suitable surrogate function is formulated based on the objective function to reduce the optimization complexity under the majorization-minimization framework. Finally, a gradient descent method is used to iteratively optimize the surrogate function, leading to a gradual process to refine the grid coordinates and the source strength distribution vector. The results of numerical simulations and experiments demonstrate that the proposed technology can circumvent the basis mismatch and thus achieve higher location and quantification accuracy, comparing to the conventional two-dimensional fix-grid approach. It can be applied to a planar array with microphones randomly distributing and does not require prior knowledge of signal-to-noise ratio (noise interference) or source sparsity. The dynamic grid compressive beamforming can provide high-resolution and low contamination imaging, allowing accurate estimation of two-dimensional DOAs and quantification of source strengths, even with a small number of microphones.

Key words: compressive beamforming, dynamic grid, planar microphone array, sound source identification, direction-of-arrival

中图分类号: