CN109633554A - 基于概率数据关联的移动声源到达时延估计方法 - Google Patents
基于概率数据关联的移动声源到达时延估计方法 Download PDFInfo
- Publication number
- CN109633554A CN109633554A CN201910056039.8A CN201910056039A CN109633554A CN 109633554 A CN109633554 A CN 109633554A CN 201910056039 A CN201910056039 A CN 201910056039A CN 109633554 A CN109633554 A CN 109633554A
- Authority
- CN
- China
- Prior art keywords
- tdoa
- time delay
- moment
- time
- data association
- 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)
- Circuit For Audible Band Transducer (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于概率数据关联的移动声源到达时延估计方法。预先在前一时刻利用卡尔曼滤波预测的声源声音信号到达麦克风的时延,广义互相关算法计算本时刻数据的互功率谱密度分布函数,取出互功率谱密度分布函数上最大的m个峰值对应的m个到达麦克风的时延,利用概率关联算法结合预测的TDOA计算本时刻估计出的各个时延可能是TDOA的概率,将这m个时延利用上述概率加权得到该时刻最终TDOA。考虑到移动声源的时延估计具有连贯性,本方法将PDA跟踪算法运用在TDOA连续估计中,把每个时刻估计出的多个时延看做目标跟踪中的多个目标量测,从而将多时刻连续TDOA估计问题转化为跟踪问题,大大增加了移动声源到达时延估计的准确性,对声阵列处理领域具有较高的价值。
Description
技术领域
本发明涉及阵列信号领域,特别是涉及麦克风阵列信号处理和到达时延估计。
背景技术
本发明的背景基于实际需要产生。近年来,在无人机迅速成为研究热点时,也带来一系列问题,比如无人机黑飞,严重影响区域安全。因此无人机防御正在成为各国政府和军方重点关注的新领域。无人机的声音具有明显的特征,可以有效侦别出空中飞行的物体。通过架设多个麦克风阵列,能有效地收集到无人机的信号,但是实际环境存在的噪声、信号在电路设备传输过程中的噪声以及一些其他的混响、多径干扰等,收集到的声音信号非常嘈杂,导致对无人机进行定位时,计算得到的信号到达麦克风的时延(TDOA)会出现比较大的误差,为了提高TDOA估计的准确性,需要对TDOA估计算法进行改进。
在目前的研究内容中,常用的改进方式是结合前一时刻的TDOA的信息,将本时刻的TDOA估计函数的取值范围限定在一个区间内,从而减小TDOA估计的波动,但这种方法对初值要求较高,且对无人机的飞行速度有要求。因此,急需一种新的方法来真正有效的估计出TDOA,同时能够保证较快的运行速度和***较低的成本。
发明内容
为了实现对无人机移动时声音信号时延的跟踪,本发明采用麦克风阵列传感器对空中无人机声音信号进行处理,可以有效估计移动中的无人机声音信号到达麦克风时延。
本发明解决其技术问题所采用的技术方案是:一种基于概率数据关联的移动声源到达时延估计方法,包括如下步骤:
(1)根据广义互相关函数计算麦克风阵列t时刻采集的两组时域声信号x1(t)、x2(t)互功率谱密度分布函数
其中为x1(t)和x2(t)的傅立叶变换结果的乘积,为x1(t)和x2(t)的频域滤波器。
(2)提取出函数中每一时刻最大的m个峰值对应的m个时延。
(3)将每一个时刻的m个时延看做到达麦克风时延TDOA的m个量测,并引入TDOA的动态方程:
Xk+1=AXk+vk,
zk+1,i=CXk+wk,i,
其中Xk是TDOA在k时刻的状态,vk是均值为0,协方差矩阵为Q的高斯白噪声。zk+1,i是在k+1时刻的第i个时延的量测,wk,i是均值为0,协方差矩阵为R的高斯白噪声,且
C=[1 0],
(4)利用卡尔曼滤波的预测部分对TDOA状态进行预测。
(5)利用概率数据关联算法计算各个时延是TDOA的概率,然后引入卡尔曼滤波的更新部分对TDOA状态进行更新。
进一步地,步骤(1)中,采用相位变换(PHAT)加权函数
进一步地,步骤(4)中,引入卡尔曼滤波的预测部分对TDOA状态进行预测:
Xk+1|k=AXk|k,
Pk+1|k=APk|kAT+Q,
其中Xk|k代表TDOA的状态估计,Pk|k代表TDOA的状态估计误差的协方差矩阵,Xk+1|k代表一步预测TDOA状态估计,Pk+1|k代表一步预测TDOA状态估计误差的协方差矩阵,AT代表A的转置。
进一步地,步骤(5)中,利用如下公式计算各个时延是TDOA的概率:
其中βi k代表在k时刻第i个时延是TDOA的概率,β0 k代表在k时刻所有时延都不是TDOA的概率;
bk=λ(2π)1/2|CPk|k-1CT+R-1|1/2,
其中λ是一个常数,表示第i个时延量测的残差,
进一步地,步骤(5)中,引入卡尔曼滤波的更新部分对TDOA状态进行更新:
Kk=Pk|k-1CT(CPk|k-1CT+R-1),
Pk|k=β0 kPk|k-1+(1-β0 k)Pc k|k+Pd k,
其中
Pc k|k=[I-KkC]Pk|k-1,
此时TDOA估计值为CXk|k。
本发明提出的基于概率数据关联的移动声源到达时延估计方法,可以估计移动中的无人机声音信号到达麦克风时延,具有成本低廉,计算迅速等特点。本发明具有以下优势:
(1)在不增加成本的条件下,直接利用同一时刻互功率谱密度分布函数提取的多个时延进行融合,提高了TDOA估计的精度。
(2)能工作在低信噪比受多径效应影响严重的环境下,拓宽了声音阵列的工作条件;
(3)本发明方法简便易行,实时性好,易实现。
附图说明
图1为到达麦克风时延图解;
图2为多时刻的最大互功率谱密度分布函数的峰值对应时延的连线;
图3为多时刻的第二大互功率谱密度分布函数的峰值对应时延的连线;
图4为将多个时延融合后的TDOA连续估计结果;
图5为声源或者无人机定位所采用的四面体阵列模型;
图6为无人机飞行定位结果。
具体实施方式
以下结合附图和具体实施例对本发明作进一步详细说明。
本发明提出的一种基于概率数据关联的移动声源到达时延估计方法,包括如下步骤:
(1)根据广义互相关函数计算麦克风阵列t时刻采集的两组时域声信号x1(t)、x2(t)互功率谱密度分布函数
其中为x1(t)和x2(t)的傅立叶变换结果的乘积,为x1(t)和x2(t)的频域滤波器。
可以采用相位变换(PHAT)加权函数
(2)提取出函数中每一时刻最大的m个峰值对应的m个时延。
(3)将每一个时刻的m个时延看做到达麦克风时延TDOA的m个量测,并引入TDOA的动态方程:
Xk+1=AXk+vk,
zk+1,i=CXk+wk,i,
其中Xk是TDOA在k时刻的状态,vk是均值为0,协方差矩阵为Q的高斯白噪声。zk+1,i是在k+1时刻的第i个时延的量测,wk,i是均值为0,协方差矩阵为R的高斯白噪声,且
C=[1 0],
(4)利用卡尔曼滤波的预测部分对TDOA状态进行预测,具体如下:
Xk+1|k=AXk|k,
Pk+1|k=APk|kAT+Q,
其中Xk|k代表TDOA的状态估计,Pk|k代表TDOA的状态估计误差的协方差矩阵,Xk+1|k代表一步预测TDOA状态估计,Pk+1|k代表一步预测TDOA状态估计误差的协方差矩阵,AT代表A的转置。
(5)利用概率数据关联算法计算各个时延是TDOA的概率,然后引入卡尔曼滤波的更新部分对TDOA状态进行更新。
计算各个时延是TDOA的概率的公式如下:
其中βi k代表在k时刻第i个时延是TDOA的概率,β0 k代表在k时刻所有时延都不是TDOA的概率;
bk=λ(2π)1/2|CPk|k-1CT+R-1|1/2,
其中λ是一个常数,表示第i个时延量测的残差,
引入卡尔曼滤波的更新部分对TDOA状态进行更新的公式如下:
Kk=Pk|k-1CT(CPk|k-1CT+R-1),
Pk|k=β0 kPk|k-1+(1-β0 k)Pc k|k+Pd k,
其中
Pc k|k=[I-KkC]Pk|k-1,
此时TDOA估计值为CXk|k。
附图1是一个麦克风到达时延的一个图解,在实际环境中,麦克风与无人机或者移动声源的距离相对于麦克风阵列阵元之间的间距要大很多,所以声源信号传播到麦克风阵列不同阵元可以看做平行入射。从图中可以明显的看到声源信号到达第一个麦克风比到达第二个麦克风要多传播的距离L’,利用这个距离除以声速就可以得到到达时延(TDOA)。
附图2是连续采集几分钟数据后,将每秒的数据利用广义互相关算法求出互功率谱密度分布函数,然后直接将互功率谱密度分布函数上最大的峰值对应的时延当作声音信号到达麦克风时延(TDOA),再将连续估计的TDOA画成连线图。从图中可以看出,如果直接取最大的峰值对应的时延当作TDOA,得到TDOA估计曲线会有比较大的波动点,这些波动点是由噪声或者多径效应产生的,
附图3是连续采集几分钟数据后,将每秒的数据利用广义互相关算法求出互功率谱密度分布函数,然后直接将互功率谱密度分布函数上第二大的峰值对应的时延当作声音信号到达麦克风时延(TDOA),再将连续估计的TDOA画成连线图。从图中可以看到,直接利用互功率谱密度分布函数第二大峰值对应时延画出来的连线的波动点相较于附图2是不同的,在附图2中是波动点的地方,附图3中却可能是平滑曲线。这说明由于噪声等干扰,TDOA不一定对应互功率谱密度分布函数的最大峰值,可能第二大,第三大,甚至第四大峰值对应的时延才是真实的TDOA。因此,我们可以引入概率数据关联(PDA)算法来计算前几个最大峰值对应时延点是真实TDOA的概率,再利用上述概率将时延加权起来,这样可以提高TDOA估计的准确性。
附图4展示了用本发明提出的时延估计算法进行TDOA估计后得到的TDOA估计连线,可以看到,相较于附图2与附图3,附图4的曲线没有任何的波动,曲线非常平滑,消除了前面由噪声和多径效应产生的干扰点。
附图5是声源或者无人机定位所采用的四面体阵列模型,图中L是阵列原点到麦克风的距离,r1是1号麦克风到达无人机的距离,S1,S2,S3,S4分别是1,2,3,4号麦克风,令(x0,y0,z0)代表无人机的空间坐标,d1i代表声源信号到达1号麦克风与到达i号麦克风的路程差,然后引入下列方程:
x0 2+y0 2+(z0-L)2=r1 2,
x0 2+(y0-L)2+z0 2=(r1+b12)2,
通过上述方程,四个方程四个未知量(x0,y0,z0,r1),可以定位无人机。
附图6展示了我们利用上述阵列进行定位后的结果,并与无人机飞行的GPS轨迹进行了对比,可以看到,在大部分区域都能较好的跟踪无人机。
上述实施例用来解释说明本发明,而不是对本发明进行限制,在本发明的精神和权利要求的保护范围内,对本发明作出的任何修改和改变,都落入本发明的保护范围。
Claims (5)
1.一种基于概率数据关联的移动声源到达时延估计方法,其特征在于,包括如下步骤:
(1)根据广义互相关函数计算麦克风阵列t时刻采集的两组时域声信号x1(t)、x2(t)互功率谱密度分布函数
其中为x1(t)和x2(t)的傅立叶变换结果的乘积,为x1(t)和x2(t)的频域滤波器。
(2)提取出函数中每一时刻最大的m个峰值对应的m个时延。
(3)将每一个时刻的m个时延看做到达麦克风时延TDOA的m个量测,并引入TDOA的动态方程:
Xk+1=AXk+vk,
zk+1,i=CXk+wk,i,
其中Xk是TDOA在k时刻的状态,vk是均值为0,协方差矩阵为Q的高斯白噪声。zk+1,i是在k+1时刻的第i个时延的量测,wk,i是均值为0,协方差矩阵为R的高斯白噪声,且
C=[1 0],
(4)利用卡尔曼滤波的预测部分对TDOA状态进行预测。
(5)利用概率数据关联算法计算各个时延是TDOA的概率,然后引入卡尔曼滤波的更新部分对TDOA状态进行更新。
2.根据权利要求1所述的一种基于概率数据关联的移动声源到达时延估计方法,其特征在于,所述步骤(1)中,采用相位变换(PHAT)加权函数
3.根据权利要求1所述的一种基于概率数据关联的移动声源到达时延估计方法,其特征在于,所述步骤(4)中,引入卡尔曼滤波的预测部分对TDOA状态进行预测:
Xk+1|k=AXk|k,
Pk+1|k=APk|kAT+Q,
其中Xk|k代表TDOA的状态估计,Pk|k代表TDOA的状态估计误差的协方差矩阵,Xk+1|k代表一步预测TDOA状态估计,Pk+1|k代表一步预测TDOA状态估计误差的协方差矩阵,AT代表A的转置。
4.根据权利要求1所述的一种基于概率数据关联的移动声源到达时延估计方法,其特征在于,所述步骤(5)中,利用如下公式计算各个时延是TDOA的概率:
其中βi k代表在k时刻第i个时延是TDOA的概率,β0 k代表在k时刻所有时延都不是TDOA的概率;
bk=λ(2π)1/2|CPk|k-1CT+R-1|1/2,
其中λ是一个常数,表示第i个时延量测的残差,
5.根据权利要求1所述的一种基于概率数据关联的移动声源到达时延估计方法,其特征在于,所述步骤(5)中,引入卡尔曼滤波的更新部分对TDOA状态进行更新:
Kk=Pk|k-1CT(CPk|k-1CT+R-1),
Pk|k=β0 kPk|k-1+(1-β0 k)Pc k|k+Pd k,
其中
Pc k|k=[I-KkC]Pk|k-1,
此时TDOA估计值为CXk|k。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910056039.8A CN109633554A (zh) | 2019-01-18 | 2019-01-18 | 基于概率数据关联的移动声源到达时延估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910056039.8A CN109633554A (zh) | 2019-01-18 | 2019-01-18 | 基于概率数据关联的移动声源到达时延估计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109633554A true CN109633554A (zh) | 2019-04-16 |
Family
ID=66062923
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910056039.8A Pending CN109633554A (zh) | 2019-01-18 | 2019-01-18 | 基于概率数据关联的移动声源到达时延估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109633554A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111474521A (zh) * | 2020-04-09 | 2020-07-31 | 南京理工大学 | 多径环境中基于麦克风阵列的声源定位方法 |
CN113640744A (zh) * | 2021-08-20 | 2021-11-12 | 歌尔科技有限公司 | 声源定位方法及音频设备 |
CN113708859A (zh) * | 2021-08-26 | 2021-11-26 | 大连工业大学 | 一种超分辨多径数量及时延联合估计方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103781095A (zh) * | 2012-10-23 | 2014-05-07 | 华为技术有限公司 | 一种tdoa测量误差的校正方法、传输点和*** |
US20140162686A1 (en) * | 2012-12-11 | 2014-06-12 | Transpacific Ip Management Group Ltd. | Femto-assisted location estimation in macro-femto heterogeneous networks |
EP3203760A1 (en) * | 2016-02-08 | 2017-08-09 | Thomson Licensing | Method and apparatus for determining the position of a number of loudspeakers in a setup of a surround sound system |
CN107102295A (zh) * | 2017-04-13 | 2017-08-29 | 杭州电子科技大学 | 基于glmb滤波的多传感器tdoa无源定位方法 |
CN107707324A (zh) * | 2017-08-28 | 2018-02-16 | 西安电子科技大学 | 一种基于相位差和最大似然估计的声信号时延估计方法 |
WO2018119910A1 (zh) * | 2016-12-29 | 2018-07-05 | 清华大学 | 一种运动声源的定位方法和定位装置 |
CN108398662A (zh) * | 2018-01-23 | 2018-08-14 | 佛山市顺德区中山大学研究院 | 一种提高空间定位精度的方法 |
-
2019
- 2019-01-18 CN CN201910056039.8A patent/CN109633554A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103781095A (zh) * | 2012-10-23 | 2014-05-07 | 华为技术有限公司 | 一种tdoa测量误差的校正方法、传输点和*** |
US20140162686A1 (en) * | 2012-12-11 | 2014-06-12 | Transpacific Ip Management Group Ltd. | Femto-assisted location estimation in macro-femto heterogeneous networks |
EP3203760A1 (en) * | 2016-02-08 | 2017-08-09 | Thomson Licensing | Method and apparatus for determining the position of a number of loudspeakers in a setup of a surround sound system |
WO2018119910A1 (zh) * | 2016-12-29 | 2018-07-05 | 清华大学 | 一种运动声源的定位方法和定位装置 |
CN107102295A (zh) * | 2017-04-13 | 2017-08-29 | 杭州电子科技大学 | 基于glmb滤波的多传感器tdoa无源定位方法 |
CN107707324A (zh) * | 2017-08-28 | 2018-02-16 | 西安电子科技大学 | 一种基于相位差和最大似然估计的声信号时延估计方法 |
CN108398662A (zh) * | 2018-01-23 | 2018-08-14 | 佛山市顺德区中山大学研究院 | 一种提高空间定位精度的方法 |
Non-Patent Citations (6)
Title |
---|
D6PISTAGE DANS UNE AMBIANCE ENCOMBR等: ""Tracking in a Cluttered Environment With Probabilistic Data Association"", 《AUTOMATICA》 * |
ILYAS POTAMITIS等: ""Tracking of Multiple Moving Speakers With Multiple Microphone Arrays"", 《IEEE TRANSACTIONS ON SPEECH AND AUDIO PROCESSING》 * |
XIANYU CHANG等: ""A Surveillance System for Drone Localization and Tracking Using Acoustic Arrays"", 《2018 IEEE 10TH SENSOR ARRAY AND MULTICHANNEL SIGNAL PROCESSING WORKSHOP (SAM)》 * |
李元实: ""面向被动式传感器网络的多目标定位跟踪研究与***设计"", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
杨露箐: "《多源信息融合理论与应用》", 28 February 2006, 北京邮电大学出版社 * |
潘泉: "《现代目标跟踪与信息融合》", 31 October 2009, 国防工业出版社 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111474521A (zh) * | 2020-04-09 | 2020-07-31 | 南京理工大学 | 多径环境中基于麦克风阵列的声源定位方法 |
CN113640744A (zh) * | 2021-08-20 | 2021-11-12 | 歌尔科技有限公司 | 声源定位方法及音频设备 |
CN113708859A (zh) * | 2021-08-26 | 2021-11-26 | 大连工业大学 | 一种超分辨多径数量及时延联合估计方法 |
CN113708859B (zh) * | 2021-08-26 | 2024-03-22 | 大连工业大学 | 一种超分辨多径数量及时延联合估计方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109188362B (zh) | 一种麦克风阵列声源定位信号处理方法 | |
CN109633554A (zh) | 基于概率数据关联的移动声源到达时延估计方法 | |
CN109541548B (zh) | 一种基于匹配场的空气声呐定位方法 | |
CN111474521B (zh) | 多径环境中基于麦克风阵列的声源定位方法 | |
CN106093864A (zh) | 一种麦克风阵列声源空间实时定位方法 | |
CN104142492A (zh) | 一种srp-phat多源空间定位方法 | |
CN103176166A (zh) | 一种用于水声被动定位的信号到达时延差跟踪算法 | |
CN109901112A (zh) | 基于多通道声获取的声学同时定位与建图方法 | |
CN111398902A (zh) | 一种直角三角形三元水听器阵被动测距测向方法 | |
CN105607042A (zh) | 用麦克风阵列时延估计定位声源的方法 | |
CN105353340B (zh) | 一种双层圆柱阵水下被动目标检测方法 | |
Xu et al. | Long-time coherent integration for radar target detection base on Radon-Fourier transform | |
CN103778288A (zh) | 基于蚁群优化的非均匀阵元噪声条件下近场声源定位方法 | |
CN115951305A (zh) | 一种基于srp-phat空间谱和gcc的声源定位方法 | |
Wang et al. | Robust direct position determination methods in the presence of array model errors | |
CN109655783A (zh) | 一种感应器列阵来波方向的估计方法 | |
Hao et al. | BFGS quasi-Newton location algorithm using TDOAs and GROAs | |
CN112731292B (zh) | 局部imf能量加权的低空飞行目标信号时延估计方法 | |
CN116008913A (zh) | 基于stm32和小型麦克风阵列的无人机检测定位方法 | |
CN112731283B (zh) | 基于多级维纳滤波器的高亚音速飞行目标声学测向方法 | |
CN115524672A (zh) | 结构化干扰与杂波下目标鲁棒智能检测方法 | |
CN109239665B (zh) | 一种基于信号子空间相似度谱和粒子滤波器的多声源连续定位方法和装置 | |
CN110736988B (zh) | 双基地pfa运动目标参数估计和成像方法 | |
CN108333575B (zh) | 基于高斯先验与区间约束的移动声源时延滤波方法 | |
Gustafsson et al. | Direction of arrival estimation in sensor arrays using local series expansion of the received signal |
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: 20190416 |
|
RJ01 | Rejection of invention patent application after publication |