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.     

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.]     

Surface decomposition method for near-field acoustic holography

Near-field acoustic holography reconstruction of the acoustic field at the surface of an arbitrarily shaped radiating structure from pressure measurements at a nearby conformal surface is obtained from the solution of a boundary integral equation. This integral equation is discretized using the equivalent source method and transformed into a matrix system that can be solved using iterative regularization methods that counteract the effect of noise on the measurements. This work considers the case when the resultant matrix system is so large that it cannot be explicitly formed and iterative methods of solution cannot be directly implemented. In this case the method of surface decomposition is proposed, where the measurement surface is divided into smaller nonoverlapping subsurfaces. Each subsurface is used to form a smaller matrix system that is solved and the result joined together to generate a global solution to the original matrix system. Numerically generated data are used to study the use of subsurface extensions to increase the continuity of the global solution, and investigate the size of the subsurfaces, as well as the distance between the measurement and the vibrating surface. Finally a vibrating ship hull structure is considered as a physical example to apply and validate the proposed methodology.     

Real-time near-field acoustic holography for continuously visualizing nonstationary acoustic fields

Near-field acoustic holography (NAH) is an effective tool for visualizing acoustic sources from pressure measurements made in the near-field of sources using a microphone array. The method involving the Fourier transform and some processing in the frequency-wavenumber domain is suitable for the study of stationary acoustic sources, providing an image of the spatial acoustic field for one frequency. When the behavior of acoustic sources fluctuates in time, NAH may not be used. Unlike time domain holography or transient method, the method proposed in the paper needs no transformation in the frequency domain or any assumption about local stationary properties. It is based on a time formulation of forward sound prediction or backward sound radiation in the time-wavenumber domain. The propagation is described by an analytic impulse response used to define a digital filter. The implementation of one filter in forward propagation and its inverse to recover the acoustic field on the source plane implies by simulations that real-time NAH is viable. Since a numerical filter is used rather than a Fourier transform of the time-signal, the emission on a point of the source may be rebuilt continuously and used for other post-processing applications.     

基于支持向量回归的patch近场声全息研究

为了解决patch近场声全息中全息数据外推问题,提出一种基于支持向量回归的一步式patch近场声全息技术。该方法首先把初始全息面上的数据当成训练样本进行学习,构造出回归函数,然后利用回归函数实现全息数据外推,最后基于统计最优近场声全息进行重建。数值仿真和实验研究的结果表明:在各个分析频率下,该方法都可以实现小孔径全息面的近场外推。从近场声全息重建结果看,即使初始全息数据受到噪声干扰,该方法也是一种有效的patch近场声全息技术。     

提高成像分辨率的远场多孔声全息方法

针对远场声全息成像分辨率较差的问题,提出多孔声全息方法。利用多个小型传声器阵列进行优化布置,可以有效地提高分辨率,达到等效大阵列的声场可视化效果。首先推导了空间多孔加窗对远场Kirchhoff衍射积分声全息带来的干涉叠加效应;然后对主瓣宽度和旁瓣高度进行了定量计算,阐述了多孔优化布置方法,实现了旁瓣抑制条件下的分辨率最优;最后利用双音箱实验,验证了多孔声全息的可行性和理论分辨率,证明了多孔声全息方法的有效性。研究表明,多孔声全息的干涉叠加方法可以有效提高远场测量的重建分辨率,多孔声全息的优化设计方法可以实现旁瓣抑制条件下的分辨率最优。     

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.     

声场分离技术及其在近场声全息中的应用

提出空间声场分离技术，突破了近场声全息（NAH）的应用局限.它们的局限在于全息面一侧的声场必须是自由声场，即要求所有的声源必须位于另一侧.利用波数域内的波场外推理论及声压的标量叠加原理，建立起声场分离技术的双全息面实现方法，利用波数域内的Euler公式及粒子振速的矢量叠加原理，建立起该技术的单全息面实现方法.该技术的一个突出优点是在具有背景噪声的全息测量情况下, 可以消除背景噪声对全息变换结果的影响.理论的推导表明该技术方法的正确性，而仿真算例和实验则显示该技术的可行性和有效性. 关键词： 声全息 波数域 声场分离 背景噪声     

循环谱密度组合切片分析在近场声全息中的应用研究

为减小循环谱密度的计算量以及提高其特征提取的准确性,提出了组合切片法:取循环自相关切片的各峰值频率作循环频率分别计算循环谱密度切片,通过组合切片分析来确定声场信号的特征频率。然后采用循环谱密度取代功率谱密度作为平面近场声全息的重建物理量。实验结果表明,此方法可针对性的提取循环平稳声场的调制特性,全息重建的结果可准确反映噪声源的位置。     

Use of nonsingular boundary integral formulation for reducing errors due to near-field measurements in the boundary element method based near-field acoustic holography

In the conformal near-field acoustic holography (NAH) using the boundary element method (BEM), the transfer matrix relating the vibro-acoustic properties of source and field depends solely on the geometrical condition of the problem. This kind of NAH is known to be very powerful in dealing with the sources having irregular shaped boundaries. When the vibro-acoustic source field is reconstructed by using this conformal NAH, one tends to position the sensors as close as possible to the source surface in order to get rich information on the nonpropagating wave components. The conventional acoustic BEM based on the Kirchhoff-Helmholtz integral equation has the singularity problem in the close near field of the source surface. This problem stems from the singular kernel of the Green function of the boundary integral equation (BIE) and the singularity can influence the reconstruction accuracy greatly. In this paper, the nonsingular BIE is introduced to the NAH calculation and the holographic BIE is reformulated. The effectiveness of nonsingular BEM has been investigated for the reduction of reconstruction error. Through interior and exterior examples, it is shown that the resolution of predicted field pressure could be improved in the close near field by employing the nonsingular BIE. Because the BEM-based NAH inevitably requires the field pressure measured in the close proximity to the source surface, the present approach is recommended for improving the resolution of the reconstructed source field.     

振动体及其辐射场的近场声全息实验研究

本文介绍我们自行研制的空气中近场声全息成像实验系统。用这一系统对单个扬声器,双扬声器及振动圆板等声源及其辐射场进行了测量,获得一些有意义的结果.这些结果与理论分析和计算机模拟良好地符合,表明近场声全息可以作为研究振动及其辐射场关系的有效手段.     

Interior near-field acoustical holography in flight

In this paper boundary element methods (BEM) are mated with near-field acoustical holography (NAH) in order to determine the normal velocity over a large area of a fuselage of a turboprop airplane from a measurement of the pressure (hologram) on a concentric surface in the interior of the aircraft. This work represents the first time NAH has been applied in situ, in-flight. The normal fuselage velocity was successfully reconstructed at the blade passage frequency (BPF) of the propeller and its first two harmonics. This reconstructed velocity reveals structure-borne and airborne sound-transmission paths from the engine to the interior space.     

Study of the comparison of the methods of equivalent sources and boundary element methods for near-field acoustic holography

Boundary element methods (BEM) based near-field acoustic holography (NAH) has been used successfully in order to reconstruct the normal velocity on an arbitrarily shaped structure surface from measurements of the pressure field on a nearby conformal surface. An alternative approach for this reconstruction on a general structure utilizes the equivalent sources method (ESM). In ESM the acoustic field is represented by a set of point sources located over a surface that is close to the structure surface. This approach is attractive mainly for its simplicity of implementation and speed. In this work ESM as an approximation of BEM based NAH is studied and the necessary conditions for the successful application of this approach in NAH is discussed. A cylindrical fuselage surface excited by a point force as an example to validate the results is used.     

基于统计最优和波叠加的联合局部近场声全息

提出了一种基于统计最优近场声全息和波叠加法的联合局部近场声全息技术。首先利用两次统计最优近场声全息的声源定位结果来指导配置等效源,其后利用波叠加法进行局部声场重构。该技术适合于中低频声场的局部重建,计算快速,重建精度高;可以在测量数据有缺失的情况下重建声场。进行了脉动球声源模型的数值仿真,并在半消声室内对电机噪声源进行了实验,仿真实验都准确地重构了声源所辐射的外部声场。该技术可以重建任意类球形声源辐射的声场。     

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

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

Planar near-field acoustical holography in a moving medium

Near-field acoustical holography (NAH) is a well-established method to study acoustic radiation near a stationary sound source in a homogeneous, stationary medium. However, the current theory of NAH is not applicable to moving sound sources, such as automobiles and trains. In this paper, the inclusion of a moving medium (i.e., moving source and receiver) is introduced in the wave equation and a new set of equations for plannar NAH is developed. Equations are developed for the acoustic pressure, particle velocity, and intensity when mean flow is either parallel or perpendicular to the hologram plane. If the source and the measurement plane are moving at the same speed, the frequency Doppler effect is absent, but a wave number Doppler effect exists. This leads to errors when reconstructing the acoustic field both towards and away from the source using static NAH. To investigate these errors, a point source is studied analytically using planar NAH with flow in one direction. The effect of the medium moving parallel to the hologram plane is noted by a shift of the radiation circle in wave number space (k-space). A k-space Green's function and a k-space filter are developed that include the effects of the moving medium.     

Compensation for source nonstationarity in multireference,scan-based near-field acoustical holography

Multireference, scan-based near-field acoustical holography is a useful measurement tool that can be applied when an insufficient number of microphones is available to make measurements on a complete hologram surface simultaneously. The scan-based procedure can be used to construct a complete hologram by joining together subholograms captured using a relatively small, roving scan array and a fixed reference array. For the procedure to be successful, the source levels must remain stationary for the time taken to record the complete hologram; that is unlikely to be the case in practice, however. Usually, the reference signal levels measured during each scan differ from each other with the result that spatial noise is added to the hologram. A procedure to suppress the effects of source level, and hence reference level, variations is proposed here. The procedure is based on a formulation that explicitly features the acoustical transfer functions between the sources and both the reference and scanning, field microphones. When it is assumed that source level changes do not affect the sources' directivity, a nonstationarity compensation procedure can be derived that is based on measured transfer functions between the reference and field microphones. It has been verified both experimentally and in numerical simulations that the proposed procedure can help suppress spatially distributed noise caused by the type of source level nonstationarity that is characteristic of realistic sources.     

Optimized microphone deployment for near-field acoustic holography: To be, or not to be random, that is the question

Arrays with sparse and random sensor deployment are known to be capable of delivering high quality far-field images without grating lobes. This raises the question of whether or not this idea can be applied to near-field imaging as well. To answer this question that has not yet been widely investigated in previous research, numerical simulations are undertaken in this paper to optimize the microphone deployment for both far-field and near-field arrays with the latter being the main focus. In the simulation, a recently introduced near-field equivalent source imaging (NESI) technique is employed for the near-field imaging. Global optimization techniques including the simulated annealing (SA) algorithm and the intra-block Monte Carlo (IBMC) algorithm are exploited to find the optimal microphone position efficiently. The combined use of the SA and the IBMC algorithms enables efficient search for satisfactory deployment with excellent beam pattern and relatively uniform distribution of microphones. In the near-field optimization, a special kind of beam pattern and cost function definition is used for the multiple-input-multiple-output (MIMO) imaging problem. As indicated by the simulation results, random deployment of microphones is necessary to avoid grating lobes in far-field imaging. In the near-field simulation, all results suggest that the optimal near-field array is the uniform rectangular array (URA) and the random deployment presents no particular benefit in near-field imaging.     

Data extrapolation method for boundary element method-based near-field acoustical holography

A new data extrapolation method for boundary element method (BEM)-based near-field acoustical holography (NAH) is proposed to reduce an error of the reconstructed result obtained from the pressure measured on an aperture small compared with the structure. The finiteness of the measurement aperture is a serious impediment to actual large-scale implementation of NAH, because NAH requires the measurement of the pressure field over a complete surface of structure. To eliminate the requirement, the wave number-space data extrapolation method for fast Fourier transform (FFT)-based NAH has been proposed. In this paper, the extension of this data extrapolation method to BEM-based NAH is presented. The effectiveness of this method is demonstrated by experiments. The experiment results confirm that the reconstruction error is sufficiently suppressed by the proposed method.     

A fault diagnosis scheme of rolling element bearing based on near-field acoustic holography and gray level co-occurrence matrix

Vibration signal analysis is the most widely used technique in condition monitoring or fault diagnosis, whereas in some cases vibration-based diagnosis is restrained because of its contact measurement. Acoustic-based diagnosis (ABD) with non-contact measurement has received little attention, although sound field may contain abundant information related to fault pattern. A new scheme of ABD for rolling element bearing fault diagnosis based on near-field acoustic holography (NAH) and gray level co-occurrence matrix (GLCM) is presented in this paper. It focuses on applying the distribution information of sound field to bearing fault diagnosis. A series of rolling element bearings with different types of fault are experimentally studied. Sound fields and corresponding acoustic images in different bearing conditions are obtained by fast Fourier transform (FFT) based NAH. GLCM features are extracted for capturing fault pattern information underlying sound fields. The optimal feature subset selected by improved F-score is fed into multi-class support vector machine (SVM) for fault pattern identification. The feasibility and effectiveness of our proposed scheme is demonstrated on the good experimental results and the comparison with the traditional ABD method. Considering test cost, the quantized level and the number of GLCM features for each characteristic frequency is suggested to be 4 and 32, respectively, with the satisfactory accuracy rate 97.5%.