CN109764956B - 一种基于组合优化正则化方法的近场声全息方法 - Google Patents
一种基于组合优化正则化方法的近场声全息方法 Download PDFInfo
- Publication number
- CN109764956B CN109764956B CN201811607626.3A CN201811607626A CN109764956B CN 109764956 B CN109764956 B CN 109764956B CN 201811607626 A CN201811607626 A CN 201811607626A CN 109764956 B CN109764956 B CN 109764956B
- Authority
- CN
- China
- Prior art keywords
- wave number
- sound pressure
- regularization
- point
- cylindrical
- 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.)
- Active
Links
Images
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Holo Graphy (AREA)
Abstract
本发明公开了一种基于组合优化正则化方法的近场声全息方法,利用在全息面上测得的各点声压,在统计最优柱面近场声全息理论方法的基础上,对选取的波数矢量采用截断奇异值和标准Tikhonov组合优化正则化方法抑制噪声和随机误差引起的波数矢量的测量误差,即组合优化正则化方法,利用组合优化正则化方法和GCV方法求取叠加系数矩阵和在全息面上测得的各点声压,将壳体结构装备表面各点声压表示为共形测量面声压的线性叠加,即可获得壳体结构装备表面声场,充分展示了其在壳体辐射噪声分析中的有效性,适用于柱形壳体结构,将柱形水下航行器射声场以可视化的方式显示出来,从而可以直观看出辐射声场的大小及分布情况,具有重要的理论意义及工程应用价值。
Description
技术领域
本发明属于机械结构声辐射信号处理领域,特别涉及一种基于组合优化正则化方法的近场声全息方法。
背景技术
水下航行器的续航力大,机动性强,独立作战能力强,但其产生的噪声波在海水中可传递至数百海里,易被地方探测,严重削弱了其声隐身性能,因此,提高其声隐身性能是保证其安全性的关键。目前,首要解决的技术问题是获得水下航行器辐射声场,以实现对其隐身性能的精确评估。但对于内部声源繁多、传播路径复杂的水下航行器壳体装备,难以对其声隐身性能进行精确评估。当前,海军强国普遍利用装备有高精度固定式测量分析***的水声试验场测试评估水下航行器声隐身性能,然而,水声试验场选址要求极高、技术复杂且成本高昂。
近场声全息(Nearfield Acoustic Holography,NAH)是一种非常有效的噪声源识别、定位与声场可视化方法,是通过近场测试获取足够多低空间频率传播波和高空间频率倏逝波成分重建出分辨率不受分析波长限制的高精度全息图像和丰富空间声场信息。然而传统NAH是通过空间二维傅立叶变换来实现解卷积运算得到重建面声压,在全息计算过程中会带来窗效应和卷绕误差,因而要求测试孔径面积至少为声源面积两倍,对于大尺寸柱形壳体结构装备,满足NAH要求的精确测试很难实现。统计最优近场声全息技术(SONAH)直接通过全息测量面复声压的线性叠加实现空间声场的重建,从根本上避免了传统NAH存在的窗效应和卷绕误差,克服了对测量尺寸的严格要求。然而,国内学者对该方法进行了大量的理论和实验研究,不过往往偏重于统计最优平面近场声全息问题,对于水下航行器这种柱形声源研究较少,而且缺少对此声学不适定性问题的研究。
发明内容
本发明的目的在于提供一种基于组合优化正则化方法的近场声全息方法,以克服现有技术的不足。
为达到上述目的,本发明采用如下技术方案:
一种基于组合优化正则化方法的近场声全息方法,包括以下步骤:
步骤1),在最大波数限定的空间波数平面上非等间隔地选取各波数矢量;
步骤2),对选取的波数矢量采用截断奇异值和标准Tikhonov组合优化正则化方法抑制噪声和随机误差引起的波数矢量的测量误差;
步骤3),利用步骤2)组合优化正则化方法和GCV方法求取叠加系数矩阵C(rS)和在全息面上测得的各点声压p(rHn),将壳体结构装备表面各点声压表示为共形测量面声压的线性叠加,即p(rS,θ,z)=C(rS)p(rHn),即可获得壳体结构装备表面声场。
进一步的,步骤1)中,在最大波数限定的空间波数平面上非等间隔地选取各波数矢量:
最优柱面近场声全息方法过程如下:
稳态声波场Helmholtz方程在柱面坐标系下的行波解为
在柱面坐标系下,令Pn(r,n,kz)为p(r,θ,z)的二维傅氏变换,为
同时给出其傅氏反变换为
令上式r=a得,
由(1)(2)(3)(4)式可解得:
式中,Pn(a,kz)为p(a,θ,z)的二维傅氏变换。
由波数矢量K=(n,kz)确定的柱面上的空间频率域单元柱面波为:
重建柱面壳体结构装备表面上任意点rS=(rS,θ,z)处波数矢量为Km的单元柱面波都可以由全息面上所有点rHn=(rH,θn,zn)处的波数矢量为Km的单元柱面波叠加得到,即
式中,rHn=(rH,θn,zn)(n=1,2,…,N)为全息柱面上N个声压测量点,M为重建柱面和全息柱面上复声压所包含的单元柱面波的数目,C(rS)为叠加系数矩阵;
由上式确定的M个线性方程所构成的线性方程组表示成矩阵形式为
b=AC(rS) (9)
得到重建面上各点声压:
式中,p(rS,θ,z)为重建面上各点声压,p(rHn)为测量面上各点声压。
进一步的,在Nyguist采样定理基础上,在波数域辐射圆内和辐射圆外采用采用不同大小的采样间隔。
进一步的,在波数域辐射圆内采样间隔小于在辐射圆外采样间隔。
进一步的,采用截断奇异值和标准Tikhonov组合优化的正则化方法:对于截断点前的较大奇异值所对应的不含噪声的低空间频率成分采用TSVD方法进行处理,不进行正则化处理;对于截断点后的小奇异值,采用Tikhonov正则化方法,控制其对噪声的过度方法,同时保留高频细节信息;
经组合正则化得到的解为
进一步的,使用基于后验准则的广义交叉验证选取截断点和正则化参数:
采用广义交叉验证方法确定截断点k,对于公式(11)中正则化参数λ的选取,以截断点k为正则化处理截断点之后奇异值的正则化参数,优化后的组合正则化方法得到的解为:
结合公式(8)、(9)和(12),求取叠加系数矩阵C(rS)。
进一步的,根据求取的叠加系数矩阵C(rS)和在全息面上测得的各点声压p(rHn),将柱形壳体结构装备表面各点声压表示为共形测量面声压的线性叠加,即p(rS,θ,z)=C(rS)p(rHn),获得柱形壳体结构装备辐射声场。
与现有技术相比,本发明具有以下有益的技术效果:
本发明公开了一种基于组合优化正则化方法的近场声全息方法,利用在全息面上测得的各点声压p(rHn),在统计最优柱面近场声全息理论方法的基础上,提出改进方法,对选取的波数矢量采用截断奇异值和标准Tikhonov组合优化正则化方法抑制噪声和随机误差引起的波数矢量的测量误差,提出了组合优化正则化方法,利用组合优化正则化方法和GCV方法求取叠加系数矩阵C(rS)和在全息面上测得的各点声压p(rHn),将壳体结构装备表面各点声压表示为共形测量面声压的线性叠加,即p(rS,θ,z)=C(rS)p(rHn),即可获得壳体结构装备表面声场,本方法充分展示了其在壳体辐射噪声分析中的有效性,适用于柱形壳体结构,将柱形水下航行器射声场以可视化的方式显示出来,从而可以直观看出辐射声场的大小及分布情况,具有重要的理论意义及工程应用价值。
本发明在统计最优柱面近场声全息及技术的基础上,在最大波数限定的空间波数平面上非等间隔地选取各波数矢量,这样选取能够有效地获取声场中低波数主要能量,又可以保证高波数域能量的获取,保证了声场信息获得的完整性,有效地减少了噪声信号对重建精度的影响。
采用本发明提供的组合优化方法,比TSVD包含更多的细节信息,同时Tikhonov正则化保留高频细节信息且抑制了小奇异值对噪声的放大作用,能够获得更高的重建精度和稳定性。
附图说明
图1为柱形壳体结构表面理论声压图;
图2为传统方法重建柱形壳体结构表面声压图;
图3为本发明重建柱形壳体结构表面声压图;
图4为本发明和传统方法在θ=0母线上声压分布曲线图;
图5为本发明中改进近场声全息方法对噪声的抑制作用图;
图6本发明中改进近场声全息方法对噪声的抑制作用图;
图7为本发明方法的重建声压沿轴向分布图;
图8为本发明与现有不同正则化方法的重建误差图;
图9为本发明与不同正则化方法的重建误差图;
图10为本发明方法对噪声的抑制作用图;
图11为柱面波数矢量图。
具体实施方式
下面结合附图对本发明做进一步详细描述:
如图1至图10所示,一种基于组合优化正则化方法的近场声全息方法,包括以下步骤:
步骤1),在最大波数限定的空间波数平面上非等间隔地选取各波数矢量:在Nyguist采样定理基础上,在波数域辐射圆内和辐射圆外采用采用不同大小的采样间隔,提高重建精度;在波数域辐射圆内采样间隔小于在辐射圆外采样间隔;
步骤2),对选取的波数矢量采用截断奇异值(TSVD)和标准Tikhonov组合优化正则化方法抑制噪声和随机误差引起的波数矢量的测量误差;
步骤3),利用步骤2)组合优化正则化方法结合GCV方法求取叠加系数矩阵C(rS)和在全息面上测得的各点声压p(rHn),将壳体结构装备表面各点声压表示为共形测量面声压的线性叠加,即p(rS,θ,z)=C(rS)p(rHn),即可获得壳体结构装备表面声场。
本发明为模拟柱形壳体结构装备,采用轴向为2.0m、半径为0.2m的柱形壳体为研究对象;使用与柱形壳体结构装备表面共形的测量面,选择合理的测量参数,布置声压传感器,获取高效的全息声压数据;
步骤1)中,在最大波数限定的空间波数平面上非等间隔地选取各波数矢量:
最优柱面近场声全息方法(SOCNAH)过程如下:
稳态声波场Helmholtz方程在柱面坐标系下的行波解为
在柱面坐标系下,令Pn(r,n,kz)为p(r,θ,z)的二维傅氏变换,为
同时给出其傅氏反变换为
令上式r=a得,
由(1)(2)(3)(4)式可解得:
式中,Pn(a,kz)为p(a,θ,z)的二维傅氏变换。
由波数矢量K=(n,kz)确定的柱面上的空间频率域单元柱面波为:
根据声场叠加原理可知,重建柱面壳体结构装备表面上任意点rS=(rS,θ,z)处波数矢量为Km的单元柱面波都可以由全息面上所有点rHn=(rH,θn,zn)处的波数矢量为Km的单元柱面波叠加得到,即
式中,rHn=(rH,θn,zn)(n=1,2,…,N)为全息柱面上N个声压测量点,M为重建柱面和全息柱面上复声压所包含的单元柱面波的数目,C(rS)为叠加系数矩阵;
由上式确定的M个线性方程所构成的线性方程组表示成矩阵形式为
b=AC(rS) (9)
所以,只要确定了式(7)中的波数矢量Km,则叠加系数矩阵C(rS)就可以唯一确定(M≥N),进而可以得到重建面上各点声压:
式中,p(rS,θ,z)为重建面上各点声压,p(rHn)为测量面上各点声压,所以问题的关键在于M波数矢量Km的确定;
然而,SOCNAH重建过程具有不适定性,即矩阵A的病态性,产生的主要原因是矩阵A的小奇异值,实际中,因为高波数倏逝波衰减较快,即使在近场范围内进行测量,也很容易被测量误差掩盖,所以当全息面的测量数据中含有噪声等误差信号时,需要采用正则化方法限制小奇异值对测量误差的放大作用;进一步,在SOCNAH中,矩阵A的奇异值分布对应的是根据采样定理确定的波数域区间内的单元柱面波;在式(8)和(9)中,可以看到,要保证解的唯一性,需要矩阵A中M≥N,当M>N时,奇异值和各项单元柱面波并不是一一对应的;实际上,较大的奇异值对应的是传播波和低波数的倏逝波,较小的奇异值对应的是高波数的倏逝波,在建立矩阵A时,使用的高波数倏逝波的项数越多,那么小奇异值的数目也越多,对测量误差的放大作用越明显;然而,如果在重建过程中完全舍弃声场中的高波数倏逝波的能量,会导致声场能量的部分丢失,无法保证重建的精度;鉴于以上分析,本发明改进了空间波数域中波数矢量数组的选取方法:
在Nyguist采样定理基础上,在波数域辐射圆内采用小的采样间隔,在辐射圆外采用大的采样间隔,提高重建精度;如图11所示该方法在最大波数限定的空间波数平面上非等间隔地选取各波数矢量,在波数域辐射圆内采用小的采样间隔,在辐射圆外采用大的采样间隔;根据近场声全息原理,声场中的能量主要集中在低波数声波;所以在选取波数矢量时,改进方法在辐射圆外的高波数区域按照满足采样定理的最大采样间隔Δl选取波数矢量,在辐射圆内低波数域以小采样间隔选取波数矢量;这样选取可以有效地获取声场中的主要能量,又可以保证高波数域能量对声场重建的贡献,同时也减少了重建过程中测量误差的放大,保证声场重建精度。
采用截断奇异值(TSVD)和标准Tikhonov组合优化的正则化方法:本方法对不同大小的奇异值区别对待,对于截断点前的较大奇异值所对应的不含噪声的低空间频率成分采用TSVD方法进行处理,不进行正则化处理;而对于截断点后的小奇异值,由于其对应包含噪声的高空间频率成分,并对测量中噪声等引起误差的信号起到放大作用,所以必须进行处理,但若采用TSVD中直接进行截断的方法,会损失高频的细节信息,降低了重建的精度,因此,对于小的奇异值,采用Tikhonov正则化方法,控制其对噪声的过度方法,同时保留高频细节信息;
综合上述分析,经组合正则化得到的解为
使用基于后验准则的广义交叉验证(GCV)选取截断点和正则化参数;
由于实际中通常无法确定声场中的噪声大小,所以再使用经验公式无法确定截断点,SNR为信噪比;故,在此处采用GCV方法确定截断点k,对于公式(11)中后面一项的中正则化参数λ的选取,直接以该截断点k为正则化处理截断点之后奇异值的正则化参数,这样既保证了组合正则化中只有一个未定的参数,又可以保证截断点后面的奇异值没有被直接赋零,所以,组合正则化比TSVD包含更多的细节信息,同时正则化抑制了小奇异值对噪声的放大作用,能够获得较高的重建精度。
基于上述分析,最终组合正则化方法得到的解为
进一步地,利用组合优化正则化去除测量误差的方法,结合公式(8)、(9)和(12),求取叠加系数矩阵C(rS);
利用步骤3)中求取的叠加系数C(rS)和在全息面上测得的各点声压p(rHn),将柱形壳体结构装备表面各点声压表示为共形测量面声压的线性叠加,即p(rS,θ,z)=C(rS)p(rHn),获得柱形壳体结构装备表面声场。
使用与柱形壳体结构装备表面共形的测量面,根据测量距离、传感器间隔、测量面积,获取待测结构外侧全息声压数据,即获取在全息面上测得的各点声压p(rHn);
为验证本发明对柱形壳体结构装备辐射声场重构的有效性,以轴向长度为2.0m、半径为0.2m的柱形面模拟柱形壳体结构装备;在壳体内部选取60个脉动球声源组成线阵声源,脉动球的间距为0.02m,振动频率为f,仿真测量用的传声器阵元间距为0.1m,全息测量尺寸为2m,测量距离为0.27m,重建距离(测量面到重建面的距离)为0.05m,为了模拟实际测量情况,在全息面上加信噪比为40dB的随机噪声。使用与柱形壳体结构装备表面共形的测量面,根据测量参数获取待测结构外侧全息声压数据;所述测量参数包括测量距离、传感器间隔和测量面积;使用传声器测量声场全息声压数据;
在全息测量面施加信噪比为40的随机噪声,在不同声源频率时,图1、图2、图3、图4分别给出了柱形壳体结构表面理论声压图、传统方法重建柱形壳体结构表面声压图、本发明重建柱形壳体结构表面声压图和两种方法沿轴向的声压分布,可以看到,改进方法相比于传统方法可以更好地实现对重构面各个位置的重建,其重建值值与理论值都有更为精确的吻合。
图5给出了改进方法对噪声的抑制作用,可以看到,在改进方法中,当信噪比由40dB减小到20dB时,重建结果的相对平均误差变化很小,当信噪比再减小到10dB时,重建误差略有增大,充分说明了改进方法对噪声的抑制作用更好。
在图5中,不同的频率下,重建的总体相对误差大多在30%以下,在频率为700Hz时达到最低,但总体相对误差仍有24%,这是因为当测量面边长大于声波长的3倍时,有着更好的重建精度,但此仿真所采用的测量面边长仅小于声波长的2倍。进一步,增大测量面边长至声波长的3.2倍,线阵声源长度由1.2m增大到2m,其他仿真条件不变,由此得到改进方法在Tikhonov正则化结合GCV正则参数选取方法下的重建误差如图6所示,由图6可以看到,在不同频带上,重建的总体相对误差均有大幅减小,其中在低频段已减小到10%以下,与图3-8比较,在信噪比同为40dB时,在低频段200~800Hz时,其重建总体相对误差都至少减小了10%,尤其在200Hz时,重建总体相对误差减小了近20%。为了直观地看到改进方法的重建效果,图7给出了改进方法的重建值与理论值的吻合程度,从中可以看到,在各个重构点上改进方法的重建值与理论值都有很好的吻合,充分验证了改进方法的有效性。
图8和图9给出了在全息尺寸分别为声波长1.7倍和3.2倍时各正则化方法重建的总体相对误差,信噪比为20dB。可以看出,组合优化方法在高频段有更好的重建性能。显示了其对高频声场重建的优势;为了进一步验证组合优化方法对噪声的抑制作用,图10给出了本发明的组合优化正则化方法对在测量面不同信噪比下重建的总体相对误差;在不同的信噪比条件下,组合方法重建的总体相对误差变化不大,尤其是在高频段,改方法的重建误差更小,充分说明了该方法对噪声有好的抑制作用。
在现有的近场声全息技术中,对于水下航行器这种柱形声源研究较少,缺少对此声学不适定性问题的研究。本发明方法在现有的统计最优近场声全息技术基础上,提出改进方法,在改进方法基础上,提出了组合优化正则化方法,有效提高了壳体结构装备表面声场重构精度和稳定性,从而为柱形壳体结构装备的声隐身性能评估提供了理论基础,具有重要的理论意义及工程应用价值。
Claims (8)
1.一种基于组合优化正则化方法的近场声全息方法,其特征在于,包括以下步骤:
步骤1),在最大波数限定的空间波数平面上非等间隔地选取各波数矢量;
步骤2),对选取的波数矢量采用截断奇异值和标准Tikhonov组合优化正则化方法抑制噪声和随机误差引起的波数矢量的测量误差;具体的:对于截断点前的较大奇异值所对应的不含噪声的低空间频率成分采用TSVD方法进行处理,不进行正则化处理;对于截断点后的小奇异值,采用Tikhonov正则化方法,控制其对噪声的过度方法,同时保留高频细节信息;
经组合正则化得到的解为
步骤3),利用步骤2)组合优化正则化方法和GCV方法求取叠加系数矩阵C(rS)和在全息面上测得的各点声压p(rHn),将壳体结构装备表面各点声压表示为共形测量面声压的线性叠加,即p(rS,θ,z)=C(rS)p(rHn),即可获得壳体结构装备表面声场。
2.根据权利要求1所述的一种基于组合优化正则化方法的近场声全息方法,其特征在于,
步骤1)中,在最大波数限定的空间波数平面上非等间隔地选取各波数矢量:
最优柱面近场声全息方法过程如下:
稳态声波场Helmholtz方程在柱面坐标系下的行波解为
在柱面坐标系下,令Pn(r,n,kz)为p(r,θ,z)的二维傅氏变换,为
同时给出其傅氏反变换为
令上式r=a得,
由(1)(2)(3)(4)式可解得:
式中,Pn(a,kz)为p(a,θ,z)的二维傅氏变换;
由波数矢量K=(n,kz)确定的柱面上的空间频率域单元柱面波为:
重建柱面壳体结构装备表面上任意点rS=(rS,θ,z)处波数矢量为Km的单元柱面波都可以由全息面上所有点rHn=(rH,θn,zn)处的波数矢量为Km的单元柱面波叠加得到,即
由上式确定的M个线性方程所构成的线性方程组表示成矩阵形式为
b=AC(rS) (9)
得到重建面上各点声压:
式中,p(rS,θ,z)为重建面上各点声压,p(rHn)为测量面上各点声压。
3.根据权利要求2所述的一种基于组合优化正则化方法的近场声全息方法,其特征在于,在Nyguist采样定理基础上,在波数域辐射圆内和辐射圆外采用采用不同大小的采样间隔。
4.根据权利要求3所述的一种基于组合优化正则化方法的近场声全息方法,其特征在于,在波数域辐射圆内采样间隔小于在辐射圆外采样间隔。
7.根据权利要求6所述的一种基于组合优化正则化方法的近场声全息方法,其特征在于,根据求取的叠加系数矩阵C(rS)和在全息面上测得的各点声压p(rHn),将柱形壳体结构装备表面各点声压表示为共形测量面声压的线性叠加,即p(rS,θ,z)=C(rS)p(rHn),获得柱形壳体结构装备辐射声场。
8.根据权利要求1所述的一种基于组合优化正则化方法的近场声全息方法,其特征在于,使用与柱形壳体结构装备表面共形的测量面,根据测量距离、传感器间隔、测量面积,获取在全息面上测得的各点声压p(rHn)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811607626.3A CN109764956B (zh) | 2018-12-27 | 2018-12-27 | 一种基于组合优化正则化方法的近场声全息方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811607626.3A CN109764956B (zh) | 2018-12-27 | 2018-12-27 | 一种基于组合优化正则化方法的近场声全息方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109764956A CN109764956A (zh) | 2019-05-17 |
CN109764956B true CN109764956B (zh) | 2020-03-31 |
Family
ID=66450899
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811607626.3A Active CN109764956B (zh) | 2018-12-27 | 2018-12-27 | 一种基于组合优化正则化方法的近场声全息方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109764956B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110780273B (zh) * | 2019-11-04 | 2022-03-04 | 电子科技大学 | 一种混合正则化方位超分辨成像方法 |
CN111174904B (zh) * | 2020-01-13 | 2021-10-29 | 中国船舶科学研究中心(中国船舶重工集团公司第七0二研究所) | 用于水下噪声源分离的声全息测试共形阵列 |
CN111707354B (zh) * | 2020-05-29 | 2021-07-13 | 西安交通大学 | 一种基于平面测试的圆柱壳体声场分步组合重构方法 |
CN111707353B (zh) * | 2020-05-29 | 2021-11-09 | 西安交通大学 | 一种基于近场声全息技术的回转曲面声场重建方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE1911231A1 (de) * | 1968-03-06 | 1970-02-12 | Holotron Corp | Ultraschall-Mammograph |
JPH09113248A (ja) * | 1995-10-24 | 1997-05-02 | Isuzu Motors Ltd | 近距離音場ホログラフィー装置 |
CN102121847B (zh) * | 2010-12-16 | 2012-09-26 | 合肥工业大学 | 一种瞬态声场重建方法 |
CN103743469B (zh) * | 2013-12-11 | 2015-12-02 | 西安交通大学 | 压电水声换能器声辐射模态测量方法及*** |
CN103743474A (zh) * | 2013-12-24 | 2014-04-23 | 广西科技大学 | 基于边界元法的近场声全息重建方法 |
CN103712684B (zh) * | 2013-12-25 | 2016-01-20 | 广西科技大学 | 声场重建方法 |
CN105844114A (zh) * | 2016-04-28 | 2016-08-10 | 广西科技大学 | 一种非共形测量的近场声全息声场重建方法 |
CN108629081A (zh) * | 2018-04-02 | 2018-10-09 | 西安交通大学 | 基于统计最优近场声全息的壳体装备表面声场重构方法 |
-
2018
- 2018-12-27 CN CN201811607626.3A patent/CN109764956B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109764956A (zh) | 2019-05-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109764956B (zh) | 一种基于组合优化正则化方法的近场声全息方法 | |
JP5496192B2 (ja) | 音響場を再構成するための方法 | |
CN107478325B (zh) | 一种非稳态平面声源的自由场还原方法 | |
CN109211399B (zh) | 一种锥面声源统计近场声全息方法 | |
US20060080418A1 (en) | Farfield analysis of noise sources | |
JP6386556B2 (ja) | 広周波数帯域音響ホログラフィ | |
Chu et al. | Deconvolution using CLEAN-SC for acoustic source identification with spherical microphone arrays | |
Yang et al. | Functional delay and sum beamforming for three-dimensional acoustic source identification with solid spherical arrays | |
CN109489796A (zh) | 一种基于单元辐射法的水下复杂结构辐射噪声源定位识别与声辐射预报方法 | |
CN111707353B (zh) | 一种基于近场声全息技术的回转曲面声场重建方法 | |
CN108629081A (zh) | 基于统计最优近场声全息的壳体装备表面声场重构方法 | |
Chen et al. | A modified Helmholtz equation least squares method for reconstructing vibroacoustic quantities on an arbitrarily shaped vibrating structure | |
CN109798975A (zh) | 一种采用声压和质点加速度测量的非稳态平面声源的自由场实时还原方法 | |
Hu et al. | Achieving high-resolution 3D acoustic imaging in a large-aspect-ratio cabin by the non-synchronous measurements | |
Thakur et al. | Sound source localization of harmonic sources in entire 3D space using just 5 acoustic signals | |
CN114252148B (zh) | 一种基于长椭球波叠加的声场重建方法 | |
Burfeindt et al. | Receive-beamforming-enhanced linear sampling method imaging | |
CN109884187B (zh) | 一种适用于板状结构的基于压缩感知的超声导波场损伤检测方法 | |
CN113465729A (zh) | 一种矩形管道内噪声源识别与声场重建方法 | |
Arcondoulis et al. | Experimental application of an acoustic beamforming array pairing method using CLEAN-SC | |
Lan-Yue et al. | Combined method of near field acoustic holography and focused beamforming for noise source identification in enclosed space | |
Jiang et al. | Sound power estimation method based on dual arrays sound intensity scaling and sound pressure superposition | |
Hald | Beamforming and wavenumber processing | |
Li et al. | Research on improving measurement accuracy of acoustic transfer function of underwater vehicle | |
Lou et al. | Inverse patch transfer function with fast iterative shrinkage-thresholding algorithm as a tool for sparse source identification |
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 |