基于概率数据关联的移动声源到达时延估计方法技术

本发明专利技术公开了一种基于概率数据关联的移动声源到达时延估计方法。预先在前一时刻利用卡尔曼滤波预测的声源声音信号到达麦克风的时延，广义互相关算法计算本时刻数据的互功率谱密度分布函数，取出互功率谱密度分布函数上最大的m个峰值对应的m个到达麦克风的时延，利用概率关联算法结合预测的TDOA计算本时刻估计出的各个时延可能是TDOA的概率，将这m个时延利用上述概率加权得到该时刻最终TDOA。考虑到移动声源的时延估计具有连贯性，本方法将PDA跟踪算法运用在TDOA连续估计中，把每个时刻估计出的多个时延看做目标跟踪中的多个目标量测，从而将多时刻连续TDOA估计问题转化为跟踪问题，大大增加了移动声源到达时延估计的准确性，对声阵列处理领域具有较高的价值。

【技术实现步骤摘要】
基于概率数据关联的移动声源到达时延估计方法
本专利技术涉及阵列信号领域，特别是涉及麦克风阵列信号处理和到达时延估计。
技术介绍
本专利技术的背景基于实际需要产生。近年来，在无人机迅速成为研究热点时，也带来一系列问题，比如无人机黑飞，严重影响区域安全。因此无人机防御正在成为各国政府和军方重点关注的新领域。无人机的声音具有明显的特征，可以有效侦别出空中飞行的物体。通过架设多个麦克风阵列，能有效地收集到无人机的信号，但是实际环境存在的噪声、信号在电路设备传输过程中的噪声以及一些其他的混响、多径干扰等，收集到的声音信号非常嘈杂，导致对无人机进行定位时，计算得到的信号到达麦克风的时延(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＝[10],(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的概率：其中βik代表在k时刻第i个时延是TDOA的概率，β0k代表在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＝β0kPk|k-1+(1-β0k)Pck|k+Pdk,其中Pck|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＝[10],(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的概率的公式如下：其中βik代表在k时刻第i个时延是TDOA的概率，β0k代表在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＝β0kPk|k-1+(1-β0k)Pck|k+Pdk,其中Pck|k＝[I-KkC]Pk|k-1,此时TDOA估计值为CXk|k。附图1是一个麦克风到达时延的一个图解，在实际环境中，麦克风与无人机或者移动声源的距离相对于麦克风阵列阵元之间的间距要大很多，所以声源信号传播到麦克风阵列不同阵元可以看做平行入射。从图中可以明显的看到声源信号到达第一个麦克风比到达第二个麦克风要多传播的距离L’，利用这个距离除以声速就可以得到到达时延(TDOA)。附图2是连续采集几分钟数据后，将每秒的数据利用广义互相关算法求出互功率谱密度分布函数，然后直接将互功率谱密度分布函数上最大的峰值对应的时延当作声音信号到达麦克风时延(TDOA)，再将连续估计的TDOA画成连线图。从图中可以看出，如果直接取最大的峰值对应的时延当作TDOA，得到TDOA估计曲线会有比较大的波动点，这些波动点是由噪声或者多径效应产生的，附图3是连续采集几分钟数据后，将每秒的数据利用广义互相关算法求出互功率谱密度分布函数，然后本文档来自技高网...

【技术保护点】
1.一种基于概率数据关联的移动声源到达时延估计方法，其特征在于，包括如下步骤：(1)根据广义互相关函数计算麦克风阵列t时刻采集的两组时域声信号x1(t)、x2(t)互功率谱密度分布函数

【技术特征摘要】
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＝[10],(4)利用卡尔曼滤波的预测部分对TDOA状态进行预测。(5)利用概率数据关联算法计算各个时延是TDOA的概率，然后引入卡尔曼滤波的更新部分对TDOA状态进行更新。2.根据权利要求1所述的一种基于概率数据关联的移动声源到达时延估计方法，其特征在于，所述步骤(1)中，采用相位变换(PHAT)加权函数3.根据权利要求1所述的一种基于概率数据关联的移动声源到达时延估计方法...

【专利技术属性】
技术研发人员：陈积明，吴泽先，史治国，吴均峰，常先宇，杨超群，
申请(专利权)人：浙江大学，
类型：发明
国别省市：浙江,33