CN106124044A - 实心球声源识别低旁瓣超高分辨率声学图像快速获取方法 - Google Patents
实心球声源识别低旁瓣超高分辨率声学图像快速获取方法 Download PDFInfo
- Publication number
- CN106124044A CN106124044A CN201610477614.8A CN201610477614A CN106124044A CN 106124044 A CN106124044 A CN 106124044A CN 201610477614 A CN201610477614 A CN 201610477614A CN 106124044 A CN106124044 A CN 106124044A
- Authority
- CN
- China
- Prior art keywords
- omega
- sound source
- microphone
- sound
- solid sphere
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 70
- 239000003814 drug Substances 0.000 title abstract 3
- 229940079593 drug Drugs 0.000 title 1
- 239000007787 solid Substances 0.000 claims description 35
- 239000011159 matrix material Substances 0.000 claims description 27
- 238000003384 imaging method Methods 0.000 claims description 20
- 238000001228 spectrum Methods 0.000 claims description 16
- 238000001514 detection method Methods 0.000 claims description 9
- 230000017105 transposition Effects 0.000 claims description 9
- 238000012546 transfer Methods 0.000 claims description 7
- 238000012935 Averaging Methods 0.000 claims description 3
- 239000013598 vector Substances 0.000 description 13
- 230000008569 process Effects 0.000 description 12
- 238000010586 diagram Methods 0.000 description 11
- 230000009977 dual effect Effects 0.000 description 6
- 230000005404 monopole Effects 0.000 description 4
- 238000012800 visualization Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 3
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000012805 post-processing Methods 0.000 description 2
- 230000001902 propagating effect Effects 0.000 description 2
- 238000013139 quantization Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000005352 clarification Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002620 method output Methods 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000008685 targeting Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Circuit For Audible Band Transducer (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开一种实心球声源识别低旁瓣超高分辨率声学图像快速获取方法。在识别三维空间内不相干声源时,本发明为快速获得清晰干净的超高分辨率声源成像图提供了有效途径,同时相对于实心球函数型延迟求和方法能更准确地量化声源,对声源识别结果的准确分析具有重要意义。
Description
技术领域
本发明涉及三维空间声源识别领域
背景技术
凭借旋转对称性好、声场记录全面、衍射作用强、记录信号信噪比高等优势,实心球形传声器阵列广泛流行于三维声源识别领域。球谐函数波束形成(Spherical HarmonicsBeamforming,SHB)是其普遍采用的声信号后处理方法。但是受球谐函数阶次截断、传声器离散采样等因素影响,该方法输出的声源成像图被大量高水平旁瓣污染,同时空间分辨率低,这些缺陷会使声源识别结果的分析具有严重不确定性。目前三维空间已有的反卷积(DAMAS)方法在SHB方法输出结果、阵列点传播函数、声源分布之间建立线性方程组,通过在反复迭代过程中引入正约束来定解该方程组,从而提取真实声源分布,可获得超高空间分辨率,但存在严重耗时问题;实心球的函数型延迟求和方法(FDAS)在实心球延迟求和(DAS)方法的基础上发展而来,能快速衰减DAS方法产生的旁瓣并能准确探测弱源,但该方法仅能轻微改善DAS方法空间分辨率,无法准确分辨彼此靠近的声源,并且实际应用中对声源的量化偏差也较大。
发明内容
本发明的目的是针对实心球阵列三维声源识别低旁瓣、超高分辨率声学图像的快速获取问题提供一种方法,同时克服FDAS方法对声源量化偏差大的不足。
为实现本发明目的而采用的技术方案如下,实心球声源识别低旁瓣超高分辨率声学图像快速获取方法:
1):构建测量系统:
参见图1,一个处于三维空间内的实心球阵列;所述实心球阵列包括一个半径为a的实心球,以及分布在所述实心球表面的Q个传声器;这些传声器的编号为q,q=1,2,…,Q;所述实心球的球心作为原点,实心球阵列所处的三维空间内任意一点的位置坐标用(r,Ω)描述,r表示所描述位置与原点间的距离,表示所描述位置的方向,θ、分别为其仰角、方位角;(a,Ωq)为q号传声器的位置坐标;
2)获取声压
各个传声器接收声压信号p(ka,Ωq);其中,处于所述三维空间内的单极子点声源,辐射声波为球面波;波数k=2πf/c,声波频率为f,声速为c;
3)通过各个传声器接收声压信号P,得到互谱矩阵C:
P为各个传声器接收的复声压信号,
B表示所有声源位置坐标组成的集合,(r0,Ω0)和(r0',Ω0')均表示声源的位置坐标,s(kr0,Ω0)表示声源的强度,t(kr0,Ω0)=[t1(kr0,Ω0),t2(kr0,Ω0),…,tQ(kr0,Ω0)]T,tq(kr0,Ω0)表示声源到q号传声器的声场传递函数,
分别为Ω0、Ωq方向的球谐函数,n、m为球谐函数的阶次,Rn(kr0,ka)为声信号传播径向函数,
其中,被减数项为球面波在自由场中传播的径向函数,减数项为实心球引起的球面波散射径向函数,为虚数单位,jn(ka)为n阶第一类球贝塞尔函数,为n阶第二类球汉克尔函数,j'n(ka)、分别为jn(ka)、的一阶导数。
上标“-”、“H”、“*”和“T”分别表示平均运算、转置共轭运算、共轭运算和转置运算;
4)确定各个频率对应的球谐函数阶次截断长度:
表示将数值向正无穷方向圆整到最近的整数;
5)将互谱矩阵C和球谐函数阶次截断长度N代入FDAS方法输出式:
其中:(rf,Ωf)为聚焦点的位置坐标,v(krf,Ωf)=[v1(krf,Ωf),v2(krf,Ωf),…,vQ(krf,Ωf)]T为聚焦列向量,vq(krf,Ωf)表示q号传声器的聚焦分量,Rn(krf,ka)为聚焦径向函数。|| ||2表示2范数,指数参数ξ推荐取值为16。
6)对(Ⅰ)构造脊探测(RD)对象,
L(Ωf)=|WF(krf,Ωf)|
“||”表示取模运算,对每个聚焦方向,可按下列步骤检测其是否为极大值脊:
步骤1.构造L(Ωf)的二阶偏导数矩阵,即黑塞矩阵:
步骤2.特征值分解矩阵Η,找出绝对值最大的特征值λ及相应的特征向量uλ,其中uλ指示最大表面曲率方向,输出值在uλ方向上若为极大值则λ为负,若为极小值则λ为正。
步骤3.计算标量ρ(Ωf)为:
其中,“sign”为符号函数,为梯度算子,“.”表示求向量内积,χ为常数,用以控制检测空间精度,本发明取为0.2。当ρ(Ωf)=1时,聚焦方向为极大值脊,保留FDAS的输出值;反之,用零取代FDAS在该聚焦方向的输出值来剔除非脊点。
7)对DAS方法结果进行反卷积(DAMAS)求取各声源的平均声压贡献:
定义声源平均声压贡献PAC(kr0,Ω0)为声源在阵列各传声器处产生声压信号自谱的均值,即为互谱矩阵C的迹,并带入DAS输出量表达式,
定义DAS的点传播函数为:
其表示(r0,Ω0)位置处的单位平均声压贡献点声源在聚焦点(rf,Ωf)处的波束形成贡献量。DAS输出量为各声源平均声压贡献与对应点传播函数的乘积的和。其中,
t(kr0,Ω0)=[t1(kr0,Ω0),t2(kr0,Ω0),…,tQ(kr0,Ω0)]T
表示声场传递函数列向量。6)中RD方法获取的脊中包含声源点,且声源点处的输出值相对很大,可默认一定动态范围外的脊和非脊处无声源,输出值为零,仅对动态范围内的脊进行DAMAS运算来获得超高空间分辨率,此时G为参与运算的脊的数目,该数值很小,可使反卷积具有极高的效率。W为G×1维列向量,对应元素为参与运算各脊处的DAS输出量。A为G×G维的点传播函数矩阵,PAC为G×1维列向量,对应元素为参与运算各脊处声源的平均声压贡献,为未知非负数。如下线性方程组成立:
W=APAC
该线性方程组采用高斯-赛德尔迭代方案求解PAC,移除点传播函数的影响,能够显著衰减旁瓣并获得超高空间分辨率。初始化基于第l次迭代计算的结果进行第l+1次迭代的具体步骤为:
步骤1.计算残差re:
其中,D为已完成l+1次迭代计算的声源点组成的集合,E为未完成l+1次迭代计算的脊点组成的集合,二者的并集包含所有参与运算脊点声源。
步骤2.确定
按照该方案依次计算每个运算脊点,即可得出
8)根据迭代所得PAC即能定位声源位置并且能绘制声源成像图实现三维声场可视化。
本发明的方法推导过程如下:
步骤1:获取各传声器接收复声压信号的互谱矩阵
图1为球阵列波束形成坐标系,原点位于阵列中心,三维空间内任意位置可用(r,Ω)描述,r表示所描述位置与原点间的距离,表示所描述位置的方向,θ、分别为其仰角、方位角。符号“●”代表嵌在实心球表面的传声器,记a为球半径,Q为传声器总数,q=1,2,…,Q为传声器索引,(a,Ωq)为q号传声器的位置坐标。假设声源为单极子点声源,辐射声波为球面波。声波频率为f,声速为c,波数k=2πf/c,(r0,Ω0)表示声源的位置坐标,tq(kr0,Ω0)表示声源到q号传声器的声场传递函数,根据散射理论,
其中,分别为Ω0、Ωq方向的球谐函数,n、m为球谐函数的阶次,上标“*”表示共轭运算,Rn(kr0,ka)为声信号传播径向函数:
其中,被减数项为球面波在自由场中传播的径向函数,减数项为实心球引起的球面波散射径向函数,为虚数单位,jn(ka)为n阶第一类球贝塞尔函数,为n阶第二类球汉克尔函数,j'n(ka)、分别为jn(ka)、的一阶导数。
用s(kr0,Ω0)表示声源的强度,B表示所有声源位置坐标组成的集合,则q号传声器接收的声压信号p(ka,Ωq)可写为:
令p=[p(ka,Ω1),p(ka,Ω2),…,p(ka,ΩQ)]T
t(kr0,Ω0)=[t1(kr0,Ω0),t2(kr0,Ω0),…,tQ(kr0,Ω0)]T分别为传声器声压信号列向量和声场传递函数列向量,上标“T”表示转置运算,对应式(3),
各传声器接收声压信号的互谱矩阵C写为:
其中,上标“-”、“H”分别表示平均运算和转置共轭运算,(r0',Ω0')也为声源位置坐标。
步骤2:获取实心球延迟求和(DAS)输出函数
过程201 定义实心球DAS输出函数
DAS输出函数为:
其中,(rf,Ωf)为聚焦点的位置坐标,
v(krf,Ωf)=[v1(krf,Ωf),v2(krf,Ωf),…,vQ(krf,Ωf)]T为聚焦列向量,|| ||2表示2范数。v(krf,Ωf)中,元素vq(krf,Ωf)表示q号传声器的聚焦分量,被定义为:
其中,为Ωf方向的球谐函数,Rn(krf,ka)为聚焦径向函数,其表达式与式(2)类同,只需将式(2)中的r0替换为rf即可,N为球谐函数阶次截断长度,需人为设定。式(6)也可写为:
其中,q'也为传声器索引,Cqq'为q号传声器接收声压信号相对于q'号传声器接收声压信号的互谱。
过程202 变形实心球DAS输出函数
定义平均声压贡献为声源在阵列各传声器处产生声压信号自谱的均值,记为PAC(kr0,Ω0),等于传声器接收声压信号互谱矩阵的迹平均。对于单声源,互谱矩阵C可写为:
C=S(kr0,Ω0)t(kr0,Ω0)tH(kr0,Ω0) (9)
其中,被定义为具有功率意义的声源强度。相应地,
其中,“tr”表示求矩阵的迹。根据式(10),用PAC(kr0,Ω0)表示S(kr0,Ω0)后代入式(9)可得:
将式(11)代入式(6)得:
psf被定义为实心球DAS声源识别方法的点传播函数,表示单位平均声压贡献单极子点声源的响应,根据式(12)该DAS方法的psf表达式为:
对单声源,该DAS方法的输出结果等于声源的平均声压贡献与相应点传播函数的乘积。
过程203 确定球谐函数阶次截断长度
球谐函数阶次截断长度按下式确定;
其中,符号表示将数值向正无穷方向圆整到最近的整数。
步骤3获取实心球函数型延迟求和(FDAS)输出函数
函数波束形成(FB)方法首先特征值分解阵列传声器接收声压信号的互谱矩阵C,得:
C=UΣUH (15)
其中,U=[u1,u2,…,uq,…,uQ]为酉矩阵,Σ=diag(σ1,σ2,…,σq,…,σQ)为对角矩阵,σq(q=1,2,…,Q)为矩阵C的特征值,满足σ1≤σ2≤…≤σQ,uq为σq对应的特征向量。然后,引入指数参数ξ,令:
最后,采用已有声源识别方法后处理并计算所得结果的ξ次方。根据步骤2,构建的FDAS方法的输出函数为:
显然,当ξ=1时,式(17)即为步骤2中DAS方法输出函数的表达式。假设三维空间内仅存在一个单极子点声源,平均声压贡献为PAC(kr0,Ω0),理论上,下列关系成立:
σ1=σ2=…=σQ-1=0 (18)
σQ=tr(C)=QPAC(kr0,Ω0) (19)
联立式(16)~(20)可得:
表明:FDAS方法对单极子点声源的响应等于该声源的平均声压贡献与其点传播函数ξ次方的乘积。
步骤4对函数型延迟求和(FDAS)结果进行脊探测(RD)
若二元函数在某位置的输出值在其最大表面曲率方向上为极值,则该位置被称为脊。特定频率及聚焦距离下,FDAS算法的输出函数WF(krf,Ωf)是关于聚焦方向Ωf=(θf,φf)的二元函数,其在声源方向的输出值远大于在其它聚焦方向的输出值,声源方向为明显的极大值脊,可通过脊检测剔除主瓣内的非脊点,以缩减主瓣宽度、提高空间分辨率。根据定义,脊也是最大表面曲率方向上一阶导数符号改变的转折点。记L(Ωf)=|WF(krf,Ωf)|为脊检测对象,“||”表示取模运算,对每个聚焦方向,可按下列过程检测其是否为极大值脊:
过程401.构造L(Ωf)的二阶偏导数矩阵,即黑塞矩阵:
过程402.特征值分解矩阵Η,找出绝对值最大的特征值λ及相应的特征向量uλ,其中uλ指示最大表面曲率方向,输出值在uλ方向上若为极大值则λ为负,若为极小值则λ为正。
过程403.计算标量ρ(Ωf)为:
其中,“sign”为符号函数,为梯度算子,“.”表示求向量内积,χ为常数,用以控制检测空间精度,本发明取为0.2。当ρ(Ωf)=1时,聚焦方向为极大值脊,保留FDAS的输出值;反之,用零取代FDAS在该聚焦方向的输出值来剔除非脊点。
步骤5对延迟求和(DAS)结果进行反卷积(DAMAS)
记G为聚焦点总数,W为G×1维列向量,对应元素为各聚焦点的延迟求和输出量,由式(12)计算得出,A为G×G维的点传播函数矩阵,由式(13)计算得出,PAC为G×1维列向量,对应元素为各聚焦点处声源的平均声压贡献,为未知非负数,对于不相干声源如下线性方程组成立:
W=APAC (24)
(24)式采用高斯-赛德尔迭代方案求解PAC,移除点传播函数的影响,能够显著衰减旁瓣并获得超高空间分辨率。初始化基于第l次迭代计算的结果进行第l+1次迭代的具体过程为:
过程501.计算残差re:
其中,D为已完成l+1次迭代计算的声源点组成的集合,E为未完成l+1次迭代计算的声源点组成的集合,二者的并集包含所有声源。
过程502.确定
按照该方案依次计算每个聚焦点,即可得出
直接用DAMAS清晰化DAS的成像结果需要考虑所有聚焦点,庞大的聚焦点数会使反卷积承受严重耗时问题。步骤4中对FDAS结果进行脊探测获取的脊中包含声源点,且声源点处的输出值相对很大,可默认一定动态范围外的脊和非脊处无声源,输出值为零,仅对动态范围内的脊进行DAMAS运算来获得超高空间分辨率,此时,G不再为聚焦点总数,而是参与运算的脊的数目,该数值很小,可使反卷积具有极高的效率。
步骤6定位声源位置并且能绘制声源成像图实现三维声场可视化。
根据步骤5求得的PAC即能定位声源并能绘制声源成像图实现三维声场可视化。
附图说明
图1是本发明球阵列波束形成坐标系图;
图2是本发明测量布置示意图;
图3是DAMAS方法3000Hz单声源识别声源成像图;
图4是FDAS方法3000Hz单声源识别声源成像图;
图5是本发明方法3000Hz单声源识别成像图;
图6是DAMAS方法3000Hz不相干双声源识别声源成像图;
图7是FDAS方法3000Hz不相干双声源识别声源成像图;
图8是本发明方法3000Hz不相干双声源识别成像图;
具体实施方式
下面结合实施例对本发明作进一步说明,但不应该理解为本发明上述主题范围仅限于下述实施例。在不脱离本发明上述技术思想的情况下,根据本领域普通技术知识和惯用手段,做出各种替换和变更,均应包括在本发明的保护范围内。
实心球声源识别低旁瓣超高分辨率声学图像快速获取方法:
1):构建测量系统:
参见图1,一个处于三维空间内的实心球阵列;所述实心球阵列包括一个半径为a的实心球,以及分布在所述实心球表面的Q个传声器;这些传声器的编号为q,q=1,2,…,Q;实施例中,Q=36;所述实心球的球心作为原点,实心球阵列所处的三维空间内任意一点的位置坐标用(r,Ω)描述,r表示所描述位置与原点间的距离,表示所描述位置的方向,θ、分别为其仰角、方位角;(a,Ωq)为q号传声器的位置坐标;作为验证,本实施例中,声源为稳态白噪声信号激励的扬声器,采用公司、半径为0.0975m、集成4958型传声器的36通道实心球阵列采样声压信号,扬声器声源到阵列中心的距离均为1m。
参见图2,为本实施例测量布置示意图。
2)获取声压
通过硬件采集系统采集各个传声器接收声压信号p(ka,Ωq);其中,波数k=2πf/c,声波频率为f,声速为c,(a,Ωq)为q号传声器的位置坐标;
3)通过各个传声器接收声压信号,得到互谱矩阵C(可采用B&K公司的pulse软件获得):
其中:B表示所有声源位置坐标组成的集合,(r0,Ω0)表示声源的位置坐标,s(kr0,Ω0)表示声源的强度,
t(kr0,Ω0)=[t1(kr0,Ω0),t2(kr0,Ω0),…,tQ(kr0,Ω0)]T,tq(kr0,Ω0)表示声源到q号传声器的声场传递函数,
分别为Ω0、Ωq方向的球谐函数,n、m为球谐函数的阶次,
上标“-”、“H”、“*”和“T”分别表示平均运算、转置共轭运算、共轭运算和转置运算;
4)确定各个频率对应的球谐函数阶次截断长度:
符号表示将数值向正无穷方向圆整到最近的整数。
5)将互谱矩阵C和球谐函数阶次截断长度N分别代入FDAS和DAS输出式:
和
其中:PAC(kr0,Ω0)为互谱矩阵C的迹,(rf,Ωf)为聚焦点的位置坐标,t(kr0,Ω0)=[t1(kr0,Ω0),t2(kr0,Ω0),…,tQ(kr0,Ω0)]T为传声器声压信号声场传递函数列向量,v(krf,Ωf)=[v1(krf,Ωf),v2(krf,Ωf),…,vQ(krf,Ωf)]T,vq(krf,Ωf)表示q号传声器的聚焦分量,ξ为指数参数,推荐取值为16,|| ||2表示2范数,上标“T”表示转置运算;
6)对WF(krf,Ωf)进行脊探测(RD)并剔除非脊点和30dB动态范围之外的脊点,然后针对剩余脊点对延迟求和结果W(krf,Ωf)进行反卷积(DAMAS)求得PAC。
7)根据PAC绘制声源成像图,实现声场可视化,准确定位声源位置。
用DAMAS方法、FDAS方法和本发明方法分别识别3000Hz单声源,成像动态范围均为30dB,聚焦声源面被设定为与阵列同心的球面,仰角从0°到180°,方位角从0°到360°,间隔均为5°,共37×73个聚焦点,聚焦距离为1m,反卷积的迭代次数均取为5000。该扬声器声源平均声压贡献的实测值为49.1649dB。
参见图3,为DAMAS方法3000Hz单声源识别声源成像图,输出主瓣各聚焦点处输出值的线性叠加为49.1866dB,与实测值相差0.0217dB,具有超高分辨率并且旁瓣很小;
参见图4,为FDAS方法3000Hz单声源识别声源成像图,输出主瓣峰值为48.2983dB,与实测值相差较大为0.8666dB且分辨率较低;
参见图5,为本发明方法3000Hz单声源识别成像图,输出主瓣各聚焦点处输出值的线性叠加为49.1649dB,与实测值相差很小为0.0275dB,同时也具有超高分辨率并且旁瓣衰减能力极好;
参见图6,为DAMAS方法3000Hz不相干双声源识别声源成像图;
参见图7,为FDAS方法3000Hz不相干双声源识别声源成像图;
参见图8,为本发明方法3000Hz不相干双声源识别成像图;
显见,本发明方法对单声源及不相干声源均能快速准确地定位声源,成像图清晰干净。相对于直接对延迟求和结果进行反卷积(DAMAS),本发明方法的计算时间仅为前者用时的千分之四,极大地提高了效率,同时与DAMAS方法具有同等的超高空间分辨及更优的旁瓣衰减能力;相对于函数型延迟求和(DAS)方法,本发明延续了前者的强旁瓣衰减能力同时分辨率获得极大提高,并且能更准确地量化声源。
Claims (1)
1.实心球声源识别低旁瓣超高分辨率声学图像快速获取方法,其特征在于:
1):构建测量系统:
一个处于三维空间内的所述实心球阵列;所述实心球阵列包括一个半径为a的实心球,以及分布在所述实心球表面的Q个所述传声器;这些传声器的编号为q,q=1,2,…,Q;所述实心球的球心作为原点,实心球阵列所处的三维空间内任意一点的位置坐标用(r,Ω)描述,r表示所描述位置与原点间的距离,表示所描述位置的方向,θ、分别为其仰角、方位角;(a,Ωq)为q号传声器的位置坐标;
2)获取传声器声压
各个传声器接收声压信号p(ka,Ωq);其中,波数k=2πf/c,声波频率为f,声速为c,(a,Ωq)为q号传声器的位置坐标;
3)通过各个传声器接收声压信号。得到互谱矩阵C:
其中:B表示所有声源位置坐标组成的集合,(r0,Ω0)表示声源的位置坐标,
s(kr0,Ω0)表示声源的强度,
t(kr0,Ω0)=[t1(kr0,Ω0),t2(kr0,Ω0),…,tQ(kr0,Ω0)]T,tq(kr0,Ω0)表示声源到q号传声器的声场传递函数,
分别为Ω0、Ωq方向的球谐函数,n、m为球谐函数的阶次,Rn(kr0,ka)为声信号传播径向函数;
上标“-”、“H”、“*”和“T”分别表示平均运算、转置共轭运算、共轭运算和转置运算;
4)确定各个频率对应的球谐函数阶次截断长度:
表示将数值向正无穷方向圆整到最近的整数;
5)将互谱矩阵C和球谐函数阶次截断长度N分别代入:
和
其中:PAC(kr0,Ω0)为互谱矩阵C的的迹,(rf,Ωf)为聚焦点的位置坐标,v(krf,Ωf)=[v1(krf,Ωf),v2(krf,Ωf),…,vQ(krf,Ωf)]T,vq(krf,Ωf)表示q号传声器的聚焦分量, 和分别为聚焦点方向和各麦克风方向的球谐函数,ξ为指数参数,||||2表示2范数;
6)对WF(krf,Ωf)进行脊探测并剔除非脊点和设定的动态范围之外的脊点,然后针对剩余脊点对延迟求和结果W(krf,Ωf)进行反卷积求得PAC°
7)根据PAC绘制声源成像图。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610477614.8A CN106124044B (zh) | 2016-06-24 | 2016-06-24 | 实心球声源识别低旁瓣超高分辨率声学图像快速获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610477614.8A CN106124044B (zh) | 2016-06-24 | 2016-06-24 | 实心球声源识别低旁瓣超高分辨率声学图像快速获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106124044A true CN106124044A (zh) | 2016-11-16 |
CN106124044B CN106124044B (zh) | 2019-05-07 |
Family
ID=57266075
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610477614.8A Active CN106124044B (zh) | 2016-06-24 | 2016-06-24 | 实心球声源识别低旁瓣超高分辨率声学图像快速获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106124044B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107389184A (zh) * | 2017-06-29 | 2017-11-24 | 江苏理工学院 | 一种二维空间窗处理声信号误差的方法 |
CN107765221A (zh) * | 2017-09-28 | 2018-03-06 | 合肥工业大学 | 适用于相干和非相干声源的反卷积声源成像算法 |
CN109631756A (zh) * | 2018-12-06 | 2019-04-16 | 重庆大学 | 一种基于混合时频域的旋转声源识别方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090028347A1 (en) * | 2007-05-24 | 2009-01-29 | University Of Maryland | Audio camera using microphone arrays for real time capture of audio images and method for jointly processing the audio images with video images |
CN101910807A (zh) * | 2008-01-18 | 2010-12-08 | 日东纺音响工程株式会社 | 声源识别测定装置、系统及方法 |
CN102901950A (zh) * | 2012-09-20 | 2013-01-30 | 浙江工业大学 | 平面阵列识别声源三维坐标的方法 |
CN103389155A (zh) * | 2013-06-26 | 2013-11-13 | 浙江工业大学 | 声品质客观参量三维空间分布数字图像生成方法 |
CN104360314A (zh) * | 2014-10-16 | 2015-02-18 | 浙江省计量科学研究院 | 一种声源识别定位系统计量校准方法及其校准装置 |
-
2016
- 2016-06-24 CN CN201610477614.8A patent/CN106124044B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090028347A1 (en) * | 2007-05-24 | 2009-01-29 | University Of Maryland | Audio camera using microphone arrays for real time capture of audio images and method for jointly processing the audio images with video images |
CN101910807A (zh) * | 2008-01-18 | 2010-12-08 | 日东纺音响工程株式会社 | 声源识别测定装置、系统及方法 |
CN102901950A (zh) * | 2012-09-20 | 2013-01-30 | 浙江工业大学 | 平面阵列识别声源三维坐标的方法 |
CN103389155A (zh) * | 2013-06-26 | 2013-11-13 | 浙江工业大学 | 声品质客观参量三维空间分布数字图像生成方法 |
CN104360314A (zh) * | 2014-10-16 | 2015-02-18 | 浙江省计量科学研究院 | 一种声源识别定位系统计量校准方法及其校准装置 |
Non-Patent Citations (1)
Title |
---|
YANG YANG ET.AL: "Functional delay and sum beamforming for three-dimensional acoustic source identification with solid spherical arrays", 《JOURNAL OF SOUND AND VIBRATION》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107389184A (zh) * | 2017-06-29 | 2017-11-24 | 江苏理工学院 | 一种二维空间窗处理声信号误差的方法 |
CN107765221A (zh) * | 2017-09-28 | 2018-03-06 | 合肥工业大学 | 适用于相干和非相干声源的反卷积声源成像算法 |
CN109631756A (zh) * | 2018-12-06 | 2019-04-16 | 重庆大学 | 一种基于混合时频域的旋转声源识别方法 |
CN109631756B (zh) * | 2018-12-06 | 2020-07-31 | 重庆大学 | 一种基于混合时频域的旋转声源识别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106124044B (zh) | 2019-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chiariotti et al. | Acoustic beamforming for noise source localization–Reviews, methodology and applications | |
AU2015292238B2 (en) | Planar sensor array | |
Chen et al. | Theory and design of compact hybrid microphone arrays on two-dimensional planes for three-dimensional soundfield analysis | |
TWI556654B (zh) | 用以推衍方向性資訊之裝置與方法和系統 | |
CN106483503B (zh) | 实心球阵列三维声源识别的快速反卷积方法 | |
US9299336B2 (en) | Computationally efficient broadband filter-and-sum array focusing | |
CN106443587B (zh) | 一种高分辨率的快速反卷积声源成像算法 | |
CN106093921B (zh) | 基于稀疏分解理论的声矢量阵宽带测向方法 | |
Chu et al. | Deconvolution for three-dimensional acoustic source identification based on spherical harmonics beamforming | |
CN112180329A (zh) | 一种基于阵元随机均匀分布球阵反卷积波束形成的汽车噪声源声成像方法 | |
Yang et al. | Functional delay and sum beamforming for three-dimensional acoustic source identification with solid spherical arrays | |
Zhang et al. | Transient nearfield acoustic holography based on an interpolated time-domain equivalent source method | |
Ginn et al. | Noise source identification techniques: simple to advanced applications | |
CN107884741A (zh) | 一种多球阵列多宽带声源快速定向方法 | |
CN103575808B (zh) | 基于多角度立体匹配的高实时定量超声检测方法 | |
CN106124044B (zh) | 实心球声源识别低旁瓣超高分辨率声学图像快速获取方法 | |
CN109343003B (zh) | 一种快速迭代收缩波束形成声源识别方法 | |
EP3025130A2 (en) | Wide-band acoustic holography | |
CN107765221A (zh) | 适用于相干和非相干声源的反卷积声源成像算法 | |
CN104655728B (zh) | 一种声学相控阵成像方法 | |
de Santana | Fundamentals of acoustic beamforming | |
Chardon et al. | A blind dereverberation method for narrowband source localization | |
Epain et al. | Super-resolution sound field imaging with sub-space pre-processing | |
CN111812581B (zh) | 基于原子范数的球面阵列声源波达方向估计方法 | |
CN105785320A (zh) | 实心球阵列三维声源识别的函数型延迟求和方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |