CN103592022B - 采用声压和质点加速度测量的实时声场分离方法 - Google Patents
采用声压和质点加速度测量的实时声场分离方法 Download PDFInfo
- Publication number
- CN103592022B CN103592022B CN201310606976.9A CN201310606976A CN103592022B CN 103592022 B CN103592022 B CN 103592022B CN 201310606976 A CN201310606976 A CN 201310606976A CN 103592022 B CN103592022 B CN 103592022B
- Authority
- CN
- China
- Prior art keywords
- time
- measurement plane
- acoustic pressure
- sound source
- formula
- 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
Links
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种采用声压和质点加速度测量的实时声场分离方法,其特征是在目标声源和干扰声源之间布置两个相互平行的测量平面H和测量平面H1,同步采集两个测量平面上的声压时域信号;利用两个测量平面上的声压时域信号,通过有限差分法获得测量平面H上的质点加速度时域信号;再利用测量平面H上的声压时域波数谱和质点加速度时域波数谱、已知的时域脉冲响应函数,实时分离出目标声源单独在测量平面H上所辐射的声压时域波数谱,进而获得目标声源单独在测量平面H上所辐射的声压时域信号。本发明方法无需进行任何求逆和正则化运算处理,计算速度快、稳定性高;本发明方法具备实时分离声场的能力,可用于噪声干扰环境下现场分析目标声源的时变辐射特性。
Description
技术领域
本发明涉及物理专业中噪声类领域声场分离方法。
背景技术
在实际工程中,目标声源位于测量平面的一侧,而在测量平面的另一侧往往存在干扰声源,这些干扰声源所产生的声场影响了对目标声源所辐射声场的准确测量,因此需要采用一定的声场分离方法将干扰声源的影响从测量结果中分离出来。到目前为止,国内外学者已提出多种声场分离方法,这些方法可大致分为五类:一是基于空间傅里叶变换法(SFT)的声场分离技术:G.V.Frisk等在1980年首次提出采用SFT法间接测量海洋底面的反射系数,并建立了基于SFT法的双面声场分离理论;M.Tamura于1990年将G.V.Frisk等提出的方法用于测量斜入射时材料的反射系数;M.T.Cheng等对M.Tamura提出的方法进行了推广,建立了柱面坐标下的双测量声场分离式,并用于实现散射声场的分离。二是基于球面波叠加的声场分离方法:1956年,J.Pachner采用球面波叠加法实现了任意波场中行波和驻波声场的分离;G.Weinreich等在1980年对J.Pachner提出的方法作了进一步改进,建立了基于双球面测量的声场分离理论;Y.Braikia等在2013年分别采用双球面声压测量的分离方法和单球面声压振速测量的分离方法实现了小空间范围内的声源识别。三是基于统计最优近场声全息(SONAH)的声场分离技术:J.Hald在对SONAH研究的基础上,提出了基于双全息面声压测量的统计最优声场分离技术;F.Jacobsen等在J.Hald提出的方法的基础上,提出了基于声压振速测量的统计最优声场分离技术;E.Fernandez-Grande等在2011年提出了基于双面振速测量的统计最优声场分离技术。四是基于边界元法(BEM)的声场分离技术:C.Langrenne等在2007年提出一种基于BEM的双面声场分离方法;随后,E.G.Williams等在2008年提出一种基于BEM和声压振速测量(Cauchy数据)的声场分离方法。五是基于等效源法(ESM)的声场分离技术:C.X.Bi提出的基于ESM的声场分离技术,适用于任意形状的测量面,且计算稳定性好、计算精度高;E.Fernandez-Grande等在2012年提出了基于ESM的双面振速测量的分离技术和单面声压振速测量的分离技术。但是,上述声场分离方法都仅仅实现了单一频率或频带下的声场分离。若要实现任意时刻下的声场分离,上述声场分离方法将不再适用。为实现任意时刻下的声场分离,X.Z.Zhang等在2012年提出一种基于双平面声压测量的时域声场分离方法。在该分离方法中,每个时刻目标声源所辐射的声压信号需要通过解卷运算得到,由于该解卷运算过程具有不适定性,需要采用奇异值分解和正则化来稳定求解过程。而正则化的采用使得该分离方法耗时长,且分离精度很大程度上依赖于正则化参数的合理选取。
发明内容
为实现任意时刻下的声场分离,本发明提供一种计算速度快、稳定性高的采用声压和质点加速度测量的实时声场分离方法,以期获得实时分离声场的能力,进而用于噪声干扰环境下现场分析目标声源的时变辐射特性。
本发明解决技术问题所采用的技术方案是:
本发明采用声压和质点加速度测量的实时声场分离方法的特点是按如下步骤进行:
步骤a、在目标声源Mo和干扰声源Md之间布置相互平行的测量平面H和测量平面H1;所述目标声源Mo和干扰声源Md均可辐射任意线性声场信号;在所述测量平面H和测量平面H1上分别均匀分布有M个测量网格点,测量平面H和测量平面H1上的测量网格大小和测量网格点位置相同;以目标声源Mo的中心为坐标原点建立笛卡尔坐标系,所述笛卡尔坐标系中的xoy坐标平面平行于测量平面H,z轴方向垂直于测量平面H;测量平面H上的测量网格点坐标为(x,y,zH),测量平面H1上的测量网格点坐标为(x,y,zH1);同步采集测量平面H上各测量网格点(x,y,zH)处的声压时域信号p(x,y,zH,t)和测量平面H1上各测量网格点(x,y,zH1)处的声压时域信号p(x,y,zH1,t);
步骤b、按照式(1)进行有限差分获得测量平面H上各测量网格点(x,y,zH)处的质点加速度时域信号a(x,y,zH,t):
式(1)中,t为时间;ρ0为介质密度;
步骤c、对于测量平面H上的声压时域信号p(x,y,zH,t)按照式(2)进行二维空间傅里叶变换获得声压时域波数谱P(kx,ky,zH,t),对测量平面H上的质点加速度时域信号a(x,y,zH,t)按照式(3)进行二维空间傅里叶变换获得质点加速度时域波数谱A(kx,ky,zH,t),
在式(2)和式(3)中,j表示虚数单位;kx、ky分别为x、y方向的波数;
步骤d、构建声压时域波数谱P(kx,ky,zH,t)、质点加速度时域波数谱A(kx,ky,zH,t)、已知的时域脉冲响应函数h(kx,ky,0,t)和目标声源Mo单独在测量平面H上所辐射的声压时域波数谱Po(kx,ky,zH,t)之间的关系如式(4)所示:
Po(kx,ky,zH,t)=0.5[P(kx,ky,zH,t)+A(kx,ky,zH,t)*h(kx,ky,0,t)] (4)
在式(4)中,“*”表示两个时间函数的卷积运算;
将式(4)中的时间t离散为tn=(n-1)Δt,其中Δt为采样时间间隔,n=1,2,...,N,N为采样点总数,式(4)的离散形式如式(5)所示:
步骤e、对于所述声压时域波数谱Po(kx,ky,zH,tn)按照式(6)进行二维空间傅里叶反变换,分离出每个时刻下目标声源Mo单独在测量平面H上所辐射的声压时域信号po(x,y,zH,tn):
本发明采用声压和质点加速度测量的实时声场分离方法的特点也在于:
所述测量平面H和测量平面H1上各测量网格点的声压时域信号p(x,y,zH,t)和p(x,y,zH1,t)是采用声压传感器阵列一次快照测量获得。
所述干扰声源Md为噪声源、反射源或散射源。
理论模型:
在辐射任意线性声场信号的目标声源Mo和辐射任意线性声场信号的干扰声源Md之间分别布置两个相互平行的测量平面H和测量平面H1。根据Euler公式,测量平面H上的质点加速度时域信号a(x,y,zH,t)与测量平面H上所测得的声压时域信号p(x,y,zH,t)的偏导数存在成正比的关系,如式(7)所示:
在式(7)中,ρ0为介质密度。
利用测量平面H上的声压时域信号p(x,y,zH,t)和测量平面H1上的声压时域信号p(x,y,zH1,t),按照式(8)进行有限差分获得测量平面H上的声压时域信号p(x,y,zH,t)的偏导数:
联合式(7)和(8)获得测量平面H上的质点加速度时域信号a(x,y,zH,t),即:
根据声波的叠加原理,在测量平面H上所测得的声压时域信号p(x,y,zH,t)等于目标声源Mo在测量平面H上所辐射的声压时域信号po(x,y,zH,t)与干扰声源Md在测量平面H所辐射的声压时域信号pd(x,y,zH,t)之和,如式(10)所示:
p(x,y,zH,t)=po(x,y,zH,t)+pd(x,y,zH,t) (10)
由于质点加速度是矢量,在测量平面H上所获得的质点加速度时域信号a(x,y,zH,t)应等于目标声源Mo在测量平面H上所辐射的质点加速度时域信号ao(x,y,zH,t)与干扰声源Md在测量平面H所辐射的质点加速度时域信号ad(x,y,zH,t)之差,如式(11)所示:
a(x,y,zH,t)=ao(x,y,zH,t)-ad(x,y,zH,t) (11)
定义函数f(x,y,z,t)关于x、y的二维空间傅里叶变换为:
式(12)中,F(kx,ky,z,t)表示f(x,y,z,t)的时域波数谱。
分别对式(10)和式(11)进行关于x、y的二维空间傅里叶变换可得:
P(kx,ky,zH,t)=Po(kx,ky,zH,t)+Pd(kx,ky,zH,t) (13)
A(kx,ky,zH,t)=Ao(kx,ky,zH,t)-Ad(kx,ky,zH,t) (14)
由声压和质点加速度在时域波数域的关系可知:
Po(kx,ky,zH,t)=Ao(kx,ky,zH,t)*h(kx,ky,0,t)(15)
Pd(kx,ky,zH,t)=Ad(kx,ky,zH,t)*h(kx,ky,0,t) (16)
式(15)和式(16)中,h(kx,ky,0,t)为声压和质点加速度之间的时域脉冲响应函数,其表达式为:
在式(17)中,c为声速,J0为阶数为0的第一类贝塞尔函数,H(t)为Heaviside函数。
将式(16)代入式(13)中得到:
P(kx,ky,zH,t)=Po(kx,ky,zH,t)+Ad(kx,ky,zH,t)*h(kx,ky,0,t) (18)
将式(14)两边分别与h(kx,ky,0,t)做卷积得到:
A(kx,ky,zH,t)*h(kx,ky,0,t)=Ao(kx,ky,zH,t)*h(kx,ky,0,t)-Ad(kx,ky,zH,t)*h(kx,ky,0,t) (19)
将式(15)代入式(19)中得到:
A(kx,ky,zH,t)*h(kx,ky,0,t)=Po(kx,ky,zH,t)-Ad(kx,ky,zH,t)*h(kx,ky,0,t) (20)
将式(18)与式(20)相加可进一步得到:
Po(kx,ky,zH,t)=0.5[P(kx,ky,zH,t)+A(kx,ky,zH,t)*h(kx,ky,0,t)] (21)
式(21)中,P(kx,ky,zH,t)和A(kx,ky,zH,t)可分别通过对测得的声压时域信号p(x,y,zH,t)和通过有限差分获得的质点加速度时域信号a(x,y,zH,t)进行二维空间傅里叶变换得到,h(kx,ky,0,t)为已知函数。利用式(21)按如下过程可实时分离出每个时刻下目标声源Mo单独在测量平面H上所辐射的声压时域波数谱Po(kx,ky,zH,t):
将式(21)中的时间t离散为tn=(n-1)Δt,其中Δt为采样时间间隔,n=1,2,...,N,N为采样点总数。
取n=1时,分离得到:
取n=2时,分离得到:
……
取n=N时,分离得到:
由式(22)到式(24)的分离过程可知,要分离出tn时刻的Po(kx,ky,zH,tn),只需要利用tn时刻的P(kx,ky,zH,tn)和ti(i=1,2,...,n)时刻的A(kx,ky,zH,ti),而P(kx,ky,zH,tn)和A(kx,ky,zH,ti)可通过对测得的声压时域信号p(x,y,zH,t)和通过有限差分获得的质点加速度时域信号a(x,y,zH,t)进行二维空间傅里叶变换得到。因此,一旦测得某个时刻下测量平面H和测量平面H1上的声压时域信号,即可分离出该时刻的声压时域波数谱,实现声场的实时分离。
对分离出的时域波数谱Po(kx,ky,zH,tn)按照式(25)进行二维空间傅里叶反变换,最终获得每个时刻下目标声源Mo单独在测量平面H上所辐射的声压时域信号po(x,y,zH,tn)
与已有技术相比,本发明有益效果体现在:
1、本发明方法的分离过程是一个简单的正向求解过程,不需要进行任何求逆和正则化运算处理,因而计算速度快、稳定性高。
2、本发明方法具备实时分离声场的能力,因而更适用于噪声干扰环境下现场分析目标声源的时变辐射特性。
附图说明
图1为本发明采用声压和质点加速度测量的实时声场分离方法示意图;
图2(a)为本发明方法的实验声源、测量平面位置分布图;
图2(b)为本发明方法的实验测量平面网格划分示意图;
图3(a)为本发明方法的实验测量点A处声压时域信号图;
图3(b)为本发明方法的实验测量点B处声压时域信号图;
图3(c)为本发明方法的实验测量点C处声压时域信号图;
图4(a)为本发明方法实验所获得的相位评价因子Ep的数值分布图;
图4(b)为本发明方法实验所获得的幅值评价因子Ea的数值分布图;
图5(a)为本发明方法t=9.38ms时刻目标声源Mo单独在测量平面H上所辐射的测量声压场pm;
图5(b)为本发明方法t=10.78ms时刻目标声源Mo单独在测量平面H上所辐射的测量声压场pm;
图5(c)为本发明方法t=9.38ms时刻的混合声压场pc;
图5(d)为本发明方法t=10.78ms时刻的混合声压场pc;
图5(e)为t=9.38ms时刻采用本发明方法分离出的声压场pe;
图5(f)为t=10.78ms时刻采用本发明方法分离出的声压场pe。
具体实施方式
参见图1,本实施例采用声压和质点加速度测量的实时声场分离方法按如下步骤分离出目标声源Mo单独在测量平面H上所辐射的声压时域信号:
步骤a、在目标声源Mo和干扰声源Md之间布置相互平行的测量平面H和测量平面H1;所述目标声源Mo和干扰声源Md均可辐射任意线性声场信号;在所述测量平面H和测量平面H1上分别均匀分布有M个测量网格点,测量平面H和测量平面H1上的测量网格大小和测量网格点位置相同;以目标声源Mo的中心为坐标原点建立笛卡尔坐标系,所述笛卡尔坐标系中的xoy坐标平面平行于测量平面H,z轴方向垂直于测量平面H;测量平面H上的测量网格点坐标为(x,y,zH),测量平面H1上的测量网格点坐标为(x,y,zH1);同步采集测量平面H上各测量网格点(x,y,zH)处的声压时域信号p(x,y,zH,t)和测量平面H1上各测量网格点(x,y,zH1)处的声压时域信号p(x,y,zH1,t)。
步骤b、按照式(26)进行有限差分获得测量平面H上各测量网格点(x,y,zH)处的质点加速度时域信号a(x,y,zH,t):
式(26)中,t为时间;ρ0为介质密度。
步骤c、对于测量平面H上的声压时域信号p(x,y,zH,t)按照式(27)进行二维空间傅里叶变换获得声压时域波数谱P(kx,ky,zH,t),对测量平面H上的质点加速度时域信号a(x,y,zH,t)按照式(28)进行二维空间傅里叶变换获得质点加速度时域波数谱A(kx,ky,zH,t),
在式(27)和式(28)中,j表示虚数单位;kx、ky分别为x、y方向的波数。
步骤d、构建声压时域波数谱P(kx,ky,zH,t)、质点加速度时域波数谱A(kx,ky,zH,t)、已知的时域脉冲响应函数h(kx,ky,0,t)和目标声源Mo单独在测量平面H上所辐射的声压时域波数谱Po(kx,ky,zH,t)之间的关系如式(29)所示:
Po(kx,ky,zH,t)=0.5[P(kx,ky,zH,t)+A(kx,ky,zH,t)*h(kx,ky,0,t)] (29)
在式(29)中,“*”表示两个时间函数的卷积运算。
将式(29)中的时间t离散为tn=(n-1)Δt,其中Δt为采样时间间隔,n=1,2,...,N,N为采样点总数。式(29)的离散形式如式(30)所示
步骤e、对于所述声压时域波数谱Po(kx,ky,zH,tn)按照式(31)进行二维空间傅里叶反变换,分离出每个时刻下目标声源Mo单独在测量平面H上所辐射的声压时域信号po(x,y,zH,tn)
测量平面H和测量平面H1上各测量网格点的声压时域信号p(x,y,zH,t)和p(x,y,zH1,t)是采用声压传感器阵列一次快照测量获得。干扰声源Md为噪声源、反射源或散射源。
方法的检验:
在测量平面H和测量平面H1的一侧布置音箱S1和音箱S3作为目标声源Mo,在另外一侧布置音箱S2作为干扰声源Md。采用本发明方法将测量平面H上目标声源Mo单独所辐射的声压时域信号分离出来,并与其测量声压时域信号进行比较。
本实验中,目标声源Mo(包含音箱S1和音箱S3)和干扰声源Md(音箱S2)都辐射余弦调制信号,其表达式为:
在式(32)中,频率f0=800Hz。
测量平面H和H1与声源之间的位置关系参见图2(a),测量平面H位于z=0m的平面上,测量平面H1位于z=-0.02m的平面上。音箱S1和音箱S3的纸盆中心分别位于(0.1m,0.45m,0.145m)和(0.55m,0.35m,0.14m)处,音箱S2的纸盆中心位于(0.4m,0.45m,-0.14m)处。测量平面H和H1的面积大小都为0.7m×0.7m,其上都均匀地分布15×15个测量点,参见图2(b)。时域信号采样频率为25600Hz,采样点数为512。
为检验本发明方法在时域内的声场分离效果,在测量平面H上选取了三个测量点,即测量点A、测量点B和测量点C,其位置分别为A(0.1m,0.45m,0m),B(0.4m,0.45m,0m)和C(0.55m,0.35m,0m),其中测量点A和测量点C分别正对着目标声源Mo的音箱S1和音箱S3,测量点B正对着干扰声源Md的音箱S2。图3(a)、图3(b)和图3(c)分别对应测量点A、测量点B和测量点C,图中曲线a表示目标声源Mo单独所辐射的测量声压时域信号,图中曲线b表示包含目标声源Mo和干扰声源Md所辐射的混合声压时域信号,图中曲线c表示采用本发明方法分离出的声压时域信号,比较图中的曲线a和曲线b可以看出,干扰声源Md对目标声源Mo在测量平面H上所辐射的声压时域信号造成了较大的干扰;比较曲线a和曲线c可以看出,采用本发明方法可以很好地消除干扰声源Md的影响,从而分离出目标声源Mo单独在测量平面H上所辐射的声压时域信号。
为了更加客观地评价本发明方法的分离效果,在此定义了两个评价因子,它们的表达式分别为:
在式(33)和式(34)中,<>表示求平均值,下标“m”表示目标声源单独辐射的测量声压值,下标“e”表示分离声压值。评价因子Ep是用来衡量目标声源单独辐射的测量声压值和分离声压值之间的相位误差,当Ep的值越靠近1时,相位误差越小。评价因子Ea是用来衡量目标声源单独辐射的测量声压值和分离声压值之间的幅值误差,当Ea的值越靠近0时,幅值误差越小。运用式(33)和(34)分别计算了测量平面H上各个测量点处的Ep和Ea。图4(a)表示所有测量点的Ep值,图中等高线值为0.9。图4(b)表示所有测量点的Ea值,图中等高线值为0.1。参见图4(a)和图4(b),在大多数测量点处,不论是相位还是幅值,测量声压值和分离声压值都吻合地较好,只是在测量平面边缘处吻合程度略低。
为检验本发明方法在空间域的声场分离效果,选取了两个时刻t=9.38ms和t=10.78ms。图5(a)和图5(b)分别为9.38ms和10.78ms时刻目标声源Mo单独在测量平面H上所辐射的测量声压场pm,图5(c)和图5(d)分别为9.38ms和10.78ms时刻的混合声压场pc(包含目标声源Mo和干扰声源Md所辐射的声压场),图5(e)和图5(f)分别为9.38ms和10.78ms时刻采用本发明方法分离出的声压场pe。比较图5(a)和图5(c)、图5(b)和图5(d)可以看出,干扰声源Md对目标声源Mo单独在测量平面H上所辐射的声压场造成了严重的干扰;比较图5(a)和图5(e)、图5(b)和图5(f)可以看出,采用本发明方法可以较好地消除干扰声源Md的影响,从而获得目标声源Mo单独在测量平面H上所辐射的声压场。
Claims (3)
1.采用声压和质点加速度测量的实时声场分离方法,其特征是按如下步骤进行:
步骤a、在目标声源Mo和干扰声源Md之间布置相互平行的测量平面H和测量平面H1;所述目标声源Mo和干扰声源Md均能够辐射任意线性声场信号;在所述测量平面H和测量平面H1上分别均匀分布有M个测量网格点,测量平面H和测量平面H1上的测量网格大小和测量网格点位置相同;以目标声源Mo的中心为坐标原点建立笛卡尔坐标系,所述笛卡尔坐标系中的xoy坐标平面平行于测量平面H,z轴方向垂直于测量平面H;测量平面H上的测量网格点坐标为(x,y,zH),测量平面H1上的测量网格点坐标为(x,y,zH1);同步采集测量平面H上各测量网格点(x,y,zH)处的声压时域信号p(x,y,zH,t)和测量平面H1上各测量网格点(x,y,zH1)处的声压时域信号p(x,y,zH1,t);
步骤b、按照式(1)进行有限差分获得测量平面H上各测量网格点(x,y,zH)处的质点加速度时域信号a(x,y,zH,t):
式(1)中,t为时间;ρ0为介质密度;
步骤c、对于测量平面H上的声压时域信号p(x,y,zH,t)按照式(2)进行二维空间傅里叶变换获得声压时域波数谱P(kx,ky,zH,t),对测量平面H上的质点加速度时域信号a(x,y,zH,t)按照式(3)进行二维空间傅里叶变换获得质点加速度时域波数谱A(kx,ky,zH,t),
在式(2)和式(3)中,j表示虚数单位;kx、ky分别为x、y方向的波数;
步骤d、构建声压时域波数谱P(kx,ky,zH,t)、质点加速度时域波数谱A(kx,ky,zH,t)、已知的时域脉冲响应函数h(kx,ky,0,t)和目标声源Mo单独在测量平面H上所辐射的声压时域波数谱Po(kx,ky,zH,t)之间的关系如式(4)所示:
Po(kx,ky,zH,t)=0.5[P(kx,ky,zH,t)+A(kx,ky,zH,t)*h(kx,ky,0,t)] (4)
在式(4)中,“*”表示两个时间函数的卷积运算;
将式(4)中的时间t离散为tn=(n-1)Δt,其中Δt为采样时间间隔,n=1,2,...,N,N为采样点总数,式(4)的离散形式如式(5)所示:
步骤e、对于所述声压时域波数谱Po(kx,ky,zH,tn)按照式(6)进行二维空间傅里叶反变换,分离出每个时刻下目标声源Mo单独在测量平面H上所辐射的声压时域信号po(x,y,zH,tn):
2.根据权利要求1所述的采用声压和质点加速度测量的实时声场分离方法,其特征是:所述测量平面H和测量平面H1上各测量网格点的声压时域信号p(x,y,zH,t)和p(x,y,zH1,t)是采用声压传感器阵列一次快照测量获得。
3.根据权利要求1所述的采用声压和质点加速度测量的实时声场分离方法,其特征是:所述干扰声源Md为噪声源、反射源或散射源。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310606976.9A CN103592022B (zh) | 2013-11-25 | 2013-11-25 | 采用声压和质点加速度测量的实时声场分离方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310606976.9A CN103592022B (zh) | 2013-11-25 | 2013-11-25 | 采用声压和质点加速度测量的实时声场分离方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103592022A CN103592022A (zh) | 2014-02-19 |
CN103592022B true CN103592022B (zh) | 2015-08-05 |
Family
ID=50082253
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310606976.9A Expired - Fee Related CN103592022B (zh) | 2013-11-25 | 2013-11-25 | 采用声压和质点加速度测量的实时声场分离方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103592022B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104596636B (zh) * | 2014-12-29 | 2018-03-27 | 广西科技大学 | 声场分离方法 |
CN105844114A (zh) * | 2016-04-28 | 2016-08-10 | 广西科技大学 | 一种非共形测量的近场声全息声场重建方法 |
CN106872019A (zh) * | 2017-01-20 | 2017-06-20 | 湖北文理学院 | 一种基于质点振速的部分场分解方法 |
CN107478325B (zh) * | 2017-09-07 | 2019-07-26 | 安徽大学 | 一种非稳态平面声源的自由场还原方法 |
CN109798975B (zh) * | 2019-03-13 | 2021-10-01 | 安徽大学 | 一种采用声压和质点加速度测量的非稳态平面声源的自由场实时还原方法 |
CN110579275B (zh) * | 2019-10-21 | 2022-03-11 | 南京南大电子智慧型服务机器人研究院有限公司 | 一种基于球形矢量传声器阵列实现声场分离的方法 |
CN110749374B (zh) * | 2019-10-22 | 2021-09-17 | 国网湖南省电力有限公司 | 一种建筑物内变压器结构传声分离方法及装置 |
CN112146751B (zh) * | 2020-09-27 | 2022-04-01 | 合肥工业大学 | 一种基于时域等效源法的实时声场分离方法 |
CN112924176B (zh) * | 2021-03-09 | 2023-03-24 | 合肥工业大学 | 一种基于时域等效源法的滚动轴承故障声学诊断方法 |
CN113063490B (zh) * | 2021-03-12 | 2022-03-01 | 北京科技大学 | 一种基于声压和质点振速双面测量的声场分离方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101403634A (zh) * | 2008-10-20 | 2009-04-08 | 合肥工业大学 | 压力速度法声场分离方法 |
CN101566495A (zh) * | 2009-06-05 | 2009-10-28 | 合肥工业大学 | 采用双面振速测量和二维空间傅立叶变换法分离声场的方法 |
CN102494755A (zh) * | 2011-12-06 | 2012-06-13 | 合肥工业大学 | 一种非稳态声场分离方法 |
CN102901559A (zh) * | 2012-09-27 | 2013-01-30 | 哈尔滨工程大学 | 一种基于单面测量和局部声全息法的分离声场方法 |
CN103217214A (zh) * | 2013-03-18 | 2013-07-24 | 哈尔滨工程大学 | 一种采用双面振速测量和局部声全息法的分离声场方法 |
CN103267571A (zh) * | 2013-05-10 | 2013-08-28 | 合肥工业大学 | 采用单面声压和振速测量的实时声场分离方法 |
-
2013
- 2013-11-25 CN CN201310606976.9A patent/CN103592022B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101403634A (zh) * | 2008-10-20 | 2009-04-08 | 合肥工业大学 | 压力速度法声场分离方法 |
CN101566495A (zh) * | 2009-06-05 | 2009-10-28 | 合肥工业大学 | 采用双面振速测量和二维空间傅立叶变换法分离声场的方法 |
CN102494755A (zh) * | 2011-12-06 | 2012-06-13 | 合肥工业大学 | 一种非稳态声场分离方法 |
CN102901559A (zh) * | 2012-09-27 | 2013-01-30 | 哈尔滨工程大学 | 一种基于单面测量和局部声全息法的分离声场方法 |
CN103217214A (zh) * | 2013-03-18 | 2013-07-24 | 哈尔滨工程大学 | 一种采用双面振速测量和局部声全息法的分离声场方法 |
CN103267571A (zh) * | 2013-05-10 | 2013-08-28 | 合肥工业大学 | 采用单面声压和振速测量的实时声场分离方法 |
Non-Patent Citations (2)
Title |
---|
基于声压-振速测量的单面声场分离技术;张永斌 等;《物理学报》;20091231(第12期) * |
声场分离技术及其在近场声全息中的应用;于飞 等;《物理学报》;20050228;第54卷(第2期);789-797 * |
Also Published As
Publication number | Publication date |
---|---|
CN103592022A (zh) | 2014-02-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103592022B (zh) | 采用声压和质点加速度测量的实时声场分离方法 | |
CN103267571B (zh) | 采用单面声压和振速测量的实时声场分离方法 | |
CN102494755B (zh) | 一种非稳态声场分离方法 | |
Sarradj et al. | A Python framework for microphone array data processing | |
Potthast | A study on orthogonality sampling | |
CN101566495B (zh) | 采用双面振速测量和二维空间傅立叶变换法分离声场的方法 | |
CN101566496B (zh) | 采用双面振速测量和等效源法分离声场的方法 | |
CN107478325B (zh) | 一种非稳态平面声源的自由场还原方法 | |
CN101403634B (zh) | 压力速度法声场分离方法 | |
CN102121847B (zh) | 一种瞬态声场重建方法 | |
CN102183298B (zh) | 不规则单全息声压测量面分离非自由声场的方法 | |
CA2873477C (en) | Systems and methods for detecting transient acoustic signals | |
CN106680376A (zh) | 一种基于三维声强阵列的隔声测量***与测量方法 | |
CN108051800B (zh) | 基于球面近场声全息重构无功声强的室内噪声源定位方法 | |
CN113063490B (zh) | 一种基于声压和质点振速双面测量的声场分离方法 | |
CN110487393A (zh) | 采用单面声压和质点振速测量的非稳态自由场还原方法 | |
CN103217214B (zh) | 一种采用双面振速测量和局部声全息法的分离声场方法 | |
CN109238444B (zh) | 一种采用稀疏测量的声场分离方法 | |
CN106324278A (zh) | 一种基于声参量阵的风速测量方法 | |
CN109798975A (zh) | 一种采用声压和质点加速度测量的非稳态平面声源的自由场实时还原方法 | |
Senet et al. | An iterative technique to determine the near surface current velocity from time series of sea surface images | |
Zhang et al. | Separation of nonstationary sound fields in the time-wavenumber domain | |
Yan et al. | Experimental Study of Sound Travel-Time Estimation Method in Stored Grain. | |
Zea et al. | Separation of rail and wheel contributions to pass-by noise with sparse regularization methods | |
Günel | On the statistical distributions of active intensity directions |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for 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: 20150805 Termination date: 20181125 |
|
CF01 | Termination of patent right due to non-payment of annual fee |