贝叶斯正则化在球面等效源法中的应用研究 (1).pdf
《贝叶斯正则化在球面等效源法中的应用研究 (1).pdf》由会员分享,可在线阅读,更多相关《贝叶斯正则化在球面等效源法中的应用研究 (1).pdf(6页珍藏版)》请在咨信网上搜索。
1、第 卷第期 年月安徽工程大学学报J o u r n a l o fA n h u iP o l y t e c h n i cU n i v e r s i t yV o l N o A p r 文章编号:()收稿日期:基金项目:国家自然科学基金资助项目();安徽省教育厅自然科学基金资助项目(K J A );安徽高校学科(专业)拔尖人才学术基金资助项目(G X B J Z D )作者简介:赵越(),女,安徽池州人,硕士研究生.通信作者:张学勇(),男,安徽合肥人,教授,博士.贝叶斯正则化在球面等效源法中的应用研究赵越,张学勇,(安徽建筑大学 数理学院,安徽 合肥 ;安徽省建筑声环境重点实验室,
2、安徽 合肥 ;安徽建筑大学 声学研究所,安徽 合肥 )摘要:为拓宽空心球阵列的应用场景,提高三维空间声场重构精度,本文开展空心球阵列结合球面等效源方法的声场重建研究.在近场条件下,由空心球阵列采集到三维空间的声场信息,结合球面等效源法,对声场重构中存在的逆问题所常用的广义交叉验证正则化参数选择方法,提出贝叶斯正则化,并利用MA T L A B对两种正则化方法的重建结果进行对比.分别探讨了两种正则化方法在不同声源距离、声源频率、重构球面半径和信噪比情况下对声场重建精度的影响.结果表明,贝叶斯正则化在基于空心球阵列的三维空间声场重建中具有明显优势.关键词:空心球阵列;等效源法;贝叶斯正则化;广义交
3、叉验证;球面近场声全息中图分类号:O 文献标志码:A声场重构是噪声源精准识别与定位的关键技术,目前声场重构技术一般是由所选择的麦克风阵列和声全息方法结合实施.空心球阵列测量成本低且与测量声场之间的相互作用小.基于球面等效源法(S p h e r i c a lE q u i v a l e n tS o u r c eM e t h o d,S E S M)的近场声全息技术不仅扩大了声场重建区域,同时适用于任何声源结构,对三维空间声场重构精度更高.但在运用S E S M时,等效源数目远大于阵列传声器数目,这导致运用等效源强信息进行声场重构时,传递矩阵条件数过大,使矩阵方程不适定,为了得到不适定
4、逆问题的精确解,通常需要加入正则化方法.经典T i k h o n o v正则化方法的核心是最优正则化参数的选取,其参数的选取方法有l c u r v e法、B a y e s i a n正则化和广义交叉验证(G e n e r a l i z e dC r o s sV a l i d a t i o n,G C V)等.近年来,李敏宗等在三维封闭空间中采用T i k h o n o v正则化方法结合G C V准则作为声场重构的最优正则化参数处理方式;平国力等,将B a y e s i a n正则化、l c u r v e曲线法、G C V和H a l d正则化公式在基于等效源法的近场声全息
5、中进行比较,总结了B a y e s i a n正则化的优越性.本文采用空心球阵列结合S E S M实现三维空间声场重构,通过MAT L A B仿真比较了G C V和B a y e s i a n正则化参数选择方法在不同声源距离、声源频率、重构球面半径和信噪比下的三维空间声场重建精度,结果表明B a y e s i a n正则化参数选择方法在三维空间声场重建性能上具有优越性.图笛卡尔坐标系与球坐标系关系图声场重建基本理论空心球阵列测量理论在笛卡尔坐标系中,空间中任意一点用(x,y,z)表示;在球坐标系中,以空心球阵列的中心为坐标原点,任意一点的位置可用(r,)表示.其中r表示该点到原点的距离,
6、()为仰角,()为方位角.具体的关系如图所示.在球坐标系下,球面上的任何压力分布都可以展开为正交球面谐波之和,即通用的球谐函数,其表达式为:Ymn(,)(n)(nm)!(nm)!Pmn(c o s)ej m,()式中,m为球谐函数的阶;n为球谐函数的次;Pmn为m阶n次的连带勒让德函数;ej m为时间因子,j.对于空心球形传声器阵列,各传声器按近似均匀分布,等距离的固定在阵列球面上,测量声源所辐射的声场信息.在三维空间中,将空心球阵列中心处于球坐标系原点处,并在模型假设中忽略阵列对声场的影响.对于(rs,s,s)处的单极子点声源所辐射的入射波,此时位于空间中任意一点处的声压可表示为,:p(r)
7、i kc QsNnjn(k r)hn(k r)nm nYmn(,)Ymn(s,s),()式中,r(r,)为空间中任一点的坐标;k为波数,即角频率与声速c之比;为介质的密度;Qs为入射波强度;jn()为第一类球贝塞尔函数;hn()为第一类球汉克尔函数;()为复共轭;N为球谐函数的截止项数,已有研究表明,基于 通道的空心球阵列,在中低频时最优截止项数取,可得到良好的重构性能.基于S E S M的近场声全息理论在近场声全息理论中,等效源法结合空心球形传声器阵列的声场重构组合称为球面等效源法.其等效源强的求解通过阵列采集的声压数据和传递函数求出.在三维空间中,假设阵列对声场不产生影响,则由球面等效源法
8、可得阵列表面上任一点的声压为:p(r)jc kYiq(ri)G(r,ri),()式中,表示空气密度;q(ri)为(re,e,e)处的第i个等效源强度;G(r,ri)为自由场格林函数.其球谐函数扩展形式为:G(r,ri)ej k|rri|rri|j knnm njn(k r)hn(k ri)Ymn()Ymn(i),()通过在阵列上布置X个采样点,在声源内部或背离声源表面放置Y个等效源,此时全息面上所测得声压的矩阵表达式为:PhG Q,()式中,Ph为阵列上X个传声器所测声压组成的列向量;Q为Y个等效源强所组成的列向量;G为传递矩阵.而在实际测量过程中,等效源数量往往多于阵列数目,即YX,这使得上
9、述方程是欠定的,故在声源重构这一逆问题中,会导致解的不适定性,为解决这一问题.重构过程中运用T i k h o n o v正则化方法,此时上述方程的正则化解可表示为:Q(GHGI)GHPh,()式中,为需要求解的正则化参数;I为单位矩阵.基于以上形式,需对传递矩阵G进行奇异值分解,即:GU(XX)S(XY)V(YY),()式中,S为包含X个奇异值的矩阵;U为X阶酉矩阵;V为Y阶酉矩阵;将其带入上式中,等效源强度的正则化解可表示为:QV(SHSI)SHUHPh,()T i k h o n o v正则化方法的核心是正则化参数的选取.对于三维空间的声场重构,B a y e s i a n正则化和G
10、C V都可用来进行正则化参数的选取,为得到最优正则化参数,需对这两种方法进行比较.正则化参数选择理论()G C V.在基于G C V的T i k h o n o v正则化方法中,假设全息面上的阵列点均不参与重建计算,此时构造的G C V函数为:G()(G(GHG I)GHI)Ph(T r(IG(GHG I)GH),()安徽工程大学学报第 卷式中,T r为矩阵的迹,G C V是一种后验误差估计方法,正则化参数在G()取最小值时达到最优.()B a y e s i a n正则化.B a y e s i a n正则化的关键思想是将正则化设想为引入先验信息,以补偿测量过程中造成的信息损失.在贝叶斯方法
11、的理论背景下,正则化参数的估计值通过极大化边际似然函数来获得.其最优正则化参数的选取方法是通过对联合概率密度函数边际化得到其最大化后验概率密度函数|p,即:|pXx(UHPh)sxX Xx(sx),()再对其取负对数,此时得到正则化参数的成本函数为:A r g m i n J()Xxl n(sx)(X)l nXXx|UHPh|sx,()正则化参数在J()取最小值时达到最优.仿真实验与分析为探究B a y e s i a n正则化在三维空间声场重构中是否具有优越性和较佳的稳健性,本节通过数值仿真,比较在近场条件下,基于空心球阵列的球面等效源法中正则化参数选择方法:G C V和B a y e s
12、i a n正则化,在不同声源距离、声源频率、重构球面半径和信噪比下的重建性能.并通过相对误差来评价G C V和B a y e s i a n正则化的重建性能.定义相对误差为:e r rPPP ,()式中,P为重建面上的理论声压;P为重建面上的重建声压.声源距离和频率对重建结果的影响()仿真条件.通过布置一个体积速度为 m/s的单极子点声源位于源面中心,设置源面中心分别距坐标原点 、m处,声源频率在 H z以 H z为步长进行仿真;声源的辐射声场采用一空心球阵列测量,该阵列有 个传声器,各传声器之间按照近似均匀分布,等距离地布置在半径为 m的球面上,并以阵列中心为坐标原点;等效源面由 个等效点源
13、组成,近似均匀地分布在源面后方 m处(避免计算出现奇异);重建面位于空心球阵列与等效源面之间,据阵列中心 m处.在实际测量过程中,可能存在由各个方面带来的噪声影响,故仿真时加入 d B的高斯白噪声.()仿真结果及分析.当声源位置分别位于 、m处时,G C V和B a y e s i a n正则化在不同声源频率下的重建误差如图所示.从仿真结果可以看出,声源频率在 H z以下时,G C V和B a y e s i a n正则化均保持良好的重建结果,相对误差均保持在以下.随着声源距离的增大、频率的增加,相对误差也随之增大,这验证了S NAH适用于小全息和中低频环境下的声场重建.总体来看,B a y
14、e s i a n正则化的重建结果较G C V波动幅度小、重建误差更小,重建性能也更优异.B a y e s i a n正则化在整个频段内重建误差均保持在 以下,并且随着频率的增大逐渐趋于稳定,可见该方法稳健性较佳.图声源位置位于 m处图声源位置位于 m处第期赵越,等:贝叶斯正则化在球面等效源法中的应用研究图声源位置位于 m处重构球面半径对重建结果的影响()仿真条件.在()的基础上,选取声源频率为 H z和 H z时;设置源面中心分别距坐标原点 m和 m处;重建面位于空心球阵列与等效源面之间,分别据阵列中心 m和 m以 m为步长进行仿真.()仿真结果及分析.在固定声源位置和声源频率下,G C
15、V和B a y e s i a n正则化在不同重构球面半径下的重建误差如图所示.从仿真结果可以看出,当声源位置相同时,声源频率越高,重建误差越大;声源频率相同时,随着重构球面半径的增大,重建误差呈现上升趋势,并在靠近声源处,呈现下降趋势.综合看来,在低频近场条件下,G C V和B a y e s i a n正则化参数选择方法的重建结果均达到最佳.相比而言,B a y e s i a n正则化在每个条件下的重建相对误差均比G C V的小,这体现出了B a y e s i a n正则化较G C V的优越性.图声源位置位于 m处,f H z图声源位置位于 m处,f H z图声源位置位于 m处,f H
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 贝叶斯正则化在球面等效源法中的应用研究 1 贝叶斯 正则 球面 等效 中的 应用 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。