CN108008355B - 基于四元正交麦克风阵列的稳健声源定位方法 - Google Patents

基于四元正交麦克风阵列的稳健声源定位方法 Download PDF

Info

Publication number
CN108008355B
CN108008355B CN201710362998.3A CN201710362998A CN108008355B CN 108008355 B CN108008355 B CN 108008355B CN 201710362998 A CN201710362998 A CN 201710362998A CN 108008355 B CN108008355 B CN 108008355B
Authority
CN
China
Prior art keywords
error
pitch angle
sound source
array
sound
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.)
Active
Application number
CN201710362998.3A
Other languages
English (en)
Other versions
CN108008355A (zh
Inventor
晏燕
陈华伟
宋宫琨琨
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201710362998.3A priority Critical patent/CN108008355B/zh
Publication of CN108008355A publication Critical patent/CN108008355A/zh
Application granted granted Critical
Publication of CN108008355B publication Critical patent/CN108008355B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S5/00Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/18Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using ultrasonic, sonic, or infrasonic waves
    • G01S5/20Position of source determined by a plurality of spaced direction-finders
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S5/00Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/18Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using ultrasonic, sonic, or infrasonic waves
    • G01S5/22Position of source determined by co-ordinating a plurality of position lines defined by path-difference measurements

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

本发明公开了基于四元正交麦克风阵列的稳健声源定位方法,包括以下步骤:(1)计算四元正交阵引入幅度误差后各声强估计方法的绝对误差;(2)计算四元正交阵引入相位误差后各声强估计方法的绝对误差;(3)确定实际应用中最佳的声源定位方法;(4)针对不同频率的信号,确定合理的阵列尺寸范围;(5)确定合理的阵列摆放位置。本发明用于确定实际应用中最佳的声源定位方法,通过对比引入失配误差后,不同阵列尺寸下的目标声源定位精度,确定合适的阵列尺寸,最后在不同信号频率下,分析俯仰角变化对方位角定位精度影响,给出实际使用时阵列安装位置的建议,对于实际的工程应用具有指导意义。

Description

基于四元正交麦克风阵列的稳健声源定位方法
技术领域
本发明涉及基于四元正交麦克风阵列的稳健声源定位方法,属于声源定位的技术领域。
背景技术
在利用正交麦克风阵列进行目标声源定位时,定位误差主要来自于两个方面,阵元失配误差和有限差分误差,而目前对于正交阵列声强估计方法的研究中,只分析了单纯由差分误差引起的方向误差,没有考虑到失配误差的影响。本发明在存在失配误差的情况下,计算各声强估计方法对应的方位角绝对误差和俯仰角绝对误差,通过分析失配误差带来的影响,给出适用于目标声源定位正交麦克风阵列的稳健设计方法,包括定位方法的合理选择,阵列尺寸的选择以及阵列位置的摆放。
发明内容
为了解决上述存在的问题,本发明公开了一种基于四元正交麦克风阵列的稳健声源定位方法,其具体技术方案如下:
基于四元正交麦克风阵列的稳健声源定位方法,包括以下步骤:
(1)计算四元正交阵引入幅度误差后各声强估计方法的绝对误差,
方位角绝对误差:
Figure GDA0003001104690000011
俯仰角绝对误差:
Figure GDA0003001104690000012
其中:
Figure GDA0003001104690000013
是声源与x轴正方向的夹角即方位角,θ是声源与z轴正方向的夹角即俯仰角,ηIox、ηIoy、ηIoz分别为x轴、y轴、z轴方向上声强的相对误差;
(2)计算四元正交阵引入相位误差后各声强估计方法的绝对误差;
(3)确定实际应用中最佳的声源定位方法,
进行目标声源方位角估计时:
振速估计采用
Figure GDA0003001104690000021
声压估计采用pa=(p1+p2+p3)/3, (4)
进行目标声源俯仰角估计时:
振速估计采用
Figure GDA0003001104690000022
声压估计采用pA=(p1+p2+p3+p4)/4,(6)
其中:j表示虚数单位,p1、p2、p3、p4分别为四个阵元接收到的信号,vx、vy、vz分别为x、y、 z轴方向上的振速,c表示声速,
Figure GDA0003001104690000023
表示波数,f0表示信号频率,ρ表示空气密度,h表示阵列间距的一半;vG表示在进行目标声源定位时,所采取的振速估计方法,pa表示在计算目标声源方位角时,所采取的声压估计方法,其为p1、p2、p3的声压平均值,pA表示在计算目标声源俯仰角时,所采取的声压估计方法,为p1、p2、p3、p4的声压平均值;
(4)针对不同频率的信号,确定合理的阵列尺寸范围;
(5)确定合理的阵列摆放位置:应避开俯仰角等于0°和180°附近的角度。
所述步骤(1)中,所述四元正交阵包含1、2、3、4号麦克风,其中1号麦克风位于原点O 处,2、3、4号麦克风依次位于x轴正方向、y轴正方向、z轴正方向,且与1号麦克风的距离为2h, x、y、z轴两两相互正交,建立三维坐标系。
声压估计有四种方法,分别是:
Figure GDA0003001104690000031
pO表示取1号麦克风上的声压作为估计值,pA表示取1、2、3、4号麦克风上的声压均值作为估计值,pW表示取1、2、3、4号麦克风上的声压加权平均值作为估计值,pT表示采用泰勒近似的方法,得到声压的估计值;
振速估计有两种方法,分别是:
Figure GDA0003001104690000032
Figure GDA0003001104690000033
得到振速和声压之后,利用I=Re{pv*}/2计算出各方向上的声强,其中*号表示求共轭,并通过声强与角度的关系,得到方位角
Figure GDA0003001104690000034
和俯仰角θ的估计,
Figure GDA0003001104690000041
Figure GDA0003001104690000042
其中,Iox、Ioy、Ioz分别为x、y、z轴方向上的声强,
振速和声压两两组合计算声强,得到八种声强估计方法,即Re{pAvG *}/2、Re{pOvG *}/2、 Re{pWvG *}/2、Re{pTvG *}/2、Re{pAvT *}/2、Re{pOvT *}/2、Re{pWvT *}/2、 Re{pTvT *}/2,八种方法依次记为A.G、O.G、W.G、T.G、A.T、O.T、W.T、T.T,前四种方法统称为方法G,后四种方法统称为方法T。
所述步骤(1)中,方法G引入幅度误差后,以通道1为基准,假设通道2,3,4与通道1之间存在的通道幅度误差分别为k12、k13、k14,通过计算分析得到方位角绝对误差和俯仰角绝对误差,方位角绝对误差和俯仰角绝对误差均在允许误差范围内;
方法T引入幅度误差后,通过计算分析得到,当方位角和俯仰角满足以下关系:θ=0°或θ=180°或
Figure GDA0003001104690000043
Figure GDA0003001104690000044
Figure GDA0003001104690000045
时,方位角绝对误差均超出允许误差,此时声源的方位角失去了定向的准确性,当方位角和俯仰角满足以下关系:
Figure GDA0003001104690000046
Figure GDA0003001104690000047
Figure GDA0003001104690000048
时,俯仰角绝对误差均超出允许误差,此时声源的俯仰角失去了定向的准确性。
所述步骤(2)中,方法G引入相位误差后,以通道1为基准,假设通道2,3,4与通道1之间存在的相位误差分别为Δφ12、Δφ13、Δφ14,通过计算分析得到,在θ=0°或θ=180°附近时,方位角绝对误差均超出允许误差,俯仰角绝对误差均在允许误差范围内;
方法T引入相位误差后,通过计算分析得到,方位角绝对误差在θ=0°或θ=180°时取值超出允许误差范围,而俯仰角绝对误差均在允许误差范围内。
所述允许误差为小于或等于5度的误差范围。
所述步骤(3)中,由于在存在失配误差的情况下,方法T在很多角度上失去了定向的准确性,而实际中失配误差很难避免,因此方法T对应的四种方法在实际应用中不应被考虑;方法G对应的四种声强估计方法,振速估计均采用同一种方法,区别在于声压估计;在振速估计方法确定的基础上,分析不同声压估计的定位性能;
方位角:在进行方位角估计时,只利用1、2、3号麦克风上的声压信息去估计方位角,整体要比把四个麦克风上的声压信息都利用上时的定位精度高,
仅利用三个麦克风的声压估计方法:
Figure GDA0003001104690000051
在Pa、Pw、Pt三种方法中,声压估计采用pa时,方位角定位精度最高;
俯仰角:在进行俯仰角估计时,声压估计采用pA时,俯仰角定位精度最高,
因此,声强估计方法是振速估计采用vG,计算方位角时,声压估计采用pa,计算俯仰角时,声压估计采用pA
所述步骤(4)中,对比不同信号频率下,阵列尺寸变化对定位精度的影响,当信号频率范围为 300Hz~1500Hz时,阵列尺寸为5cm~9cm,信号频率范围为1500Hz~2500Hz时,阵列尺寸为3cm~6cm,信号频率范围为2500Hz~3500Hz时,阵列尺寸为1cm~3cm。
所述步骤(5)中,通过分析,由于相位误差的影响,在俯仰角等于0°或180°附近时,方位角失去了定向的准确性,并且随着相位误差和信号频率的变化,无法定向的区域范围也在发生变化;避开无法定向的区域,选择合理的俯仰角范围,选择阵列尺寸为4cm,当相位误差为0.5°,信号频率为 300Hz~1500Hz时,若要求方位角的定位精度在5°以内,则要避开俯仰角为0°~10°及170°~180°度的角度,即俯仰角的合理范围是10°~170°;当信号频率为1500Hz~2500Hz时,俯仰角的合理范围是3°~ 177°,当信号频率范围为2500Hz~3500Hz时,俯仰角的合理范围是1°~179°。
本发明的有益效果是:
本发明分析了失配误差对定位精度的影响,给出一种基于四元正交麦克风阵列的稳健声源定位方法,包括定位方法的合理选择,阵列尺寸的选择,阵列位置的安装建议,对于实际的工程应用具有指导意义。
附图说明
图1为阵列结构和坐标定义示意图,
图2(a)为方法G引入幅度误差后方位角绝对误差的全局表现图,
图2(b)为方法G引入幅度误差后俯仰角绝对误差的全局表现图,
图3(a)为方法T引入幅度误差后方位角绝对误差大于30度的俯视图,
图3(b)为方法T引入幅度误差后俯仰角绝对误差大于30度的俯视图,
图4(a)为方法G引入相位误差后方位角绝对误差的全局表现图,
图4(b)为方法G引入相位误差后俯仰角绝对误差的全局表现图,
图5(a)为方法T引入相位误差后方位角绝对误差的全局表现图,
图5(b)为方法T引入相位误差后俯仰角绝对误差的全局表现图,
图6(a)为方位角误差对比图,
图6(b)为俯仰角误差对比图,
图7为方位角误差和俯仰角误差随阵列尺寸改变的变化图,
图8为方位角误差随俯仰角改变的变化图。
具体实施方式
下面结合附图和具体实施方式,进一步阐明本发明。应理解下述具体实施方式仅用于说明本发明而不用于限制本发明的范围。
步骤1:计算四元正交阵引入幅度误差后各声强估计方法的绝对误差。
本发明是在如图1所示的麦克风阵列基础上进行的。四个全向麦克风等间隔的分布在图1所示的三棱锥上。通过计算,四个阵元接收到的信号分别为:
Figure GDA0003001104690000061
Figure GDA0003001104690000062
Figure GDA0003001104690000063
Figure GDA0003001104690000064
为了获得目标声源定向的绝对误差,我们引入:
Figure GDA0003001104690000071
Figure GDA0003001104690000072
Figure GDA0003001104690000073
其中R为声源到麦克风中心O点的距离,
Figure GDA0003001104690000079
分别为x轴、y轴、z轴方向上声强的相对误差,Iox、Ioy、Ioz分别表示x轴、y轴、z轴方向上声强的真实值,
Figure GDA0003001104690000074
Figure GDA0003001104690000075
则表示对应的估计值。
由于方法G对应的四种声强估计方法,目标声源定向绝对误差的整体表现类似,为简化计算,我们选择当声压估计为pO的情况进行计算和分析。最终通过公式(1)、(2)计算得到,加入幅度误差后,方法G方位角和俯仰角的绝对误差分别为:
Figure GDA0003001104690000076
Figure GDA0003001104690000077
其中:
Figure GDA0003001104690000078
分析:
1)方位角绝对误差:分析式(13),当θ=0°或θ=180°时,虽然使得分母为0,但与此同时分子也为零,所以不会导致此时误差很大,因此整体不会出现误差很大的点。
2)俯仰角绝对误差:式(14)中的各项均取值较小,因此整体也不会有误差较大的点。
我们选择波速尺寸积kh为0.2,方位角变化范围为-180°~180°,俯仰角变化范围为0°~180°,设阵元间的相位特性一致,各阵元间的幅度误差扰动服从0~0.03的均匀分布,在全空间内O.G方法方位角绝对误差和俯仰角绝对误差表现情况如图2(a),图2(b)所示。能够看出,表现结果与上文理论分析一致。
由于方法T对应的四种估计方法,目标声源定向绝对误差的整体表现类似,为简化计算,我们同样选择声压估计为pO的情况进行计算和分析。最终通过公式(1)、(2)计算得到,加入幅度误差后,方法T的方位角和俯仰角的绝对误差分别为:
Figure GDA0003001104690000081
Figure GDA0003001104690000082
其中:
Figure GDA0003001104690000091
分析:
1)方位角绝对误差:分析式(16),对于整个式子来说,其分母有sinθ这一项,如果取值趋于0,就会导致整个公式趋于无穷大。能够求出此时的θ为0°或者180°。对式(16)中的每一项进行分析:若θxz=π/4,那么大括号内就有一项分母为0,则会导致公式结果较大,根据θxz与方位角和俯仰角的关系,我们能够求出此时
Figure GDA0003001104690000092
同样的,若θyz=π/4,也会导致公式结果变大,此时
Figure GDA0003001104690000093
当θxy=π/4,也会使公式结果变大,此时
Figure GDA0003001104690000094
也即
Figure GDA0003001104690000095
Figure GDA0003001104690000096
因此当方位角和俯仰角满足以下关系:θ=0°或θ=180°,
Figure GDA0003001104690000097
方位角绝对误差较大,此时声源的方位角失去了定向的准确性。
2)俯仰角绝对误差:分析式(17),同方位角误差分析类似,若θxz=π/4,那么大括号内就有一项分母为0,则会导致公式结果较大,此时
Figure GDA0003001104690000098
同样的,若θyz=π/4,也会导致公式结果变大,此时
Figure GDA0003001104690000099
当θxy=π/4,也会导致公式结果变大,此时
Figure GDA00030011046900000910
也即
Figure GDA00030011046900000911
Figure GDA00030011046900000912
因此,当方位角和俯仰角满足以下关系:
Figure GDA00030011046900000913
俯仰角绝对误差较大,此时声源的俯仰角失去了定向的准确性。
我们选择波速尺寸积kh为0.2,方位角变化范围为-180°~180°,俯仰角变化范围为0°~180°,设阵元间的相位特性一致,各阵元间的幅度误差扰动服从0~0.03的均匀分布,在全空间内O.T方法方位角误差和俯仰角误差变化的俯视图如图4(a),图4(b)所示。能够看出,表现结果与上文理论分析一致。
步骤2:计算四元正交阵引入相位误差后各声强估计方法的绝对误差。
方法G加入通道相位误差后,通过计算得到:
Figure GDA00030011046900000914
Figure GDA0003001104690000101
分析:
1)方位角绝对误差:分析式(19),当θ=0°或θ=180°时,分母为0,此时由于相位误差的存在导致分子不为零,因此方位角绝对误差在θ=0°或θ=180°的附近取值较大。
2)俯仰角绝对误差:式(20)中的各项均取值较小,因此整体不会有误差取值较大的点。
我们选择波速尺寸积kh为0.2,方位角变化范围为-180°~180°,俯仰角变化范围为0°~180°,设阵元间的幅度特性一致,各阵元间的相位误差扰动服从0°~1°的均匀分布,在全空间内O.G方法方位角误差和俯仰角误差变化如图3(a),图3(b)所示。能够看出,表现结果与上文理论分析一致。
方法T加入通道相位误差后,通过计算得到
Figure GDA0003001104690000102
Figure GDA0003001104690000111
分析:
1)方位角绝对误差:分析式(21),能够看出当sinθ=0时方位角绝对误差取值较大,又由于Δφ很小,导致[sin(E-Δφ12)-sin(E)],[sin(F-Δφ13)-sin(F)],[sin(Q-Δφ14)-sin(Q)] 很小,那么大括号中导致误差较大的项起不到主导作用,因而不会出现取值较大的情况。
2)俯仰角绝对误差:分析式(22),各项结果的取值都较小,不会出现使俯仰角绝对误差取值过大的情况。
我们选择波速尺寸积kh为0.2,方位角变化范围为-180°~180°,俯仰角变化范围为0°~180°,设阵元间的幅度特性一致,各阵元间的相位误差扰动服从0°~1°的均匀分布,在全空间内O.T方法各方位角误差和俯仰角误差变化如图5(a),图5(b)所示。能够看出,表现结果与上文理论分析一致。
步骤3:确定实际应用中最佳的声源定位方法
通过分析及附图所示结果,能够得到,方法T在存在失配误差的情况下,在很多角度上失去了定向的准确性,而失配误差在实际应用中又很难避免,因此方法T对应的四种估计方法在实际应用中不应被考虑。
对比方法G四种估计方法的性能优劣:我们选择波速尺寸积kh的变化范围为0.1~1.2,方位角变化范围为-180°~180°,俯仰角变化范围为0°~180°,各阵元间的幅度误差扰动服从0~0.03的均匀分布,各阵元间的相位误差扰动服从0°~1°的均匀分布,进行100次蒙特卡洛实验。对这次结果求取均值和方差, 利用均方根误差作为评价指标。图6(a)表示,进行方位角估计,对比声压采用pa,pw,pt,pO, pA,pW,pT这七种方法性能的优劣,能够看出,整体上,只利用1,2,3号麦克风的声压信息,比 4个麦克风声压信息全部利用性能要好。且精度较高的是声压估计采用pa。图6(b)表示,进行俯仰角估计,对比声压采用pO,pA,pW,pT这四种方法性能的优劣,能够看出,精度较高的是声压估计采用pA
因此,本发明给出最终的声强估计方案,首先振速估计采用vG,在进行方位角估计时声压估计采用pa,进行俯仰角估计时,声压估计采用pA
步骤4:针对不同频率的信号,确定合理的阵列尺寸范围。
下面给出可听音范围内,阵列尺寸的合理选择。频率范围为20Hz~20000Hz,方位角变化范围为 -180°~180°,俯仰角变化范围为0°~180°,各阵元间的相位误差扰动服从0°~1°的均匀分布,各阵元间的幅度误差扰动服从0~0.03的均匀分布。图7表示频率在1900Hz~2300Hz这个范围时,方位角均方根误差和俯仰角均方根误差随阵列尺寸改变的变化情况。能够看出,在这个频段内,方位角均方根误差在阵列尺寸D为3cm~7cm时,误差较小,俯仰角均方根误差在阵列尺寸D为2cm~4cm时,误差较小。因此对于1900Hz~2300Hz这个频率范围,我们最终选择的阵列尺寸为3cm~4cm,其他频段类似分析,得到下表:
Figure GDA0003001104690000121
步骤5:确定合理的阵列摆放位置
频率变化范围为20Hz~20000Hz,各阵元间的幅度误差扰动服从0~0.03的均匀分布,相位误差为0.5°, 1°,1.5°,2°时,当阵列尺寸为4cm时,如果要求方位角定位精度在5°以内,俯仰角的合理范围如下表所示。图8表示,当各阵元间相位误差为2°,信号频率范围为700Hz~1100Hz时,方位角绝对误差随俯仰角变化的表现情况,能够看出,如果要求定位精度在5°以内,俯仰角的合理范围是28°~152°。其他频段类似分析,得到以下各表:
1)当相位误差为0.5°时
Figure GDA0003001104690000131
2)当相位误差为1°时
Figure GDA0003001104690000132
3)当相位误差为1.5°时
Figure GDA0003001104690000133
Figure GDA0003001104690000141
4)当相位误差为2°时
Figure GDA0003001104690000142
本发明方案所公开的技术手段不仅限于上述技术手段所公开的技术手段,还包括由以上技术特征任意组合所组成的技术方案。
以上述依据本发明的理想实施例为启示,通过上述的说明内容,相关工作人员完全能够在不偏离本项发明技术思想的范围内,进行多样的变更以及修改。本项发明的技术性范围并不局限于说明书上的内容,必须要根据权利要求范围来确定其技术性范围。

Claims (9)

1.基于四元正交麦克风阵列的稳健声源定位方法,其特征在于,包括以下步骤:
(1)计算四元正交阵引入幅度误差后各声强估计方法的绝对误差,
方位角绝对误差:
Figure FDA0003001104680000011
俯仰角绝对误差:
Figure FDA0003001104680000012
其中:
Figure FDA0003001104680000013
是声源与x轴正方向的夹角即方位角,θ是声源与z轴正方向的夹角即俯仰角,
Figure FDA0003001104680000014
Figure FDA0003001104680000015
分别为x轴、y轴、z轴方向上声强的相对误差;
(2)计算四元正交阵引入相位误差后各声强估计方法的绝对误差;
(3)确定实际应用中最佳的声源定位方法,
进行目标声源方位角估计时:
振速估计采用
Figure FDA0003001104680000016
声压估计采用pa=(p1+p2+p3)/3,
进行目标声源俯仰角估计时:
振速估计采用
Figure FDA0003001104680000017
声压估计采用pA=(p1+p2+p3+p4)/4,
其中:j表示虚数单位,p1、p2、p3、p4分别为四个阵元接收到的信号,vx、vy、vz分别为x、y、z轴方向上的振速,c表示声速,
Figure FDA0003001104680000021
表示波数,f0表示信号频率,ρ表示空气密度,h表示阵列间距的一半;vG表示在进行目标声源定位时,所采取的振速估计方法,pa表示在计算目标声源方位角时,所采取的声压估计方法,其为p1、p2、p3的声压平均值,pA表示在计算目标声源俯仰角时,所采取的声压估计方法,为p1、p2、p3、p4的声压平均值;
(4)针对不同频率的信号,确定合理的阵列尺寸范围;
(5)确定合理的阵列摆放位置:避开俯仰角等于0°和180°附近的角度。
2.根据权利要求1所述的基于四元正交麦克风阵列的稳健声源定位方法,其特征在于,所述步骤(1)中,所述四元正交阵包含1、2、3、4号麦克风,其中1号麦克风位于原点O处,2、3、4号麦克风依次位于x轴正方向、y轴正方向、z轴正方向,且与1号麦克风的距离为2h,x、y、z轴两两相互正交,建立三维坐标系。
3.根据权利要求2所述的基于四元正交麦克风阵列的稳健声源定位方法,其特征在于,声压估计有四种方法,分别是:
pO=p1
pA=(p1+p2+p3+p4)/4
pW=(3p1+p2+p3+p4)/6
pT=(p2+p3+p4-p1)/2,
pO表示取1号麦克风上的声压作为估计值,pA表示取1、2、3、4号麦克风上的声压均值作为估计值,pW表示取1、2、3、4号麦克风上的声压加权平均值作为估计值,pT表示采用泰勒近似的方法,得到声压的估计值;
振速估计有两种方法,分别是:
Figure FDA0003001104680000031
Figure FDA0003001104680000032
得到振速和声压之后,利用I=Re{pv*}/2计算出各方向上的声强,其中*号表示求共轭,并通过声强与角度的关系,得到方位角
Figure FDA0003001104680000033
和俯仰角θ的估计,
Figure FDA0003001104680000034
Figure FDA0003001104680000035
其中,Iox、Ioy、Ioz分别为x、y、z轴方向上的声强,
振速和声压两两组合计算声强,得到八种声强估计方法,即Re{pAvG *}/2、Re{pOvG *}/2、Re{pWvG *}/2、Re{pTvG *}/2、Re{pAvT *}/2、Re{pOvT *}/2、Re{pWvT *}/2、Re{pTvT *}/2,八种方法依次记为A.G、O.G、W.G、T.G、A.T、O.T、W.T、T.T,A.G、O.G、W.G、T.G统称为方法G,A.T、O.T、W.T、T.T统称为方法T。
4.根据权利要求3所述的基于四元正交麦克风阵列的稳健声源定位方法,其特征在于,所述步骤(1)中,方法G引入幅度误差后,以通道1为基准,假设通道2,3,4与通道1之间存在的通道幅度误差分别为k12、k13、k14,通过计算分析得到方位角绝对误差和俯仰角绝对误差,方位角绝对误差和俯仰角绝对误差均在允许误差范围内;
方法T引入幅度误差后,通过计算分析得到,当方位角和俯仰角满足以下关系:θ=0°或θ=180°或
Figure FDA0003001104680000041
Figure FDA0003001104680000042
Figure FDA0003001104680000043
时,方位角绝对误差均超出允许误差,此时声源的方位角失去了定向的准确性,当方位角和俯仰角满足以下关系:
Figure FDA0003001104680000044
Figure FDA0003001104680000045
Figure FDA0003001104680000046
时,俯仰角绝对误差均超出允许误差,此时声源的俯仰角失去了定向的准确性。
5.根据权利要求4所述的基于四元正交麦克风阵列的稳健声源定位方法,其特征在于,所述步骤(2)中,方法G引入相位误差后,以通道1为基准,假设通道2,3,4与通道1之间存在的相位误差分别为Δφ12、Δφ13、Δφ14,通过计算分析得到,在θ=0°或θ=180°附近时,方位角绝对误差均超出允许误差,俯仰角绝对误差均在允许误差范围内;
方法T引入相位误差后,通过计算分析得到,方位角绝对误差在θ=0°或θ=180°时取值超出允许误差范围,而俯仰角绝对误差均在允许误差范围内。
6.根据权利要求4或5所述的基于四元正交麦克风阵列的稳健声源定位方法,其特征在于,所述允许误差为小于或等于5度的误差范围。
7.根据权利要求4所述的基于四元正交麦克风阵列的稳健声源定位方法,其特征在于,所述步骤(3)中,由于在存在失配误差的情况下,方法T在很多角度上失去了定向的准确性,而实际中失配误差很难避免,因此方法T对应的四种方法在实际应用中不应被考虑;方法G对应的四种声强估计方法,振速估计均采用同一种方法,区别在于声压估计;在振速估计方法确定的基础上,分析不同声压估计的定位性能;
方位角:在进行方位角估计时,只利用1、2、3号麦克风上的声压信息去估计方位角,整体要比把四个麦克风上的声压信息都利用上时的定位精度高,
仅利用三个麦克风的声压估计方法:
pa=(p1+p2+p3)/3
pw=(2p1+p2+p3)/4
pt=(p2+p3)/2,
在Pa、Pw、Pt三种方法中,声压估计采用pa时,方位角定位精度最高;
俯仰角:在进行俯仰角估计时,声压估计采用pA时,俯仰角定位精度最高,
因此,声强估计方法是振速估计采用vG,计算方位角时,声压估计采用pa,计算俯仰角时,声压估计采用pA
8.根据权利要求4所述的基于四元正交麦克风阵列的稳健声源定位方法,其特征在于,所述步骤(4)中,对比不同信号频率下,阵列尺寸变化对定位精度的影响,当信号频率范围为300Hz~1500Hz时,阵列尺寸为5cm~9cm,信号频率范围为1500Hz~2500Hz时,阵列尺寸为3cm~6cm,信号频率范围为2500Hz~3500Hz时,阵列尺寸为1cm~3cm。
9.根据权利要求2所述的基于四元正交麦克风阵列的稳健声源定位方法,其特征在于,所述步骤(5)中,通过分析,由于相位误差的影响,在俯仰角等于0°或180°附近时,方位角失去了定向的准确性,并且随着相位误差和信号频率的变化,无法定向的区域范围也在发生变化;避开无法定向的区域,选择合理的俯仰角范围,选择阵列尺寸为4cm,当相位误差为0.5°,信号频率为300Hz~1500Hz时,若要求方位角的定位精度在5°以内,则要避开俯仰角为0°~10°及170°~180°度的角度,即俯仰角的合理范围是10°~170°;当信号频率为1500Hz~2500Hz时,俯仰角的合理范围是3°~177°,当信号频率范围为2500Hz~3500Hz时,俯仰角的合理范围是1°~179°。
CN201710362998.3A 2017-05-18 2017-05-18 基于四元正交麦克风阵列的稳健声源定位方法 Active CN108008355B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710362998.3A CN108008355B (zh) 2017-05-18 2017-05-18 基于四元正交麦克风阵列的稳健声源定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710362998.3A CN108008355B (zh) 2017-05-18 2017-05-18 基于四元正交麦克风阵列的稳健声源定位方法

Publications (2)

Publication Number Publication Date
CN108008355A CN108008355A (zh) 2018-05-08
CN108008355B true CN108008355B (zh) 2021-06-15

Family

ID=62048819

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710362998.3A Active CN108008355B (zh) 2017-05-18 2017-05-18 基于四元正交麦克风阵列的稳健声源定位方法

Country Status (1)

Country Link
CN (1) CN108008355B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108844741B (zh) * 2018-07-05 2019-10-18 安徽大学 一种基于麦克风均匀面阵滤波的高速列车轴承故障诊断方法
CN109164416B (zh) * 2018-09-12 2022-10-14 南京信息工程大学 一种三平面五元传声器阵列的声源定位方法
CN110554358B (zh) * 2019-09-25 2022-12-13 哈尔滨工程大学 一种基于虚拟球阵列扩展技术的噪声源定位识别方法
CN111157951B (zh) * 2020-01-13 2022-02-25 东北大学秦皇岛分校 一种基于差分麦克风阵列的三维声源定位方法
CN111239256B (zh) * 2020-03-16 2021-10-08 吉林大学 适用于三维结构的声发射源定位方法
CN112129402B (zh) * 2020-08-21 2021-07-13 东风汽车集团有限公司 一种异响源探测装置
CN114488009A (zh) * 2021-12-01 2022-05-13 南京惠斯通智能科技有限责任公司 一种基于矢量传感器的声源定位方法
CN114151737B (zh) * 2021-12-08 2024-05-24 天津大学 用于管道泄漏检测的球形内检测器及定位系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101644765A (zh) * 2009-04-23 2010-02-10 中国科学院声学研究所 一种用于水声换能器线阵的幅相误差校正方法
CN101650220A (zh) * 2009-09-04 2010-02-17 合肥工业大学 互谱声强有限差分误差修正方法
RU2515179C1 (ru) * 2012-11-13 2014-05-10 Открытое Акционерное Общество "НИИ гидросвязи "Штиль" Способ определения направления на гидроакустический маяк-ответчик в условиях многолучевого распространения навигационного сигнала
CN104076331A (zh) * 2014-06-18 2014-10-01 南京信息工程大学 一种七元麦克风阵列的声源定位方法
CN104898086A (zh) * 2015-05-19 2015-09-09 南京航空航天大学 适用于微型麦克风阵列的声强估计声源定向方法
CN106501770A (zh) * 2016-10-26 2017-03-15 黑龙江大学 基于幅相误差阵列的远近场宽带混合源中近场源定位方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101644765A (zh) * 2009-04-23 2010-02-10 中国科学院声学研究所 一种用于水声换能器线阵的幅相误差校正方法
CN101650220A (zh) * 2009-09-04 2010-02-17 合肥工业大学 互谱声强有限差分误差修正方法
RU2515179C1 (ru) * 2012-11-13 2014-05-10 Открытое Акционерное Общество "НИИ гидросвязи "Штиль" Способ определения направления на гидроакустический маяк-ответчик в условиях многолучевого распространения навигационного сигнала
CN104076331A (zh) * 2014-06-18 2014-10-01 南京信息工程大学 一种七元麦克风阵列的声源定位方法
CN104898086A (zh) * 2015-05-19 2015-09-09 南京航空航天大学 适用于微型麦克风阵列的声强估计声源定向方法
CN106501770A (zh) * 2016-10-26 2017-03-15 黑龙江大学 基于幅相误差阵列的远近场宽带混合源中近场源定位方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
An Improved Method for 2D DOA Estimation Based On a Vector Sound Intensity Array;Jie Shi 等;《2006 8th international Conference on Signal Processing》;20070410;正文全文 *
Comparison of methods for processing acoustic intensity from orthogonal multimicrophone probes;Curtis P. Wiederhold 等;《Journal of the Acoustical Society of America》;20140313;第2841-2852页 *
三维压差式矢量水听器定向性能分析;石敏;《船舰科学技术》;20080831;第30卷(第4期);第65-68、75页 *
低空目标声测无源定向理论与算法研究;陈华伟;《中国优秀博硕士学位论文全文数据库 (博士) 工程科技Ⅱ辑》;20061115(第11期);正文全文 *
声强测量中相位失配误差及校正方法;董为民;《西南林学院学报》;19990630;第19卷(第2期);第131-134页 *

Also Published As

Publication number Publication date
CN108008355A (zh) 2018-05-08

Similar Documents

Publication Publication Date Title
CN108008355B (zh) 基于四元正交麦克风阵列的稳健声源定位方法
CN103308889B (zh) 复杂环境下被动声源二维doa估计方法
CN102804809B (zh) 音频源定位
CN111123192B (zh) 一种基于圆形阵列和虚拟扩展的二维doa定位方法
CN106501773B (zh) 基于差分阵列的声源方向定位方法
CN104360315A (zh) 基于LabVIEW的麦克风阵列声源定位方法及装置
Xu et al. Optimal sensor deployment for 3D AOA target localization
CN104237849A (zh) 一种双五元十字阵被动声定位融合方法
CN109932689A (zh) 一种适用于特定定位场景的任意阵列优化方法
CN104076331A (zh) 一种七元麦克风阵列的声源定位方法
CN104330768B (zh) 一种基于声矢量传感器的机动声源方位估计方法
CN107884741A (zh) 一种多球阵列多宽带声源快速定向方法
Ren et al. A novel multiple sparse source localization using triangular pyramid microphone array
CN106093866A (zh) 一种适用于空心球阵列的声源定位方法
CN112285647B (zh) 一种基于稀疏表示与重构的信号方位高分辨估计方法
JP2008249702A (ja) 音響測定装置及び音響測定方法
CN109164416B (zh) 一种三平面五元传声器阵列的声源定位方法
CN107202975B (zh) 一种二维矢量阵阵元姿态误差校正方法
CN113702903A (zh) 基于目标水下极低频矢量电磁场的阵列被动定位跟踪方法
CN109375197B (zh) 一种小尺寸矢量阵低频散射校正方法
CN108595758A (zh) 一种任意形式传感器阵列最优宽带波束图综合方法
Suo et al. Single acoustic vector sensor DOA enhanced by unsaturated bistable stochastic resonance with linear amplitude response constrained
CN108919176B (zh) 一种单矢量传感器增秩music测向技术
CN114488009A (zh) 一种基于矢量传感器的声源定位方法
CN111736111A (zh) 一种基于同心均匀分布双圆阵的相干源doa的估计方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant