舱室空间内部声源定位是声学定位研究中的一项难题，同时也是一项应用需求广泛的技术，如飞机舱室噪声源定位，会议室智能话筒传音控制系统，车载电话系统，以及各种室内声学跟踪或监控系统。封闭空间的主要特点是结构复杂、混响程度较高，小尺度空间常常还有声波干涉、衍射效应，这对声源定位技术的精度和稳定性都提出了更高的要求。
目前，常用的封闭空间声源定位方法有基于麦克风阵列的方法和基于双耳听觉机理的双传声器方法。前者主要包括可控波束形成，高分辨率谱估计以及到达时间差。这些方法及其不断更新的算法与特定的传声器阵列结合，可获得良好的定位精度[1, 2, 3]，因而已在大区域开放式环境中得到了广泛的应用。但是这些方法存在算法复杂、阵型需专门设计、经济成本较高以及在封闭空间中容易受到噪声和混响的影响等问题[4, 5, 6, 7, 8]，因此往往难以直接应用于中小尺度的舱室空间。
针对上述问题，本文提出了一种新的封闭空间声源定位方法，即基于双通道空间脉冲响应匹配的声源定位方法，并通过实际房间中的定位实验，对其可行性及主要影响因素进行分析和研究。
1 基于空间脉冲响应匹配的声源定位方法 1.1 基本原理声场空间脉冲响应(SIR)是指声场中接收位置收到的由脉冲声源辐射的信号序列。同一房间内，当声源与接收点的位置固定时，可以认为二者之间的传播通道是唯一的，即脉冲响应是唯一的确定函数，它包含了封闭声场的所有声学特性，可以准确地反映声学环境(方向、距离、空间感等)[9]。SIR由直达声、前次反射声和混响声组成[10]，其形成过程和结果分别如图 1和图 2所示。目前对空间脉冲响应的研究主要在脉冲响应的测量方法，室内和双耳脉冲响应的模拟，脉冲响应各阶段的分析，音乐厅等的音质评估和设计等方面。
如果事先能够通过计算或测量获得该声场中任意声源位置到接收点的SIR，则可通过双传声器接收的信号计算得到声源至2个传声器通道响应之间的相互关系，再与数据库里的各个采样声源对应的参数进行比对，即可确定真实声源的位置。
|
| 图 1 空间脉冲响应的形成过程 |
|
| 图 2 空间脉冲响应 |
基于上述原理，具体算法设计如下：
设未知声源S发射的信号s(t)为单频信号，hα、hβ分别为未知声源S到双传声器的空间脉冲响应，忽略本底噪声的影响，则双传声器接收到的信号分别为
对(1)式和(2)式进行傅里叶变换得
f为声源信号的频率。由于声源信号s(t)未知，为消去其影响，用(3)式除以(4)式，得
从(5)式可以看出，实际定位时并不需要知道声源信号s(t)，而只需知道双通道传递函数的比例关系。
因此，可以对实际房间中声源可能分布的区域进行网格划分，并测量所有网格中心(假定声源位置)至固定传声器位置的SIR，将它们作傅里叶变换后得G(f)=((G1α(f),G1β(f),(G2α(f),G2β(f),…,(Gnα(f),Gnβ(f))),令
为划分的网格数。
将Φ,Ψ作互相关运算，即
式中，corrcoef表示皮尔逊相关系数运算，
表示平均值。
将基于实测声信号计算得到的Φ与数据库中的所有Ψ作互相关，计算出R1,R2,…,Rn，找出相关系数最大的1组数据Ri，其中位置编号
定义为：
认为声源位置S位于该组数据所对应的网格中心点，从而实现未知声源的定位。
选取一个普通房间(如图 3所示)。
|
| 图 3 实验房间 |
建立数据库时声源分布网格如图 4所示，所有网格处于相同的高度。脉冲响应采样装置和声源定位测试装置连接图如图 5所示。在房间内均匀划分m×m=n个网格，网格间距为d。房间大小为8.3 m×5.5 m×3.8 m。信噪比大约为5 dB。两传声器分别位于Rα(4.2,4,3.6)，Rβ(5.3,2.8,3.6)处。采样时用图 5里的采样装置，声源高度固定，利用DIRAC测量系统得到声源到双传声器的脉冲响应并存入数据库。定位时用图 5里的实测装置，测试声源采用时长约2 s的信号，在每个网格中心进行5次重复定位，结果出现最多的位置即认为是未知声源的位置，定位流程如图 6所示。
|
| 图 4 实验房间示意图 |
|
| 图 5 实验装置图 |
|
| 图 6 定位流程图 |
在图 3所示的房间中按图 6中的实验步骤进行实验验证。n取100，分别采用直径为50 cm的球形声源，直径为14 cm的扬声器和直径为7 cm的扬声器作为发声装置，受到声源尺寸的限制，网格划分间隔(对应于空间分辨率)也有所不同，统计所得的100个网格的定位实验结果见表 1。
由表 1可以看出，发声装置的尺寸越小，空间分辨率越高，因为发声装置本身占据一定的空间，当其尺寸足够小，可认为是点声源时，空间分辨率可达到最高。不过，几种情况下的定位准确率都维持在很高的水平。
3.2 混响环境下的定位性能研究为了进一步探讨不同混响程度下的定位准确率，分别在混响时间T60=0.438 s、1.769 s、4.373 s的消声室、1号普通房间以及2号普通房间中进行定位实验，所有房间划分网格数n均取25，网格间距d均取20 cm，发声装置采用直径14 cm的扬声器，对25个位置重复定位125次，统计得到的定位准确率均为100%。这表明，在不同的混响环境中，该方法均具有较高的定位准确率，即受混响影响小。
3.3 不同信噪比条件下的定位性能研究网格数n均取25，网格间距d均取20 cm，发声装置采用直径14 cm的扬声器，实验房间如图 6所示。通过在传声器附近播放白噪声和改变测试声源的音量来获得不同的信噪比SNR，分别选取SNR约为0 dB、1 dB、2 dB、3 dB、4 dB、5 dB、6 dB时按图 6所示的实验步骤进行实验，得到不同信噪比下25个网格的定位准确率如图 7所示。
|
| 图 7 不同信噪比下的定位准确率 |
从图 7可以看出，背景噪声对该方法影响显著，当信噪比低于2 dB时，因不满足辐射声能的要求，定位准确率不到10%，信噪比高于2 dB后，定位准确率均超过90%。这表明，该方法在信噪比满足一定的数值时可以维持很高的定位精度。
3.4 缓慢移动声源的定位研究室内环境中的声源一般移动缓慢，可以忽略多普勒效应。仍然在图 3所示的房间内进行实验，声源以50 cm/s的速度按不同路线移动，其中1条移动轨迹如图 8所示。表 2给出了按2种不同路线移动时统计的定位结果。
|
| 图 8 采样点及声源移动路线示意图 |
| 组别 | 移动方向 | 定位结果 |
| 移动方向1 | (1-6-11-16-21) | (21-16-11-6-1) |
| 1,6,11,16,21 | 21,16,11,6,1 | |
| 移动方向2 | (21-22-23-24-25) | (25-24-23-22-21) |
| 21,22,23,24,25 | 25,24,23,22,21 |
结合图 8和表 2结果可以看出，当声源缓慢移动时，跟踪定位结果准确率为100%，说明该方法在声源缓慢移动的情况下也可以获得很高的定位准确率。
4 结 论本文基于空间中声源到接收点脉冲响应的唯一性，提出了一种新的舱室声源定位方法，即基于双通道空间脉冲响应匹配的声源定位方法，并研究了该方法的空间分辨率以及不同的混响、信噪比和声源移动情况下的定位准确率。该方法计算量小，兼容性高，应用前景广泛。
| [1] | Huang F, Sheng W, Ma X. Modified Projection Approach for Robust Adaptive Array Beamforming[J]. Signal Processing, 2012, 10:1758-1763 |
| Click to display the text | |
| [2] |
王永良,陈辉,彭应宁,等. 空间谱估计理论与算法[M]. 北京:清华大学出版社, 2004 Wang Yongliang, Chen Hui, Peng Yingning, et al. Spatial Spectral Estimation Theory and Algorithm[M]. Beijing:Tsinghua University Press, 2004(in Chinese) |
| Cited By in Cnki (1) | |
| [3] | El Korso M N, Boyer R, Renaux A, et al. Conditional and Unconditional Cramer-Rao Bounds for Near-Field Source Localization[J]. IEEE Trans on Signal Processing, 2010, 58(5):2901-2907 |
| Click to display the text | |
| [4] |
朱铮宇. 基于传声器阵列的声源定位技术研究[D]. 西安:陕西师范大学, 2011 Zhu Zhengyu. Study on Sound Source Positioning Technology Based on Microphone Array[D]. Xi'an, Shaanxi Normal University, 2011(in Chinese) |
| Cited By in Cnki (46) | |
| [5] | Yuan Q, Chen Q, Sawaya K. Accurate DOA Estimation Using Array Antenna with Arbitrary Geometry[J]. IEEE Trans on Antennas and Propagation,2005, 53(4):1352-1357 |
| Click to display the text | |
| [6] | Tarik Yardibi, Nikolas S Zawodny, Chris Bahr, et al. Comparison of Microphone Array Processing Techniques for Aeroacoustic Measurements[J]. International Journal of Aeroacoustics, 2010, 9(6):733-762 |
| Click to display the text | |
| [7] |
马巍. 声矢量传感器稳健空间谱估计技术研究[D]. 哈尔滨:哈尔滨工程大学,2013 Ma Wei. Research on Acoustic Vector Sensor Robust SpaceSpectrum Estimation[D]. Harbin, Harbin Engineering University,2013(in Chinese) |
| [8] | Roman N, Wang D L. Binaural Tracking of Multiple Moving Sources[J]. IEEE Trans on Audio, Speech, and Language Processing, 2008, 16(4):728-739 |
| Click to display the text | |
| [9] |
杨春花. 运用m序列测量房间脉冲响应的技术研究[D]. 天津:天津大学,2007 Yang Chunhua. Research on Room Impulse Response Measure Using M-Sequence[D]. Tianjin, Tianjin University, 2007(in Chinese) |
| Cited By in Cnki (4) | |
| [10] | Heinz R. Binaural Room Simulation Based on an Image Source Model with Addition of Statistical Methods to Include the Diffuse Sound Scattering of Walls and to Predict the Reverberant Tail[J]. Applied Acoustics, 1993, 38(2):145-159 |
| Click to display the text |
