一种适用于空心球阵列的声源定位方法
【专利摘要】本发明公开了一种适用于空心球阵列的声源定位方法。该方法结合空心球阵列和最大似然估计的优点，在不影响空间声场的情况下，使用最大似然估计策略，准确定位空间声源方向。本发明的显著优点在于：(1)将传声器接收到的信号变换到球谐域，可降低计算复杂度；(2)使用最大似然估计策略，能够获得较高的空间分辨率和定位准确度；(3)避免了球贝塞尔函数的除法，有效解决了空心球阵列在球贝塞尔函数零点对应频率处鲁棒性弱的问题；(4)可以直接扩展到宽带，降低了球阵列宽带定位算法扩展的复杂度，适用于混响声场；(5)除了适用于空心球阵列以外，该方法同样适用于刚性球阵列以及使用单指向传声器的球阵列。
【专利说明】
一种适用于空心球阵列的声源定位方法 一、
技术领域
[0001] 本发明提出了一种适用于空心球阵列的声源定位方法，结合空心球阵列和最大似 然估计的优点，在不影响空间声场的情况下，使用最大似然估计策略，处理球谐域信号，避 免球贝塞尔函数除法，准确定位空间声源方向，定位方法可有效扩展到宽带信号的处理场 景。 二、
【背景技术】
[0002] 球形传声器阵列具有旋转对称结构，能够设计指向空间任意方向的波束，被广泛 应用在波达方向估计、房间脉冲响应分析、声源定位以及房间结构预测中，为近年研究的热 门话题。空心球阵列能够在不影响测量声场的情况下，获得空间信息，因而更适用于实际应 用场景。
[0003] 文献[I.Balmages&B.Rafaely，Open_sph ere designs for spherical microphone arrays]提出了一种使用双层空心球阵列的结构，在声源定位中，能够避免球 贝塞尔函数零点处空心球阵列鲁棒性较弱的问题，但是需要两倍的传声器数目。文献 [B.Rafaely，The spherical-shell microphone array]提出了一种球体阵列结构，能够避 免球贝塞尔函数零点处空心球阵列鲁棒性较弱的问题，但是传声器单元的半径不固定，实 现车交为复杂。文南犬[B · Rafaely，Bessel nulls recovery in spherical microphone arrays for time-limited signals]在球阵列宽带声源定位时，跳过球贝塞尔函数零点对 应的频率，鲁棒性较高，但是方法灵活度较低。专利CN 102866385A和专利CN 103592628A分 别提出了一种基于球谐域子空间分解和一种基于球谐域实值权重波束形成的声源定位方 法，可以应用于空心球阵列，并且获得较高的空间分辨率，但是上述两种方法并未扩展到宽 带，不适用于混响声场。
[0004] 本专利针对现有球阵列声源定位技术的弱点，提出一种适用于空心球阵列的声源 定位方法。 三、
【发明内容】
[0005] 发明目的：本发明提出来一种适用于空心球阵列的声源定位方法，使用空心球阵 列，在不改变测量声场的情况下，使用最大似然策略，准确定位声源位置。
[0006] 技术方案:为达到上述发明的目的，本发明所提出的技术方案如下:将球形传声器 阵列接收到的时域信号，经过快速傅里叶变换和球谐变换，将信号变换到球谐域，使用最大 似然估计策略，避免球贝塞尔函数的除法，最终准确估计空间声源的方向。
[0007] 本发明中信号在球谐域的表达式卩1?(1〇 = [?(),()(10，？1,-1(10，？1,()(10，？1,1(10，...， PN,N(k) ]T中每一项可以表不为：
[0009]其中，Ωq=(θq，φq)为第q个传声器的位置，θ和φ分别是空间仰角和方位角， p(k， Ω。)为频域信号，Yn,m为η阶m次球谐函数，k为波数，Q为传声器数目，N为球阵列阶数，（· r表 示共辄。
[0010]本发明使用最大似然估计策略，进行单频空间声源定位，计算方法如下：
[0012] 其中，作表示伪逆，cU(k，W)为方向Ψ处的球谐域导向向量，可以表示为：
[0013] dr?(k,W)=B(k)y*(W) (3)
[0014] B(k) =diag{bo(k)，bi(k)，bi(k)，bi(k)，· · ·，bN(k)} (4)
[0015] γ(Ψ) = [Υ0,0(Ψ),Υι,-ι(Ψ),Υι,0(Ψ),Υι,ι(Ψ), . . . ,Υν,ν(Ψ)]Τ (5)
[0016] 对于空心球阵列1^(1〇=4对％(1?1)，九(1?1)为n阶球贝塞尔函数，a为球阵列半径。 公式(2)为单频声源定位方法，回避了球贝塞尔函数匕(1〇的除法，解决的匕(1〇零点对应频 率处鲁棒性弱的问题。
[0017] 在混响声场中，宽带空间声源定位往往会有更好的准确性，本发明的定位方法可 以很容易地扩展到宽带场景，定位方法如下：
[0019] 式中，kmi4Pkmax分别表示最低频率和最高频率对应的波数。
[0020] 本发明提出的声源定位方法，适用于多种阵列形式，如空心球阵列、刚性球阵列和 使用单指向传声器的球阵列。 四、
【附图说明】
[0021] 图1是本发明提出的声源定位方法的流程图
[0022] 图2是本发明在消声室中声源定位实验结果(a)单频(b)宽带
[0023] 图3是本发明在视听室中声源定位实验结果(a)单频(b)宽带 五、
【具体实施方式】
[0024] 附图1是使用本发明提出的声源定位方法的流程图，对传声器接收到的时域信号 进行分帧，对每帧信号做快速傅里叶变换，获得频域信号P(k，Q q)，对频域信号进行球谐变 换，获得球谐域信号p?(k)，球谐变换可以表示为：
[0026] 其中，Ωq=(θq，φq)为第q个传声器的位置，θ和φ分别是空间仰角和方位角， p(k， Ω。)为频域信号，Yn,m为η阶m次球谐函数，k为波数，Q为传声器数目，N为球阵列阶数，（· Γ表 示共辄。
[0027] 使用最大似然策略，进行单频空间声源定位，计算方法如下：
[0029] 其中，印表示伪逆，cU(k，W)为方向Ψ处的球谐域导向向量，可以表示为：
[0030] cU(k，W)=B(k)y*(W) (9)
[0031] B(k) =diag{bo(k)，bi(k)，bi(k)，bi(k)，· · ·，bN(k)} (10)
[0032] γ(Ψ) = [Υ0,0(Ψ),Υι,-ι(Ψ),Υι,0(Ψ),Υι,ι(Ψ), . . . ,Υν,ν(Ψ)]Τ (11)
[0033] 对于空心球阵列1^(1〇=4对％(1?1)，九(1?1)为n阶球贝塞尔函数，a为球阵列半径。
[0034] 对于混响声场，使用单频声源定位误差较大，本发明提出的声源定位方法可以直 接扩展到宽带，使用最大似然估计策略，进行宽带空间声源定位，计算方法如下：
[0036]式中，kmi4Pkmax分别表示最低频率和最高频率对应的波数。
[0037] 使用本发明提出的声源定位方法，分别在消声室和视听室中进行了实验，实验中 使用一个8单元均勾分布空心球阵列，球阵列半径为13.9cm，球谐分解展开阶数N=1，每帧 时域信号帧长K = 1024,采样频率为fs = 16kHz，单频信号频率满足ka = l，宽带信号为白噪 声，宽带声源定位时使用的频率范围kae [0. 5，1 ]。图2给出了消声室中声源定位实验结果， 图中黑色圆圈表示实际声源位置，可以看出消声室中，使用单频和宽带定位方法，均可以获 得较好的声源定位结果，对于一个l〇s的宽带信号，定位的均方根误差为〇.58°。图3给出了 试听室中声源定位实验结果，图中可以看出试听室中，由于混响的存在，使用单频定位方法 无法准确定位声源的方向，而使用本发明提出的宽带声源定位方法，可以准确定位声源的 方向。对于一个10S的宽带信号，定位的均方根误差为4.43°。
[0038] 本实施例只给出了空心球阵列的相关结果，而本发明提出的声源定位方法同样适 用于刚性球阵列和使用单指向传声器的球阵列。
【主权项】
1. 一种适用于空屯、球阵列的声源定位方法，其处理方法包含如下步骤：（1)根据空屯、球 阵列中传声器的位置，W及传声器接收到的时域信号，通过快速傅里叶变换和球谐变换，获 得信号在球谐域的表达式Pn"化）；第二步、根据信号在球谐域的表达式Pn"(k)，W及空间任 意一个方向Ψ处的球谐域导向向量CU化，Ψ )，使用最大似然估计策略，获得波数k对应频率 (k=2村/C，C表示声速，f表示频率）下的定位结果；第Ξ步、根据设定的波数上下限kmin和 kmax，将算法扩展到宽带，获得宽带下的定位结果。2. 如权利要求1所述的一种适用于空屯、球阵列的最大似然声源定位方法，其特征在于： 信号在球谐域的表达式Pnm化）= [P〇,日化），P1,-1化），P1,日化），P1,1化），...，PN,^k)]T中的每一 项可W表示为：(1) 其中，Ω q= (Θ。，Φ q)为第q个传声器的位置，Θ和Φ分别是空间仰角和方位角，P化，Ω q) 为频域信号，Yn,m为η阶m次球谐函数，k为波数，Q为传声器数目，N为球阵列阶数，（· Γ表示 共辆。3. 如权利要求1所述的一种适用于空屯、球阵列的声源定位方法，其特征在于:使用最大 似然估计策略，避免球贝塞尔函数的除法。4. 如权利要求1所述的一种适用于空屯、球阵列的声源定位方法，其特征在于:单频信号 声源定位计算方法如下：(2) 其中，表示伪逆，dnm化，Ψ)为方向Ψ处的球谐域导向向量，可W表示为： dnm化，Ψ)=Β(1〇/(Ψ) (3) B(k)=diag{bo(k)，bi(k)，bi(k)，bi(k)，. . .，bN(k)} (4) 7(Ψ) = [Υο'ο(Ψ)，Υι'-ι(Ψ)，Υι'ο(Ψ)，Υι'ι(Ψ)，...，化'Ν(Ψ)]τ (5) 对于空屯、球阵列bn化）=仙門。化a), jn化a)为η阶球贝塞尔函数，a为球阵列半径。5. 如权利要求1所述的一种适用于空屯、球阵列的声源定位方法，其特征在于：宽带信号 声源定位计算方法如下：(6) 式中，kmin和kmax分别表示最低频率和最高频率对应的波数。 6 .如权利要求1所述的一种适用于空屯、球阵列的声源定位方法，其特征在于:该方法适 用于多种阵列形式，如空屯、球阵列、刚性球阵列和使用单指向传声器的球阵列。
【文档编号】G01S5/20GK106093866SQ201610367282
【公开日】2016年11月9日
【申请日】2016年5月27日
【发明人】胡玉祥, 卢晶
【申请人】南京大学