基于巴特沃斯滤波器的近场声全息滤波方法.pdf
《基于巴特沃斯滤波器的近场声全息滤波方法.pdf》由会员分享,可在线阅读,更多相关《基于巴特沃斯滤波器的近场声全息滤波方法.pdf(6页珍藏版)》请在咨信网上搜索。
1、第45卷第2 4期2023年1 2 月舰船科学技术SHIP SCIENCEANDTECHNOLOGYVol.45,No.24Dec.,2023基于巴特沃斯滤波器的近场声全息滤波方法卢锦芳1,楼京俊,刘树勇1,毛为民3(1.海军工程大学动力工程学院,湖北武汉430 0 33;2.海军工程大学舰船与海洋学院,湖北武汉430033;3.中国人民解放军91 0 0 7 部队,浙江舟山31 6 0 0 0)摘要:针对近场声全息研究中声场重构时,指数滤波器无法兼顾通带衰减率和高频带衰减率的局限性,采用巴特沃斯滤波器作为滤波函数以消除噪声干扰,对2 个同向振动的脉动球源的重构声场进行仿真研究,给出巴特沃斯滤
2、波器随参数变化的特性曲线。基于重构误差最小原则,给出该滤波器的最优参数,并与指数滤波器效果进行对比。结果表明,巴特沃斯滤波器可在声场重构过程中进行有效滤波,且在高频段,巴特沃斯滤波器重建误差明显小于指数滤波器重建误差,滤波效果更好,重建精度更高。关键词:近场声全息;波数域滤波;巴特沃斯滤波器;指数滤波器中图分类号:TB532文章编号:1 6 7 2-7 6 49(2 0 2 3)2 4-0 1 6 0-0 6Near-field acoustic holography filtering method based on Butterworth filter(1.College of Power
3、 Engineering,Naval University of Engineering,Wuhan 430033,China;2.College of Naval Architecture and Ocean Engineering,Naval University of Engineering,Wuhan 430033,China;3.No.91007 Unit of PLA Liberation Army,Zhoushan 316000,China)Abstract:In view of the limitation that the exponential filter cannot
4、take into account both passband attenuation rateand high frequency attenuation rate when the sound field is reconstructed in the study of near-field acoustic holography,But-terworth filter is adopted as the filter function to eliminate the noise.The reconstructed sound field of two pulsating spheric
5、alsources with the same direction of vibration is simulated.The characteristic curve of Butterworth filter changing with theparameters is given.Based on the principle of minimum reconstruction error,the optimal parameters of the filter are given,and the filtering effect is compared with that of the
6、exponential filter.The results show that Butterworth filter can effectivelyfilter in the process of sound field reconstruction,and the reconstruction error of Butterworth filter is significantly smallerthan that of exponential filter in high frequency band and its filtering effect is better and the
7、reconstruction accuracy is higher.Key words:near-field acoustic holography;k-space filtering;Butterworth filter;exponential filter0引言近场声全息技术(NearfieldAcousticHolography,NAH)由Maynard等在2 0 世纪8 0 年代提出 1-2 1,是声学研究中一种具有高分辨率重建精度的噪声源定位和声场可视化方法,在声场重构以及振动噪声控制等领域有着广泛应用。近场声全息通过多个传感器组成的测量面在被测声源物体表面的近场区域获取全息测量数据
8、,建立重建面与测量面之间的空间变换关系,进收稿日期:2 0 2 2-1 1-30基金项目:湖北省自然科学基金面上项目(2 0 2 2 CFB405);国家自然科学基金青年科学基金资助项目(52 2 0 1 38 9)作者简介:卢锦芳(1 999),女,博士研究生,研究方向为近场声全息。文献标识码:ALU Jin-fang,LOU Jing-jun,LIU Shu-yong,MAO Wei-mindoi:10.3404/j.issn.1672-7649.2023.24.029而重建三维空间的声压场、声强场、质点振速场或者预测远场指向性,从而为声源识别定位、有效控制结构振动噪声提供技术参考 3。在
9、声场逆向重构过程中,位于高波数区的逝波成分将与噪声误差一起被逆传递因子按指数规律放大,导致重建结果产生很大误差,严重影响声场重建的准确性 4。针对这一问题,在重建过程中,必须对高波数区噪声误差进行滤波处理,以确保重建精度,为此国内外学者进行了广泛研究。Veronesi等 5 提出一第45卷种在波数域加低通滤波窗的方法来抑制高波数成分波对噪声的放大效应,该方法采取指数滤波器对波数域噪声信号进行处理,其算法简单,可操作性强,但该方法没有考虑声波频率、信噪比以及测量面与源面间距离等因素的影响,使用中常会出现因截止波数取值不当以致滤波效果不理想的情况。赵晨等 6 对指数滤波器中窗函数参数的合理设置进行
10、分析,针对截止波数和窗函数陡度系数2 个关键参数,总结了在不同声源频率下取得最优重建效果的截止波数和陡度系数的最优取值,对窗函数的设计优化具有重要意义。张德俊7 在系统存在噪声情况下,对指数滤波器的截止波数进行修正,同时还考虑到指数滤波器仅在测量有限距离时有效,随着测量距离增大,滤波窗收效甚微,提出一种带约束条件的最小二乘法滤波函数。罗禹贡等 8 总结了候逝波随距离变化的规律并设计了相应的空间频域滤波器。由于指数滤波器仅对被逆传递函数按指数规律放大的噪声误差有效,对其他区域噪声误差过滤效果不佳,何其伟等 9 提出一种波数域循环加窗逼近法,通过循环加窗筛选最接近理论声压幅值的数据作为重建数据,从
11、而确保声场重构的准确性。莫登沅等 1 0-1 2 将空间域滤波和波数域滤波相结合,提出一种空间域即波数域联合迭代滤波窗,该滤波窗可对重建面边缘误差进行有效消除,并且经过多次选代,明显提高了声场重建精度。在指数滤波器实际使用中发现,窗函数陡度系数的取值变化决定着窗函数的衰减速率,陡度系数取值越小,指数窗函数衰减越快,相较于截止波数的取值变化,陡度系数取值变化对反演效果影响较小,故其对指数滤波器衰减程度的调整也有限。此外,指数滤波器的陡度系数会同时影响通带内和高频带内滤波曲线的衰减速度,且不能时时兼顾两者之间的平衡,甚至会导致有价值的候逝波被剔除,影响滤波效果,进而影响声场重建精度。本文针对上述问
12、题,选取2 个脉动球源作为研究对象,采用巴特沃斯数字滤波器(Butterworthfilter)作为滤波函数,对脉动球源声场反演过程进行数值仿真研究。基于重构误差最小原则,得到了使声场重建误差最小的巴特沃斯滤波器最优参数,并比较了该参数条件下,巴特沃斯滤波器与常用指数滤波器的效果。1基于Fourier空间声场变换的近场声全息在均匀理想流体介质中,小振幅声波的三维波动方程为:卢锦芳,等:基于巴特沃斯滤波器的近场声全息滤波方法卷积转化为波数域乘积可得:P(kx,ky,ZH)=P(kx,ky,zs)GD(kx,ky,zH-zs)。式中:P(kx,ky,zH)为测量面声压p(xH,yH,zH)的空间F
13、ourier变换;P(kx,ky,zs)为重建面声压p(xs,ys,zs)的空间Fourier变换;Gp(k,ky,zH-zs)为Green函数gD(x,y,z)的空间Fourier变换。其解析表达式为:Gp(kr,ky,ZH-zs)=ei(cH-s)k.。式中:zZH为测量面z轴坐标,zs为声源面z轴坐标,kx、k y 分别为x轴和y轴方向波数,k,取值存在以下2种情况(+k=k2被称之为辐射圆,辐射圆区分了声波2 种不同的传播模式)1 3。当+k0z0的空间为自由声场时,由Green公式可得任意平面上的声压为:p(XH,H,ZH)=JJ,p(xs,ys,zs)gD.(XH-Xs,yH-ys
14、,ZH-zs)dxsdyso式中:S为用于积分的无穷大平面;P(xH,YH,ZH)为测量面声压;p(xs,ys,zs)为重建面声压;gD(xH-xs,yH-y,zH-zs)为满足Dirichlet边界条件的格林函数。其表达式为:8D(a,2)=k2时:k,=i/(k2+k)-k2。当kz=-i/(k+磅)-k2时,其对应的解没有物理意义,在无穷远处趋向于无穷大,故此时对k,取值为kz=i/(+)-k2。在此情况下,k为纯虚数,GD(ks,k,ZH-zs)为波数空间的高频成分,在声波传播过程中,其幅值在传播方向随距离变化按指数规律衰减,称其为候逝波,为了获取到候逝波中包含的重要信息,需在声源近场
15、对其进行测量。当已知测量面上声压分布p(xH,YH,zH)时,对式(5)进行空间二维Fourier逆变换,可得重建面上声压为:p(xs,ys,zs)=F-IP(kx,ky,zH)GB(k,ky,ZH-zs)。(9)式中,G(kx,ky,zH-zs)为声压重建的逆向传递因子。基于空间声场变换的平面近场声全息逆向重构流程,如图1 所示。测量面数据分析测量面声压二维空间测量面角谱数据采集傅里叶通P(k.k,zn)变换G(k,k,ZH-z)重建面声压二维空间加窗后重建面角谱加窗滤波p(a,y,z)博里叶过P(k,k,zm)f(c,M变换图1 基于空间声场变换的平面近场声全息逆向重构流程图Fig.1 R
16、everse reconstruction flow chart of planar near-field acous-tical holograph based on spatial sound field transformation2滤波器在实际声场反演过程中,位于高波数区的候逝波成分将与噪声误差一起被逆传递函数按指数规律放大,若不对噪声信号进行低通滤波,将导致重建结果产生很大误差,严重影响声场重建的准确性,因此必须在声场逆向重构过程中进行滤波处理2.1指数滤波器指数滤波器是声场重构过程中常用的滤波方法,其算法简单、可操作性强、应用广泛。截止波数和陡度系数是该方法的2 个关键系数,截止波
17、数对参与声场重构过程的测量面声压角谱范围做出了限定,其取值受声波频率、信噪比以及测量面与源面间距离等因素的影响。截止波数取值过大,则达不到滤波效果,重建误差大;截止波数取值偏小,将导致过滤波,影响重建精度。陡度系数决定着窗函数的衰减速率,其取值越小,窗函数在截止波数处截止得越陡峭,且陡度系数会同时影响通带和高频带内滤波曲线的衰减速度,有时无法兼顾两者之间的平衡,甚至会导致有价舰船科学技术值的逝波被剔除,影响滤波效果,进而影响声场重建准确性。(8)指数滤波器计算公式为:exp(kr/kc-1)/a,k,kc,2V(Gex(1-/k)/al,k,ke.1式中:kr=V(k+);k e 为截止波数;
18、为窗函数陡度系数。2.2巴特沃斯滤波器巴特沃斯滤波器被称为最大平坦滤波器,其在通带的频率响应曲线符合最大限度平坦特性,是连续衰减的滤波器,该滤波器能够保持增益的平坦特性,尤为适合设计低通滤波器进行滤波处理,消除高频噪声干扰。巴特沃斯滤波器过渡带的陡峭程度与阶次成正比,随着阶次的增加,幅频特性曲线越陡峭,巴特沃斯滤波器更接近于理想低通状态,抑制噪声的能力变强。随着截止频率增加,系统响应速度变快,但系统稳定性也随之下降 4。巴特沃斯滤波器计算公式为:1重建面角谱f(C,M)=-P(k.,k,z)1+C2式中:C=(100.1p-1)为巴特沃斯滤波器截止频率;p为通带最大允许衰减系数;N为巴特沃斯滤
19、波器阶数;;k=/&+)。3仿真研究与参数分析3.1两脉动球源声场模型如图2 所示,两脉动球源的中心坐标分别为C1=(0.175,0,0)和C2=(-0.175,0,0),半径ro为0.0 1 m,两脉动球同相振动,表面振动速度Ua=1m/s,频率f-1000Hz。重建面zs=0.01m,测量面zH=0.05m,重建距离为0.0 4m,采样间距为0.0 5m,测点数量为2 52 5,在测量面上等间距分布。脉动球源辐射声压表达式:r式中:r为单个声源到测点距离;二绝对值为声压振幅。A=-i.kpcoQ0。4元式中:p为空气密度;co为空气中声速;Q=4元rUa第45卷(10)(11)2N。krA
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 巴特沃斯 滤波器 近场 全息 滤波 方法
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。