统计最优近场声全息重建精度和计算速度优化方法

分析了统计最优近场声全息(Statistically optimized nearfield acousticholography,简称SONAH)的重建过程发现采用SONAH重建噪声源表面法向振速时误差较大的原因在于正则化参数选取较小。在此基础上提出一种单元平面波优化选择方法,该方法保留了单元平面波中的全部传播波和部分倏逝波,去除了一些较高波数的倏逝波成分,保证了重建过程中正则化参数的准确选取。另外,采用单元平面波优化选择方法还可以降低SONAH中传递矩阵的阶数,从而解决SONAH的计算速度随着测量点数目的增加急剧变慢的问题。通过数值仿真和实验对所提出的单元平面波优化选择方法的有效性进行了验证,结果表明采用该方法后SONAH的计算速度和法向振速的重建精度都得到了较大提高。      

Optimal secondary source position in exterior spherical acoustical holophony

Exterior spherical acoustical holophony is a branch of spatial audio reproduction that deals with the rendering of a given free-field radiation pattern (the primary field) by using a compact spherical loudspeaker array (the secondary source). More precisely, the primary field is known on a spherical surface surrounding the primary and secondary sources and, since the acoustic fields are described in spherical coordinates, they are naturally subjected to spherical harmonic analysis. Besides, the inverse problem of deriving optimal driving signals from a known primary field is ill-posed because the secondary source cannot radiate high-order spherical harmonics efficiently, especially in the low-frequency range. As a consequence, a standard least-squares solution will overload the transducers if the primary field contains such harmonics. Here, this is avoided by discarding the strongly decaying spherical waves, which are identified through inspection of the radiation efficiency curves of the secondary source. However, such an unavoidable regularization procedure increases the least-squares error, which also depends on the position of the secondary source.This paper deals with the above-mentioned questions in the context of far-field directivity reproduction at low and medium frequencies. In particular, an optimal secondary source position is sought, which leads to the lowest reproduction error in the least-squares sense without overloading the transducers. In order to address this issue, a regularization quality factor is introduced to evaluate the amount of regularization required. It is shown that the optimal position improves significantly the holophonic reconstruction and maximizes the regularization quality factor (minimizes the amount of regularization), which is the main general contribution of this paper. Therefore, this factor can also be used as a cost function to obtain the optimal secondary source position.     

单层传声器阵列重建相干声场的近场声全息方法研究

针对传声器阵列两侧存在相干声源的非自由声场重建问题,提出基于球面谐波函数扩展近场声全息理论的相干声场重建方法。该方法在已知测量面两侧声源几何位置时,使用单层传声器阵列获取测量面处的声压分布,通过最小二乘法获得与目标声源和干扰噪声源响应对应的最优球波函数扩展项数和最优系数向量,结合测点位置的空间坐标进行声波分解,并分别重建出各声源在测量面上的声压分布。为了验证方法的有效性,分别给出了相干噪声源为球形声源和非球形声源的仿真验证,并在全消声室内对双扬声器产生的相干声场的重建进行了实验验证。结果表明:该方法对球形声源和非球形声源干扰下的声场重建都具有较好的效果,球形声源干扰下的重建精度更高。      

HELS法在循环平稳声场全息重建中的理论与实验研究

Helmholtz 方程最小二乘法利用一组球面波基函数拟合声源产生的声场,根据重建和实际声压的误差最小原则,利用最小二乘法确定基函数展开的项数以及对应的权重系数,该方法具有计算效率高和需要测点少的优点,在实际工程中有很大的实用性.Helmholtz 方程最小二乘法和其他近场声全息方法一样都是针对平稳声场,对非平稳声场的分析很少.对于实际工程中经常遇到的一类特殊非平稳声场——循环平稳声场,现有的技术多以单通道信号分析为主,其高阶统计量在故障诊断领域应用较广.分析了循环平稳声场中Helmholtz方程最小二乘 关键词： 声全息 循环平稳 Helmholtz 方程 球面波     

Theoretical and experimental investigation of holographic reconstruction and prediction of multi-source hybrid acoustic field

Recently, NAH has been developed as an efficient tool for noise identification, noise localization and acoustic field visualization. Compared with the traditional acoustic ra-diation calculation problem, the solving problem in the NAH is an inverse acoustic problem. By measuring partial acoustic field information, such as complex pressures or particle velocities on the hologram surface, but not the surface normal velocities of the vibrating body, the surface information can be reconstructed,…     

Spatial resolution limits for the reconstruction of acoustic source strength by inverse methods

One method for deducing the strength of an acoustic source distribution from measurement of the radiated field involves the inversion of the matrix of frequency response functions relating the field measurement points to the strengths of a number of point sources used to represent the source distribution. In practice, the frequency response function matrix to be inverted may very often be ill-conditioned. This ill-conditioning will also often result in an ill-posed problem and thus regularization algorithms are used to produce reasonable solutions. For this purpose, Tikhonov regularization has been applied, and generalized cross-validation (GCV) has been introduced as an effective method for determining the proper amount of regularization without prior knowledge of either the source distribution or the contaminating errors. In the present work, the emphasis is placed on the relationship between the spatial resolution of the reconstructed source distribution and the small singular values of the frequency response function matrix to be inverted. However, the use of Tikhonov regularization often suppresses the effect of small singular values and these are in turn often associated with high spatial frequencies of the source distribution. Thus, the process of regularization produces a useful estimate of the acoustic source strength distribution but with a limited spatial resolution. Furthermore, in the field of Fourier acoustics, the spatial resolution of the reconstructed source distribution is usually limited by the wavelength of the radiation. This paper expresses the relationship between estimation accuracy, spatial resolution, noise-level and source/sensor geometry, when a range of inverse sound radiation problems are regularised using Tikhonov regularization with GCV. The results presented form the basis of guidelines that enable the reconstruction of acoustic source strength with a resolution that is finer than the intrinsic half-wavelength limit.     

Hybrid near-field acoustic holography

Hybrid near-field acoustical holography (NAH) is developed for reconstructing acoustic radiation from an arbitrary object in a cost-effective manner. This hybrid NAH is derived from a modified Helmholtz equation least squares (HELS) formula that expands the acoustic pressure in terms of outgoing and incoming waves. The expansion coefficients are determined by solving an overdetermined linear system of equations obtained by matching the assumed-form solution to measured acoustic pressures through the least squares. Measurements are taken over a conformal surface around a source at close range so that the evanescent waves can be captured. Next, the modified HELS is utilized to regenerate as much acoustic pressures on the conformal surface as necessary and take them as input to the Helmholtz integral formulation implemented numerically by boundary element method (BEM). The acoustic pressures and normal velocities on the source surface are reconstructed by using a modified Tikhnov regularization (TR) with its regularization parameter determined by generalized cross validation (GCV) method. Results demonstrate that this hybrid NAH combines the advantages of HELS and inverse BEM. This is because a majority of the input data are regenerated but not measured, thus the efficiency of reconstruction is greatly enhanced. Meanwhile, the accuracy of reconstruction is ensured by the Helmholtz integral theory and modified TR together with GCV method, provided that HELS converges fast enough on the measurement surface. Numerical examples of reconstructing acoustic quantities on the surface of a simplified engine block are demonstrated. [Work supported by NSF.]     

Reconstruction of the unsteady rotating forces of fan’s blade from far-field sound pressure

The aim of this paper is to investigate an inverse method used to evaluate the unsteady rotating forces acting on the fluid by the fan’s blade. A simple model based on the tonal noise produced by an axial fan and validated with a directivity experience is used to derive a discretized form of the direct problem and to simulate acoustic pressures at known spatial positions in the radiated field. The inversion of this direct problem is ill-conditioned and requires a regularization technique to stabilize the solution for small perturbations in the measured acoustic pressures. The reconstruction reveals that the conditioning of the inverse problem depends on the aeroacoustic source and the far-field microphones number as well as on the studied frequency. Tikhonov regularization can provide an appropriate regularization parameter leading to an accurate reconstruction of imposed unsteady rotating forces even with the presence of measurement noise.     

声辐射模态和球谐波函数在球源声场重建中的比较

声辐射模态和球谐波函数在球源声场重建中均有应用,两者具有很多相似性,都具有相同的多极子辐射模式和"群组"分布特性。有观点认为球形声源的声辐射模态向量即为球谐波函数,但两者的内在关系一直缺乏有效的证明。为分析两种基函数的相似性及在球源声场重建领域的差异,在理论上建立了两者的内在联系,在应用上比较研究了基于两种基函数的声场重建方法的有效性和可靠性。数值算例结果表明,在中高频段时基于两种函数的声场重建效果相当,但在低频段尤其是存在干扰源的非自由场条件下的声场重建应用中,基于声辐射模态的声场重建方法能得到更高的重建精度。可见,声辐射模态具有较为明显的声场重建优越性,这为利用声辐射模态和球谐波函数进行球源声场重建提供了理论参考。      

Regularization for improving the deconvolution in real-time near-field acoustic holography

Near-field acoustic holography is a measuring process for locating and characterizing stationary sound sources from measurements made by a microphone array in the near-field of the acoustic source plane. A technique called real-time near-field acoustic holography (RT-NAH) has been introduced to extend this method in the case of nonstationary sources. This technique is based on a formulation which describes the propagation of time-dependent sound pressure signals on a forward plane using a convolution product with an impulse response in the time-wavenumber domain. Thus the backward propagation of the pressure field is obtained by deconvolution. Taking the evanescent waves into account in RT-NAH improves the spatial resolution of the solution but makes the deconvolution problem "ill-posed" and often yields inappropriate solutions. The purpose of this paper is to focus on solving this deconvolution problem. Two deconvolution methods are compared: one uses a singular value decomposition and a standard Tikhonov regularization and the other one is based on optimum Wiener filtering. A simulation involving monopoles driven by nonstationary signals demonstrates, by means of objective indicators, the accuracy of the time-dependent reconstructed sound field. The results highlight the advantage of using regularization and particularly in the presence of measurement noise.     

利用声辐射模态重构任意目标的散射声场

水下目标散射声场的重构可以作为水下目标散射特性的研究基础。本文主要利用声辐射模态对水下目标进行散射声场重构研究。首先,在借助声传递矩阵给出的任意结构声辐射模态的流体域求解方法基础上,通过理论证明了目标的散射声压与声辐射模态具有函数关系。其次,借助声场分布模态的概念,同时考虑到声场分布模态病态及声压测量易受噪声污染,提出基于声辐射模态的正则化散射声场重构算法。仿真结果表明,波数越低,重构所需声辐射模态阶数越少,在较高波数时仅需总模态数的大约20%即可对声场进行重构。与基于边界元的声场重构算法相比,计算量减小了至少80%,且克服了赫姆霍兹积分方程最小二乘法仅对球壳结构的重构效果较好而不适用于长条形结构重构的缺陷。     

Visualization of acoustic radiation from a vibrating bowling ball.

This paper presents visualization of acoustic radiation from a vibrating bowling ball using the Helmholtz equation least squares (HELS) method. In conducting the experiments, the ball is excited by a vibration shaker using stationary random signals. The radiated acoustic pressures are measured using two microphones and taken as input to the HELS formulations. The reconstructed acoustic pressures on the bowling ball surface are compared with those measured at the same locations. Also shown are comparisons of the reconstructed and measured acoustic pressure spectra at various locations on the bowling ball surface. Results demonstrate that the accuracy of reconstruction based on measurements over a conformal surface is much higher than that over a finite planar surface. This is because the latter often extends beyond the near-field region, making the accuracy of measurements inconsistent. Nevertheless, satisfactory reconstruction of acoustic pressure fields over the entire bowling ball surface can still be obtained based on the measurements taken over a finite planar surface on one side of the source. In a similar manner, the normal component of the surface velocity is reconstructed. Once these acoustic quantities are determined, the time-averaged acoustic intensity is calculated. Also presented are the formulations for estimating a priori the numbers of expansion functions and measurements required by the HELS method and the guidelines for determining the reconstruction error and optimum measurement locations, given the overall dimensions of the source and the highest frequency of interest in reconstruction.     

Krylov subspace iterative methods for boundary element method based near-field acoustic holography

The reconstruction of the acoustic field for general surfaces is obtained from the solution of a matrix system that results from a boundary integral equation discretized using boundary element methods. The solution to the resultant matrix system is obtained using iterative regularization methods that counteract the effect of noise on the measurements. These methods will not require the calculation of the singular value decomposition, which can be expensive when the matrix system is considerably large. Krylov subspace methods are iterative methods that have the phenomena known as "semi-convergence," i.e., the optimal regularization solution is obtained after a few iterations. If the iteration is not stopped, the method converges to a solution that generally is totally corrupted by errors on the measurements. For these methods the number of iterations play the role of the regularization parameter. We will focus our attention to the study of the regularizing properties from the Krylov subspace methods like conjugate gradients, least squares QR and the recently proposed Hybrid method. A discussion and comparison of the available stopping rules will be included. A vibrating plate is considered as an example to validate our results.     

Combined Helmholtz equation-least squares method for reconstructing acoustic radiation from arbitrarily shaped objects

A combined Helmholtz equation-least squares (CHELS) method is developed for reconstructing acoustic radiation from an arbitrary object. This method combines the advantages of both the HELS method and the Helmholtz integral theory based near-field acoustic holography (NAH). As such it allows for reconstruction of the acoustic field radiated from an arbitrary object with relatively few measurements, thus significantly enhancing the reconstruction efficiency. The first step in the CHELS method is to establish the HELS formulations based on a finite number of acoustic pressure measurements taken on or beyond a hypothetical spherical surface that encloses the object under consideration. Next enough field acoustic pressures are generated using the HELS formulations and taken as the input to the Helmholtz integral formulations implemented through the boundary element method (BEM). The acoustic pressure and normal component of the velocity at the discretized nodes on the surface are then determined by solving two matrix equations using singular value decomposition (SVD) and regularization techniques. Also presented are in-depth analyses of the advantages and limitations of the CHELS method. Examples of reconstructing acoustic radiation from separable and nonseparable surfaces are demonstrated.     

Extracting the Green's function from measurements of the energy flux

Existing methods for Green's function extraction give the Green's function from the correlation of field fluctuations recorded at those points. In this work it is shown that the Green's function for acoustic waves can be retrieved from measurements of the integrated energy flux through a closed surface taken from three experiments where two time-harmonic sources first operate separately, and then simultaneously. This makes it possible to infer the Green's function in acoustics from measurements of the energy flux through an arbitrary closed surface surrounding both sources. The theory is also applicable to quantum mechanics where the Green's function can be retrieved from measurement of the flux of scattered particles through a closed surface.     

Modified adaptive algebraic tomographic reconstruction of gas distribution from incomplete projection by a two-wavelength absorption scheme

A modified adaptive algebraic reconstruction technique (MAART) with an auto-adjustment relaxation parameter and smoothness regularization is developed to reveal the tomographic reconstruction of H 2 O distribution in combustion from incomplete projections. Determinations of relaxation parameter and regularization factor are discussed with regard to the consideration of improvement in reconstruction and reduction of computational burden. A two-wavelength scheme from tunable diode laser absorption measurement, 7444.36 and 7185.59 cm 1 , is used to simplify the nonlinear solution problem for obtaining the tomographic distributions of concentration and temperature simultaneously. Results of calculations demonstrate that MAART can perform the reconstruction results more efficiently with little complex modification and much lower iterations as compared with the traditional algebraic reconstruction technique (ART) or simultaneous iterative reconstruction technique (SIRT) methods. The stability of the algorithm is studied by reconstruction from projections with random noise at different levels to demonstrate the dependence of reconstruction results on precise line-of-sight measurements.     

Effect of sensor proximity over the non-conformal hologram plane in the near-field acoustical holography based on the inverse boundary element method

Near-field acoustical holography (NAH) is a useful tool for the identification and visualization of vibro-acoustic sources. In particular, NAH can be applied to many practical sources having irregular shape if the inverse boundary element method (BEM) is employed. Once the relation between the source and the radiated field is defined in the transfer matrix modeled by the BEM, the reconstruction of acoustic parameters on the source surface can be conducted by the multiplication of the inverse transfer matrix and the field data measured over the hologram surface. The usual practical way to measure the field data radiated from an irregular shaped source is to adopt a regular hologram surface, which can be a flat, cylindrical, or spherical shape, for the measurement ease. Then, the hologram surface is not conformal to the source surface and the resulting transfer matrix becomes further ill-posed than the conformal case. To investigate the effect of sensor proximity and distance variability on the reconstruction error, simulation and measurement were conducted for an interior problem comprised of a parallelepiped rigid box with a vibrating end plate. Flat, tilted, and randomized hologram surfaces were adopted in the test. It was shown that the reconstruction error is greatly affected by the conditioning of transfer matrix which is related to the positioning of sensor. To improve the reconstruction accuracy, an investigation was conducted to find proper field points among a large number of overdetermined field points. The number of field data was reduced gradually under various reduction schemes using condition number, effective independence value, and sensor distance. It was demonstrated that the quality of reconstruction result given by the non-conformal measurement can be improved by removing some field points that contribute to the ill-conditioning of the inverse problem. A small improvement of the reconstruction accuracy was observed by reducing the field points in the overdetermined situation. However, further reduction of the field points, becoming an underdetermined situation, yielded a drastic improvement of the reconstruction accuracy.     

Regularization methods for near-field acoustical holography.

The reconstruction of the pressure and normal surface velocity provided by near-field acoustical holography (NAH) from pressure measurements made near a vibrating structure is a linear, ill-posed inverse problem due to the existence of strongly decaying, evanescentlike waves. Regularization provides a technique of overcoming the ill-posedness and generates a solution to the linear problem in an automated way. We present four robust methods for regularization; the standard Tikhonov procedure along with a novel improved version, Landweber iteration, and the conjugate gradient approach. Each of these approaches can be applied to all forms of interior or exterior NAH problems; planar, cylindrical, spherical, and conformal. We also study two parameter selection procedures, the Morozov discrepancy principle and the generalized cross validation, which are crucial to any regularization theory. In particular, we concentrate here on planar and cylindrical holography. These forms of NAH which rely on the discrete Fourier transform are important due to their popularity and to their tremendous computational speed. In order to use regularization theory for the separable geometry problems we reformulate the equations of planar, cylindrical, and spherical NAH into an eigenvalue problem. The resulting eigenvalues and eigenvectors couple easily to regularization theory, which can be incorporated into the NAH software with little sacrifice in computational speed. The resulting complete automation of the NAH algorithm for both separable and nonseparable geometries overcomes the last significant hurdle for NAH.     

三维复杂温度场高精度声学测量方法

声学温度场检测技术通过多路径声波传播时间数据,反演被测区域的温度分布。提供了一种高精度的三维复杂温度场的声学测量方法。首先从射线声学角度给出了三维非均匀温度场中声波传播路径的数学模型。在此基础上,将三维温度场的重建问题转化为声波传播路径的求解和温度场的反演问题,建立了基于多项式修正径向基函数(RBF-PR)和改进的Tikhonov正则化三维温度场重建算法。采用两种典型的炉膛三维温度场模型,在信噪比SNR=35 dB下进行了数值模拟,分析了声波传播路径在非均匀温度场中的弯曲特性、算法的重建质量和抗噪性,同时进行了实际炉膛内二维温度场的重建。结果表明了提出的考虑声线弯曲的温度场重建算法具有精度高,抗噪性强、适用性好的特点,为实现高精度的复杂温度场的声学测量提供了有效方法。      

半自由声场的全息重建和预测实验研究

半自由声场环境下的声源重建和声场预测研究对声全息技术走向实际应用具有非常重要的意义.在提出基于分布源边界点法的半自由声场全息重建和预测方法的基础上，对此展开了实验研究.并将重建和预测的结果与常规方法重建和预测的结果进行了比较和讨论，说明了重建预测过程中反射声压的影响和考虑反射声压的必要性，证明了所提出方法在解决半自由声场环境下存在地面反射时的声源重建和声场预测时的有效性和准确性.还提出了采用奇异值截断滤波和Tikhonov正则化方法来削弱测量误差的影响，从而进一步优化了重建结果，提高了全息成像的可信度. 关键词： 声全息 半自由场 边界点 声辐射 反射声