CN109238444A - 一种采用稀疏测量的声场分离方法 - Google Patents

一种采用稀疏测量的声场分离方法 Download PDF

Info

Publication number
CN109238444A
CN109238444A CN201810916431.0A CN201810916431A CN109238444A CN 109238444 A CN109238444 A CN 109238444A CN 201810916431 A CN201810916431 A CN 201810916431A CN 109238444 A CN109238444 A CN 109238444A
Authority
CN
China
Prior art keywords
sound
source
measuring surface
sound source
sound field
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
CN201810916431.0A
Other languages
English (en)
Other versions
CN109238444B (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.)
Shanghai University of Engineering Science
Original Assignee
Shanghai University of Engineering Science
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 Shanghai University of Engineering Science filed Critical Shanghai University of Engineering Science
Priority to CN201810916431.0A priority Critical patent/CN109238444B/zh
Publication of CN109238444A publication Critical patent/CN109238444A/zh
Application granted granted Critical
Publication of CN109238444B publication Critical patent/CN109238444B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H17/00Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明涉及一种采用稀疏测量的声场分离方法,属于声学技术领域。其特征在于:在多个待分离声场的声源之间设置两个测量面,使用随机稀疏阵列测得测量面上的声压信号,在测量面两侧各布置一个等效声源面,等效声源面上布置有点声源,结合等效源原理与波叠加原理将声压信号与点声源进行拟合,形成拟合声场;根据压缩感知理论通过奇异值分解获取拟合声场声压向量的稀疏基,并建立声场分离模型;将测量到的声压信号和拟合声场声压向量的稀疏基代入声场分离模型,通过稀疏正则化方法求得权重向量的稀疏解,从而分别得出测量面两侧声源产生的声场,实现声场分离。本发明只需较少的测量点即可对声源进行声场分离,有较高的分离精度和空间分辨率。

Description

一种采用稀疏测量的声场分离方法
技术领域
本发明涉及一种采用稀疏测量的声场分离方法,属于声学技术领域。
背景技术
声场分离技术在散射声计算、反射系数测量以及非自由声场中的声场重建等领域都具有广泛的应用。目前声场分离技术主要有:于飞等提出的基于空间二维傅里叶变换的声场分离技术、J.Hald等提出的基于的统计最优声场分离技术、F.Jacobsen等提出的镜像法声场分离技术等。但这些声场分离方法大都需要两个测量面或是单个测量面上的声压和质点振速,且其分辨率受到采样定理限制,测量成本较高。宋玉来等和毛锦等分别提出的单测量面声压测量的声场分离技术在一定程度上减少了测点,但这两种方法都需要预知声源的大概位置,影响了其使用范围。
发明内容
本发明基于压缩感知理论提出一种采用稀疏测量的声场分离方法,该方法用一组稀疏正交基表示叠加的声场,并通过稀疏正则化方法获取权重向量的稀疏解,可以在较少测点的基础上实现声场的分离。
技术方案如下:
一种采用稀疏测量的声场分离方法,在多个待分离声场的声源之间设置两个测量面,使用随机稀疏阵列测得测量面上的声压信号,在测量面两侧各布置一个等效声源面,等效声源面上布置有点声源,结合等效源原理与波叠加原理将声压信号与点声源进行拟合,形成拟合声场;根据压缩感知理论通过奇异值分解获取拟合声场声压向量的稀疏基,并建立声场分离模型;将测量到的声压信号和拟合声场声压向量的稀疏基代入声场分离模型,通过稀疏正则化方法求得权重向量的稀疏解,从而分别得出测量面两侧声源产生的声场,实现声场分离。
进一步地,声压信号与点声源进行拟合是将同在一侧的等效声源面与测量面之间定义一个虚拟边界面,虚拟边界面与等效声源面之间的距离为δ,通过将测量面上的声压信号与等效声源的强度相除得到等效声源面与虚拟边界面之间的传递矩阵。
更进一步地,拟合声场的传递矩阵进行奇异值分解后,声源在虚拟边界面上产生的声压能用传递矩阵的左奇异矩阵的列向量张成的空间表示,传递矩阵的左奇异矩阵即为声源在虚拟边界面上产生的声压向量的一组稀疏基。
进一步地,稀疏测量的方法采用如下步骤:
步骤1:在声源附近设置测量面H1和H2,采集两个测量面H1和H2上的声压信号,测量面H1和H2上的声压可表示为PH1=P11+P21、 PH2=P12+P22;其中P11和P21分别表示声源1和声源2在测量面H1上产生的声压;P12和P22分别表示声源1和声源2在测量面H2上产生的声压。
步骤2:在测量面H1和H2的两侧布置等效声源面Q1和Q2,在测量面 H1与等效声源面Q1之间定义一个虚拟边界面B1、在测量面H2与等效声源面 Q2之间定义一个虚拟边界面B2,将测量面两侧声源产生的声场分别用Q1和 Q2面上的点声源进行拟合,则声源1在虚拟边界面B1上产生的声压PB11=GB11q、声源强度q=(GB11)-1PB11、等效声源面Q1与虚拟边界面B1之间的传递矩阵 GB11=exp(ik|rm-rn|)/(4π|rm-rn|) ;其中i为虚数单位、k为波数、rm和rn分别表示第m个空间场点和第n个等效声源的位置坐标;
步骤3:对传递矩阵GB11进行奇异值分解,则有GB11=U11S11V11 H,令PB11=U11S11V11 Hq=U11w1,w1=S11V11 Hq即权重向量;其中U11和V11分别为传递矩阵GB11的左奇异矩阵和右奇异矩阵,S11为包含奇异值的对角矩阵,w1为权重向量,上标“H”表示共轭转置,此时U11也为声源1在虚拟边界面B1上产生的声压向量的一组稀疏基;用同样方法获得声源1在虚拟边界面B2上的一组稀疏基U12、声源2在虚拟边界面B1上的一组稀疏基U21、声源2在虚拟边界面B2上的一组稀疏基U22;;
步骤4:建立声场分离模型,将声源1在测量面H1上产生的声压P11表示为:P11=GH11q=GH11(GB11)-1PB11=GH11(GB11)-1U11w1,同理得P21=GH21(GB21)-1U21w2、P12=GH12(GB12)-1U12w1、P22=GH22(GB22)-1U22w2,测量面H1和H2上的声压向量则表示为 PH1=Φ11w121w2、PH2=Φ12w122w2;其中Φ11=GH11(GB11)-1U11、Φ12=GH12(GB12)-1U12、Φ21=GH21(GB21)-1U21、Φ22=GH22(GB22)-1U22
步骤5:求解稀疏基,将步骤4中的声压向量PH1、PH2的表达式写成矩阵形式:PH=Φw,通过稀疏正则化的方法求解权重向量w的稀疏解 Jemp=||PH-Φw||2+η||w||1;其中PH=[PH1 PH2]Tw=[w1 w2]T,Jemp表示待最小化的目标函数,η表示正则化参数。
进一步地,步骤2中所述的虚拟边界面B1和B2上离散点的个数相等、等效声源面Q1和Q2上点声源的个数相等。
进一步地,待分离声场的声源为任意形状、稀疏或非稀疏分布。
有益效果:
1)本发明通过稀疏正则化的方法求解拟合声场稀疏基的权重向量,只需要较少的测量点即可实现声场分离,且能保持较高的分离精度和空间分辨率。
2)使用随机稀疏阵列对测量面上的声压信号,测量点个数远少于现有声场分离技术所需的测点数。
3)本发明适用于对任意形状、任意分布的声源进行声场分离。
附图说明
图1为本发明等效声源面、虚拟边界面与测量面位置分布图;
图2为本发明采用的传声器位置分布示意图;
图3为1000Hz时测量面H1上的声压测量值、声源1实部产生的声压理论值以及本发明方法分离的声压值;
图4为1000Hz时测量面H1上的声压测量值、声源1虚部产生的声压理论值以及本发明方法分离的声压值;
图5为在不同的频率下本发明方法与常规等效源分离方法的分离误差。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明:
如图1所示一种采用稀疏测量的声场分离方法,在多个待分离声场的声源之间设置两个测量面,使用随机稀疏阵列测得测量面上的声压信号,在测量面两侧各布置一个等效声源面,等效声源面上布置有点声源,结合等效源原理与波叠加原理将声压信号与点声源进行拟合,形成拟合声场;根据压缩感知理论通过奇异值分解获取拟合声场声压向量的稀疏基,并建立声场分离模型;将测量到的声压信号和拟合声场声压向量的稀疏基代入声场分离模型,通过稀疏正则化方法求得权重向量的稀疏解,从而分别得出测量面两侧声源产生的声场,实现声场分离。
声压信号与点声源进行拟合是将同在一侧的等效声源面与测量面之间定义一个虚拟边界面,虚拟边界面与等效声源面之间的距离为δ,通过将测量面上的声压信号与等效声源的强度相除得到等效声源面与虚拟边界面之间的传递矩阵。
拟合声场的传递矩阵进行奇异值分解后,声源在虚拟边界面上产生的声压能用传递矩阵的左奇异矩阵的列向量张成的空间表示,传递矩阵的左奇异矩阵即为声源在虚拟边界面上产生的声压向量的一组稀疏基。
稀疏测量的方法采用如下步骤:
步骤1:在声源附近设置测量面H1和H2,采集两个测量面H1和H2上的声压信号,测量面H1和H2上的声压可表示为PH1=P11+P21、PH2=P12+P22;其中P11和P21分别表示声源1和声源2在测量面H1上产生的声压;P12和P22分别表示声源1和声源2在测量面H2上产生的声压。
步骤2:在测量面H1和H2的两侧布置等效声源面Q1和Q2,在测量面H1与等效声源面Q1之间定义一个虚拟边界面B1、在测量面H2与等效声源面Q2之间定义一个虚拟边界面B2,令虚拟边界面B1和B2上离散点的个数相等、等效声源面Q1和Q2上点声源的个数相等;将测量面两侧声源产生的声场分别用Q1和 Q2面上的点声源进行拟合;则声源1在虚拟边界面B1上产生的声压 PB11=GB11q,源强q=(GB11)-1PB11,等效声源面Q1与虚拟边界面B1之间的传递矩阵 GB11=exp(ik|rm-rn|)/(4π|rm-rn|);其中i为虚数单位、k为波数、rm和rn分别表示第m个空间场点和第n个等效声源的位置坐标;
步骤3:对传递矩阵GB11进行奇异值分解,则有GB11=U11S11V11 H,令 PB11=U11S11V11 Hq=U11w1,w1=S11V11 Hq即权重向量;其中U11和V11分别为传递矩阵GB11的左奇异矩阵和右奇异矩阵,S11为包含奇异值的对角矩阵,w1为权重向量,上标“H”表示共轭转置,此时U11也为声源1在虚拟边界面B1上产生的声压向量的一组稀疏基;用同样方法获得声源1在虚拟边界面B2上的一组稀疏基U12、声源2在虚拟边界面B1上的一组稀疏基U21、声源2在虚拟边界面B2上的一组稀疏基U22;;
步骤4:建立声场分离模型,将声源1在测量面H1上产生的声压P11表示为:P11=GH11q=GH11(GB11)-1PB11=GH11(GB11)-1U11w1,同理得P21=GH21(GB21)-1U21w2、P12=GH12(GB12)-1U12w1、 P22=GH22(GB22)-1U22w2,测量面H1和H2上的声压向量则表示为 PH1=Φ11w121w2、PH2=Φ12w122w2;其中Φ11=GH11(GB11)-1U11、Φ12=GH12(GB12)-1U12、Φ21=GH21(GB21)-1U21、Φ22=GH22(GB22)-1U22
步骤5:求解稀疏基,将步骤4中的声压向量PH1、PH2的表达式写成矩阵形式:PH=Φw,通过稀疏正则化的方法求解权重向量w的稀疏解Jemp=||PH-Φw||2+η||w||1;其中PH=[PH1 PH2]Tw=[w1 w2]T,Jemp表示待最小化的目标函数,η表示正则化参数。
步骤2中所述的虚拟边界面B1和B2上离散点的个数相等、等效声源面Q1和Q2上点声源的个数相等。
待分离声场的声源为任意形状、稀疏或非稀疏分布。
实施例:采用两个半径为0.01m的脉动球作为声源,其表面振速为1m/s,记位于(0.25,0.25,-0.4)m处的声源为声源1,位于(0.25,0.25,0.4)m处的声源为声源2。测量面位于zH1=-0.025m和zH2=0.025m,在每个测量面上布置36 个测点,测点分布如图2所示,两个等效声源面与虚拟边界面之间的距离均为0.075m。在两个测量面上加入高斯白噪声且信噪比为30dB,通过将声源 1和声源2的声压信号与点声源进行叠加为拟合声场,获取拟合声场声压向量的稀疏基,根据声压信号和稀疏基求解权重向量w的稀疏解,根据权重系数分别得出测量面两侧声源产生的声场,实现声场分离。
图3和图4显示的是频率为1000Hz时测量面H1上的声压测量值、声源1 产生的声压理论值以及本发明方法分离的声压值。可以看出,由于两个声源同时作用在测量面上,因此测量面上的声压值与理论值的声压存在较大差别,而本发明方法分离的声压值和理论值的实部和虚部都吻合地很好。从图5中可以看出,此时的分离误差只有4.2%,证明了本发明方法的有效性。
以上所述仅为本发明的较佳实施例而已,并不用于限制本发明,凡在本发明的原则和精神之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (6)

1.一种采用稀疏测量的声场分离方法,其特征在于:在多个待分离声场的声源之间设置两个测量面,使用随机稀疏阵列测得测量面上的声压信号,在测量面两侧各布置一个等效声源面,等效声源面上布置有点声源,结合等效源原理与波叠加原理将声压信号与点声源进行拟合,形成拟合声场;根据压缩感知理论通过奇异值分解获取拟合声场声压向量的稀疏基,并建立声场分离模型;将测量到的声压信号和拟合声场声压向量的稀疏基代入声场分离模型,通过稀疏正则化方法求得权重向量的稀疏解,从而分别得出测量面两侧声源产生的声场,实现声场分离。
2.如权利要求1所述的声场分离方法,其特征在于:所述声压信号与点声源进行拟合是将同在一侧的等效声源面与测量面之间定义一个虚拟边界面,虚拟边界面与等效声源面之间的距离为δ,通过将测量面上的声压信号与等效声源的强度相除得到等效声源面与虚拟边界面之间的传递矩阵。
3.如权利要求2所述的声场分离方法,其特征在于:所述拟合声场的传递矩阵进行奇异值分解后,声源在虚拟边界面上产生的声压能用传递矩阵的左奇异矩阵的列向量张成的空间表示,传递矩阵的左奇异矩阵即为声源在虚拟边界面上产生的声压向量的一组稀疏基。
4.如权利要求1或2或3所述的声场分离方法,其特征在于:所述的稀疏测量的方法采用如下步骤:
步骤1:在声源附近设置测量面H1和H2,采集两个测量面H1和H2上的声压信号,测量面H1和H2上的声压可表示为PH1=P11+P21、PH2=P12+P22;其中P11和P21分别表示声源1和声源2在测量面H1上产生的声压;P12和P22分别表示声源1和声源2在测量面H2上产生的声压。
步骤2:在测量面H1和H2的两侧布置等效声源面Q1和Q2,在测量面H1与等效声源面Q1之间定义一个虚拟边界面B1、在测量面H2与等效声源面Q2之间定义一个虚拟边界面B2,将测量面两侧声源产生的声场分别用Q1和Q2面上的点声源进行拟合,则声源1在虚拟边界面B1上产生的声压PB11=GB11q、声源强度q=(GB11)-1PB11、等效声源面Q1与虚拟边界面B1之间的传递矩阵
GB11=exp(ik|rm-rn|)/(4π|rm-rn|);
其中i为虚数单位、k为波数、rm和rn分别表示第m个空间场点和第n个等效声源的位置坐标;
步骤3:对传递矩阵GB11进行奇异值分解,则有GB11=U11S11V11 H,令PB11=U11S11V11 Hq=U11w1,w1=S11V11 Hq即权重向量;其中U11和V11分别为传递矩阵GB11的左奇异矩阵和右奇异矩阵,S11为包含奇异值的对角矩阵,w1为权重向量,上标“H”表示共轭转置,此时U11也为声源1在虚拟边界面B1上产生的声压向量的一组稀疏基;用同样方法获得声源1在虚拟边界面B2上的一组稀疏基U12、声源2在虚拟边界面B1上的一组稀疏基U21、声源2在虚拟边界面B2上的一组稀疏基U22;;
步骤4:建立声场分离模型,将声源1在测量面H1上产生的声压P11表示为:P11=GH11q=GH11(GB11)-1PB11=GH11(GB11)-1U11w1,同理得P21=GH21(GB21)-1U21w2、P12=GH12(GB12)-1U12w1、P22=GH22(GB22)-1U22w2,测量面H1和H2上的声压向量则表示为PH1=Φ11w121w2、PH2=Φ12w122w2;其中Φ11=GH11(GB11)-1U11、Φ12=GH12(GB12)-1U12、Φ21=GH21(GB21)-1U21、Φ22=GH22(GB22)- 1U22
步骤5:求解稀疏基,将步骤4中的声压向量PH1、PH2的表达式写成矩阵形式:PH=Φw,通过稀疏正则化的方法求解权重向量w的稀疏解Jemp=||PH-Φw||2+η||w||1;其中PH=[PH1PH2]Tw=[w1 w2]T,Jemp表示待最小化的目标函数,η表示正则化参数。
5.如权利要求4所述的声场分离方法,其特征在于:步骤2中所述的虚拟边界面B1和B2上离散点的个数相等、等效声源面Q1和Q2上点声源的个数相等。
6.如权利要求1所述的声场分离方法,其特征在于:所述的待分离声场的声源为任意形状、稀疏或非稀疏分布。
CN201810916431.0A 2018-08-13 2018-08-13 一种采用稀疏测量的声场分离方法 Active CN109238444B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810916431.0A CN109238444B (zh) 2018-08-13 2018-08-13 一种采用稀疏测量的声场分离方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810916431.0A CN109238444B (zh) 2018-08-13 2018-08-13 一种采用稀疏测量的声场分离方法

Publications (2)

Publication Number Publication Date
CN109238444A true CN109238444A (zh) 2019-01-18
CN109238444B CN109238444B (zh) 2021-02-05

Family

ID=65070153

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810916431.0A Active CN109238444B (zh) 2018-08-13 2018-08-13 一种采用稀疏测量的声场分离方法

Country Status (1)

Country Link
CN (1) CN109238444B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111337793A (zh) * 2020-04-22 2020-06-26 华北电力大学 一种适用于配电网压缩感知故障定位技术的稀疏测点配置方法
CN111537058A (zh) * 2020-04-16 2020-08-14 哈尔滨工程大学 一种基于Helmholtz方程最小二乘法的声场分离方法
CN113063490A (zh) * 2021-03-12 2021-07-02 北京科技大学 一种基于声压和质点振速双面测量的声场分离方法
CN113707172A (zh) * 2021-06-02 2021-11-26 西安电子科技大学 稀疏正交网络的单通道语音分离方法、系统、计算机设备

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0810450B1 (en) * 1996-05-29 2000-05-24 Ge Yokogawa Medical Systems, Ltd. Method and apparatus for ultrasound imaging
CN1952627A (zh) * 2006-10-27 2007-04-25 合肥工业大学 一种声场分离方法
CN101566496A (zh) * 2009-06-05 2009-10-28 合肥工业大学 采用双面振速测量和等效源法分离声场的方法
EP2738762A1 (en) * 2012-11-30 2014-06-04 Aalto-Korkeakoulusäätiö Method for spatial filtering of at least one first sound signal, computer readable storage medium and spatial filtering system based on cross-pattern coherence
CN105490682A (zh) * 2015-11-27 2016-04-13 安庆师范学院 基于svd和omp的压缩传感信号恢复算法
CN106548780A (zh) * 2016-10-28 2017-03-29 南京邮电大学 一种语音信号的压缩感知重构方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0810450B1 (en) * 1996-05-29 2000-05-24 Ge Yokogawa Medical Systems, Ltd. Method and apparatus for ultrasound imaging
CN1952627A (zh) * 2006-10-27 2007-04-25 合肥工业大学 一种声场分离方法
CN101566496A (zh) * 2009-06-05 2009-10-28 合肥工业大学 采用双面振速测量和等效源法分离声场的方法
CN101566496B (zh) * 2009-06-05 2010-09-01 合肥工业大学 采用双面振速测量和等效源法分离声场的方法
EP2738762A1 (en) * 2012-11-30 2014-06-04 Aalto-Korkeakoulusäätiö Method for spatial filtering of at least one first sound signal, computer readable storage medium and spatial filtering system based on cross-pattern coherence
CN105490682A (zh) * 2015-11-27 2016-04-13 安庆师范学院 基于svd和omp的压缩传感信号恢复算法
CN106548780A (zh) * 2016-10-28 2017-03-29 南京邮电大学 一种语音信号的压缩感知重构方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
CHUAN-XING BI: "Sound field reconstruction using compressed modal equivalent point source method", 《THE JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA》 *
GILLES CHARDON: "Near-field acoustic holography using sparse regularization and compressive sampling principles", 《THE JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA》 *
于飞: "双全息面分离声场技术及其在声全息中的应用", 《声学学报》 *
冯振: "CS-MRI中稀疏信号支撑集混合检测方法", 《计算机工程》 *
平国力: "等效源法近场声全息求解算法研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 *
贺岩松: "基于等效源法的双测量面相干声源分离方法", 《仪器仪表学报》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111537058A (zh) * 2020-04-16 2020-08-14 哈尔滨工程大学 一种基于Helmholtz方程最小二乘法的声场分离方法
CN111337793A (zh) * 2020-04-22 2020-06-26 华北电力大学 一种适用于配电网压缩感知故障定位技术的稀疏测点配置方法
CN111337793B (zh) * 2020-04-22 2021-05-25 华北电力大学 一种适用于配电网压缩感知故障定位技术的稀疏测点配置方法
CN113063490A (zh) * 2021-03-12 2021-07-02 北京科技大学 一种基于声压和质点振速双面测量的声场分离方法
CN113707172A (zh) * 2021-06-02 2021-11-26 西安电子科技大学 稀疏正交网络的单通道语音分离方法、系统、计算机设备
CN113707172B (zh) * 2021-06-02 2024-02-09 西安电子科技大学 稀疏正交网络的单通道语音分离方法、系统、计算机设备

Also Published As

Publication number Publication date
CN109238444B (zh) 2021-02-05

Similar Documents

Publication Publication Date Title
CN109238444A (zh) 一种采用稀疏测量的声场分离方法
Barth et al. Acoustic tomographic imaging of temperature and flow fields in air
Dougherty Functional beamforming
CN107478325B (zh) 一种非稳态平面声源的自由场还原方法
CN104237849B (zh) 一种双五元十字阵被动声定位融合方法
Holstein et al. Acoustic tomography on the basis of travel-time measurement
CN104596636B (zh) 声场分离方法
CN107247251A (zh) 基于压缩感知的三维声源定位方法
Turhan et al. On characterization canal surfaces around timelike horizontal biharmonic curves in Lorentzian Heisenberg group Heis3
CN106483503A (zh) 实心球阵列三维声源识别的快速反卷积方法
CN107543601A (zh) 一种复杂外形声源表面瞬态加速度重建方法
CN102944872A (zh) 雷达散射截面近场到近场的变换方法
CN109764956A (zh) 一种基于组合优化正则化方法的近场声全息方法
Zhang et al. Iterative regularization method in generalized inverse beamforming
Langrenne et al. Measurement of confined acoustic sources using near-field acoustic holography
Zhou et al. Three-dimensional identification of flow-induced noise sources with a tunnel-shaped array of MEMS microphones
Shen et al. Periodic boundary based FFT-FISTA for sound source identification
CN109798975A (zh) 一种采用声压和质点加速度测量的非稳态平面声源的自由场实时还原方法
CN106124044B (zh) 实心球声源识别低旁瓣超高分辨率声学图像快速获取方法
Cheng et al. Conical statistical optimal near-field acoustic holography with combined regularization
Fernandez Comesana et al. Measuring musical instruments directivity patterns with scanning techniques
Huthwaite et al. On the convergence of finite element scattering simulations
CN102914354B (zh) 一种三维组合式水听器
Zhang et al. Separation of nonstationary sound fields in the time-wavenumber domain
CN202974420U (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