CN108333575A - 基于高斯先验与区间约束的移动声源时延滤波方法 - Google Patents

基于高斯先验与区间约束的移动声源时延滤波方法 Download PDF

Info

Publication number
CN108333575A
CN108333575A CN201810105261.8A CN201810105261A CN108333575A CN 108333575 A CN108333575 A CN 108333575A CN 201810105261 A CN201810105261 A CN 201810105261A CN 108333575 A CN108333575 A CN 108333575A
Authority
CN
China
Prior art keywords
time delay
time
estimated result
cross
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.)
Granted
Application number
CN201810105261.8A
Other languages
English (en)
Other versions
CN108333575B (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201810105261.8A priority Critical patent/CN108333575B/zh
Publication of CN108333575A publication Critical patent/CN108333575A/zh
Application granted granted Critical
Publication of CN108333575B publication Critical patent/CN108333575B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/539—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section

Landscapes

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

Abstract

本发明公开了一种基于高斯先验与区间约束的移动声源时延滤波方法。该方法利用移动声源的时延估计结果具有连贯性的特点，结合其先验信息以及实际声阵列尺寸的约束条件，优化时延估计的结果。首先利用阵列具体尺寸限定时延的可行区间，利用广义互相关函数计算出互功率谱密度函数，然后利用上一时刻时延估计结果构造高斯先验概率函数，最后更新预测得到当前时刻的时延估计互功率谱密度函数，求解其峰值位置对应的时刻，并根据其时延估计结果是否为零进一步判断得到当前时刻的时延估计结果。该方法对于解决移动声源目标在低信噪比情况下的时延估计具有重要的应用价值，能够极大提高移动声源的时延估计结果。

Description

基于高斯先验与区间约束的移动声源时延滤波方法
技术领域
本发明涉及移动声源目标时延信息估计领域，特别涉及低信噪比且多径效应严重的情况下的移动声源时延信息估计。
背景技术
近年来，无人机市场迅速发展，无人机的数量呈现井喷式的增长；与此同时，无人机“黑飞”、“滥飞”的事件层出不穷，对公共安全和个人隐私带来了严重的影响，因此，无人机的监管与防御成了各国学术界和工业界研究的重点。而声音传感器可以实现对入侵无人机的检测定位，虽然声音检测的距离较短，但是可以通过增加传感的数量和覆盖面积，增加其有效的检测定位距离，因此利用声音传感器对无人机进行检测定位得到了广泛的关注。
无人机作为一个移动的声源目标，尤其是在外界噪声干扰和低信噪比的情况下，给采用声音传感器进行定位带来了巨大的挑战。利用广义互相关函数进行时延估计，然后进行声源目标定位是广泛采用的方法，其中时延估计结果的准确性直接影响定位的精度，因此也是该问题的核心所在。但是，在采用广义互相关函数计算移动声源时延信息时，由于环境中物体的发射，会出现声音传播的多径效应，使得广义互相关函数中出现虚假的峰值，由此引发时延估计的不准确；同时，当声源的信噪比较低时，信号之间的相关性较低，在零值附近出现尖锐的峰值，因而进一步加重了时延估计的误差。
在目前的研究内容中，大多采用波束形成的方法进行信号增强，从而提升有效信号的信噪比，增加时延估计的准确性。但是这样会大大增加声音传感器的数量，同时计算的复杂度也会提升，使得系统的实时性也难以保证。因此，急需一种新的方法来增加时延估计的准确性，同时能够保证较快的运行速度和系统较低的成本。而本发明在不改变声阵列形式和声音传感器数量的基础上，结合时延估计的先验概率和约束条件来增加其准确性，同时算法复杂度较低，具有很好的适用性。
发明内容
本发明的目的在于针对现有技术的不足，提供一种基于高斯先验与区间约束的移动声源时延滤波方法。该方法在不改变声阵列形式的基础上，利用时延估计的先验信息和阵列形式的约束条件，进行时延曲线的滤波估计，以提供更加准确的无人机声源时延估计结果。
本发明的目的是通过以下技术方案来实现的：一种基于高斯先验与区间约束的移动声源时延滤波方法，该方法包括以下步骤：
(1)在声音传感器覆盖的范围内，仅有一个移动声源目标，且其声音信号强度能够被声音传感器接收感知；保证进行互相关函数计算的声音传感器之间时间同步；根据广义互相关函数计算两组时域声信号x1(t)、x2(t)互功率谱密度分布函数Ry1y2(τ)：
其中为x1(t)和x2(t)的傅立叶变换结果，为x1(t)和x2(t)的频域滤波器同时令表示当前时刻的互功率谱密度分布函数，表示上一时刻的互功率谱密度分布函数；
(2)利用得到的上一时刻的时延估计结果tpast计算先验概率密度分布函数GPDF(t)；
(3)计算当前时刻互功率谱密度分布函数其中t满足实际声阵列模型下的时延约束条件；
(4)利用GCC(t)峰值和峰值联合估计出最终的时延结果。
进一步地，根据实际声阵列模型，写出时间约束条件τ∈[-d/c,d/c]，其中τ表示声源到达两个声传感器之间的时延，d表示两个声传感器之间的距离，c表示声速；将τ作为GCC(t)的约束区间，减少时延估计的搜索区间，降低噪点的影响。
进一步地，所述步骤(1)中，采用相位变换(PHAT)加权函数
进一步地，所述步骤(2)中，利用得到的时延估计结果tpast计算高斯分布的先验概率密度分布函数GPDF(t)：
其中μ＝tpast，B表示上一时刻所对应的值左右两侧零点所限定的区间范围；构造的高斯分布的先验概率密度分布函数，其均值为上一时刻的时延估计结果，方差根据上一时刻的互功率谱密度分布函数实时更新。
进一步地，所述步骤(4)中，计算GCC(t)峰值所对应的时刻τnow，如果τnow≠0，那么取τnow作为最终的时延估计结果；如果τnow＝0，那么计算峰值所对应的时刻作为最终的时延估计结果。
本发明提出的基于高斯先验与区间约束的移动声源时延滤波方法，可以更加精确地计算出移动声源的时延结果，从而为移动声源目标定位奠定基础，本发明具有以下优势：
1、不改变声阵列系统的具体形式，不增加系统的成本和复杂度，利用声源移动的特点构造先验概率密度函数，联合预测时延估计结果，大大提高了准确度；
2、当移动声源的信噪比较低，同时环境中存在较为严重的多径效应时，利用该方法能够取得很好的时延估计结果，能够有效抑制时延估计中的噪点；
3、算法简便易行，实时性好，易实现。
附图说明
图1为实验所用声阵列模型图；
图2为时延估计的互功率谱密度函数；
图3为广义互相关法时延估计结果；
图4为多径效应；
图5为利用高斯先验预测当前时刻的互功率谱密度函数；
图6为滤波后时延估计结果。
具体实施方式
以下结合附图对本发明的实施作如下详述：
附图1是实验所用声阵列装置图，该设备部署在浙江大学行政楼的楼顶。声音传感器采用声传科技的CHZ-213，同时配置前置放大器。首先声阵列系统前端由2个四面体阵型的声音阵列组成，每个声阵列安装4个声音传感器，每个声音传感器距离四面体底部中心点的距离为1m；声音信号通过NI-9234四通道数据采集卡进行采集，每个采集卡保证单个四面体阵列上的4个声音传感器信号采集同步，采集频率为25600Hz。
附图2时延估计的互功率谱密度函数，其峰值所对应的横坐标的时刻，即为时延估计的结果。具体的实现方式为：将两个声音传感器接收到的时域信号x1(t)，x2(t)，每0.5秒截取一段数据，对两路数据进行傅里叶变换得到频域结果X1(w)，X2(w)，预滤波器采用相位变换加权(PHAT)，那么具体的互功率谱密度函数可以由下式计算：
其中X2 *(w)表示X2(w)的共轭，然后在时间尺度上对R进行变换，将其横坐标与时间对应即可得到然后搜索峰值所对应的横坐标，即为广义互相关时延估计的结果。
附图3是采用广义互相关法对单个四面体阵列估计的时延估计结果，由于有4个声音传感器，两两之间做互相关计算，可以得到6组时延估计结果，4个声音传感器分别标号为1，2，3，4；tmn表示两个声音传感器之间的额时延估计结果，m＝1,2,3,4，n＝1,2,3,4且m≠n；在本图中，只取其中的三组时延估计结果进行观察，从时延估计曲线可以看出，时延估计存在较大的跳变，时延结果的估计非常不准确，存在很多噪点。
附图4是声音信号传输的多径效应，它是将多个时刻的互相关功率谱密度函数以强度图的形式表现出来，从图中可以看到有多杂乱的轨迹，这是由于声音通过地面或者其他物体反射造成，正因为多径效应的存在，造成了时延估计巨大的误差。
附图5展示了利用高斯先验预测当前时刻的互功率谱密度函数，其具体的实施步骤为：根据上一时刻的时延估计结果以及互功率谱密度分布函数，计算高斯分布先验概率密度分布函数：
其中μ＝tpast，tpast表示上一时刻的时延估计的最终结果，B表示上一时刻互功率谱密度分布函数所对应的值左右两侧零点限定的区间范围；根据计算当前时刻的互功率谱密度分布函数，其中t满足t∈[-d/c,d/c]约束条件，表示当前时刻的互功率谱密度函数；为了防止时延估计结果陷入到零点附近，需要做以下处理：计算互功率谱密度分布函数GCC(t)峰值所对应的时间τnow，如果τnow≠0，那么取τnow作为最终的时延估计结果；如果τnow＝0，那么计算峰值所对应的时间作为最终的时延估计结果。
附图6为滤波后得时延估计结果图，通过采用本文提出的基于高斯先验与区间约束的移动声源时延滤波方法可以有效地去除时延估计曲线中的噪点，大大提高了时延估计的准确度。

Claims (5)

1.一种基于高斯先验与区间约束的移动声源时延滤波方法，其特征在于，包括以下步骤：
(1)根据广义互相关函数计算两组时域声信号x1(t)、x2(t)互功率谱密度分布函数
其中为x1(t)和x2(t)的傅立叶变换结果，为x1(t)和x2(t)的频域滤波器，同时令表示当前时刻的互功率谱密度分布函数，表示上一时刻的互功率谱密度分布函数；
(2)利用得到的上一时刻的时延估计结果tpast计算先验概率密度分布函数GPDF(t)；
(3)计算当前时刻互功率谱密度分布函数其中t满足实际声阵列模型下的时延约束条件；
(4)利用GCC(t)峰值和峰值联合估计出最终的时延结果。
2.根据权利要求1中所述的一种基于高斯先验与区间约束的移动声源时延滤波方法，其特征在于，根据实际声阵列模型，写出时间约束条件τ∈[-d/c,d/c]，其中τ表示声源到达两个声传感器之间的时延，d表示两个声传感器之间的距离，c表示声速；将τ作为GCC(t)的约束区间，减少时延估计的搜索区间，降低噪点的影响。
3.根据权利要求1中所述的一种基于高斯先验与区间约束的移动声源时延滤波方法，其特征在于，所述步骤(1)中，采用相位变换(PHAT)加权函数
4.根据权利要求1中所述的一种基于高斯先验与区间约束的移动声源时延滤波方法，其特征在于，所述步骤(2)中，利用得到的上一时刻的时延估计结果tpast计算高斯分布的先验概率密度分布函数GPDF(t)：
其中μ＝tpast，B表示上一时刻所对应的值左右两侧零点所限定的区间范围；构造的高斯分布的先验概率密度分布函数，其均值为上一时刻的时延估计结果，方差根据上一时刻的互功率谱密度分布函数实时更新。
5.根据权利要求1中所述的一种基于高斯先验与区间约束的移动声源时延滤波方法，其特征在于，所述步骤(4)中，计算GCC(t)峰值所对应的时刻τnow，如果τnow≠0，那么取τnow作为最终的时延估计结果；如果τnow＝0，那么计算峰值所对应的时刻作为最终的时延估计结果。
CN201810105261.8A 2018-02-02 2018-02-02 基于高斯先验与区间约束的移动声源时延滤波方法 Active CN108333575B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810105261.8A CN108333575B (zh) 2018-02-02 2018-02-02 基于高斯先验与区间约束的移动声源时延滤波方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810105261.8A CN108333575B (zh) 2018-02-02 2018-02-02 基于高斯先验与区间约束的移动声源时延滤波方法

Publications (2)

Publication Number Publication Date
CN108333575A true CN108333575A (zh) 2018-07-27
CN108333575B CN108333575B (zh) 2020-10-20

Family

ID=62927711

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810105261.8A Active CN108333575B (zh) 2018-02-02 2018-02-02 基于高斯先验与区间约束的移动声源时延滤波方法

Country Status (1)

Country Link
CN (1) CN108333575B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110007276A (zh) * 2019-04-18 2019-07-12 太原理工大学 一种声源定位方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020048376A1 (en) * 2000-08-24 2002-04-25 Masakazu Ukita Signal processing apparatus and signal processing method
WO2003093868A1 (en) * 2002-05-01 2003-11-13 Quester Tangent Corporation Compensation of sonar image data primarily for seabed classification
CN102438189A (zh) * 2011-08-30 2012-05-02 东南大学 基于双通路声信号的声源定位方法
CN103439688A (zh) * 2013-08-27 2013-12-11 大连理工大学 一种用于分布式麦克风阵列的声源定位系统及定位方法
CN106683122A (zh) * 2016-12-16 2017-05-17 华南理工大学 一种基于高斯混合模型和变分贝叶斯的粒子滤波方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020048376A1 (en) * 2000-08-24 2002-04-25 Masakazu Ukita Signal processing apparatus and signal processing method
WO2003093868A1 (en) * 2002-05-01 2003-11-13 Quester Tangent Corporation Compensation of sonar image data primarily for seabed classification
CN102438189A (zh) * 2011-08-30 2012-05-02 东南大学 基于双通路声信号的声源定位方法
CN103439688A (zh) * 2013-08-27 2013-12-11 大连理工大学 一种用于分布式麦克风阵列的声源定位系统及定位方法
CN106683122A (zh) * 2016-12-16 2017-05-17 华南理工大学 一种基于高斯混合模型和变分贝叶斯的粒子滤波方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Y. OUALIL等: "A multiple hypothesis gaussian mixture filter for acoustic source localization and tracking", 《PROCEEDINGS OF INTERNATIONAL WORKSHOP ON ACOUSTIC SIGNAL ENHANCEMENT》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110007276A (zh) * 2019-04-18 2019-07-12 太原理工大学 一种声源定位方法及系统

Also Published As

Publication number Publication date
CN108333575B (zh) 2020-10-20

Similar Documents

Publication Publication Date Title
CN106772246B (zh) 基于声阵列的无人机实时检测和定位系统及方法
Bimpas et al. Detection of water leaks in supply pipes using continuous wave sensor operating at 2.45 GHz
CN104469154A (zh) 一种基于麦克风阵列的摄像头引导装置及引导方法
CN105589066B (zh) 一种利用垂直矢量阵估计水下匀速运动航行器参数的方法
CN107272005B (zh) 基于可靠声路径下目标回波到达时延和到达角度的主动定位方法
CN109444896B (zh) 水声阵列定位系统及其定位方法
CN105323772B (zh) 一种基于智能手机的传感器网络节点自定位方法
CN109541548A (zh) 一种基于匹配场的空气声呐定位方法
CN108828566A (zh) 基于拖曳线列阵的水下脉冲信号识别方法
CN103235294A (zh) 一种基于外辐射源定位的微弱信号分离估计方法
CN109444864A (zh) 一种深海微弱多目标深度长时累积估计方法
CN108680910A (zh) 基于波形认知的调频广播外辐射源雷达目标检测方法
CN113281727B (zh) 一种基于水平线列阵的输出增强的波束形成方法及其系统
CN107707324A (zh) 一种基于相位差和最大似然估计的声信号时延估计方法
CN104076404A (zh) 运用多通道相干抑制地磁背景噪声的磁异常探测方法
CN108828501A (zh) 在室内声场环境中对移动声源进行实时跟踪定位的方法
CN116609725B (zh) 一种利用深海垂直阵的窄带线谱目标深度估计方法及系统
CN109597021A (zh) 一种波达方向估计方法及装置
CN106443792A (zh) 一种水下多缆声学网络定位方法
CN105388470B (zh) 一种估计目标运动参数的方法
CN105676169A (zh) 一种脉冲声目标的定向方法及装置
CN110488253A (zh) 基于互相关函数延迟求和的多途时延差估计方法
CN108333575A (zh) 基于高斯先验与区间约束的移动声源时延滤波方法
CN103901422A (zh) 一种水下目标回波几何亮点结构特征提取方法
CN103376443B (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
GR01 Patent grant
GR01 Patent grant