统计最优平面近场声全息原理与声场分离技术.docxVIP

统计最优平面近场声全息原理与声场分离 统计最优平面近场声全息原理与声场分离技术 3 李卫兵? 陈 剑 于 飞 毕传兴 陈心昭 ( 合肥工业大学动态测试中心 ,合肥 230009) (2004 年 3 月 24 日收到 ;2004 年 6 月 17 日收到修改稿) 测量孔径尺寸的有限性 ,在基于空间傅里叶变换的平面近场声全息中会带来窗效应和卷绕误差 . 为了克服窗 效应和卷绕误差 ,引入了统计最优平面近场声全息技术 . 运用声场叠加原理 ,证明了统计最优平面近场声全息的 理论公式 . 通过在空间波数域限定 kx , ky 的取值范围 ,并离散其确定的空间波数面的途径 ,提出了一种确定波数矢 量的方法 . 为了克服常规统计最优平面近场声全息技术的应用局限性 ———全息面一侧的声场必须为自由声场 ,提 出了适用于统计最优平面近场声全息的 、基于双全息面测量的空间声场分离技术 . 通过实验和数值仿真对理论推 导的正确性进行了验证 . 关键词 : 统计最优 , 平面近场声全息 , 波数矢量 , 声场分离 PACC : 4320 , 4363 空间域中全息面上复声压的线性叠加 ,来计算重建 面上的复声压和表面粒子振速. 不存在卷积运算 , 所以它可以从根本上解决窗效应和卷绕误差7 . 本 文运用声场叠加原理 ,从理论上证明了 SOPNAH 的 全息重建与预测公式. 以 Nyguist 采样定理为基础 , 采取在空间波数域限定 kx , ky 的最大 、最小值 ,并离 散其确定的空间波数面的途径 ,提出了一种确定波 数矢量的新方法 . 为了克服常规 SOPNAH 的应用局 限性 ———全息面一侧的声场必须为自由声场 ,推导 了适用于 SOPNAH 的 、基于双全息面测量的空间声 场分离公式 ,并且进行了一系列数值仿真和实验来 验证理论推导的正确性 . 11 引 言 在平 面 近 场 声 全 息 中 , 根 据 Helmholtz 方 程 和 Kirchiholff 公式知 : 与源面平行的平面上的复声压 , 可由源面复声压与 Dirichlet 边界条件下的格林函数 二维空间卷积得到. 而基于空间傅里叶变换的平面 近场声全息 ( FFT- based planar NAH) 是通过空间傅里 叶变换来实现解卷积运算的 ,它把源面复声压与格 林函数二维空间卷积转化为空间波数 (频率) 域的简 单乘法 ,从而使基于空间傅里叶变换的平面近场声 全息技术成为解决声学逆问题的一种快速有效的方 法1 —5 . 由于实际全息面复声压是在有限测量孔径 21 理论模型的推导 (有限的测量面积) 上离散进行的 ,并且在数值重建 过程中 ,是用离散傅里叶变换来实现解卷积运算的 , 这样就会在全息计算过程中带来窗 效 应 和 卷 绕 误 差 . 在测量孔径面积足够大 ( 至少是声源尺寸的两 倍) ,以至其边缘区域的声压幅值可以近似为零的情 况下 ,可以通过补零的方式来减小窗效应和卷绕误 差 ;否则 ,测量孔径边缘区域不连续的声压经傅里叶 变换后 会 带 来 更 大 的 误 差6 . 对 于 面 积 较 大 的 声 源 ,要满足这样的测量孔径要求很困难 . 统计最优平面近场声全息 ( SOPNAH) ,直接通过 设 p ( r ) = p ( x , y , z ) 为 空 间 声 场 中 任 意 点 ( x , y , z) 处的复声压 ,其频率为 f , k = 2πfΠc 为声波 数 , c 为空气中声波的传播速度. 根据基于空间傅 里叶变换的平面近场声全息的基本理论可知 : 当声 源全位于 z ≤0 的负半空间 ,且 z 0 的正半空间为 自由声场时 ,近场范围内 z 0 的任何平面上的声压 可以看成是无数个空间波数域的平面传播波和倏逝 波的叠加 ,即有 3 国家自然科学基金 ( 批准号 和高等学校博士学科点专项科研基金 ( 批准号 :20020359005) 资助的课题. ? 通讯联系人. E- mail : hf - lwb @sina . com 1254物理学报54 卷MMp ( x , y , z)= ∑∑P ( kx , ky )= c ( r ) ∑P ( Km) exp [ - i ( Km ·r ) ]∑ n SH nkx ky×exp [ - i ( kx x +n = 1Nm = 1ky y + kz z) ] , (1)= ∑cn ( rS 1254 物 理 学 报 54 卷 M M p ( x , y , z) = ∑∑P ( kx , ky ) = c ( r ) ∑P ( Km) exp [ - i ( Km