CN102494755A - 一种非稳态声场分离方法 - Google Patents
一种非稳态声场分离方法 Download PDFInfo
- Publication number
- CN102494755A CN102494755A CN2011104022519A CN201110402251A CN102494755A CN 102494755 A CN102494755 A CN 102494755A CN 2011104022519 A CN2011104022519 A CN 2011104022519A CN 201110402251 A CN201110402251 A CN 201110402251A CN 102494755 A CN102494755 A CN 102494755A
- Authority
- CN
- China
- Prior art keywords
- measuring surface
- sound source
- acoustic pressure
- time
- infin
- 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
Images
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种非稳态声场分离方法,其特征是对于同时存在目标声源Mo和干扰声源Md,并且均辐射非稳态的声压信号的被测声场,在被测声场中布置一个与测量面S1平行、且相隔距离为Δz的辅助测量面S2,同步测量辅助测量面S2与测量面S1上的声压时域信号;并将声压时域信号分别进行二维空间傅里叶变换获得其声压波数谱;利用声压波数谱和时域脉冲响应函数,并通过解卷运算分离出只有目标声源Mo在测量面S1上所辐射的声压波数谱;对分离出的声压波数谱进行二维空间反傅里叶变换,最终获得测量面S1上只有目标声源Mo所辐射的声压时域信号。本发明方法既可以实现声场在空间域的分离,也可以实现声场在时域的分离,特别适用于非自由声场环境下非稳态声场的分离。
Description
技术领域
本发明涉及物理专业中噪声类领域声场分离方法。
背景技术
在实际工程中,目标声源位于测量面的一侧,而在测量面的另一侧往往存在干扰声源,这些干扰声源所产生的声场影响了对目标声源所辐射声场的准确测量,因此需要采用一定的声场分离方法将干扰声源的影响从测量结果中分离出来。到目前为止,国内外学者已提出多种声场分离方法,这些方法可大致分为五类:一是基于空间傅里叶变换法(SFT)的声场分离技术:G.V.Frisk等在1980年首次提出采用SFT法间接测量海洋底面的反射系数,并建立了基于SFT法的双面声场分离理论;M.Tamura于1990年将G.V.Frisk等提出的方法用于测量斜入射时材料的反射系数;M.T.Cheng等对M.Tamura提出的方法进行了推广,建立了柱面坐标下的双测量面声场分离公式,并用于实现散射声场的分离。二是基于统计最优近场声全息(SONAH)的声场分离技术:J.Hald在对SONAH研究的基础上,提出了基于双全息面声压测量的统计最优声场分离技术;F.Jacobsen等在J.Hald提出的方法的基础上,提出了基于声压和速度测量的统计最优声场分离技术。三是基于球面波叠加的声场分离方法:1956年,J.Pachner采用球面波叠加法实现了任意波场中行波和驻波声场的分离;G.Weinreich等在1980年对J.Pachner提出的方法作了进一步改进,建立了基于双球面测量的声场分离理论。四是基于边界元法(BEM)的声场分离技术:C.Langrenne等在2007年提出一种基于边界元法的双面声场分离方法;随后,E.G.Williams等在2008年提出一种基于逆边界元法和声压速度测量(Cauchy数据)的声场分离方法。五是基于等效源法的声场分离技术:C.X.Bi提出的基于等效源法的声场分离技术,适用于任意形状测量面,且计算稳定性好、计算精度高。但是,上述声场分离方法都仅仅考虑了稳态声场,仅实现了声场在空间域的分离。当目标声源和干扰声源辐射的声场是非稳态时,既需要实现声场在空间域的分离,也需要实现声场在时域的分离,因此上述声场分离方法将不再适用。
发明内容
本发明是为避免上述现有技术所存在的不足之处,提供一种非稳态声场分离方法,既实现声场在空间域的分离,又实现声场在时域的分离。
本发明解决技术问题所采用的技术方案是:
本发明非稳态声场分离方法的特点是对于同时存在目标声源Mo和干扰声源Md,目标声源Mo和干扰声源Md均辐射非稳态的声压信号的被测声场,在所述被测声场中布置一个与测量面S1平行、且相隔距离为Δz的辅助测量面S2,同步测量辅助测量面S2与测量面S1上的声压时域信号;将测量面S1和辅助测量面S2上的声压时域信号分别进行二维空间傅里叶变换获得其声压波数谱;利用所述声压波数谱和时域脉冲响应函数,并通过解卷运算分离出只有目标声源Mo在测量面S1上所辐射的声压波数谱;对分离出的声压波数谱进行二维空间反傅里叶变换,从而获得测量面S1上只有目标声源Mo所辐射的时域信号。
本发明方法的特点按如下步骤进行:
步骤a、同步测量测量面S1和辅助测量面S2上的声压时域信号
在由目标声源Mo和干扰声源Md构成的非稳态被测声场中,位于目标声源Mo和干扰声源Md之间有测量面S1,在测量面S1与干扰声源Md之间设置一个与测量面S1平行、且相隔距离为Δz的辅助测量面S2;在测量面S1与辅助测量面S2上分别均匀分布M个测量网格点,同步测量辅助测量面S2与测量面S1上各个网格点处的声压时域信号;所述Δz的取值大于零,且不大于测量网格点的间隔。
步骤b、对测量面S1与辅助测量面S2上的声压时域信号分别进行二维空间傅里叶变换获得声压波数谱,变换过程为:
在式(11)和式(12)中:
j表示虚数单位;kx、ky分别为x、y方向的波数分量;
z1、z2分别为测量面S1、辅助测量面S2的z方向的空间坐标;
p(x,y,z1,t)为测量面S1上的声压时域信号;
p(x,y,z2,t)为辅助测量面S2上的声压时域信号;
P(kx,ky,z1,t)为测量面S1上的声压波数谱;
P(kx,ky,z2,t)为辅助测量面S2上的声压波数谱;
步骤c、建立测量面S1与辅助测量面S2上的声压波数谱、时域脉冲响应函数、目标声源Mo在测量面S1上所辐射的声压波数谱之间的传递关系
P(kx,ky,z1,t)-P(kx,ky,z2,t)*h(kx,ky,Δz,t)
=Po(kx,ky,z1,t)*[δ(t)-h(kx,ky,Δz,t)*h(kx,ky,Δz,t)] (13)
式(13)中
“*”表示卷积运算;
δ(t)为已知的Dirac delta函数;
h(kx,ky,Δz,t)为已知的时域脉冲响应函数;
Po(kx,ky,z1,t)为目标声源Mo在测量面S1上所辐射的声压波数谱,为待求量;
步骤d、通过解卷运算分离出只有目标声源Mo在测量面S1上所辐射的声压波数谱Po(kx,ky,z1,t),其解卷过程如下:
Po=B-1A (14)
式(14)中
A=[A(kx,ky,t1) A(kx,ky,t2) Λ A(kx,ky,tN)]T (15)
Po=[Po(kx,ky,z1,t1) Po(kx,ky,z1,t2) Λ Po(kx,ky,z1,tN)]T (16)
在式(15)至式(17)中,tn(n=1,2,Λ N)表示各离散时间点,N表示时间离散点数;
在式(15)中,A(kx,ky,tn)表示当t=tn时A(kx,ky,t)的取值,其中
A(kx,ky,t)=P(kx,ky,z1,t)-P(kx,ky,z2,t)*h(kx,ky,Δz,t) (18)
在式(16)中,Po(kx,ky,z1,tn)表示当t=tn时Po(kx,ky,z1,t)的取值;
在式17中,B(kx,ky,tn)表示当t=tn时B(kx,ky,t)的取值,其中
B(kx,ky,t)=δ(t)-h(kx,ky,Δz,t)*h(kx,ky,Δz,t) (19)
步骤e、对由步骤d所分离出的波数谱Po(kx,ky,z1,t)进行二维空间反傅里叶变换,获得测量面S1上只有目标声源Mo所辐射的时域信号,其变换过程如下:
本发明非稳态声场分离方法的特点也在于:
所述各网格点上的声压时域信号是采用双传声器阵列在两测量面上一次快照获得。
所述的测量面S1和辅助测量面S2均为平面。
所述的干扰声源为噪声源、反射源或散射源。
理论模型:
在测量面S1的两侧分别布置着目标声源Mo和干扰声源Md,在测量面S1与干扰声源Md之间设置一个与测量面S1平行、且相隔距离为Δz的辅助测量面S2。
根据声波的叠加原理,在测量面S1上所测得的声压等于目标声源Mo在测量面S1上所辐射的声压与干扰声源Md在测量面S1所辐射的声压之和,即
p(x,y,z1,t)=po(x,y,z1,t)+pd(x,y,z1,t) (21)
式(21)中,下标“o”表示目标声源Mo,下标“d”表示干扰声源Md。
同理,在辅助测量面S2上所测得的声压等于目标声源Mo在辅助测量面S2上所辐射的声压与干扰声源Md在辅助测量面S2所辐射的声压之和,即
p(x,y,z2,t)=po(x,y,z2,t)+pd(x,y,z2,t) (21)
定义声压p(x,y,z,t)关于x、y的二维空间傅里叶变换为
P(kx,ky,z,t)表示声压波数谱。
分别对式(21)和式(22)进行关于x、y的二维空间傅里叶变换可得
P(kx,ky,z1,t)=Po(kx,ky,z1,t)+Pd(kx,ky,z1,t) (24)
P(kx,ky,z2,t)=Po(kx,ky,z2,t)+Pd(kx,ky,z2,t) (25)
由声压在时域波数域的传播方程可知
Po(kx,ky,z2,t)=Po(kx,ky,z1,t)*h(kx,ky,Δz,t) (26)
Pd(kx,ky,z1,t)=Pd(kx,ky,z2,t)*h(kx,ky,Δz,t) (27)
式(26)和(27)中,h(kx,ky,Δz,t)为时域脉冲响应函数,其表达式为
式(28)中,c为声速,J1表示第一类贝塞尔函数,阶数为1,H(t)表示Heaviside函数。
将式(26)代入式(25)中得
P(kx,ky,z2,t)=Po(kx,ky,z1,t)*h(kx,ky,Δz,t)+Pd(kx,ky,z2,t) (29)
将式(27)代入式(24)中得
P(kx,ky,z1,t)=Po(kx,ky,z1,t)+Pd(kx,ky,z2,t)*h(kx,ky,Δz,t) (30)
将式(29)两边分别与时域脉冲响应函数h(kx,ky,Δz,t)进行卷积,然后与式(30)相减可进一步得到
P(kx,ky,z1,t)-P(kx,ky,z2,t)*h(kx,ky,Δz,t)
=Po(kx,ky,z1,t)-Po(kx,ky,z1,t)*h(kx,ky,Δz,t)*h(kx,ky,Δz,t)
=Po(kx,ky,z1,t)*[δ(t)-h(kx,ky,Δz,t)*h(kx,ky,Δz,t)] (31)
式(31)中,P(kx,ky,z1,t)和P(kx,ky,z2,t)可分别通过对测得的声压p(x,y,z1,t)和p(x,y,z2,t)进行二维空间傅里叶变换得到,δ(t)和h(kx,ky,Δz,t)均为已知函数,通过对式(31)进行解卷运算,求得目标声源Mo在测量面S1上所辐射的声压波数谱Po(kx,ky,z1,t)。
解卷过程:
定义
A(kx,ky,t)=P(kx,ky,z1,t)-P(kx,ky,z2,t)*h(kx,ky,Δz,t) (32)
B(kx,ky,t)=δ(t)-h(kx,ky,Δz,t)*h(kx,ky,Δz,t) (33)
则式(31)重写成为:
A(kx,ky,t)=Po(kx,ky,z1,t)*B(kx,ky,t) (34)
在实际计算过程中,需要将时间t离散为t1,t2,Λ tN。根据离散卷积方程,式(34)的离散形式可表示为:
对于每一个波数对(kx,ky),可将式(35)写成矩阵形式
A=BPo (36)
式(36)中,
A=[A(kx,ky,t1) A(kx,ky,t2) Λ A(kx,ky,tN)]T (37)
Po=[Po(kx,ky,z1,t1) Po(kx,ky,z1,t2) Λ Po(kx,ky,z1,tN)]T (38)
对式(36)直接求逆可求得
Po=B-1A (40)
一旦求得目标声源Mo在测量面S1上所辐射的声压波数谱Po(kx,ky,z1,t),再对其进行二维空间反傅里叶变换,即可获得测量面S1上只有目标声源Mo所辐射的时域信号,其变换过程为:
在上述求逆过程中,为保证求解过程的稳定性,采用Tikhonov正则化方法,正则化参数用L曲线法来进行选取。
通过上述方法,既实现了信号在空间域的分离,也实现了信号在时域的分离。
与已有技术相比,本发明的有益效果:
1、本发明考虑了声场在时域内的分离,因而与传统的方法相比,本发明方法特别适用于对非稳态声场进行分离。
2、本发明方法也可以用于对稳态声场进行分析。
附图说明
图1为本发明双测量面非稳态声场分离方法示意图;
图2(a)为本发明方法声源、测量面位置分布图;
图2(b)为本发明方法测量面网格划分示意图:图中测量面S1上的测量点A和测量点C点分别正对着目标声源Mo和干扰声源Md,测量点B位于测量面S1的中心,测量点D靠近测量面S1的边缘;
图3(a)为本发明方法测量点A处声压时域信号图:图中虚线表示测量面S1上的直接测量值(包含目标声源Mo和干扰声源Md所辐射的声压),实线表示目标声源Mo在测量面S1上所辐射的声压理论值,点线表示采用本发明方法分离出的声压值;
图3(b)为本发明方法测量点B处声压时域信号图:图中虚线表示测量面S1上的直接测量值(包含目标声源Mo和干扰声源Md所辐射的声压),实线表示目标声源Mo在测量面S1上所辐射的声压理论值,点线表示采用本发明方法分离出的声压值;
图3(c)为本发明方法测量点C处声压时域信号图:图中虚线表示测量面S1上的直接测量值(包含目标声源Mo和干扰声源Md所辐射的声压),实线表示目标声源Mo在测量面S1上所辐射的声压理论值,点线表示采用本发明方法分离出的声压值;
图3(d)为本发明方法测量点D处声压时域信号图:图中虚线表示测量面S1上的直接测量值(包含目标声源Mo和干扰声源Md所辐射的声压),实线表示目标声源Mo在测量面S1上所辐射的声压理论值,点线表示采用本发明方法分离出的声压值;
图4(a)为本发明方法所获得的时域评价因子T1的数值分布图,图中等高线值为0.99;
图4(b)为本发明方法所获得的时域评价因子T2的数值分布图,图中等高线值为0.05;
图5(a)为本发明方法t1=4.2ms时刻目标声源Mo在测量面S1上所辐射的理论声压场pt;
图5(b)为本发明方法t1=4.2ms时刻直接测量的声压场pm(包含目标声源Mo和干扰声源Md所辐射的声压场);
图5(c)为t1=4.2ms时刻采用本发明方法分离出的声压场pc;
图5(d)为本发明方法t2=5.0ms时刻目标声源Mo在测量面S1上所辐射的理论声压场pt;
图5(e)为本发明方法t2=5.0ms时刻直接测量的声压场pm(包含目标声源Mo和干扰声源Md所辐射的声压场);
图5(f)为t2=5.0ms时刻采用本发明方法分离出的声压场pc;
具体实施方式
参见图1,本实施例中,在测量面S1的两侧分别布置着目标声源Mo和干扰声源Md,干扰声源Md可以为噪声源、反射或散射源。在测量面S1与干扰声源Md之间设置一与测量面S1平行、且相隔距离为Δz的辅助测量面S2。在测量面S1与辅助测量面S2上分别均匀分布有M个测量网格点。Δz的取值大于零,且不大于测量网格点的间隔。
具体实施步骤为:
a、采用双传声器阵列一次快照获得测量面S1与辅助测量面S2上的声压时域信号;
b、对测量面S1和辅助测量面S2上的声压时域信号分别进行二维空间傅里叶变换获得其声压波数谱,其变换过程如下:
在式(42)和式(43)中:
j表示虚数单位;kx、ky分别为x、y方向的波数分量;
z1、z2分别为测量面S1、辅助测量面S2的z方向的空间坐标;
p(x,y,z1,t)为测量面S1上的声压时域信号;
p(x,y,z2,t)为辅助测量面S2上的声压时域信号;
P(kx,ky,z1,t)为测量面S1上的声压波数谱;
P(kx,ky,z2,t)为辅助测量面S2上的声压波数谱;
c、建立测量面S1与辅助测量面S2上的声压波数谱、时域脉冲响应函数、目标声源Mo在测量面S1上所辐射的声压波数谱之间的传递关系
P(kx,ky,z1,t)-P(kx,ky,z2,t)*h(kx,ky,Δz,t)
=Po(kx,ky,z1,t)*[δ(t)-h(kx,ky,Δz,t)*h(kx,ky,Δz,t)] (44)
在式(44)中:
“*”表示卷积运算;
δ(t)为已知的Dirac delta函数;
h(kx,ky,Δz,t)为已知的时域脉冲响应函数;
Po(kx,ky,z1,t)为目标声源Mo在测量面S1上所辐射的声压波数谱,为待求量;
d、通过解卷运算分离出只有目标声源Mo在测量面S1上所辐射的声压波数谱Po(kx,ky,z1,t),其解卷过程如下:
Po=B-1A (45)
式(45)中
A=[A(kx,ky,t1) A(kx,ky,t2) Λ A(kx,ky,tN)]T (46)
Po=[Po(kx,ky,z1,t1) Po(kx,ky,z1,t2) Λ Po(kx,ky,z1,tN)]T (47)
在式(46)至式(48)中,tn(n=1,2,Λ N)表示各离散时间点,N表示时间离散点数;
在式(46)中,A(kx,ky,tn)表示当t=tn时A(kx,ky,t)的取值,其中
A(kx,ky,t)=P(kx,ky,z1,t)-P(kx,ky,z2,t)*h(kx,ky,Δz,t) (49)
在式(47)中,Po(kx,ky,z1,tn)表示当t=tn时Po(kx,ky,z1,t)的取值;
在式(48)中,B(kx,ky,tn)表示当t=tn时B(kx,ky,t)的取值,其中
B(kx,ky,t)=δ(t)-h(kx,ky,Δz,t)*h(kx,ky,Δz,t) (50)
e、对由步骤d所分离出的波数谱Po(kx,ky,z1,t)进行二维空间反傅里叶变换,获得测量面S1上只有目标声源Mo所辐射的时域信号,其变换过程如下:
方法的检验:
在测量面S1两侧分别布置一个单极子,两个单极子均辐射非稳态信号。采用本发明的声场分离方法将测量面S1上目标声源Mo所辐射的声压分离出来,并与其理论声压值进行比较。
本实施例中,作为目标声源Mo的单极子辐射一线性频率调制(频率调制范围为[200,1800]Hz)、高斯幅值调制信号;作为干扰声源Md的单极子辐射一Morlet小波信号,其表达式为
式中,f0=800Hz。
测量面S1与辅助测量面S2的位置关系参见图2(a),测量面S1位于z1=0m的平面上,辅助测量面S2位于z2=0.03m的平面上,因此测量面S1与辅助测量面S2之间的距离Δz为0.03m。测量面S1与辅助测量面S2的大小均为0.648m×0.648m,其上均匀地分布19×19个测量点,参见图2(b)。目标声源Mo位于(0.252,0.324,-0.15)m处,干扰声源Md位于(0.432,0.324,0.18)m处。时域信号采样频率为34400Hz,采样点数为256。
在本实施例的数据处理过程中,为进一步减少截断误差对分离结果的影响,将测量面矩阵维数由19×19通过补零扩展为37×37。
为检验本发明方法在时域内的声场分离效果,在测量面S1上选取了四个测量点,即测量点A、测量点B、测量点C和测量点D,其位置分别为A(0.252,0.324,0)m、B(0.324,0.324,0)m、C(0.432,0.324,0)m、D(0.540,0.324,0)m,其中测量点A和测量点C分别正对着目标声源Mo和干扰声源Md,测量点B位于测量面S1的中心,测量点D靠近测量面S1的边缘。参见图3(其中,图3(a)、3(b)、3(c)、3(d)分别对应测量点A、测量点B、测量点C和测量点D,图中虚线表示测量面S1上的直接测量值(包含目标声源Mo和干扰声源Md所辐射的声压),图中实线表示目标声源Mo在测量面S1上所辐射的声压理论值,图中点线表示采用本发明方法分离出的声压值),比较图中的虚线和实线可以看出,干扰声源Md对目标声源Mo在测量面S1上所辐射的声压造成了较大的干扰;比较图中的实线和点线可以看出,采用本发明方法可以很好地消除干扰声源Md的影响,从而分离出只有目标声源Mo在测量面S1上所辐射的声压值。
为了更加客观地评价本发明方法在时域内的分离效果,在此定义了两个评价因子,它们的表达式分别为
式中,<>表示求平均值,下标“t”表示声压理论值,下标“c”表示声压分离值,上标“rms”表示求均方根值。评价因子T1是用来衡量声压理论值和声压分离值之间的相位误差,当T1的值越靠近1时,相位误差越小。评价因子T2是用来衡量声压理论值和声压分离值之间的幅值误差,当T2的值越靠近0时,幅值误差越小。运用公式(53)、(54)分别计算了测量面S1上各个测量点处的T1和T2。参见图4(a)和图4(b),在大多数测量点处,不论是相位还是幅值,声压理论值和声压分离值都吻合地较好,只是在测量面边缘处吻合程度略低。
为检验本发明方法在空间域的声场分离效果,选取了两个时刻t1=4.2ms和t2=5.0ms。参见图5(图5(a)和5(d)分别为t1和t2时刻目标声源Mo在测量面S1上所辐射的理论声压场pt,图5(b)和5(e)分别为t1和t2时刻直接测量的声压场pm(包含目标声源Mo和干扰声源Md所辐射的声压场),图5(c)和5(f)分别为t1和t2时刻采用本发明方法分离出的声压场pc),比较图5(a)和图5(b)、图5(d)和图5(e)可以看出,干扰声源Md对目标声源Mo在测量面S1上所辐射的理论声压场造成了严重的干扰;比较图5(a)和图5(c)、图5(d)和图5(f)可以看出,采样本发明方法可以很好地消除干扰声源Md的影响,从而还原只有目标声源Mo在测量面S1上所辐射的声压场。
为了更加客观地评价本发明方法在空间域的分离效果,在此定义一个空间误差衡量标准,其表达式为
该标准反映了目标声源Mo在测量面S1上所辐射的理论声压场与采用本发明方法分离出的声压场之间的相对误差。运用公式(55)计算了每个时刻的相对误差。参见图6,在大多数时刻,目标声源Mo在测量面S1上所辐射的理论声压场与采用本发明方法分离出的声压场之间的相对误差都较小;较大的误差出现在信号的两端,这是因为在信号的两端理论声压值较小(即公式(55)的分母较小),从而造成误差变大。
上述例子表明不论是在时域还是在空间域,采用本发明方法都可以很好地将干扰声源Md的影响从测量声压中消除,从而分离出目标声源Mo在测量面S1上所辐射的声压场。
Claims (4)
1.一种非稳态声场分离方法,其特征是对于同时存在目标声源Mo和干扰声源Md,目标声源Mo和干扰声源Md均辐射非稳态的声压信号的被测声场,在所述被测声场中布置一个与测量面S1平行、且相隔距离为Δz的辅助测量面S2,同步测量辅助测量面S2与测量面S1上的声压时域信号;将测量面S1和辅助测量面S2上的声压时域信号分别进行二维空间傅里叶变换获得其声压波数谱;利用所述声压波数谱和时域脉冲响应函数,并通过解卷运算分离出只有目标声源Mo在测量面S1上所辐射的声压波数谱;对分离出的声压波数谱进行二维空间反傅里叶变换,从而获得测量面S1上只有目标声源Mo所辐射的时域信号。
2.根据权利要求1所述的非稳态声场分离方法,其特征是按如下步骤进行:
步骤a、同步测量测量面S1和辅助测量面S2上的声压时域信号
在由目标声源Mo和干扰声源Md构成的非稳态被测声场中,位于目标声源Mo和干扰声源Md之间有测量面S1,在测量面S1与干扰声源Md之间设置一个与测量面S1平行、且相隔距离为Δz的辅助测量面S2;在测量面S1与辅助测量面S2上分别均匀分布M个测量网格点,同步测量辅助测量面S2与测量面S1上各个网格点处的声压时域信号;所述Δz的取值大于零,且不大于测量网格点的间隔;
步骤b、对测量面S1与辅助测量面S2上的声压时域信号分别进行二维空间傅里叶变换获得声压波数谱,变换过程为:
在式(1)和式(2)中:
j表示虚数单位;kx、ky分别为x、y方向的波数分量;
z1、z2分别为测量面S1、辅助测量面S2的z方向的空间坐标;
p(x,y,z1,t)为测量面S1上的声压时域信号;
p(x,y,z2,t)为辅助测量面S2上的声压时域信号;
P(kx,ky,z1,t)为测量面S1上的声压波数谱;
P(kx,ky,z2,t)为辅助测量面S2上的声压波数谱;
步骤c、建立测量面S1与辅助测量面S2上的声压波数谱、时域脉冲响应函数、目标声源Mo在测量面S1上所辐射的声压波数谱之间的传递关系
P(kx,ky,z1,t)-P(kx,ky,z2,t)*h(kx,ky,Δz,t)
=Po(kx,ky,z1,t)*[δ(t)-h(kx,ky,Δz,t)*h(kx,ky,Δz,t)] (3)
式(3)中
“*”表示卷积运算;
δ(t)为已知的Dirac delta函数;
h(kx,ky,Δz,t)为已知的时域脉冲响应函数;
Po(kx,ky,z1,t)为目标声源Mo在测量面S1上所辐射的声压波数谱,为待求量;
步骤d、通过解卷运算分离出只有目标声源Mo在测量面S1上所辐射的声压波数谱Po(kx,ky,z1,t),其解卷过程如下:
Po=B-1A (4)
式(4)中
A=[A(kx,ky,t1) A(kx,ky,t2) Λ A(kx,ky,tN)]T (5)
Po=[Po(kx,ky,z1,t1) Po(kx,ky,z1,t2) Λ Po(kx,ky,z1,tN)]T (6)
在式(5)至式(7)中,tn(n=1,2,Λ N)表示各离散时间点,N表示时间离散点数;在式(5)中,A(kx,ky,tn)表示当t=tn时A(kx,ky,t)的取值,其中
A(kx,ky,t)=P(kx,ky,z1,t)-P(kx,ky,z2,t)*h(kx,ky,Δz,t) (8)
在式(6)中,Po(kx,ky,z1,tn)表示当t=tn时Po(kx,ky,z1,t)的取值;
在式(7)中,B(kx,ky,tn)表示当t=tn时B(kx,ky,t)的取值,其中
B(kx,ky,t)=δ(t)-h(kx,ky,Δz,t)*h(kx,ky,Δz,t) (9)
步骤e、对由步骤d所分离出的波数谱Po(kx,ky,z1,t)进行二维空间反傅里叶变换,获得测量面S1上只有目标声源Mo所辐射的时域信号,其变换过程如下:
2.根据权利要求1所述的非稳态声场分离方法,其特征是所述各网格点上的声压时域信号是采用双传声器阵列在两测量面上一次快照获得。
3.根据权利要求1所述的非稳态声场分离方法,其特征是所述测量面S1和辅助测量面S2均为平面。
4.根据权利要求1所述的非稳态声场分离方法,其特征是所述干扰声源Md为噪声源、反射源或散射源。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110402251 CN102494755B (zh) | 2011-12-06 | 2011-12-06 | 一种非稳态声场分离方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110402251 CN102494755B (zh) | 2011-12-06 | 2011-12-06 | 一种非稳态声场分离方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102494755A true CN102494755A (zh) | 2012-06-13 |
CN102494755B CN102494755B (zh) | 2013-11-06 |
Family
ID=46186599
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110402251 Expired - Fee Related CN102494755B (zh) | 2011-12-06 | 2011-12-06 | 一种非稳态声场分离方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102494755B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103267571A (zh) * | 2013-05-10 | 2013-08-28 | 合肥工业大学 | 采用单面声压和振速测量的实时声场分离方法 |
CN103592022B (zh) * | 2013-11-25 | 2015-08-05 | 合肥工业大学 | 采用声压和质点加速度测量的实时声场分离方法 |
CN105181121A (zh) * | 2015-05-29 | 2015-12-23 | 合肥工业大学 | 采用加权迭代等效源法的高精度近场声全息算法 |
CN105588637A (zh) * | 2016-01-31 | 2016-05-18 | 广西科技大学 | 一种复杂稳定声场声压测试装置 |
CN106289505A (zh) * | 2016-07-21 | 2017-01-04 | 合肥工业大学 | 一种分离静止声源辐射声场与旋转声源辐射声场的方法 |
CN107478325A (zh) * | 2017-09-07 | 2017-12-15 | 安徽大学 | 一种非稳态平面声源的自由场还原方法 |
CN109798975A (zh) * | 2019-03-13 | 2019-05-24 | 安徽大学 | 一种采用声压和质点加速度测量的非稳态平面声源的自由场实时还原方法 |
CN111537058A (zh) * | 2020-04-16 | 2020-08-14 | 哈尔滨工程大学 | 一种基于Helmholtz方程最小二乘法的声场分离方法 |
CN111707353A (zh) * | 2020-05-29 | 2020-09-25 | 西安交通大学 | 一种基于近场声全息技术的回转曲面声场重建方法 |
CN111561991B (zh) * | 2020-06-08 | 2021-03-30 | 中国海洋大学 | 基于边缘填充和傅里叶变换的近场声全息方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1952627A (zh) * | 2006-10-27 | 2007-04-25 | 合肥工业大学 | 一种声场分离方法 |
CN101566495A (zh) * | 2009-06-05 | 2009-10-28 | 合肥工业大学 | 采用双面振速测量和二维空间傅立叶变换法分离声场的方法 |
-
2011
- 2011-12-06 CN CN 201110402251 patent/CN102494755B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1952627A (zh) * | 2006-10-27 | 2007-04-25 | 合肥工业大学 | 一种声场分离方法 |
CN101566495A (zh) * | 2009-06-05 | 2009-10-28 | 合肥工业大学 | 采用双面振速测量和二维空间傅立叶变换法分离声场的方法 |
Non-Patent Citations (2)
Title |
---|
CHUAN-XING BI ET AL.: "Sound field separation technique based on equivalent source method and its application in nearfield acoustic holography", 《ACOUSTICAL SOCIETY OF AMERICA》 * |
于飞等: "双全息面分离声场技术及其在声全息中的应用", 《声学学报》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103267571A (zh) * | 2013-05-10 | 2013-08-28 | 合肥工业大学 | 采用单面声压和振速测量的实时声场分离方法 |
CN103592022B (zh) * | 2013-11-25 | 2015-08-05 | 合肥工业大学 | 采用声压和质点加速度测量的实时声场分离方法 |
CN105181121B (zh) * | 2015-05-29 | 2018-02-06 | 合肥工业大学 | 采用加权迭代等效源法的高精度近场声全息方法 |
CN105181121A (zh) * | 2015-05-29 | 2015-12-23 | 合肥工业大学 | 采用加权迭代等效源法的高精度近场声全息算法 |
CN105588637A (zh) * | 2016-01-31 | 2016-05-18 | 广西科技大学 | 一种复杂稳定声场声压测试装置 |
CN105588637B (zh) * | 2016-01-31 | 2018-07-20 | 柳州市展虹科技有限公司 | 一种复杂稳定声场声压测试装置 |
CN106289505A (zh) * | 2016-07-21 | 2017-01-04 | 合肥工业大学 | 一种分离静止声源辐射声场与旋转声源辐射声场的方法 |
CN106289505B (zh) * | 2016-07-21 | 2019-01-11 | 合肥工业大学 | 一种分离静止声源辐射声场与旋转声源辐射声场的方法 |
CN107478325A (zh) * | 2017-09-07 | 2017-12-15 | 安徽大学 | 一种非稳态平面声源的自由场还原方法 |
CN109798975A (zh) * | 2019-03-13 | 2019-05-24 | 安徽大学 | 一种采用声压和质点加速度测量的非稳态平面声源的自由场实时还原方法 |
CN111537058A (zh) * | 2020-04-16 | 2020-08-14 | 哈尔滨工程大学 | 一种基于Helmholtz方程最小二乘法的声场分离方法 |
CN111707353A (zh) * | 2020-05-29 | 2020-09-25 | 西安交通大学 | 一种基于近场声全息技术的回转曲面声场重建方法 |
CN111707353B (zh) * | 2020-05-29 | 2021-11-09 | 西安交通大学 | 一种基于近场声全息技术的回转曲面声场重建方法 |
CN111561991B (zh) * | 2020-06-08 | 2021-03-30 | 中国海洋大学 | 基于边缘填充和傅里叶变换的近场声全息方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102494755B (zh) | 2013-11-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102494755B (zh) | 一种非稳态声场分离方法 | |
CN103267571B (zh) | 采用单面声压和振速测量的实时声场分离方法 | |
CN103592022B (zh) | 采用声压和质点加速度测量的实时声场分离方法 | |
EP2685288B1 (en) | Methods and Systems for Reconstruction of Low Frequency Particle Velocity Wavefields and Deghosting of Seismic Streamer Data | |
CN101566496B (zh) | 采用双面振速测量和等效源法分离声场的方法 | |
CN101566495B (zh) | 采用双面振速测量和二维空间傅立叶变换法分离声场的方法 | |
CN102680071B (zh) | 采用振速测量和局部近场声全息法的噪声源识别方法 | |
CN102047146B (zh) | 对低频地震数据测量结果的重构 | |
CN101403634B (zh) | 压力速度法声场分离方法 | |
US20150139444A1 (en) | Systems and methods for detecting transient acoustic signals | |
US11327195B2 (en) | Correction of source motion effects in seismic data recorded in a marine survey using a moving source | |
CN104678384B (zh) | 一种波束域的声压差互相关谱分析水下目标速度估计方法 | |
CN102736063A (zh) | 近场声源定位方法 | |
CN102183298A (zh) | 不规则单全息声压测量面分离非自由声场的方法 | |
Bleistein et al. | Migration/inversion: think image point coordinates, process in acquisition surface coordinates | |
CN103217214B (zh) | 一种采用双面振速测量和局部声全息法的分离声场方法 | |
CN109211399A (zh) | 一种锥面声源统计近场声全息方法 | |
US11650343B2 (en) | Directional designature of marine seismic survey data | |
CN109764956A (zh) | 一种基于组合优化正则化方法的近场声全息方法 | |
CN103308157A (zh) | 振动分布欠采样条件下的结构低频辐射声功率的估计方法 | |
Sarkar et al. | Information and linearity of time-domain complex demodulated amplitude and phase data in shallow water | |
US10317543B2 (en) | Estimation of a far field signature in a second direction from a far field signature in a first direction | |
Wapenaar et al. | On the retrieval of the directional scattering matrix from directional noise | |
US20160061979A1 (en) | Methods and systems to remove particle-motion-sensor noise from vertical-velocity data | |
Minakov et al. | Acoustic waveform inversion for ocean turbulence |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20131106 Termination date: 20161206 |
|
CF01 | Termination of patent right due to non-payment of annual fee |