CN107247251B - 基于压缩感知的三维声源定位方法 - Google Patents

基于压缩感知的三维声源定位方法 Download PDF

Info

Publication number
CN107247251B
CN107247251B CN201710493279.5A CN201710493279A CN107247251B CN 107247251 B CN107247251 B CN 107247251B CN 201710493279 A CN201710493279 A CN 201710493279A CN 107247251 B CN107247251 B CN 107247251B
Authority
CN
China
Prior art keywords
sound source
microphone array
dimensional
matrix
expression model
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.)
Expired - Fee Related
Application number
CN201710493279.5A
Other languages
English (en)
Other versions
CN107247251A (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.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
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 Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN201710493279.5A priority Critical patent/CN107247251B/zh
Publication of CN107247251A publication Critical patent/CN107247251A/zh
Application granted granted Critical
Publication of CN107247251B publication Critical patent/CN107247251B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/20—Position of source determined by a plurality of spaced direction-finders

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Circuit For Audible Band Transducer (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

本发明公开了一种基于压缩感知的三维声源定位方法，用于解决现有三维声源定位方法抗噪性能差的技术问题。技术方案是通过麦克风阵列获得声源信号的测量值，对选定的三维声源区域进行均匀的网格划分，将每一个网格节点作为潜在的声源位置。进而根据自由场格林函数的Helmholtz方程建立网格节点与麦克风阵列之间的测量矩阵，获得麦克风阵列测量值与未知声源信号之间的三维窄带声源定位稀疏表达模型。通过对稀疏表达模型中的麦克风阵列测量值进行奇异值分解，获得变形后的声源定位稀疏表达模型。最后采用压缩感知OMP算法对变形后的表达模型进行迭代求解，获得声源区域各网格节点的声源强度，对声源进行定位。提高了声源定位的抗噪性能。

Description

基于压缩感知的三维声源定位方法
技术领域
本发明涉及一种三维声源定位方法，特别涉及一种基于压缩感知的三维声源定位方法。
背景技术
基于压缩感知的声源定位方法，能够在不满足奈奎斯特采样定理的情况下对声源进行定位，而且能够获得超分辨率的定位结果。
文献“Ning F,Wei J,Qiu L,et al.Three-dimensional acoustic imaging withplanar microphone arrays and compressive sensing[J].Journal of Sound andVibration,2016,380:112-128.”公开了一种基于压缩感知正交匹配追踪(orthogonalmatching pursuit，OMP)算法的三维声源定位方法。该方法能够基于麦克风阵列测量技术与压缩感知OMP算法对三维区域内的声源进行定位。首先对所选三维声源区域进行网格划分，建立网格节点与麦克风阵列之间的测量矩阵，获得声源与麦克风阵列之间的三维窄带声源定位稀疏表达模型。进而通过OMP算法迭代求解各网格节点的声源强度，达到声源定位的目的。虽然文献所述方法能够获得超分辨率的结果，但是当信噪比(Signal to NoiseRatio,SNR)低至5dB时，该方法不能对声源进行定位，即抗噪性能差。
综上所述，文献中采用压缩感知OMP算法的三维声源定位方法抗噪性能差。
发明内容
为了克服现有三维声源定位方法抗噪性能差的不足，本发明提供一种基于压缩感知的三维声源定位方法。该方法通过麦克风阵列获得声源信号的测量值，经过信号处理得到所需窄带频率点的测量数据。对选定的三维声源区域进行均匀的网格划分，将每一个网格节点作为潜在的声源位置。进而根据自由场格林函数的Helmholtz方程建立网格节点与麦克风阵列之间的测量矩阵，获得麦克风阵列测量值与未知声源信号之间的三维窄带声源定位稀疏表达模型。通过对稀疏表达模型中的麦克风阵列测量值进行奇异值分解，获得变形后的声源定位稀疏表达模型。最后采用压缩感知OMP算法对变形后的表达模型进行迭代求解，获得声源区域各网格节点的声源强度，对声源进行定位。由于采用了压缩感知OMP算法，并通过对奇异值分解的使用，提高了声源定位方法的抗噪性能。
本发明解决其技术问题所采用的技术方案：一种基于压缩感知的三维声源定位方法，其特点是包括以下步骤：
步骤一、通过由M个麦克风传感器组成的麦克风阵列对三维空间的声源信号进行采集，获得时域的麦克风阵列测量数据。
步骤二、对所采集的麦克风阵列测量数据进行分块、加窗以及快速傅立叶变换处理，获得麦克风阵列在频域的测量数据。
步骤三、在每个麦克风传感器所对应的数据块内选择所需窄带频率点的测量数据y，维度为M×B，B为划分的数据块总数。
步骤四、建立三维窄带声源定位的稀疏表达模型，包括三维声源区域网格的划分和测量矩阵的建立。
对选定声源的三维空间区域进行网格划分。以每个网格节点作为潜在声源，构建未知声源信号x，维度为M×N，N为声源区域划分网格节点的数目。未知声源信号x由网格节点处的声源强度组成。
根据自由场格林函数的Helmholtz方程建立网格节点与麦克风阵列之间的测量矩阵A，维度为M×N。公式如下：
Figure BDA0001331603790000021
式中，Cmn为归一化因子，j为虚数单位，f为窄带频率，c为声速，dmn为第m个麦克风传感器与第n个网格节点间的距离。
三维窄带声源定位的稀疏表达模型即为：
y＝Ax+e， (2)
式中，e为麦克风阵列测量数据中包含的噪声项。
步骤五、对稀疏表达模型中的麦克风阵列测量值y进行奇异值分解处理，获得变形的三维声源定位稀疏表达模型。
对稀疏表达模型中的麦克风阵列频域测量值y进行奇异值分解，获得测量值y，未知声源信号x和噪声项e的信号子空间，
ySV＝yVDK，xSV＝xVDK，eSV＝eVDK， (3)
DK＝[IK0]T， (4)
式中，V为对测量值y奇异值分解后的酉矩阵，IK为K×K阶单位矩阵，0为K×(B-K)阶零矩阵，K为声源信号的稀疏度。
为便于后续压缩感知OMP算法的使用，将已建立的表达模型变形为
Figure BDA0001331603790000031
式中，
Figure BDA0001331603790000032
和
Figure BDA0001331603790000033
分别由ySV，xSV和eSV中每一列首尾相接堆叠而来。由N个矩阵A为对角线元素组成新的矩阵
Figure BDA0001331603790000034
步骤六、通过压缩感知OMP算法对表达模型进行迭代求解，获得声源区域内各网格节点处的声源强度。OMP算法的求解过程如下：
(a)初始化残差
Figure BDA0001331603790000035
索引集合
Figure BDA0001331603790000036
迭代计数值t＝1。
(b)通过解决简单的优化问题寻找
Figure BDA0001331603790000037
中对应列的索引λt，
Figure BDA0001331603790000038
(c)增加索引集合和所选原子的矩阵
Figure BDA0001331603790000039
φ0为一个空矩阵。
(d)使用最小二乘获得新的信号估计值
Figure BDA00013316037900000310
(e)计算获得新的近似并更新残差值
Figure BDA00013316037900000311
(f)如果新的残差值rt小于OMP算法迭代的误差限ε，迭代计数值t加1，并返回到步骤(b)继续迭代求解。OMP算法迭代的误差限ε通过求解麦克风阵列测量值中所含噪声的l2范数获得。
最后，将所求解的向量
Figure BDA00013316037900000312
转变为矩阵xSV，对xSV进行列向量的求和平均，获得声源区域内各网格节点处的声源强度。
本发明的有益效果是：该方法通过麦克风阵列获得声源信号的测量值，经过信号处理得到所需窄带频率点的测量数据。对选定的三维声源区域进行均匀的网格划分，将每一个网格节点作为潜在的声源位置。进而根据自由场格林函数的Helmholtz方程建立网格节点与麦克风阵列之间的测量矩阵，获得麦克风阵列测量值与未知声源信号之间的三维窄带声源定位稀疏表达模型。通过对稀疏表达模型中的麦克风阵列测量值进行奇异值分解，获得变形后的声源定位稀疏表达模型。最后采用压缩感知OMP算法对变形后的表达模型进行迭代求解，获得声源区域各网格节点的声源强度，对声源进行定位。由于采用了压缩感知OMP算法，并通过对奇异值分解的使用，提高了声源定位方法的抗噪性能。经测试，当信噪比低至5dB时，本发明方法仍然能够对声源进行定位。
下面结合附图和具体实施方式对本发明作详细说明。
附图说明
图1是本发明基于压缩感知的三维声源定位方法的流程图。
图2是本发明方法中三维声源定位的传播模型。
图3是图1中压缩感知正交匹配追踪算法的流程图。
图4是本发明方法在信噪比为5dB下的定位结果。
图5是背景技术方法在信噪比为5dB下的定位结果。
图中，方块表示仿真中所采用的真实声源，圆点表示通过声源定位方法所获得的重构声源。
具体实施方式
参照图1-5。本实施例采用位于S1(-0.2,0,0.43)，S2(0.2,0,0.43)、S3(-0.2,0,0.74)和S4(-0.2,0,0.74)的四个单极子声源作为仿真中的声源信号。基于压缩感知的三维声源定位方法具体步骤如下：
步骤一、首先采用阵列口径为1m，包含60个麦克风传感器的阵列对仿真信号进行采集，获得每个麦克风的时域测量数据。
步骤二、对所获得的麦克风阵列测量数据进行分块、加窗以及快速傅里叶变换，获得麦克风阵列的频域测量数据。
在该过程中所分数据块的长度均为1024，对每个数据块加汉宁窗处理，并通过快速傅里叶变换得到麦克风阵列频域测量数据。为了探究所提出方法的抗噪性能，仿真中对麦克风阵列测量数据加入了高斯白噪声。选取窄带频率为5000Hz，SNR为5dB。
步骤三、建立三维窄带声源定位稀疏表达模型，包括对三维声源区域的网格划分和测量矩阵的建立。
首先选定0.8m×0.8m×0.68m的三维声源区域，三维区域与阵列平面最小距离为0.4m。对所选声源区域进行网格划分，获得60×60×51个网格节点。然后根据自由场格林函数的Helmholtz方程建立网格节点与麦克风阵列之间的测量矩阵A，其维度为60×183600。具体公式如下：
Figure BDA0001331603790000051
式中，Cmn为归一化因子，j为虚数单位，f为窄带频率，c为声速，dmn为第m个麦克风传感器与第n个网格节点间的距离。
三维窄带声源定位的稀疏表达模型即为
y＝Ax+e， (2)
式中e为麦克风阵列测量数据中包含的噪声项。
步骤四、对稀疏表达模型中的麦克风阵列测量值y进行奇异值分解处理，获得变形的三维声源定位稀疏表达模型。
对稀疏表达模型中的麦克风阵列频域测量值y进行奇异值分解，获得测量值y，未知声源信号x和噪声项e的信号子空间，
ySV＝yVDK，xSV＝xVDK，eSV＝eVDK， (3)
DK＝[IK0]T， (4)
式中，V为对测量值y奇异值分解后的酉矩阵，IK为K×K阶单位矩阵，0为K×(B-K)阶零矩阵，K为声源信号的稀疏度。由于实际情况中，声源信号的稀疏度往往是未知的。研究中发现，稀疏度K的取值对声源定位的结果影响很小，因而本实施例选取K＝1。
为便于后续压缩感知OMP算法的使用，将已建立的表达模型变形为
Figure BDA0001331603790000052
式中
Figure BDA0001331603790000053
和
Figure BDA0001331603790000054
分别由ySV，xSV和eSV中每一列首尾相接堆叠而来。由N个矩阵A为对角线元素组成新的矩阵
Figure BDA0001331603790000055
步骤五、采用压缩感知OMP算法对变形后的表达模型进行迭代求解。OMP算法的迭代终止条件选取误差限，误差限ε通过计算仿真中所加入高斯白噪声的能量而获得。具体步骤如下：
(a)初始化残差
Figure BDA0001331603790000056
索引集合
Figure BDA0001331603790000057
迭代计数值t＝1。
(b)通过解决简单的优化问题寻找
Figure BDA0001331603790000058
中对应列的索引λt，
Figure BDA0001331603790000061
(c)增加索引集合和所选原子的矩阵
Figure BDA0001331603790000062
φ0为一个空矩阵。
(d)使用最小二乘获得新的信号估计值
Figure BDA0001331603790000063
(e)计算获得新的近似并更新残差值
Figure BDA0001331603790000064
(f)如果新的残差值rt小于OMP算法迭代的误差限ε，迭代计数值t加1，并返回到步骤(b)继续迭代求解。OMP算法迭代的误差限ε通过求解麦克风阵列测量值中所含噪声的l2范数获得。
最后，将所求解的向量
Figure BDA0001331603790000065
转变为矩阵xSV，对xSV进行列向量的求和平均，获得声源区域内各网格节点处的声源强度。
从图4和图5可以发现，当信噪比等于5dB时，背景技术方法对声源定位误差过大，本发明方法能够对每个声源进行定位，提高了抗噪性能。

Claims (1)

1.一种基于压缩感知的三维声源定位方法，其特征在于包括以下步骤：
步骤一、通过由M个麦克风传感器组成的麦克风阵列对三维空间的声源信号进行采集，获得时域的麦克风阵列测量数据；
步骤二、对所采集的麦克风阵列测量数据进行分块、加窗以及快速傅立叶变换处理，获得麦克风阵列在频域的测量数据；
步骤三、在每个麦克风传感器所对应的数据块内选择所需窄带频率点的测量值y，维度为M×B，B为划分的数据块总数；
步骤四、建立三维窄带声源定位的稀疏表达模型，包括三维声源区域网格的划分和测量矩阵的建立；
对选定声源的三维空间区域进行网格划分；以每个网格节点作为潜在声源，构建未知声源信号x，维度为M×N，N为声源区域划分网格节点的数目；未知声源信号x由网格节点处的声源强度组成；
根据自由场格林函数的Helmholtz方程建立网格节点与麦克风阵列之间的测量矩阵A，维度为M×N；公式如下：
Figure FDA0002270735250000011
式中，Cmn为归一化因子，j为虚数单位，f为窄带频率，c为声速，dmn为第m个麦克风传感器与第n个网格节点间的距离；
三维窄带声源定位的稀疏表达模型即为：
y＝Ax+e， (2)
式中，e为麦克风阵列测量数据中包含的噪声项；
步骤五、对稀疏表达模型中的麦克风阵列测量值y进行奇异值分解处理，获得变形的三维声源定位稀疏表达模型；
对稀疏表达模型中的麦克风阵列测量值y进行奇异值分解，获得测量值y，未知声源信号x和噪声项e的信号子空间，
ySV＝yVDK，xSV＝xVDK，eSV＝eVDK， (3)
DK＝[IK0]T， (4)
式中，V为对测量值y奇异值分解后的酉矩阵，IK为K×K阶单位矩阵，0为K×(B-K)阶零矩阵，K为声源信号的稀疏度；
为便于后续压缩感知OMP算法的使用，将已建立的表达模型变形为
Figure FDA0002270735250000021
式中，
Figure FDA0002270735250000022
和
Figure FDA0002270735250000023
分别由ySV，xSV和eSV中每一列首尾相接堆叠而来；由N个矩阵A为对角线元素组成新的矩阵
Figure FDA0002270735250000024
步骤六、通过压缩感知OMP算法对表达模型进行迭代求解，获得声源区域内各网格节点处的声源强度；OMP算法的求解过程如下：
(a)初始化残差
Figure FDA0002270735250000025
索引集合
Figure FDA0002270735250000026
迭代计数值t＝1；
(b)通过解决简单的优化问题寻找
Figure FDA0002270735250000027
中对应列的索引λt，
Figure FDA0002270735250000028
(c)增加索引集合和所选原子的矩阵
Figure FDA0002270735250000029
φ0为一个空矩阵；
(d)使用最小二乘获得新的信号估计值
Figure FDA00022707352500000210
(e)计算获得新的近似并更新残差值
Figure FDA00022707352500000211
(f)如果新的残差值rt小于OMP算法迭代的误差限ε，迭代计数值t加1，并返回到步骤(b)继续迭代求解；OMP算法迭代的误差限ε通过求解麦克风阵列测量值中所含噪声的l2范数获得；
最后，将所求解的向量
Figure FDA00022707352500000212
转变为矩阵xSV，对xSV进行列向量的求和平均，获得声源区域内各网格节点处的声源强度。
CN201710493279.5A 2017-06-20 2017-06-20 基于压缩感知的三维声源定位方法 Expired - Fee Related CN107247251B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710493279.5A CN107247251B (zh) 2017-06-20 2017-06-20 基于压缩感知的三维声源定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710493279.5A CN107247251B (zh) 2017-06-20 2017-06-20 基于压缩感知的三维声源定位方法

Publications (2)

Publication Number Publication Date
CN107247251A CN107247251A (zh) 2017-10-13
CN107247251B true CN107247251B (zh) 2020-04-03

Family

ID=60018395

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710493279.5A Expired - Fee Related CN107247251B (zh) 2017-06-20 2017-06-20 基于压缩感知的三维声源定位方法

Country Status (1)

Country Link
CN (1) CN107247251B (zh)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110133595B (zh) * 2018-02-09 2023-05-23 北京搜狗科技发展有限公司 一种声源测向方法、装置和用于声源测向的装置
CN108828525A (zh) * 2018-06-25 2018-11-16 大连大学 混响房间内分布式麦克风阵列多声源定位系统
CN109001680A (zh) * 2018-06-25 2018-12-14 大连大学 声源定位中的块稀疏优化算法
CN108828526A (zh) * 2018-06-25 2018-11-16 大连大学 多声源定位中二维网格中的信号s在空-频域的表达具有块稀疏性的方法
CN109001681A (zh) * 2018-06-25 2018-12-14 大连大学 多声源定位中构造压缩观测矩阵的方法
CN109932689A (zh) * 2019-02-24 2019-06-25 华东交通大学 一种适用于特定定位场景的任意阵列优化方法
CN110221249A (zh) * 2019-05-16 2019-09-10 西北工业大学 基于压缩感知的宽带声源定位方法
CN111664932A (zh) * 2020-05-22 2020-09-15 重庆大学 一种基于贝叶斯压缩感知的声源识别方法
CN111640188A (zh) * 2020-05-29 2020-09-08 中国地质大学（武汉） 基于Mumford-Shah算法框架的抗噪声三维网格优化方法
CN112147571B (zh) * 2020-08-31 2023-10-31 南京理工大学 基于正则正交匹配追踪和蝙蝠算法的声源方位角估计方法
CN112255590B (zh) * 2020-10-26 2024-03-01 中国电子科技集团公司第三研究所 基于模糊函数匹配的低空声源反演定位方法及装置
CN112526451B (zh) * 2020-10-29 2024-04-19 清华大学 基于麦克风阵列成像的压缩波束形成及系统
CN113721194B (zh) * 2021-07-30 2024-04-09 南京师范大学 一种基于mwcs的近场语音信号三维定位方法
CN113655440B (zh) * 2021-08-09 2023-05-30 西南科技大学 一种自适应折中预白化的声源定位方法
CN115113139B (zh) * 2022-05-12 2024-02-02 苏州清听声学科技有限公司 基于传声器阵列的声源识别方法、装置及电子设备

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103400349A (zh) * 2013-07-20 2013-11-20 西安电子科技大学 基于盲压缩感知的图像重构方法
CN105068048A (zh) * 2015-08-14 2015-11-18 南京信息工程大学 基于空间稀疏性的分布式麦克风阵列声源定位方法
CN105490682A (zh) * 2015-11-27 2016-04-13 安庆师范学院 基于svd和omp的压缩传感信号恢复算法
CN105974358A (zh) * 2016-05-25 2016-09-28 天津商业大学 基于压缩感知的智能天线doa估计方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2943980C (en) * 2014-05-27 2018-11-27 Halliburton Energy Services, Inc. Acoustic deblurring for downwell sensors

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103400349A (zh) * 2013-07-20 2013-11-20 西安电子科技大学 基于盲压缩感知的图像重构方法
CN105068048A (zh) * 2015-08-14 2015-11-18 南京信息工程大学 基于空间稀疏性的分布式麦克风阵列声源定位方法
CN105490682A (zh) * 2015-11-27 2016-04-13 安庆师范学院 基于svd和omp的压缩传感信号恢复算法
CN105974358A (zh) * 2016-05-25 2016-09-28 天津商业大学 基于压缩感知的智能天线doa估计方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
A robust super-resolution approach with sparsity constraint in acoustic imaging;Ning Chu,et al;《Applied Acoustics 》;20130907;p197-208 *
three-dimensional acoustic imaging with planar microphone arrays and compressive sensing;Fangli Ning,et al;《Journal of Sound and Vibration 》;20160624;p112-128 *
压缩感知声源定位方法研究;宁方立等;《机械工程学报》;20161031;第42-52页 *
基于压缩感知的信号参数估计;刘姗姗;《中国优秀硕士学位论文全文数据库 信息科技辑》;20170215;全文 *

Also Published As

Publication number Publication date
CN107247251A (zh) 2017-10-13

Similar Documents

Publication Publication Date Title
CN107247251B (zh) 基于压缩感知的三维声源定位方法
JP5496192B2 (ja) 音響場を再構成するための方法
CN109375171B (zh) 一种基于正交匹配追踪算法的声源定位方法
CN109655799B (zh) 基于iaa的协方差矩阵向量化的非均匀稀疏阵列测向方法
CN108802683B (zh) 一种基于稀疏贝叶斯学习的源定位方法
EP3025130B1 (en) Wide-band acoustic holography
CN107576931B (zh) 一种基于协方差低维度迭代稀疏重构的相关/相干信号波达方向估计方法
CN107544051A (zh) 嵌套阵列基于k‑r子空间的波达方向估计方法
CN107493106B (zh) 一种基于压缩感知的频率和角度联合估计的方法
CN109884592B (zh) 一种面向低频高斯噪声源的声源定位仿真方法
CN110687528B (zh) 自适应波束形成器生成方法及系统
CN111812581B (zh) 基于原子范数的球面阵列声源波达方向估计方法
CN110927669A (zh) 一种用于无线声传感器网络的cs多声源定位方法及系统
CN109343003B (zh) 一种快速迭代收缩波束形成声源识别方法
CN112147571A (zh) 基于正则正交匹配追踪和蝙蝠算法的声源方位角估计方法
CN111664932A (zh) 一种基于贝叶斯压缩感知的声源识别方法
CN112285647A (zh) 一种基于稀疏表示与重构的信号方位高分辨估计方法
CN110907892B (zh) 一种球麦克风阵列语音信号到达角估计方法
CN113381793B (zh) 一种面向相干信源估计的无网格波达方向估计方法
Nascimento et al. Acoustic imaging using the Kronecker array transform
CN113640891B (zh) 一种基于奇异谱分析的瞬变电磁探测数据噪声滤除方法
Bahr et al. Wavespace-based coherent deconvolution
CN114648041A (zh) 一种基于平行稀疏阵列的二维欠定doa估计算法
CN115267673A (zh) 考虑重建网格偏移的稀疏声源成像方法、系统
CN112558009A (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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200403