CN102033223A - 使用麦克风阵列对声源定位的方法 - Google Patents
使用麦克风阵列对声源定位的方法 Download PDFInfo
- Publication number
- CN102033223A CN102033223A CN2010106110878A CN201010611087A CN102033223A CN 102033223 A CN102033223 A CN 102033223A CN 2010106110878 A CN2010106110878 A CN 2010106110878A CN 201010611087 A CN201010611087 A CN 201010611087A CN 102033223 A CN102033223 A CN 102033223A
- Authority
- CN
- China
- Prior art keywords
- microphone
- sound source
- distance
- position angle
- mistiming
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
提供一种使用麦克风阵列对声源定位的方法,三个麦克风构成正三角形,该方法包括:建立坐标系,其原点与正三角形的重心重合,第一麦克风位于坐标系的纵轴上;将正三角形的重心与正三角形的三个顶点连接并延长,将全圆周分为6个相等的区间;计算声源分别到达第一至第三麦克风的时间差;将时间差乘以音速,得到声源到第二麦克风的距离与声源到第一麦克风的距离之间的第一距离差以及声源到第三麦克风的距离与声源到第一麦克风的距离之间的第二距离差;根据第一和第二距离差的符号以及第一、第二距离差的绝对值大小比较关系,初步确定声源相对于正三角形重心的方位角Φ所落入的区间;采用0.618法对声源与正三角形重心的距离R以及方位角Φ逼近优化点。
Description
技术领域
本发明涉及声源的定位,更具体地讲,涉及一种使用麦克风(MIC)阵列来对声源定位的方法。
背景技术
阵列信号处理已经广泛应用于通信、雷达、声纳、医学和航天航空等诸多领域。近年来,随着反恐战的开展,各国已经研制了各种探测***点或枪炮射击点的位置的设备。在实际应用中,由于视觉定位往往不易实现,从而发展了一些声学定位的***,例如,美国的回力棒***、以色列拉斐尔公司的反狙击手声探侧***,加拿大麦克唐纳·迪特维利公司和加拿大防务研究和发展委员会也共同开发了一种名为“雪貂”(Ferret)的小型武器探测和定位***。Microsoft公司开发的新的Vista操作***集成了对MIC阵列的支持,Intel公司的HDAudoi规范可以实现16个MIC、32KHz采样,这些为声波信号阵列处理提供了良好的支持。
目前,使用麦克风阵列来对声源进行定位的方法总体来说可以分为两类:定向波速形成和时间延迟估计。定向波速形成是基于最大输出功率可控波束方法的一种形式,该方法对麦克风阵列接收到的语音信号进行滤波延迟并加权求和形成特定指向的波束,然后计算各个指向上的波束能量,其中搜索空间中功率最大的方向就是对声源位置的一个最大近似。在基于时间延迟估计的方法中,由于声源发出的信号到达各个麦克风存在时间差,该方法根据各个通道信号间的相关关系对这个时间延迟进行估计,从而来对声源位置进行估计。
定向波速形成的方法为了搜索能够使波束功率最大的位置,需要对整个空间进行扫描,需要进行大量计算,因此很难满足实时的要求。
基于时间延迟估计的方法计算速度快,能达到实时的要求。这类方法可以适用于如同语音这样的非平稳宽带信源,但是在只有单个信源时能较好地工作,定位的精度与到达时间差(TDOA)的估计精度有关。
发明内容
本发明的目的在于提供一种使用麦克风阵列对声源定位的方法,该方法能够快速对声源定位。
根据本发明的一方面,提供一种使用麦克风阵列对声源定位的方法,所述麦克风阵列为构成正三角形的三个麦克风,所述方法包括:建立坐标系,所述坐标系的原点与所述正三角形的重心重合,三个麦克风中的第一麦克风位于坐标系的纵轴上;将正三角形的重心与正三角形的三个顶点连接并延长,从而将以所述正三角形的重心为圆心的全圆周分为6个相等的区间;计算声源分别到达三个麦克风中的第一麦克风、第二麦克风、第三麦克风的时间差;将计算的时间差乘以音速,得到声源到第二麦克风的距离与声源到第一麦克风的距离之间的第一距离差以及声源到第三麦克风的距离与声源到第一麦克风的距离之间的第二距离差;根据第一距离差和第二距离差的符号以及第一距离差、第二距离差的绝对值大小比较关系,初步确定声源相对于正三角形重心的方位角Φ所落入的区间;采用0.618法对声源与正三角形重心的距离R以及方位角Φ逼近优化点,获得最终估计的距离R和方位角Φ。
根据本发明的另一方面,提供一种使用麦克风阵列对声源定位的方法,所述麦克风阵列为构成正三角形的三个麦克风,所述方法包括:建立坐标系,所述坐标系的原点与所述正三角形的重心重合,三个麦克风中的第一麦克风位于坐标系的纵轴上;计算声源分别到达三个麦克风中的第一麦克风、第二麦克风、第三麦克风的时间差;将计算的时间差乘以音速,得到声源分别到达第一麦克风、第二麦克风、第三麦克风的距离差;根据三角形边角关系建立关于声源到各个麦克风的距离差的联立方程组;对联立方程组中的任意两个方程进行联立求解,得到声源与正三角形重心的距离R以及声源相对于正三角形重心的方位角Φ。
附图说明
通过结合附图,从下面的实施例的描述中,本发明这些和/或其它方面及优点将会变得清楚,并且更易于理解,其中:
图1是根据本发明的采用正三角形麦克风阵列对声源定位的示意图;
图2是根据本发明第一实施例的采用正三角形麦克风阵列对声源定位的方法的流程图;
图3是根据本发明的采用时间延迟特性进行方位分割的示意图;
图4是根据本发明第二实施例的采用正三角形麦克风阵列对声源定位的方法的流程图;
具体实施方式
本发明针对基于时间延迟定位声源的方法进行了改进和创新。在本发明中,采用了准L1相关技术提高时间延迟估计的速度和抗野值干扰能力,进一步采用时间延迟特性进行方位分割以提高计算速度和去除局部极值点,同时可采用0.618法加快逼近速度,从而可以用普通微型计算机在数十毫秒级解决声源定位问题,满足实用性的速度和精度要求。
以下,参照附图来详细说明本发明的原理以及实施例。在本发明中,采用由三个麦克风构成的麦克风阵列来对声源定位,三个麦克风被布置为正三角形,即,三个麦克风分别位于正三角形的三个顶点上。
图1是根据本发明的采用正三角形麦克风阵列对声源定位的示意图,图2是根据本发明第一实施例的采用正三角形麦克风阵列对声源定位的方法的流程图。
参照图1和图2,在步骤201,建立坐标系,坐标系的原点O(0,0)与第一麦克风a、第二麦克风b、第三麦克风c构成的正三角形的重心重合,三个麦克风中的一个麦克风(例如第一麦克风a)位于坐标系的纵轴上。
在步骤202,计算声源p分别到达第一麦克风a、第二麦克风b、第三麦克风c的时间差。
在步骤203,将计算的时间差乘以声音在空气中传播的速度(音速),得到声源p分别到达第一麦克风a、第二麦克风b、第三麦克风c的距离差。
在计算所述时间差的步骤中,可采用准L1算法估计所述时间差。普通相关函数可以认为是2范数下的测度,L1相关是1范数下的测度,L1相关与L2相关有几乎相同的各种性质。但是,L1相关比L2相关容易计算,且抗干扰能力更强。
两个麦克风处的两个声音信号x(n)和y(n)的L1互相关函数定义如下:
其中,n表示声音信号的采样时间点,N表示声音信号的总采样点数,m表示时间偏移量。m是采样间隔的整数倍。
由于上式的第一项对平稳随机过程是常数,因此选择最大的互相关值Rxy(m),等效于选择最小的第二项,这就是准L1互相关函数。
对于本发明,对声音信号进行采样,已知两个麦克风处的两个声音信号x(n)和y(n)的波形基本相同但时间延迟不同,针对上述互相关函数改变m值,使得两个声音信号x(n)和y(n)的互相关值Rxy(m)最大,即,调整x(n)与y(n)错位对齐,此时的m值即是两个声音信号x(n)和y(n)的时间差。另外,本领域技术人员应该理解,在上述互相关函数中,abs()表示绝对值函数,x(n)和y(n)表示声音信号的幅度。
互相关比对求时间差的方法(即,准L1算法)虽然可以应用所有采样点的延迟信息,但因为数字化,所以比对只是在整采样点处进行,而最佳延迟量一般不在整采样点处。对此问题可以用提高采样率的方法提高估计精度。采样率提高,导致采样间隔较密,可以更准确地估计延迟时间差。但采样率受硬件限制,存储设备也受硬件限制。因此,在没有提高采样率的情况下,在本发明中还可以采用三点二次曲线插值的方法估计两个采样点之间的值,再进行互相关法估计,从而提高了估计精度。
如图1所示,正三角形的重心为坐标系的原点O(0,0)。正三角形的边长为L,声源p到各个麦克风a、b、c的距离差为:
D1=R2-R1(1)
D2=R3-R1(2)
D3=R3-R2(3)
声源p的极坐标为(R,Φ),R是声源p与正三角形的重心的距离,Φ是声源p相对于正三角形的重心的方位角,R和Φ是待求未知数。
三个麦克风a、b、c的顶点极坐标分别为麦克风a:(r,π/2),麦克风b:(r,7π/6),麦克风c:(r,-π/6)。
另外,参照图2,各个麦克风a、b、c与正三角形的重心的距离r有下述关系:
在步骤204,根据三角形边角关系建立关于声源p到各个麦克风a、b、c的距离差的联立方程组。
参照图2,根据三角形边角关系可以得出:
将方程(4)-(7)代入方程(1)-(3),可以得到声源p到各个麦克风a、b、c的距离差的联立方程组如下:
在步骤205,对联立方程组中的任意两个方程进行联立求解,得到声源p的极坐标,从而完成对声源S的定位。
例如,对三个方程(8)-(10)中的任意两个方程进行联立求解,均可以解出方位角Φ和距离R的值,得到声源p的极坐标(R,Φ)。
由于实际测量误差,还可以进一步采用方程(8)(9)、(9)(10)、(8)(10)分别联立求解,分别求出方位角Φ和距离R的三组解,再取三组解的平均值作为方位角Φ和距离R的估计值。
图3是根据本发明的采用时间延迟特性进行方位分割的示意图,图4是根据本发明第二实施例的采用正三角形麦克风阵列对声源定位的方法的流程图。
在步骤401,建立坐标系,坐标系的原点O(0,0)与第一麦克风a、第二麦克风b、第三麦克风c构成的正三角形的重心重合,三个麦克风中的一个麦克风(例如第一麦克风a)位于坐标系的纵轴上。
为了防止逼近算法进入局部最小值,考虑时间差的符号和大小,在步骤402,将正三角形的重心与正三角形的三个顶点连接并延长,从而将以正三角形的重心为圆心的全圆周分为6个相等的区间1-6,如图3所示。
在步骤403,计算圆周上任一声源p分别到达第一麦克风a、第二麦克风b、第三麦克风c的时间差。
在第一实施例中描述的准L1算法以及三点二次曲线插值的方法同样可应用于第二实施例中,这里不再重复描述。
在步骤404,将计算的时间差乘以声音在空气中传播的速度(音速),得到声源p到第二麦克风b的距离与声源p到第一麦克风a的距离之间的第一距离差D1以及声源p到第三麦克风c的距离与声源p到第一麦克风a的距离之间的第二距离差D2。
在步骤405,根据第一距离差D1和第二距离差D2的符号以及D1、D2的绝对值大小比较关系,初步确定方位角Φ所落入的区间(即,方位角Φ落入区间1-6中的哪个区间中)。
参照图3和图4,假设D1为圆周上任一声源p到麦克风b的距离pb减去声源p到麦克风a的距离pa所得的距离差,即D1=pb-pa;假设D2为圆周上声源p到麦克风c的距离pc减去声源p到麦克风a的距离pa所得的距离差,即D2=pc-pa,则有如下表所示的关系:
区间 | Φ | D1 | D2 | 性质 |
1 | 0°-60° | + | + | |D1|<|D2| |
2 | 60°-120° | - | + | |D1|<|D2| |
3 | 120°-180° | - | - | |D1|<|D2| |
4 | 180°-240° | - | - | |D1|>|D2| |
5 | 240°-300° | + | - | |D1|>|D2| |
6 | 300°-360° | + | + | |D1|>|D2| |
因此,根据D1、D2的符号以及D1、D2的绝对值大小比较关系,可以初步确定方位角所落入的区间,从而使每一个取值范围内具有单一的优化点,避免了局部最小值。同时减少估计的值域范围,提高估计速度。
在步骤406,用0.618法对双变量(即,距离R和方位角Φ)逼近优化点,获得最终的估计结果。
0.618法也称为黄金分割法,0.618为黄金分割点。可对二维单调唯一最小点的进行0.618法选优。
在本实施例中,先取距离R的初始值。距离R的初始值在50-400米的范围内选取,优选为200米。如果距离R小于50米,则声源借助视力可被容易确定。如果距离R大于400米,则接收的数据差比已经很小,包含的距离信息较少,不便于识别定位。因为将针对距离R的初始值进行0.618法逼近,因此,距离R的初始值可认为是对声源p定位的测量范围。
然后,固定距离R的初始值,在初步确定的方位角Φ的范围内,采用0.618法缩小优化点所在的区间,并逼近误差最小的方位角Φ,直到两次逼近的方位角Φ的误差变化小于预定阈值(例如相对误差为0.01)为止。接着,固定该逼近的方位角Φ,采用0.618法缩小距离R的范围,直到两次逼近的距离R的误差变化小于另一预定阈值为止,得到逼近的距离R。因此,对方位角Φ及距离R进行逼近,直到二者的误差变化达到指定的阈值以内,得到估计的结果。
具体地讲,对于方位角Φ,如上所述,根据D1、D2的符号以及D1、D2的绝对值大小比较关系,可以初步确定方位角Φ落入区间1-6中的哪个区间中。
接着,可按照下面的方式逼近误差最小的方位角Φ:
(a1)取方位角Φ为确定的区间上限S1乘以0.618;
(a2)根据三角形的边角关系,利用所述方位角Φ以及距离R的初始值求出声源p到各个麦克风的距离R1、R2、R3;
(a3)将计算出的距离R1、R2、R3分别除以音速,可以求得声源p到达各个麦克风的时间;
(a4)计算声源p到达第二麦克风b与声源p到达第一麦克风a的时间差t′ba1,并计算声源p到达第三麦克风c与声源p到达第一麦克风a的时间差t′ca1;
(a5)利用先前通过准L1互相关算法(互相关值)得到的声源p到达第二麦克风b与声源p到达第一麦克风a的时间差tba以及声源p到达第三麦克风c与声源p到达第一麦克风a的时间差tca,计算误差error1=|tba-t′ba1|+|tca-t′ca1|;
(a6)取方位角Φ为上限S1乘以(1-0.618),按照步骤(a2)-(a5)计算误差error2=|tba-t′ba2|+|tca-t′ca2|;这里,t′ba2和t′ca2分别是方位角Φ为上限S1乘以(1-0.618)时计算的声源p到达第二麦克风b与声源p到达第一麦克风a的时间差以及计算的声源p到达第三麦克风c与声源p到达第一麦克风a的时间差;
(a7)如果|error1-error2|大于预定阈值(例如可以是0.01),并且error1>error2,则确定方位角Φ的范围在确定的区间的下限S2与上限S1乘以0.618之间;如果|error1-error2|大于所述预定阈值,并且error1<error2,则确定方位角Φ的范围在上限S1乘以(1-0.618)与上限S1之间;
(a8)根据再次确定的方位角Φ所落入的区间,重复上述步骤(a1)-(a7),直到|error1-error2|小于所述预定阈值。由于最后所取的两个方位角Φ之间相差很小,因此最后所取的两个方位角Φ之一或二者的平均值为最终估计的方位角Φ。
例如,假设初步确定方位角落入区间1中(0到60°),则取方位角Φ为60°×0.618,为37.08°。因为在准L1相关算法中已经得到声源p到达各个麦克风的时间差tba、tca,因此可根据三角形的边角关系,将方位角Φ(37.08°)以及距离R的初始值代入方程(5)-(7),可求出声源p到各个麦克风的距离R1、R2、R3。将计算出的距离R1、R2、R3分别除以音速,可以求得声源p到达各个麦克风的时间,进而可以求得声源p到达第二麦克风b与声源p到达第一麦克风a的时间差t′ba1,并求得声源p到达第三麦克风c与声源p到达第一麦克风a的时间差t′ca1。然后,计算误差error1=|tba-t′ba1|+|tca-t′ca1|。
接着,取方位角Φ为60°×(1-0.618),为22.92°。按照以上同样的方式求出误差error2=|tba-t′ba2|+|tca-t′ca2|。
当|error1-error2|大于预定阈值(例如0.01)时,继续按照0.618法搜索。当error1>error2时,所求方位角Φ的范围是在0~37.08°之间,取方位角Φ为37.08°×0.618,计算error1,之后取方位角Φ为37.08°×(1-0.618),计算error2;反之,当error1<error2时,所求方位角Φ的范围是在22.92°~60°之间,取方位角Φ为22.92°+(60°-22.92°)×0.618,计算error1,之后取方位角Φ为22.92°+(60°-22.92°)×(1-0.618),计算error2。按照上述相同的方式继续求解。直到|error1-error2|小于所述预定阈值为止,最后所取的两个方位角Φ之一或二者的平均值为最终估计的方位角Φ。
在最终估计了方位角Φ之后,可对距离R进行逼近估计。
具体地讲,距离R的初始值以及零值确定了一区间(0,初始值),可按照下面的方式逼近误差最小的距离R:
(b1)取距离R为确定的区间的上限S3乘以0.618;
(b2)根据三角形的边角关系,利用最终估计的方位角Φ以及距离R求出声源p到各个麦克风的距离R1、R2、R3;
(b3)将计算出的距离R1、R2、R3分别除以音速,可以求得声源p到达各个麦克风的时间;
(b4)计算声源p到达第二麦克风b与声源p到达第一麦克风a的时间差t′ba3,并计算声源p到达第三麦克风c与声源p到达第一麦克风a的时间差t′ca3;
(b5)利用先前通过准L1互相关算法(互相关值)得到的声源p到达第二麦克风b与声源p到达第一麦克风a的时间差tba以及声源p到达第三麦克风c与声源p到达第一麦克风a的时间差tca,计算误差error3=|tba-t′ba3|+|tca-t′ca3|;
(b6)取距离R为上限S3乘以(1-0.618),按照步骤(b2)-(b5)计算误差error4=|tba-t′ba4|+|tca-t′ca4|;这里,t′ba4和t′ca4分别是距离R为上限S3乘以(1-0.618)时计算的声源p到达第二麦克风b与声源p到达第一麦克风a的时间差以及计算的声源p到达第三麦克风c与声源p到达第一麦克风a的时间差;
(b7)如果|error3-error4|大于另一预定阈值(例如可以是0.01),并且error3>error4,则确定距离R的范围在确定的区间的下限S4与上限S3乘以0.618之间;如果|error3-error4|大于所述另一预定阈值,并且error3<error4,则确定距离R的范围在上限S3乘以(1-0.618)与上限S3之间;
(b8)根据再次确定的距离R所落入的区间,重复上述步骤(b1)-(b7),直到|error3-error4|小于所述另一预定阈值。由于最后所取的两个距离R之间相差很小,因此最后所取的两个距离R之一或二者的平均值为最终估计的距离R。
例如,如果距离R的初始值为200米,则取距离R为200×0.618,为123.6米。因为在准L1相关算法中已经得到声源p到达各个麦克风的时间差tba、tca,因此可根据三角形的边角关系,将距离R(123.6米)以及前面已经求得的最终估计的方位角Φ代入方程(5)-(7),可求出声源p到各个麦克风的距离R1、R2、R3。将计算出的距离R1、R2、R3分别除以音速,可以求得声源p到达各个麦克风的时间,进而可以求得声源p到达第二麦克风b与声源p到达第一麦克风a的时间差t′ba3,并求得声源p到达第三麦克风c与声源p到达第一麦克风a的时间差t′ca3。然后,计算误差error3=|tba-t′ba3|+|tca-t′ca3|。
接着,取距离R为200×(1-0.618),为76.4米。按照以上同样的方式求出误差error4=|tba-t′ba4|+|tca-t′ca4|。
当|error3-error4|大于预定阈值(例如0.01)时,继续按照0.618法搜索。当error3>error4时,所求距离R的范围是在0~123.6米之间,取距离R为123.6×0.618,计算error3,之后取距离R为123.6×(1-0.618),计算error4;反之,当error3<error4时,所求距离R的范围是在76.4~200米之间,取距离R为76.4+(200-76.4)×0.618,计算error3,之后取距离R为76.4+(200-76.4)×(1-0.618),计算error4。按照上述相同的方式继续求解。直到|error3-error4|小于所述预定阈值为止,最后所取的两个距离R之一或二者的平均值为最终估计的距离R。
此外,由于在实际应用中,距离R的估计不如方位角Φ的估计重要,因此也可以在初步估计方位角Φ之后,采用固定的一些R值(例如,50米、100米、150米、200米、300米)进行估计,找出误差最小的值即可。
0.618法比二分之一法能够更快地得到逼近结果。
根据本发明,分析了延时特性,提出了方位分割的方法,采用0.618法加快逼近速度,从而减少计算的复杂度,实现了快速估计的算法。对延时的估计,提出了准L1相关的算法,减少了乘法运算,提高了估计的精度。另外,采用了插值算法进一步提高延时的估计精度,从而得到了快速准确的算法。
另外,在麦克风布局上,采用正三角形布局对各方向的敏感度一致性好。
虽然针对由3个麦克风构成的正三角形阵列来描述本发明的麦克风定位方法,但是也可采用由4个麦克风构成的正方形阵列来实现本发明的麦克风定位方法,只需要调整相应的几何关系即可。同样,采用正方形的麦克风布局对各方向的敏感度一致性好。
虽然本发明是参照其示例性的实施例被具体描述和显示的,但是本领域的普通技术人员应该理解,在不脱离由权利要求限定的本发明的精神和范围的情况下,可以对其进行形式和细节的各种改变。
Claims (8)
1.一种使用麦克风阵列对声源定位的方法,所述麦克风阵列为构成正三角形的三个麦克风,所述方法包括:
建立坐标系,所述坐标系的原点与所述正三角形的重心重合,三个麦克风中的第一麦克风位于坐标系的纵轴上;
将正三角形的重心与正三角形的三个顶点连接并延长,从而将以所述正三角形的重心为圆心的全圆周分为6个相等的区间;
计算声源分别到达三个麦克风中的第一麦克风、第二麦克风、第三麦克风的时间差;
将计算的时间差乘以音速,得到声源到第二麦克风的距离与声源到第一麦克风的距离之间的第一距离差以及声源到第三麦克风的距离与声源到第一麦克风的距离之间的第二距离差;
根据第一距离差和第二距离差的符号以及第一距离差、第二距离差的绝对值大小比较关系,初步确定声源相对于正三角形重心的方位角Φ所落入的区间;
采用0.618法对声源与正三角形重心的距离R以及方位角Φ逼近优化点,获得最终估计的距离R和方位角Φ。
3.根据权利要求2所述的方法,其中,采用0.618法逼近优化点的步骤包括:
对距离R取初始值,固定距离R的所述初始值,在初步确定的方位角Φ的范围内,采用0.618法缩小优化点所在的区间,并逼近误差最小的方位角Φ,直到两次逼近的方位角Φ的误差变化小于预定阈值;
固定逼近的方位角Φ,采用0.618法缩小距离R的范围,直到两次逼近的距离R的误差变化小于另一预定阈值,得到逼近的距离R,由此获得最终估计的距离R和方位角Φ。
4.根据权利要求3所述的方法,其中,逼近的方位角Φ的步骤包括:
(a1)取方位角Φ为确定的区间的上限S1乘以0.618;
(a2)根据三角形的边角关系,利用方位角Φ以及距离R的初始值求出声源到第一麦克风、第二麦克风、第三麦克风的距离;
(a3)将计算出的距离分别除以音速,计算声源到达第一麦克风、第二麦克风、第三麦克风的时间;
(a4)计算声源到达第二麦克风与声源到达第一麦克风的时间差t′ba1,并计算声源到达第三麦克风与声源到达第一麦克风的时间差t′ca1;
(a5)利用先前通过声音信号的互相关值得到的声源到达第二麦克风与声源到达第一麦克风的时间差tba以及声源到达第三麦克风与声源到达第一麦克风的时间差tca,计算误差error1=|tba-t′ba1|+|tca-t′ca1|;
(a6)取方位角Φ为上限S1乘以(1-0.618),按照步骤(a2)至步骤(a5)计算误差error2=|tba-t′ba2|+|tca-t′ca2|,其中,t′ba2和t′ca2分别是方位角Φ为上限S1乘以(1-0.618)时计算的声源到达第二麦克风与声源到达第一麦克风的时间差以及计算的声源到达第三麦克风与声源到达第一麦克风的时间差;
(a7)如果|error1-error2|大于预定阈值,并且error1>error2,则确定方位角Φ的范围在确定的区间的下限S2与上限S1乘以0.618之间;如果|error1-error2|大于所述预定阈值,并且error1<error2,则确定方位角Φ的范围在上限S1乘以(1-0.618)与上限S1之间;
(a8)根据再次确定的方位角Φ所落入的区间,重复上述步骤(a1)至步骤(a7),直到|error1-error2|小于所述预定阈值,最后所取的两个方位角Φ之一或二者的平均值为最终估计的方位角Φ。
5.根据权利要求4所述的方法,其中,距离R的初始值以及零值确定一区间,逼近的距离R的步骤包括:
(b1)取距离R为确定的区间上限S3乘以0.618;
(b2)根据三角形的边角关系,利用最终估计的方位角Φ以及距离R求出声源到第一麦克风、第二麦克风、第三麦克风的距离R1、R2、R3;
(b3)将计算出的距离分别除以音速,可以求得声源到达第一麦克风、第二麦克风、第三麦克风的时间;
(b4)计算声源到达第二麦克风与声源到达第一麦克风的时间差t′ba3,并计算声源到达第三麦克风与声源到达第一麦克风的时间差t′ca3;
(b5)利用先前通过声音信号的互相关值得到的声源到达第二麦克风与声源到达第一麦克风的时间差tba以及声源到达第三麦克风与声源到达第一麦克风的时间差tca,计算误差error3=|tba-t′ba3|+|tca-t′ca3|;
(b6)取距离R为上限S3乘以(1-0.618),按照步骤(b2)至步骤(b5)计算误差error4=|tba-t′ba4|+|tca-t′ca4|,其中,t′ba4和t′ca4分别是距离R为上限S3乘以(1-0.618)时计算的声源到达第二麦克风与声源到达第一麦克风的时间差以及计算的声源到达第三麦克风与声源到达第一麦克风的时间差;
(b7)如果|error3-error4|大于另一预定阈值,并且error3>error4,则确定距离R的范围在确定的区间的下限S4与上限S3乘以0.618之间;如果|error3-error4|大于所述另一预定阈值,并且error3<error4,则确定距离R的范围在上限S3乘以(1-0.618)与上限S3之间;
(b8)根据再次确定的距离R所落入的区间,重复上述步骤(b1)至步骤(b7),直到|error3-error4|小于所述另一预定阈值,最后所取的两个距离R之一或二者的平均值为最终估计的距离R。
6.一种使用麦克风阵列对声源定位的方法,所述麦克风阵列为构成正三角形的三个麦克风,所述方法包括:
建立坐标系,所述坐标系的原点与所述正三角形的重心重合,三个麦克风中的第一麦克风位于坐标系的纵轴上;
计算声源分别到达三个麦克风中的第一麦克风、第二麦克风、第三麦克风的时间差;
将计算的时间差乘以音速,得到声源分别到达第一麦克风、第二麦克风、第三麦克风的距离差;
根据三角形边角关系建立关于声源到各个麦克风的距离差的联立方程组;
对联立方程组中的任意两个方程进行联立求解,得到声源与正三角形重心的距离R以及声源相对于正三角形重心的方位角Φ。
7.根据权利要求6所述的方法,还包括:对联立方程组中的三组两个不同的方程进行联立求解,分别求出方位角Φ和距离R的三组解,再取三组解的平均值作为方位角Φ和距离R的估计值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201010611087A CN102033223B (zh) | 2010-12-29 | 2010-12-29 | 使用麦克风阵列对声源定位的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201010611087A CN102033223B (zh) | 2010-12-29 | 2010-12-29 | 使用麦克风阵列对声源定位的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102033223A true CN102033223A (zh) | 2011-04-27 |
CN102033223B CN102033223B (zh) | 2012-10-03 |
Family
ID=43886372
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201010611087A Expired - Fee Related CN102033223B (zh) | 2010-12-29 | 2010-12-29 | 使用麦克风阵列对声源定位的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102033223B (zh) |
Cited By (32)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103054669A (zh) * | 2012-12-25 | 2013-04-24 | 汪勇 | 一种耳塞式声音辅助感知装置 |
CN103310786A (zh) * | 2013-06-17 | 2013-09-18 | 西北工业大学 | 智能食堂打卡机的三维语音识别方法 |
CN103439689A (zh) * | 2013-08-21 | 2013-12-11 | 大连理工大学 | 一种分布式麦克风阵列中麦克风位置估计*** |
CN103473864A (zh) * | 2013-06-17 | 2013-12-25 | 西北工业大学 | 智能食堂打卡机的语音识别及指纹结算方法 |
CN103714551A (zh) * | 2014-01-20 | 2014-04-09 | 哈尔滨工业大学 | 基于区间消去的汽车仪表自动检测视觉***的相机外参数标定方法 |
CN103760479A (zh) * | 2014-01-10 | 2014-04-30 | 三峡大学 | 便携式电晕定位仪 |
CN103995252A (zh) * | 2014-05-13 | 2014-08-20 | 南京信息工程大学 | 一种三维空间声源定位方法 |
CN104035065A (zh) * | 2014-06-23 | 2014-09-10 | 河北工业大学 | 基于主动旋转的声源定向装置及其应用方法 |
CN104200805A (zh) * | 2014-08-30 | 2014-12-10 | 长城汽车股份有限公司 | 汽车驾驶员语音助手 |
CN104422922A (zh) * | 2013-08-19 | 2015-03-18 | 中兴通讯股份有限公司 | 一种移动终端实现声源定位的方法及装置 |
CN104459625A (zh) * | 2014-12-14 | 2015-03-25 | 南京理工大学 | 基于轨道移动双麦克风阵列的声源定位装置及方法 |
CN104502893A (zh) * | 2014-12-10 | 2015-04-08 | 北京智谷睿拓技术服务有限公司 | 定位方法、定位装置以及用户设备 |
CN104572009A (zh) * | 2015-01-28 | 2015-04-29 | 合肥联宝信息技术有限公司 | 一种自适应外界环境的音频控制方法及装置 |
CN104848478A (zh) * | 2015-04-28 | 2015-08-19 | 广东美的制冷设备有限公司 | 空调器中语音指令识别阈值的调整方法和*** |
CN104977560A (zh) * | 2014-04-04 | 2015-10-14 | 株式会社理光 | 移动设备、定位方法和定位*** |
CN104991573A (zh) * | 2015-06-25 | 2015-10-21 | 北京品创汇通科技有限公司 | 一种基于声源阵列的定位跟踪方法及其装置 |
CN105611014A (zh) * | 2014-11-21 | 2016-05-25 | 中兴通讯股份有限公司 | 移动终端通话语音降噪方法及装置 |
CN105629202A (zh) * | 2016-03-23 | 2016-06-01 | 浙江工业大学 | 伞形可重组三维传声器阵列识别声源三维坐标的方法 |
CN106535128A (zh) * | 2016-11-22 | 2017-03-22 | 四川九洲电器集团有限责任公司 | 一种多点定位方法及装置 |
CN106911985A (zh) * | 2017-03-20 | 2017-06-30 | 广州视源电子科技股份有限公司 | 时间差确定的方法和装置 |
CN107154256A (zh) * | 2017-06-27 | 2017-09-12 | 山东省计算中心(国家超级计算济南中心) | 基于声源定位的声掩蔽***及自适应调整方法 |
CN107592600A (zh) * | 2016-07-06 | 2018-01-16 | 深圳市三诺声智联股份有限公司 | 一种基于分布式麦克风的拾音筛选方法及拾音装置 |
CN107767137A (zh) * | 2016-08-23 | 2018-03-06 | ***通信有限公司研究院 | 一种信息处理方法、装置及终端 |
CN107800967A (zh) * | 2017-10-30 | 2018-03-13 | 维沃移动通信有限公司 | 一种拍摄方法及移动终端 |
CN108416996A (zh) * | 2018-04-28 | 2018-08-17 | 上海市城市建设设计研究总院(集团)有限公司 | 基于城市道路附属设施的声源定位***及其使用方法 |
CN109019349A (zh) * | 2018-06-19 | 2018-12-18 | 上海振华重工(集团)股份有限公司 | 故障检测方法、检测器、计算机存储介质及起重机 |
CN109068235A (zh) * | 2017-06-12 | 2018-12-21 | 田中良 | 用于精确计算声音在麦克风阵列处的到达方向的方法 |
CN110390947A (zh) * | 2018-04-23 | 2019-10-29 | 北京京东尚科信息技术有限公司 | 声源位置的确定方法、***、设备和存储介质 |
CN110620836A (zh) * | 2019-09-06 | 2019-12-27 | 中国民航大学 | 一种手机听声定位技术 |
CN112964256A (zh) * | 2019-12-13 | 2021-06-15 | 佛山市云米电器科技有限公司 | 室内定位方法、智能家电设备及计算机可读存储介质 |
CN114740428A (zh) * | 2022-04-13 | 2022-07-12 | 厦门大学 | 一种基于改进三角形法的水下被动定位方法与*** |
CN115150712A (zh) * | 2022-06-07 | 2022-10-04 | 中国第一汽车股份有限公司 | 一种车载麦克风***及汽车 |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8995228B2 (en) | 2012-10-22 | 2015-03-31 | Symbol Technologies, Inc. | Ultrasonic locationing using only time difference of arrival measurements |
US11795032B2 (en) | 2018-11-13 | 2023-10-24 | Otis Elevator Company | Monitoring system |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10227849A (ja) * | 1997-02-14 | 1998-08-25 | Fuji Xerox Co Ltd | 音源位置計測装置、カメラ撮影制御装置、音源位置記録装置及び音源位置計測方法 |
EP1331490A1 (en) * | 2000-10-02 | 2003-07-30 | Chubu Electric Power Co., Inc. | Sound source probing system |
CN1952684A (zh) * | 2005-10-20 | 2007-04-25 | 松下电器产业株式会社 | 利用麦克风定位声源的方法和装置 |
CN101201399A (zh) * | 2007-12-18 | 2008-06-18 | 北京中星微电子有限公司 | 一种声源定位方法及*** |
CN101295015A (zh) * | 2007-04-23 | 2008-10-29 | 财团法人工业技术研究院 | 声源定位***及声源定位方法 |
-
2010
- 2010-12-29 CN CN201010611087A patent/CN102033223B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10227849A (ja) * | 1997-02-14 | 1998-08-25 | Fuji Xerox Co Ltd | 音源位置計測装置、カメラ撮影制御装置、音源位置記録装置及び音源位置計測方法 |
EP1331490A1 (en) * | 2000-10-02 | 2003-07-30 | Chubu Electric Power Co., Inc. | Sound source probing system |
CN1952684A (zh) * | 2005-10-20 | 2007-04-25 | 松下电器产业株式会社 | 利用麦克风定位声源的方法和装置 |
CN101295015A (zh) * | 2007-04-23 | 2008-10-29 | 财团法人工业技术研究院 | 声源定位***及声源定位方法 |
CN101201399A (zh) * | 2007-12-18 | 2008-06-18 | 北京中星微电子有限公司 | 一种声源定位方法及*** |
Non-Patent Citations (3)
Title |
---|
***等: "基于麦克风阵列的语音信号实时时延估计", 《吉林大学学报(信息科学版)》 * |
范京等: "平面麦克风阵列估计声源的方向和距离", 《全国第4届信号和智能信息处理与应用学术会议论文集》 * |
邵怀宗等: "基于麦克风阵列的声源定位研究", 《云南民族大学学报(自然科学版)》 * |
Cited By (52)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103054669A (zh) * | 2012-12-25 | 2013-04-24 | 汪勇 | 一种耳塞式声音辅助感知装置 |
CN103473864B (zh) * | 2013-06-17 | 2016-02-03 | 西北工业大学 | 智能食堂打卡机的语音识别及指纹结算方法 |
CN103473864A (zh) * | 2013-06-17 | 2013-12-25 | 西北工业大学 | 智能食堂打卡机的语音识别及指纹结算方法 |
CN103310786A (zh) * | 2013-06-17 | 2013-09-18 | 西北工业大学 | 智能食堂打卡机的三维语音识别方法 |
CN103310786B (zh) * | 2013-06-17 | 2015-10-21 | 西北工业大学 | 智能食堂打卡机的三维语音识别方法 |
CN104422922A (zh) * | 2013-08-19 | 2015-03-18 | 中兴通讯股份有限公司 | 一种移动终端实现声源定位的方法及装置 |
CN103439689B (zh) * | 2013-08-21 | 2015-12-23 | 大连理工大学 | 一种分布式麦克风阵列中麦克风位置估计*** |
CN103439689A (zh) * | 2013-08-21 | 2013-12-11 | 大连理工大学 | 一种分布式麦克风阵列中麦克风位置估计*** |
CN103760479B (zh) * | 2014-01-10 | 2016-11-02 | 三峡大学 | 便携式电晕定位仪 |
CN103760479A (zh) * | 2014-01-10 | 2014-04-30 | 三峡大学 | 便携式电晕定位仪 |
CN103714551B (zh) * | 2014-01-20 | 2017-02-08 | 哈尔滨工业大学 | 基于区间消去的汽车仪表自动检测视觉***的相机外参数标定方法 |
CN103714551A (zh) * | 2014-01-20 | 2014-04-09 | 哈尔滨工业大学 | 基于区间消去的汽车仪表自动检测视觉***的相机外参数标定方法 |
CN104977560A (zh) * | 2014-04-04 | 2015-10-14 | 株式会社理光 | 移动设备、定位方法和定位*** |
CN104977560B (zh) * | 2014-04-04 | 2017-12-05 | 株式会社理光 | 移动设备、定位方法和定位*** |
CN103995252A (zh) * | 2014-05-13 | 2014-08-20 | 南京信息工程大学 | 一种三维空间声源定位方法 |
CN104035065A (zh) * | 2014-06-23 | 2014-09-10 | 河北工业大学 | 基于主动旋转的声源定向装置及其应用方法 |
CN104200805B (zh) * | 2014-08-30 | 2018-01-19 | 长城汽车股份有限公司 | 汽车驾驶员语音助手 |
CN104200805A (zh) * | 2014-08-30 | 2014-12-10 | 长城汽车股份有限公司 | 汽车驾驶员语音助手 |
CN105611014A (zh) * | 2014-11-21 | 2016-05-25 | 中兴通讯股份有限公司 | 移动终端通话语音降噪方法及装置 |
WO2016078369A1 (zh) * | 2014-11-21 | 2016-05-26 | 中兴通讯股份有限公司 | 移动终端通话语音降噪方法及装置、存储介质 |
CN104502893B (zh) * | 2014-12-10 | 2017-05-10 | 北京智谷睿拓技术服务有限公司 | 定位方法、定位装置以及用户设备 |
CN104502893A (zh) * | 2014-12-10 | 2015-04-08 | 北京智谷睿拓技术服务有限公司 | 定位方法、定位装置以及用户设备 |
CN104459625A (zh) * | 2014-12-14 | 2015-03-25 | 南京理工大学 | 基于轨道移动双麦克风阵列的声源定位装置及方法 |
CN104459625B (zh) * | 2014-12-14 | 2017-07-21 | 南京理工大学 | 基于轨道移动双麦克风阵列的声源定位装置及方法 |
CN104572009A (zh) * | 2015-01-28 | 2015-04-29 | 合肥联宝信息技术有限公司 | 一种自适应外界环境的音频控制方法及装置 |
CN104848478A (zh) * | 2015-04-28 | 2015-08-19 | 广东美的制冷设备有限公司 | 空调器中语音指令识别阈值的调整方法和*** |
CN104848478B (zh) * | 2015-04-28 | 2017-11-10 | 广东美的制冷设备有限公司 | 空调器中语音指令识别阈值的调整方法和*** |
CN104991573A (zh) * | 2015-06-25 | 2015-10-21 | 北京品创汇通科技有限公司 | 一种基于声源阵列的定位跟踪方法及其装置 |
CN105629202A (zh) * | 2016-03-23 | 2016-06-01 | 浙江工业大学 | 伞形可重组三维传声器阵列识别声源三维坐标的方法 |
CN105629202B (zh) * | 2016-03-23 | 2017-12-29 | 浙江工业大学 | 伞形可重组三维传声器阵列识别声源三维坐标的方法 |
CN107592600B (zh) * | 2016-07-06 | 2024-04-02 | 深圳市三诺声智联股份有限公司 | 一种基于分布式麦克风的拾音筛选方法及拾音装置 |
CN107592600A (zh) * | 2016-07-06 | 2018-01-16 | 深圳市三诺声智联股份有限公司 | 一种基于分布式麦克风的拾音筛选方法及拾音装置 |
CN107767137A (zh) * | 2016-08-23 | 2018-03-06 | ***通信有限公司研究院 | 一种信息处理方法、装置及终端 |
CN106535128B (zh) * | 2016-11-22 | 2020-04-03 | 四川九洲电器集团有限责任公司 | 一种多点定位方法及装置 |
CN106535128A (zh) * | 2016-11-22 | 2017-03-22 | 四川九洲电器集团有限责任公司 | 一种多点定位方法及装置 |
CN106911985B (zh) * | 2017-03-20 | 2018-09-18 | 广州视源电子科技股份有限公司 | 时间差确定的方法和装置 |
WO2018171151A1 (zh) * | 2017-03-20 | 2018-09-27 | 广州视源电子科技股份有限公司 | 时间差确定的方法和装置 |
CN106911985A (zh) * | 2017-03-20 | 2017-06-30 | 广州视源电子科技股份有限公司 | 时间差确定的方法和装置 |
CN109068235A (zh) * | 2017-06-12 | 2018-12-21 | 田中良 | 用于精确计算声音在麦克风阵列处的到达方向的方法 |
CN107154256B (zh) * | 2017-06-27 | 2023-11-24 | 山东省计算中心(国家超级计算济南中心) | 基于声源定位的声掩蔽***及自适应调整方法 |
CN107154256A (zh) * | 2017-06-27 | 2017-09-12 | 山东省计算中心(国家超级计算济南中心) | 基于声源定位的声掩蔽***及自适应调整方法 |
CN107800967A (zh) * | 2017-10-30 | 2018-03-13 | 维沃移动通信有限公司 | 一种拍摄方法及移动终端 |
CN110390947A (zh) * | 2018-04-23 | 2019-10-29 | 北京京东尚科信息技术有限公司 | 声源位置的确定方法、***、设备和存储介质 |
CN110390947B (zh) * | 2018-04-23 | 2024-04-05 | 北京京东尚科信息技术有限公司 | 声源位置的确定方法、***、设备和存储介质 |
CN108416996A (zh) * | 2018-04-28 | 2018-08-17 | 上海市城市建设设计研究总院(集团)有限公司 | 基于城市道路附属设施的声源定位***及其使用方法 |
CN109019349A (zh) * | 2018-06-19 | 2018-12-18 | 上海振华重工(集团)股份有限公司 | 故障检测方法、检测器、计算机存储介质及起重机 |
CN110620836A (zh) * | 2019-09-06 | 2019-12-27 | 中国民航大学 | 一种手机听声定位技术 |
CN110620836B (zh) * | 2019-09-06 | 2021-01-15 | 中国民航大学 | 一种手机听声定位方法 |
CN112964256B (zh) * | 2019-12-13 | 2024-02-27 | 佛山市云米电器科技有限公司 | 室内定位方法、智能家电设备及计算机可读存储介质 |
CN112964256A (zh) * | 2019-12-13 | 2021-06-15 | 佛山市云米电器科技有限公司 | 室内定位方法、智能家电设备及计算机可读存储介质 |
CN114740428A (zh) * | 2022-04-13 | 2022-07-12 | 厦门大学 | 一种基于改进三角形法的水下被动定位方法与*** |
CN115150712A (zh) * | 2022-06-07 | 2022-10-04 | 中国第一汽车股份有限公司 | 一种车载麦克风***及汽车 |
Also Published As
Publication number | Publication date |
---|---|
CN102033223B (zh) | 2012-10-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102033223B (zh) | 使用麦克风阵列对声源定位的方法 | |
Lo et al. | Broadband passive acoustic technique for target motion parameter estimation | |
CN106054134A (zh) | 一种基于tdoa的快速定位方法 | |
CN104035065A (zh) | 基于主动旋转的声源定向装置及其应用方法 | |
EP2081050B1 (en) | Multiple object localisation with a network of receivers | |
CN111487589B (zh) | 一种基于多源传感器网络的目标落点定位方法 | |
Lo et al. | Aircraft flight parameter estimation using acoustic multipath delays | |
CN110703202B (zh) | 基于多声学波浪滑翔机和水面无人艇的水下脉冲声定位*** | |
CN102833016B (zh) | 一种在多径环境下针对同源信号时差的测定方法 | |
CN109188362A (zh) | 一种麦克风阵列声源定位信号处理方法 | |
CN105353351A (zh) | 一种基于多信标到达时间差改进型定位方法 | |
Sun et al. | Array geometry calibration for underwater compact arrays | |
Valente et al. | Real-time TDOA measurements of an underwater acoustic source | |
Zhang et al. | A passive acoustic positioning algorithm based on virtual long baseline matrix window | |
EP3882639B1 (en) | Acoustic air data system with radially paired receivers | |
CN109254265A (zh) | 一种基于麦克风阵列的鸣笛车辆定位方法 | |
Luo et al. | Accurate underwater localization through phase difference | |
Lebarbenchon et al. | Evaluation of an open-source implementation of the SRP-PHAT algorithm within the 2018 LOCATA challenge | |
Lo et al. | Flight path estimation using frequency measurements from a wide aperture acoustic array | |
Flückiger et al. | Optimization of receiver arrangements for passive emitter localization methods | |
Majid et al. | Lightweight audio source localization for swarm robots | |
RU2620925C1 (ru) | Способ определения направления и дальности до источника сигналов | |
Nonsakhoo et al. | Angle of arrival estimation by using stereo ultrasonic technique for local positioning system | |
CN105917244B (zh) | 将脉冲源定位在色散介质中的方法 | |
RU2625094C1 (ru) | Способ определения пеленга и дальности до источника сигналов |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20121003 Termination date: 20121229 |