CN108680901A - 一种新型的声源方位定位方法 - Google Patents
一种新型的声源方位定位方法 Download PDFInfo
- Publication number
- CN108680901A CN108680901A CN201810330568.8A CN201810330568A CN108680901A CN 108680901 A CN108680901 A CN 108680901A CN 201810330568 A CN201810330568 A CN 201810330568A CN 108680901 A CN108680901 A CN 108680901A
- Authority
- CN
- China
- Prior art keywords
- sound source
- microphone
- sound
- sound bearing
- localization method
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/18—Position-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/22—Position 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)
- Circuit For Audible Band Transducer (AREA)
Abstract
本发明公开了一种新型的声源方位定位方法,利用声源方位估计算法进行声源方位定位。建立四元传声器阵列,将声源方位估计算法引入到阵列模型中,在正演数据的基础上进行声源方位数据反演。仿真结果表明:该新型算法定位精度较高,误差较低,分别比传统的直三棱锥式声源方位定位算法的水平偏角和仰角误差降低了大约10%、1%,具有较好的效果。
Description
技术领域
本发明涉及一种声源方位估计方法,特别涉及一种新型的声源方位定位方法,属于声源定位领域。
背景技术
声源方位定位是指对空间内的声源进行位置和角度数据的确定。声源方位估计算法是根据时延估计值得到的关于声源位置的定位方程组,该定位方程组为非线性方程组。声源定位技术是一种通过釆集声源声音信号并对其进行处理,再根据相关算法确定被接受的信号相对于接受传感器来自什么方向和距离是多少,最后得到所接收声源的三维位置的被动声定位技术,有着较为广阔的应用前景。
具有代表性的是中国电子科技集团第三研究所研制并开发出的声源定位***,该***是由麦克风阵列和信号处理单元组成,其中信号处理单元是一个独立的包装,可以对麦克风阵列接收的信息进行相关算法的分析,测量声源的方位角、俯仰角和距离等信息。该***的测量范围为0°到360°,有效作用距离1000米,方位角误差±3°,俯仰角误差±6°,距离误差10%到30%。
我国的声源定位***研究还处于萌芽阶段,各种理论算法不够完善,软硬件的配置也不够高级,这还需要各领域的学者们进一步的学习和研究,使得我国的声源定位***研究能够快速发展。
发明内容
为解决传统的基于直三棱锥式声源方位定位精度较低的问题,本发明提供一种新型的声源方位定位方法,将声源方位估计算法运用到声源方位定位中,建立传声器阵列来模拟空间中声源的存在进行定位研究,以提高声源方位定位的水平。
本发明为解决上述技术问题采用以下技术方案:
本发明提供一种新型的声源方位定位方法,该方法通过构建四元传声器阵列,反演声源位置;所述四元传声器阵列包括第一至第四传声器,以四元传声器阵列的阵元中心为坐标原点,设各传声器到坐标原点的距离均为第一至第四传声器的坐标分别为 反演得到声源的坐标(x,y,z):z=rcosθ,为水平偏角,θ为仰角,其中,t21为声源传播到第二、第一传声器的相对时延值,t31为声源传播到第三、第一传声器的相对时延值,t41为声源传播到第四、第一传声器的相对时延值,c为声源传播速度。
作为本发明的进一步技术方案,声源到阵元中心的距离
作为本发明的进一步技术方案,t21=t2-t1,t1、t2分别为声源传播到第一、第二传声器的时间。
作为本发明的进一步技术方案,t41=t4-t1,t1、t4分别为声源传播到第一、第四传声器的时间。
作为本发明的进一步技术方案,t31=t3-t1,t1、t3分别为声源传播到第一、第三传声器的时间。
本发明采用以上技术方案与现有技术相比,具有以下技术效果:
1)弥补传统的基于直三棱锥式的声源方位定位算法的缺陷,通过利用声源方位估计算法运算量较低、定位精度较高的优点,有效解决了基于直三棱锥式声源方位定位算法定位精度较低的问题,增强了声源方位定位的时效性与精确性;
2)采用声源到达传声器的相对时间差,从算法上大大减小声音在传播过程中由于衰减和外界干扰对声源方位定位精度的影响;
3)新型的基于四元传声器阵列的声源方位定位算法能有效地对声源进行方位定位,达到了运算量较低、定位精度较高的效果,足以表明该声源方位估计算法是非常有效的。
附图说明
图1是四元传声器阵列模型;
图2是传统的直三棱锥式阵列模型。
具体实施方式
下面结合附图对本发明的技术方案做进一步的详细说明:
本发明的一种新型的声源方位定位方法,建立四元传声器阵列,利用声源方位估计算法,在正演数据的基础上进行声源方位数据反演,进行声源方位定位。
本发明通过构建如图1所示的四元传声器阵列,反演声源位置。声源传播速度为c,以四元传声器阵列的阵元中心为坐标原点(0,0,0),设各传声器到坐标原点的距离均为(d的取值根据四元传声器阵列的实际构建确定),四个传声器A、B、C、D的坐标分别为 假设声源N的空间坐标为(x,y,z),水平偏角为仰角为θ,r为声源到四元传声器阵列的阵元中心的距离。声源N传播到传声器A、B、C、D的时间分别为t1、t2、t3、t4,r1为声源到传声器A的距离。设定3组相对时延值:t21=t2-t1,t31=t3-t1,t41=t4-t1,其对应的声程差为:d21=ct21、d31=ct31、d41=ct41。
声源N到原点、传声器A、B、C、D的距离分别用式(1)至(5)表示:
x2+y2+z2=r2 (1)
式(3)、(4)、(5)分别减去式(2),可得出:
xd-yd=2r1d21+d21 2 (6)
2xd=2r1d31+d31 2 (7)
xd+yd=2r1d41+d41 2 (8)
根据式(7)可得出:
式(8)减去式(6)可得出:
式(6)加上式(8)再减去式(7)可得出:
将直角坐标系下声源N的位置坐标(x,y,z)用球面坐标系下的位置坐标表示,可得出:
z=rcosθ (14)
式(13)除以式(12)可得出:
由于ri>>di1,d≈di1,i=1,2,3,所以:
式(1)减去式(2)并将式(9)代入,可得出:
由式(17)可得出:
根据图1,可得出:
将式(9)、(10)、(11)代入式(19),可得出:
用d21=ct21,d31=ct31,d41=ct41对式(16)、(18)、(20)进行代换,可得出:
如图2所示的传统直三棱锥式阵列模型,声源传播速度为c,假设声源M的空间位置为(x',y',z'),S为声源M到原点(0,0,0,)的距离。直三棱锥式阵列中的四个传声器A’、B’、C’、D’的坐标分别为(0,0,0)、(d',0,0)、(0,d',0)、(0,0,d'),声源M传播到传声器A’、B’、C’、D’的时间分别为t0'、t1'、t2'、t3',设定3组相对时延值:Δt1=t1'-t0',Δt2=t2'-t0',Δt3=t3'-t0'。将声源M与X轴正半轴的夹角定义为水平偏角将声源M与X-Y-O平面的夹角定义为仰角θ'。
令S1=cΔt1,S2=cΔt2,S3=cΔt3,表示出声源M到传声器A、B、C、D距离的几何等式,可得出:
x'2+y'2+z'2=S2 (24)
(x'-d')2+y'2+z'2=(S+S1)2 (25)
x'2+(y'-d')2+z'2=(S+S2)2 (26)
x'2+y'2+(z'-d')2=(S+S3)2 (27)
式(24)减去式(25),可得出:
其中,xx'为声源坐标正演值参数x'对应的仿真值。
式(24)减去式(26),可得出:
其中,yy'为声源坐标正演值参数y'对应的仿真值。
式(24)减去式(27),可得出:
其中,zz'为声源坐标正演值参数z'对应的仿真值。
将直角坐标系下声源M的位置坐标(x',y',z')用球面坐标系下的位置坐标表示有:
zz'=S·cosθ' (33)
式(32)除以式(31)可得出:
其中,为水平偏角正演值对应的仿真值。
由于Si≈d,i=1,2,3,可以近似得出:
其中,θ”为仰角正演值θ'对应的仿真值。
此外,下面通过仿真比较初始声源位置与仿真所得位置,可以分析误差来评价本发明方法对声源定位的效果,如表1和表2所示。仿真结果表明:本发明的新型算法定位精度较高,误差较低,分别比传统的直三棱锥式声源方位定位算法的水平偏角和仰角误差降低了大约10%、1%,具有较好的效果。
表1新型的基于四元传声器阵列的声源方位定位算法仿真结果表
表2传统的基于直三棱锥式阵列的声源方位定位算法仿真结果表
本发明首先基于传声器阵列进行建模,以模拟空间声源的存在;继而利用声源方位估计算法结合模型进行公式推导;通过将反演结果与起始声源数据进行对比,最终进行误差分析。仿真数据表明,与传统的基于直三棱锥式声源方位定位算法相比,其误差大大降低,能有效的解决传统的基于直三棱锥式的声源方位定位算法定位精度较低的问题,提高了声源方位定位的精确性与时效性。
以上所述,仅为本发明中的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉该技术的人在本发明所揭露的技术范围内,可理解想到的变换或替换,都应涵盖在本发明的包含范围之内,因此,本发明的保护范围应该以权利要求书的保护范围为准。
Claims (5)
1.一种新型的声源方位定位方法,其特征在于,该方法通过构建四元传声器阵列,反演声源位置;所述四元传声器阵列包括第一至第四传声器,以四元传声器阵列的阵元中心为坐标原点,设各传声器到坐标原点的距离均为则第一至第四传声器的坐标分别为 反演得到声源的坐标(x,y,z):z=rcosθ,为水平偏角,θ为仰角,其中,t21为声源传播到第二、第一传声器的相对时延值,t31为声源传播到第三、第一传声器的相对时延值,t41为声源传播到第四、第一传声器的相对时延值,c为声源传播速度,。
2.根据权利要求1所述的一种新型的声源方位定位方法,其特征在于,声源到阵元中心的距离
3.根据权利要求1所述的一种新型的声源方位定位方法,其特征在于,t21=t2-t1,t1、t2分别为声源传播到第一、第二传声器的时间。
4.根据权利要求1所述的一种新型的声源方位定位方法,其特征在于,t41=t4-t1,t1、t4分别为声源传播到第一、第四传声器的时间。
5.根据权利要求1所述的一种新型的声源方位定位方法,其特征在于,t31=t3-t1,t1、t3分别为声源传播到第一、第三传声器的时间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810330568.8A CN108680901A (zh) | 2018-04-13 | 2018-04-13 | 一种新型的声源方位定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810330568.8A CN108680901A (zh) | 2018-04-13 | 2018-04-13 | 一种新型的声源方位定位方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108680901A true CN108680901A (zh) | 2018-10-19 |
Family
ID=63799448
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810330568.8A Pending CN108680901A (zh) | 2018-04-13 | 2018-04-13 | 一种新型的声源方位定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108680901A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110716179A (zh) * | 2019-09-29 | 2020-01-21 | 浙江海洋大学 | 一种基于鸣声的鸟类定位***及方法 |
CN111060872A (zh) * | 2020-03-17 | 2020-04-24 | 深圳市友杰智新科技有限公司 | 基于麦克风阵列的声源定位方法、装置和计算机设备 |
CN111323746A (zh) * | 2020-03-19 | 2020-06-23 | 哈尔滨工程大学 | 一种双圆阵的方位-等效时延差被动定位方法 |
CN111460362A (zh) * | 2020-03-30 | 2020-07-28 | 南京信息工程大学 | 一种基于四元传声器阵列群的声源定位数据互补方法 |
WO2021017256A1 (zh) * | 2019-07-31 | 2021-02-04 | 苏州协昌环保科技股份有限公司 | 声波感知智能电磁脉冲阀以及电磁脉冲阀异常检测方法 |
CN116299181A (zh) * | 2023-03-17 | 2023-06-23 | 成都理工大学 | 一种声源三维空间定位*** |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106679651A (zh) * | 2017-02-08 | 2017-05-17 | 北京地平线信息技术有限公司 | 声源定位方法、装置和电子设备 |
-
2018
- 2018-04-13 CN CN201810330568.8A patent/CN108680901A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106679651A (zh) * | 2017-02-08 | 2017-05-17 | 北京地平线信息技术有限公司 | 声源定位方法、装置和电子设备 |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2021017256A1 (zh) * | 2019-07-31 | 2021-02-04 | 苏州协昌环保科技股份有限公司 | 声波感知智能电磁脉冲阀以及电磁脉冲阀异常检测方法 |
CN110716179A (zh) * | 2019-09-29 | 2020-01-21 | 浙江海洋大学 | 一种基于鸣声的鸟类定位***及方法 |
CN111060872A (zh) * | 2020-03-17 | 2020-04-24 | 深圳市友杰智新科技有限公司 | 基于麦克风阵列的声源定位方法、装置和计算机设备 |
CN111060872B (zh) * | 2020-03-17 | 2020-06-23 | 深圳市友杰智新科技有限公司 | 基于麦克风阵列的声源定位方法、装置和计算机设备 |
CN111323746A (zh) * | 2020-03-19 | 2020-06-23 | 哈尔滨工程大学 | 一种双圆阵的方位-等效时延差被动定位方法 |
CN111323746B (zh) * | 2020-03-19 | 2023-05-05 | 哈尔滨工程大学 | 一种双圆阵的方位-等效时延差被动定位方法 |
CN111460362A (zh) * | 2020-03-30 | 2020-07-28 | 南京信息工程大学 | 一种基于四元传声器阵列群的声源定位数据互补方法 |
CN116299181A (zh) * | 2023-03-17 | 2023-06-23 | 成都理工大学 | 一种声源三维空间定位*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108680901A (zh) | 一种新型的声源方位定位方法 | |
CN104237849B (zh) | 一种双五元十字阵被动声定位融合方法 | |
CN104360315A (zh) | 基于LabVIEW的麦克风阵列声源定位方法及装置 | |
CN104181506B (zh) | 一种基于改进phat加权时延估计的声源定位方法及其实现*** | |
CN103926579B (zh) | 一种开放式阵列信号的仿真方法 | |
CN106646645B (zh) | 一种重力正演加速方法 | |
CN106842128A (zh) | 运动目标的声学跟踪方法及装置 | |
CN105589100A (zh) | 一种微地震震源位置和速度模型同时反演方法 | |
CN103792513B (zh) | 一种雷声定位***及方法 | |
CN104035065A (zh) | 基于主动旋转的声源定向装置及其应用方法 | |
CN109635494A (zh) | 一种飞行试验与地面仿真气动力数据综合建模方法 | |
CN104330768B (zh) | 一种基于声矢量传感器的机动声源方位估计方法 | |
CN103994765A (zh) | 一种惯性传感器的定位方法 | |
CN102609940A (zh) | 利用地面激光扫描技术进行测量对象表面重建时点云配准误差处理方法 | |
CN105300382B (zh) | 一种基于惯性测量单元的大角度船体变形测量方法 | |
CN105066996B (zh) | 自适应矩阵卡尔曼滤波姿态估计方法 | |
CN104181513B (zh) | 一种雷达天线阵元位置的校正方法 | |
CN105353351A (zh) | 一种基于多信标到达时间差改进型定位方法 | |
CN110058299A (zh) | 地震定位方法、装置及终端设备 | |
CN115952691B (zh) | 多站无源时差交叉联合定位***的优化布站方法及装置 | |
CN109100685A (zh) | 一种双面四元十字阵的被动声定位融合算法 | |
Zhu et al. | A calibration method of USBL installation error based on attitude determination | |
CN111401180B (zh) | 神经网络识别模型训练方法、装置、服务器及存储介质 | |
CN104360384B (zh) | 基于纵横波能量自动扫描的微地震事件定位方法及装置 | |
CN106339081A (zh) | 一种基于商用设备的免携带设备手掌定位人机交互方法 |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181019 |
|
RJ01 | Rejection of invention patent application after publication |