基于反卷积波束形成的噪声源识别方法研究进展

丁 浩，傅栋林，张锐之

（浙江工业大学 特种装备制造与先进加工技术教育部/浙江省重点实验室，浙江 杭州310014）

噪声水平是衡量机械产品的重要环保指标，欧美国家建立了大量的相关技术标准作为市场准入的环保条件，这对国产林用机械产品的出口形成了无形的壁垒。为降低机械噪声水平，目前主要是通过声源识别定位方法进行噪声源识别定位，然后采取针对性手段降噪。波束形成法具有计算快、测量方便等优点，对静止或运动的中高频、远距离声源都有很好的识别能力［1］。传统波束形成法（delay and sum beamforming，DAS）不仅在真实声源位置输出具有一定宽度的主瓣，还在非声源位置输出旁瓣，为防止过大的旁瓣导致混淆或湮灭主瓣［2］，如何控制旁瓣水平是当前国内外研究热点［3］。2004年，美国BROOKS等［4］首先提出了反卷积波束形成法（deconvolution approach for the mapping of acoustic sources method，DAMAS），此法显著控制了旁瓣水平，克服了传统波束形成法的局限性，近10 a取得了快速发展。本文从反卷积波束形成法理论原理及现有反卷积波束形成法的对比这2个方面来概述。

1 反卷积算法原理

反卷积波束形成法是一种逆问题求解算法。该算法以传统波束形成作为预处理，得到所有扫描点的均方声压值。假设对流层和剪切层折射并不影响噪声传播到麦克风阵列，可以得出第n个假设声源点传播到第m个麦克风的声压信号矩阵：

式（1）中：m∶n表示第n个声源点到第m个麦克风，e为导向因子。Qn为网格上第n个假设声源点的声压平方矩阵。首先计算：

式（2）中：m′为除m以外的任意麦克风。可得第n个假设声源点的互谱矩阵：

式（4）中：N为总的假设声源点数。可得扫描点的信号功率谱为：

式（5）中： 导向向量e＾=［e1e2… em0］T。 将式（3）和式（4）代入式（5）， 可得：

式（6）中： ［］n′是式（3）中的内容，n′表示聚焦的网格点，n表示假设声源点。简化后，可得：

采用高斯-赛德尔迭代法求解此线性方程得出声源分布。在求解线性方程组时，DAMAS算法能更好地考虑不同网格位置的相互影响，使旁瓣变小甚至消失，呈现更直观的结果［3］，能更准确地计算出声源位置和声源强度。即等于波束形成算法得出的均方声压值，可得：

2 反卷积算法的比较

在算法的旁瓣抑制能力方面，反卷积波束形成法相比于传统波束形成法具有更好的旁瓣抑制能力。目前，声源识别领域各主要算法都已经进行了大量研究对比［5,7］，研究结果表明：频率的提高可有效改善算法的旁瓣抑制能力，且传统DAS算法的旁瓣抑制能力弱于其他算法。当声源远离中心时，反卷积波束形成法2（DAMAS2）和基于傅里叶变化的非负最小二乘算法（FFT-NNLS）算法定位声源失效。这主要由于DAMAS2和FFT-NNLS都基于阵列点扩散函数位移不变性假设，当声源远离中心时，其声源位置已位于有效区域外，因此DAMAS2和FFT-NNLS不能对该类远离中心的声源点进行精准定位。XENAKI等［8］对这类算法的有效区域进行了深入研究，有效区域一般为Z轴18度角范围内。在相干声源情况下，基于空间相干的洁净算法（CLEAN-SC）不能准确识别声源，这是由于CLEAN-SC在迭代中删除与波峰相干的声源部分，当声源较靠近时，其中一个声源被当作旁瓣而被消除。杨洋等［6］也在研究中指出CLEAN-SC不适用于相干声源。

在算法的定位精度方面，迭代算法随着迭代次数的增加，算法的精度都有不同程度的提升，分析算法迭代次数与标准差的关系［7］。在低迭代次数时洁净算法（CLEAN）和CLEAN-SC这2种算法的精度较高，具有较高的计算精度稳定性，而DAMAS2，非负最小二乘算法（NNLS）和FFT-NNLS具有较大偏差，且对迭代次数不敏感，迭代次数的增加不能大幅提升其计算精度，DAMAS在高迭代次数下声源识别性能快速迫近CLEAN和CLEAN-SC，对于迭代次数较敏感。空间分辨率的优劣也是算法精度的重要指标,算法分辨率越高，区分较近声源的能力就越高。通过对风机的噪声分布研究表明［9］：DAMAS具有很高的分辨率，可以获得出色的噪声源分布图，即使声源距离较近，也能够很好地将它们区分，CLEAN-SC也具有很高的分辨率，但由于其算法的理论原因，较近声源，只能识别其中一个，而DAS算法分辨率较差，无法区分较近声源，定位精度降低。

计算效率是算法能否被广泛运用的关键点，表1是多种典型算法的计算时间，可得出DAS，DAMAS2，和CLEAN-SC具有较高的计算效率和实时性，其次为DAMAS、稀疏约束反卷积波束形成法（SC-DAMAS），但协方差矩阵拟合法（CMF）、稀疏约束的超分辨率反卷积波束形成算法（SC-RDAMAS）计算时间较长，不适用于实时测试。

表1 算法计算时间Table 1 Computation time of algorithms

综上所述，主要代表性反卷积算法的特点如表2所示。

表2 主要代表性反卷积算法Table 2 Main representative algorithm

3 总结与展望

当前，反卷积波束形成算法在高频下都具有优秀的旁瓣抑制能力，且定位精度随着迭代次数的上升而提高，但这些算法也存在着各自的局限性，如在计算效率，实时性，扫描范围，适用声源类型等方面，因此在现有算法的基础上，需进一步提出更加全面的反卷积波束形成改进算法，使算法具有更好的普遍适用价值。同时，要提高声源识别的准确性，除了选择合适的算法，开发声源识别性能更优的传声器阵列也是一研究热点。布置形式更加优化、更加合理的传声器，使传声器阵列声源识别的空间分辨率更高，最大旁瓣水平更低，有效动态范围更大［16］。

［1］ 陈心昭.噪声源识别技术的进展［J］.合肥工业大学学报（自然科学版），2009，32（5）：609-614.CHEN Xinzhao.Progress of techniques for noise source identification ［J］.J Hefei Univ Technol Nat Sci,2009,32（5）：609-614.

［2］ YARDIBI T,BAHR C,ZAWODNY N,et al.Uncertainty analysis of the standard delay-and-sum beamformer and array calibration ［J］.J Sound Vibr,2010,329（13）：2654-2682.

［3］ 杨洋，褚志刚，江洪，等.反卷积DAMAS2波束形成声源识别研究［J］.仪器仪表学报，2013，34（8）：1779-1786.YANG Yang,CHU Zhigang,JIANG Hong,et al.Research on DAMAS2 beamforming sound source identification ［J］.Chin J Sci Instrum,2013,34（8）：1779-1786.

［4］ BROOKS T F,HUMPHREYS W M.A deconvolution approach for the mapping of acoustic sources （DAMAS） determined from phased microphone arrays ［J］.J Sound Vibr,2006,294（4）：856-879.

［5］ YARDIBI T,ZAWODNY N S,BAHR C,et al.Comparison of microphone array processing techniques for aeroacoustic measurements ［J］.Int J Aeroacoust,2010,9（6）：733-762.

［6］ 杨洋，褚志刚.基于CLEAN-SC清晰化波束形成的汽车前围板隔声薄弱部位识别［J］.声学技术，2015，34（5）： 449-456.YANG Yang,CHU Zhigang.Weak position identification of sound insulation for car dash panel based on CLEAN-SC clearness beamforming ［J］.Tech Acoust,2015,34（5）：449-456.

［7］ CHU Zhigang,YANG Yang.Comparison of deconvolution methods for the visualization of acoustic sources based on cross-spectral imaging function beamforming ［J］.Mech Syst Sign Proc,2014,48（1/2）：404-422.doi：org/10.1016/j.ymssp.2014.03.012.

［8］ XENAKI A,JACOBSEN F,TIANA-ROIG E,et al.Improving the resolution of beamforming measurements on wind turbines ［C］//ICA.Proceedings of 20th International Congress on Acoustics.Sydney：ICA,2010：23-27.

［9］ RAMACHANDRAN R C,RAMAN G,DOUGHERTY R P.Wind turbine noise measurement using a compact microphone array with advanced deconvolution algorithms ［J］.J Sound Vibr,2014,333（14）：3058-3080.

［10］ CHU Ning,PICHERAL J,MOHAMMAD-DJAFARI A.A robust super-resolution approach with sparsity constraint for near-field wideband acoustic imaging ［C］//IEEE Computer Socirty.ISSPIT’11 Proceedings of the 2011 IEEE International Symposium on Signal Processing and Information Technology.Washington D C：IEEE Computer Society,2011：310-315.doi：10.1109/ISSPIT.2011.6151579.

［11］ DOUGHERTY R.Extensions of DAMAS and benefits and limitations of deconvolution in beamforming ［C］//AIAA Association.11th AIAA/CEAS Aeroacoustics Conference.Monterey：AIAA Association,2005.doi：org/10.2514/6.2005-2961.

［12］ EHRENFRIED K,KOOP L.Comparison of iterative deconvolution algorithms for the mapping of acoustic sources［J］.AIAA J,2007,45（7）：1584-1595.doi：org/10.2514/1.26320.

［13］ YARDIBI T,LI Jian,STOICA P,et al.Sparsity constrained deconvolution approaches for acoustic source mapping［J］.J Acoust Soc Am,2008,123（5）：2631-2642.

［14］ SIJTSMA P.CLEAN based on spatial source coherence ［J］.Int J Aeroacoust,2007,6（4）：357-374.

［15］ RÜBSAMEN M,PESAVENTO M.Steering vector non-identifiability in covariance matrix fitting based beamforming［C］//IEEE Computer Socirty.Proceedings of the IEEE 7th Sensor Array and Multichannel Signal Processing Workshop （SAM）.Hoboken：IEEE Computer Socirty,2012：433-436.

［16］ 刘鹏，刘志红，李贺，等.传声器阵列可视化效果影响因素研究［J］.测试技术学报,2016,30（3）：231-235.LIU Peng,LIU Zhihong,LI He,et al.Study on the influencing factors of microphone array visualization ［J］.J Test Measyr Technol,2016,30（3）：231-235.