贝叶斯正则化在球面等效源法中的应用研究 (1).pdf

第 卷第期 年月安徽工程大学学报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 )作者简介:赵越(),女,安徽池州人,硕士研究生.通信作者:张学勇(),男,安徽合肥人,教授,博士.贝叶斯正则化在球面等效源法中的应用研究赵越,张学勇,(安徽建筑大学 数理学院,安徽 合肥 ;安徽省建筑声环境重点实验室,安徽 合肥 ;安徽建筑大学 声学研究所,安徽 合肥 )摘要:为拓宽空心球阵列的应用场景,提高三维空间声场重构精度,本文开展空心球阵列结合球面等效源方法的声场重建研究.在近场条件下,由空心球阵列采集到三维空间的声场信息,结合球面等效源法,对声场重构中存在的逆问题所常用的广义交叉验证正则化参数选择方法,提出贝叶斯正则化,并利用MA T L A B对两种正则化方法的重建结果进行对比.分别探讨了两种正则化方法在不同声源距离、声源频率、重构球面半径和信噪比情况下对声场重建精度的影响.结果表明,贝叶斯正则化在基于空心球阵列的三维空间声场重建中具有明显优势.关键词:空心球阵列;等效源法;贝叶斯正则化;广义交叉验证;球面近场声全息中图分类号: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时,等效源数目远大于阵列传声器数目,这导致运用等效源强信息进行声场重构时,传递矩阵条件数过大,使矩阵方程不适定,为了得到不适定逆问题的精确解,通常需要加入正则化方法.经典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正则化公式在基于等效源法的近场声全息中进行比较,总结了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表示该点到原点的距离,()为仰角,()为方位角.具体的关系如图所示.在球坐标系下,球面上的任何压力分布都可以展开为正交球面谐波之和,即通用的球谐函数,其表达式为:Ymn(,)(n)(nm)!(nm)!Pmn(c o s)ej m,()式中,m为球谐函数的阶;n为球谐函数的次;Pmn为m阶n次的连带勒让德函数;ej m为时间因子,j.对于空心球形传声器阵列,各传声器按近似均匀分布,等距离的固定在阵列球面上,测量声源所辐射的声场信息.在三维空间中,将空心球阵列中心处于球坐标系原点处,并在模型假设中忽略阵列对声场的影响.对于(rs,s,s)处的单极子点声源所辐射的入射波,此时位于空间中任意一点处的声压可表示为,:p(r)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的近场声全息理论在近场声全息理论中,等效源法结合空心球形传声器阵列的声场重构组合称为球面等效源法.其等效源强的求解通过阵列采集的声压数据和传递函数求出.在三维空间中,假设阵列对声场不产生影响,则由球面等效源法可得阵列表面上任一点的声压为: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,这使得上述方程是欠定的,故在声源重构这一逆问题中,会导致解的不适定性,为解决这一问题.重构过程中运用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 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正则化的关键思想是将正则化设想为引入先验信息,以补偿测量过程中造成的信息损失.在贝叶斯方法的理论背景下,正则化参数的估计值通过极大化边际似然函数来获得.其最优正则化参数的选取方法是通过对联合概率密度函数边际化得到其最大化后验概率密度函数|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 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的球面上,并以阵列中心为坐标原点;等效源面由 个等效点源组成,近似均匀地分布在源面后方 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 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 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 z图声源位置位于 m处,f H z信噪比对重建结果的影响()仿真条件.在()和()的基础上,选取声源频率为 H z时;设置源面中心距坐标原点 m处;重建面据阵列中心 m处.信噪比设置为 d B以d B为步长增长,并选取信噪比分别为 d B和 d B时,理论与重建结果进行比较.()仿真结果及分析.不同信噪比下两种参数选择方法相对误差如表所示.由表可以看出,在低频近场条件下,随着信噪比的升高,重建误差逐渐减小,这是因为信噪比越高时,抑制噪声的能力就越强,重建效果也就越好.在 d B和 d B时理论声压级与两种方法重建声压级分别见图、和图 ,可见在高信噪比条件下,两种正则化参数选择方法的重建结果均良好;而在低信噪比情况下,G C V重建的性能并不是很稳定,在信噪比越低的时候,效果越不好,B a y e s i a n正则化较G C V仍保持一个良好的重建性能.因此,B a y e s i a n正则化参数选择方法具有明显优势.安徽工程大学学报第 卷表不同信噪比下两种参数选择方法相对误差表(f H z)方法信噪比 d B d B d B d B d B d B d B d B d BG C