能量法估计全息面信噪比研究

代 振，何其伟，万海波

(海军工程大学 动力工程学院，湖北 武汉 430033)

能量法估计全息面信噪比研究

代 振，何其伟，万海波

(海军工程大学 动力工程学院，湖北 武汉 430033)

在近场声全息技术中，全息面信噪比是一个非常重要的参数，但无法直接测量得到。为了解决该问题，给出了基于能量法的波数域内信噪比估值公式，并根据该公式分析估值误差的来源，指出全息面采样间隔对估值具有重要影响。为了判别估值是否准确，提出了选取全息面不同区域的数据进行重复估计的判别方法。当估值不准确时，指出可以通过减小全息面采样间隔来降低估值误差。仿真给出满足不同误差条件的不同采样间隔所对应的最大信噪比估值，可作为实际应用中的参考。

平面近场声全息；信噪比；估值误差

0 引 言

近年来，近场声全息技术[1–2]（NAH）在噪声识别与定位、振动和声辐射特性研究等领域发展迅速，尤其是 Williams等[3]在 20 世纪 80 年代初提出的基于 FFT的平面 NAH 技术，因其理论简单，测量方便，易于计算和分析，在工程实际中应用广泛。

在应用平面 NAH 技术进行声场重建时，由于背景噪声分布在整个波数域范围内，逆向重构时位于高波数区域内的背景噪声将会被急剧放大，导致全息重构结果出现较大误差。因此，需要对全息面复声压数据进行波数域滤波[4]。目前常用的滤波窗函数主要有维纳（Wiener）滤波窗[5]以及带约束条件的最小二乘滤波窗[6]。在上述 2 种窗函数中，信噪比都是重要参数，能够直接决定滤波截止波数的选定，严重影响波数域滤波和声场全息重构效果。

但是，信噪比无法在测量过程中直接得到，需要根据全息面测量数据进行估计。辛雨等[7]提出以辐射圆外的波数域成分数据的均方根值近似代替全息面的噪声信号均方根值，采用能量法进行信噪比估计，但其估值严重依赖全息面数据的选取，而且无法判定估值是否准确。本文在其研究基础上给出波数域内的信噪比估值公式，提出一种判别估值是否精确的方法，并指出可以减小全息面采样间隔来降低估值误差。本文以 2 个点声源为例进行仿真，给出在不同采样间隔下满足不同误差条件的最大信噪比估值，可作为实际应用中的参考。

1 平面 NAH 基本原理[8]

设全息面测量得到的复声压数据为 pH(x,y,zH)，则有：

式中： Ps(x′,y′,Zs) 为声源面声压； gD(x,y,z) 为 Dirichlet 边界条件下的格林函数。

定义二维傅里叶变换为：

对式（1）两边取二维傅里叶变换，并由卷积定理可得：

式中 GD(kx,ky,z) 为 gD(x,y,z) 的二维傅里叶变换。其表达式为：

当k2x+k2y≤k2时，

当k2x+k2y＞k2时，

式中kx,ky分别对应 x、y 方向的空间波数。

将式（5）代入式（4），可得：

2 信噪比定义与噪声能量计算

信噪比即 SNR（Signal to Noise Ratio）是指信号中有用信号和噪声信号的比值，其大小通常用分贝数表示。信噪比可以使用信号或噪声的有效功率或能量定义，也可以使用信号和噪声声压的幅值定义。本文使用能量定义信噪比，公式为：

式中：ES和 EN分别为全息面信号和噪声的能量；E 为全息面总能量。

由式（10）可知，只要获得全息面噪声信号的能量，就可以估计出全息面信噪比。记全息面实际测量声压信号为 p( x,y,zH)， 噪声信号为 pN(x,y,zH)，且噪声信号为平稳随机分布，其均值为 0。M，N 分别为全息面上 x，y 方向的总测量点数，分别表示数学期望和方差，有

由式（11）可得：

式（12）本质上是噪声信号在空间域内的能量计算公式，只要估算出噪声信号在空间域内的方差，就可以估计出全息面内噪声的能量。但是全息面记录的声压信号是理论声压信号和噪声信号的混合，无法将噪声信号分离出来估计其方差。所以需要在波数域中进行讨论。根据采样定理可知，波数，记 kr=将波数域划分为 3 个区域，如图 1 所示。

其中 Ω1对应 kr＜k 的波数区域，即辐射圆以内的波数成分，记 GD为重建面到全息面的传递算子，则Ω1区域内 GD幅值恒为 1，对应在传播过程中幅值不随距离衰减的“平面波”，该区域通常包含较大能量； Ω2区域对应 kr＜k＜kmax的波数成分，此时 GD为一负指数函数，对应在传播过程中幅值随距离衰减的“倏逝波”； Ω3区域对应 kmax＜kr的波数成分，同样为“倏逝波”，而且具有更高的波数成分，在传播过程中衰减的更为迅速。一般而言， Ω3区域内的声压信号几乎全被噪声信号所淹没，无法再应用于重构，因此可以利用 Ω3区域内的波数成分估计全息面噪声信号的能量。

考虑到噪声信号在空间域内是平稳随机分布的，则其在波数域内也必然是平稳随机分布的，因此全息面噪声信号在整个波数域内的方差可以用 Ω3区域内的方差近似估计。即式中： PN(kx,ky,zH) 为 pN(x,y,zH) 的二维傅里叶变换；PN(Ω3) 为全息面声压角谱中位于 Ω3区域内的波数成分。而由 Parseval（巴塞伐）定理知，噪声信号在空间域中的总能量等于其波数域内的总能量，即

式中 ‖•‖2为矩阵的 2–范数； b 为 Ω3区域内的总测量点数。由式（13）和式（14）可知，全息面中噪声信号的总能量 E 可以估计为：

全息面测量声压总能量为：

将式（16）代入式（10）可得：

3 信噪比估值判别

由于噪声信号平稳随机分布，当全息面的位置和大小确定时，实际的全息面信噪比也随之确定。为便于仿真分析，会向全息面声压数据中添加一定信噪比下的噪声。以添加的信噪比值作为理论信噪比值，研究利用式（17）进行信噪比估计的估值误差。

假设某一全息面下理论信噪比值为 SNR，信噪比估值为 SNR’，记估值误差为 ε，其表达式如下

由式可知， ε 出现的根本原因是将 Ω3区域内的波数成分全部视为噪声信号的波数成分，而事实上该区域内的波数成分依然是声压信号和噪声信号波数的混合。因此， ε 的 大小与 Ω3区域内的声压信号衰减程度密切相关。从图 1 可以看出，kmax越大， Ω3区域内的有用声压信号衰减程度就越高，ε 值就越小，估值也就越准确。由采样定理，波数域内不发生混迭的最高波数成分 kmax=π/Δ ， Δ 越小，估值越准确。如果使采样间隔 Δ 不变，只在一定范围内改变全息面的大小，则信噪比估值应该近似。但是在实际测量中，为节约成本，一般只测量 1 次，所以，可以选取该全息面内不同区域的数据进行重复估计。

下面以 2 个点声源为例进行仿真：2 点源中心坐标分别取（0, 0.25, 0）和（0, –0.25, 0），点源强度为0.001 m3/s；分析频率 f = 1 000 Hz，声速取 340 m/s，分析波长 λ =0.34m ；取 Δ =0.0625m ， L=2m，全息面 ZH= 0.07 m，以全息面坐标原点为中心，从全息面中选取 9 个正方形区域 Li×Li(i=1,2···9)，其中 Li由下式确定：

由式（9）可得 L9=1m，并得到 9 个信噪比估值SNR′i(i=1,2···9)， 分别找到其最大估值 SN Rmax和最小估值 SN Rmin，仿真结果如图 2 所示。

由图 2 可发现，当理论信噪比小于 25 dB 时，理论值和估值极为接近，而且此时 SN Rmax和 SN Rmin也非常相近；当理论信噪比大于 25 dB 时，理论值与估值的差值越来越大，估值误差迅速增加。同时，SNRmax和 SN Rmin的距离也逐渐增大。在理论信噪比为 25 dB 时， SNRmax为 24.62 dB， SNRmin为 24.07 dB，估值误差 ε 仅为 1.52%。假定对估值误差的要求为εr， 令 εr=3%，则可认为，在 25 dB 以下的理论信噪比，其估值都精确。也就是说，在 εr=3% 的条件下，采样间隔 Δ =0.0625m 所对应的最大信噪比估值是 25 dB。

事实上，理论信噪比未知，只能得到估计值SNRmax和 SN Rmin，因此需要对估值误差 ε 进行修正。记修正后的估值误差为 ε′，其表达式为：

如果 ε′＜εr，就认为估值精确。理论信噪比为 25 dB时，有 ε′=2.35%＜3%，满足精度要求。

在同样仿真条件下，修正前后的误差对比如图 3所示。

由图 3 可看出，修正后的估值误差 ε′不会随着信噪比的增加而一直增大，而是分为 3 个阶段，第 1 阶段信噪比区间为 10～25 dB， ε′在 3% 左右波动；第 2阶段信噪比区间为 25～35 dB， ε′逐渐增大到 10% 左右；第 3 阶段信噪比区间为 30～50 dB， ε′基本稳定在10%。

在对信噪比估值进行判别时，如果取 εr=15%，可以发现在理论信噪比为 50 dB 时，信噪比估值只有32 dB，真实的估值误差 ε 高达 36%，但修正后的估值误差 ε′只有 10%，这显然不合理。所以，在进行估值判别时要合理选择 εr。考虑到信噪比越高，噪声干扰越小，对全息重构的影响也就越小， εr可以随着信噪比的增加而适当增加。通常情况下，信噪比小于 25 dB时，一般要求 εr在 3% 左右，最高不超过 10%。

综合以上分析，信噪比估值的判别方法如下：

1）保持采样间隔 Δ 不变，以全息面坐标原点为中心，选择合适的 n 值（可以令 Ln=L/2），从全息面中选择 n 个正方形区域 Li×Li(i=1,2···n)；

2）对所取区域的声压数据分别进行信噪比估计，得到一系列估值 SN R′i，找到其最大估值 SN Rmax和最小估值 SN Rmin，并计算 ε′；

3）根据实际测量情况，选择合适的 εr，如果ε′＜εr， 认为估值精确，并将 SN Rmax作为理论信噪比值。

4 减小信噪比估值误差的方法

如果判别后发现信噪比估值不准确，可以通过减少全息面采样间隔 Δ 来提高估值的准确性以降低估值误差。但减小全息面采样间隔进行第 2 次测试时，全息面大小能够保持不变，但是其位置可能会发生细微变化。因此选择 2 个不同位置的全息面，在不同采样间隔下进行仿真，结果如图 4 所示。

从图 4 可看出，采样间隔一定时， ε 值随着信噪比的增加急剧增加；而信噪比一定时，尤其是在高信噪比区域， ε 值随着采样间隔的减小迅速降低。而且，在 Δ ≤0.0625m 时，只要 Δ 不变，即使全息面的位置发生变化， ε 值也几乎不变，这为测试工作提供了方便。

实际测试中，如果全息面大小不变，减小采样间隔就意味着测量点数会增加，相应的工作量以及测量

成本也会随之增加。为方便测量工作，表 1 给出了 εr不同时，不同采样间隔所对应的最大信噪比估值。实际测量中，可以根据测试环境、声源可能的分布状况等综合考虑，参照表 1，选择合适的全息面采样间隔，在使信噪比估值误差满足要求时尽可能减少测量点数。

表 1 不同采样间隔下的最大信噪比估值Tab. 1 The maximum SNR estimation value in different sampling interval

5 结 语

本文通过分析波数域内噪声信号的分布情况，提出波数域内的全息面信噪比估计公式，并根据该估值公式分析了信噪比估值误差的来源，指出全息面采样间隔对估值误差具有重要影响。提出了一种选择不同区域的全息数据进行重复估计的方法对信噪比估值进行判断，并指出可以通过减小采样间隔来减降低信噪比估值误差。以两点源为例进行仿真，给出了在不同采样间隔下满足不同精度条件的最大信噪比估值，为实际应用提供参考。

[1]朱海潮.近场声全息技术的新进展[J].海军工程大学学报, 2013,25(6):81–87. ZHU Hai-chao.New progress in near-field acoustic holography [J].Journal of Naval University of Engineering,2013,25(6): 81–87.

[2]陈心昭, 毕传兴. 近场声全息技术及其应用[M]. 北京: 科学出版社, 2013.

[3]WILLIAMS E.G, MAYNARD J.D, SKUDRZYK E. Sound reconstruction using a microphone array [J]. Acoust. Soc. Am, 1980, 68(1): 340–344.

[4]李卫兵, 陈剑, 毕传兴, 等. 波数域滤波迭代近场声全息[J]. 机械工程学报,2004, 40(12): 14–19. LI Wei-bin,CHEN Jian,BI Chuan-xing,et al. Iterative near-field acoustic holography with k-space filter [J] .Chinese Journal of Mechanical Engineering,2004,40(12):14–19.

[5]ZHANG De-jun. Underwater acoustical holographic imaging by a square array system[J]. Acoustical Imaging, 1986, (15): 647–658.

[6]张德俊. 近场声全息对振动体及其辐射的成像[J]. 物理学进展, 1996, 16(3/4): 614–623. ZHANG De-jun. Imaging for vibration mode and radiation field of vibrating object useing NAH[J]. Progress in Physics, 1996, 16(3/4): 614–623.

[7]辛雨, 张永斌, 毕传兴, 等. 基于空间傅里叶变换的平面近场声全息中信噪比估计方法研究[J]. 计量学报, 2010, 31(6): 537–542. XIN Yu, ZHANG Yong-bin, BI Chuan-xing, et al. The estimation of the signal-to-noise in the planar nearfield acoustic holography based on the spatial fourier transform[J]. ACAT Metrologica Sinica, 2010, 31(6): 537–542.

[8]于飞, 陈剑, 周广林, 等. 噪声源识别的近场声全息方法和数值仿真分析[J]. 振动工程学报, 2003, 16(3): 339–343. YU Fei,CHEN Jian,ZHOU Guang-lin,et al. A study on noise sources identification useing NAH method and its numerical simulation analysis[J]. Journal of Vibration Engineering, 2003, 16(3): 339–343.

Estimation of the signal to noise ration in the planar acoustical holography based on energy method

DAI Zhen, HE Qi-wei, WAN Hai-bo
(Naval Engineering University, Wuhan 430033, China)

The signal to noise ratio (SNR) is a very important parameter in theplanar acoustical holography. However, it can’t be determined directly. To solve the problem, a formula in ‘-space’ is given based on theenergy method. According to the formula, the sampling interval has a great influence on the SNR estimation value. To determine the accuracy of the estimation value, the measurement data on the different areas of the holography surface is chosen to estimate. When the result is not accurate, the estimation value can be improved by decreasing the sampling interval. The maximum estimation value is given in different estimation precision and different sampling interval to provide a reference for practical applications.

planar near-field acoustic holography；signal to noise ratio；reconstruction error

O322

A

1672–7619(2017)03–0150–05

10.3404/j.issn.1672–7619.2017.03.031

2016–05–27；

2016–06–17

代振(1991–)，男，硕士研究生，研究方向为近场声全息。