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.     

Compensation for discarded singular values in vibro-acoustic inverse methods

The paper addresses the inverse problem where source strengths are back-calculated from a sound pressure field sampled at several points. Regularization techniques, such as singular value discarding or Tikhonov regularization, are commonly used to improve estimates of source strength in such situations. However, over-regularization can result in even worse errors. A simple procedure is proposed here to compensate for errors of over-regularization. The basis is to constrain the solution such that the spatial mean of the measured and reconstructed sound pressure are equal. In other words, to set the overall sound power of the equivalent (calculated) sources equal to that of the real source. It is argued that the overall sound power is the most stable and reliable quantity on which to base source strength estimates. Examples of both singular value discarding and Tikhonov regularization are given.     

On the choice of expansion functions in the Helmholtz equation least-squares method

This paper examines the performance of Helmholtz equation least-squares (HELS) method in reconstructing acoustic radiation from an arbitrary source by using three different expansions, namely, localized spherical waves (LSW), distributed spherical waves (DSW), and distributed point sources (DPS), under the same set of measurements. The reconstructed acoustic pressures are validated against the benchmark data measured at the same locations as reconstruction points for frequencies up to 3275 Hz. Reconstruction is obtained by using Tikhonov regularization or its modification with the regularization parameter selected by error-free parameter-choice methods. The impact of the number of measurement points on the resultant reconstruction accuracy under different expansion functions is investigated. Results demonstrate that DSW leads to a better-conditioned transfer matrix, yields more accurate reconstruction than both LSW and DPS, and is not affected as much by the change in measurement points. Also, it is possible to obtain optimal locations of the auxiliary sources for DSW, LSW, and DPS by taking an independent layer of measurements. Use of these auxiliary sources and an optimal combination of regularization and error-free parameter choice methods can yield a satisfactory reconstruction of acoustic quantities on the source surfaces as well as in the field in the most cost-effective manner.     

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

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

Acoustical inverse problems regularization: Direct definition of filter factors using Signal-to-Noise Ratio

Acoustic imaging aims at localization and characterization of sound sources using microphone arrays. In this paper a new regularization method for acoustic imaging by inverse approach is proposed. The method first relies on the singular value decomposition of the plant matrix and on the projection of the measured data on the corresponding singular vectors. In place of regularization using classical methods such as truncated singular value decomposition and Tikhonov regularization, the proposed method involves the direct definition of the filter factors on the basis of a thresholding operation, defined from the estimated measurement noise. The thresholding operation is achieved using modified filter functions. The originality of the approach is to propose the definition of a filter factor which provides more damping to the singular components dominated by noise than that given by the Tikhonov filter. This has the advantage of potentially simplifying the selection of the best regularization amount in inverse problems. Theoretical results show that this method is comparatively more accurate than Tikhonov regularization and truncated singular value decomposition.     

Sound source reconstruction using inverse boundary element calculations

Whereas standard boundary element calculations focus on the forward problem of computing the radiated acoustic field from a vibrating structure, the aim in this work is to reverse the process, i.e., to determine vibration from acoustic field data. This inverse problem is brought on a form suited for solution by means of an inverse boundary element method. Since the numerical treatment of the inverse source reconstruction results in a discrete ill-posed problem, regularization is imposed to avoid unstable solutions dominated by errors. In the present work the emphasis is on Tikhonov regularization and parameter-choice methods not requiring an error-norm estimate for choosing the right amount of regularization. Several parameter-choice strategies have been presented lately, but it still remains to be seen how well these can handle industrial applications with real measurement data. In the present work it is demonstrated that the L-curve criterion is robust with respect to the errors in a real measurement situation. In particular, it is shown that the L-curve criterion is superior to the more conventional generalized cross-validation (GCV) approach for the present tire noise studies.     

炉膛三维温度场重建中Tikhonov正则化和截断奇异值分解算法比较

在煤粉锅炉诊断中火焰辐射能图像扮演着越来越重要的角色, 通过电荷耦合器件(CCD)获得的辐射能图像可以重建出炉内火焰三维温度场, CCD 用于获取视场角内的辐射能图像. 温度场重建的矩阵方程是一个严重病态的方程, 本文使用两种算法(Tikhonov正则化算法和截断奇异值分解(TSVD)算法)来重建温度场. 应用广义交叉检验算法来选取正确的正则化参数. 数值模拟的环境为一个10 m×10 m×10 m的三维炉膛, 系统被划分为10×10×10的1000个网格, 每个网格单元都是边长为1 m的立方体. 在正问题求解所得到的CCD接受信号基础上加上不同随机误差以模拟测量时的CCD接受信号. 研究两种算法重建后的温度重建误差、两者的重建时间, 以及最高温度的重建效果. 初步的研究结果显示, 一般情况下基于Tikhonov算法重建的温度场比基于TSVD算法重建的温度场误差要小, 计算所需时间短, 最高温度重建更准确.     

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,…     

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.     

利用源强密度声辐射模态重建声场

为了利用声场中少量测点声压数据精确重建复杂结构的辐射声场,提出了源强密度声辐射模态分析理论和声场重建公式.在结构表面定义的空间上,利用以源强密度分布函数为参量的结构辐射声功率泛函表达式定义了一个线性自伴正辐射算子,该算子的特征函数为结构的源强密度声辐射模态.然后通过对矩形平板和带有半球帽的圆柱体的源强密度声辐射模态的分析,证明了源强密度声辐射模态具有空间滤波特性,并利用该性质建立了声场重建公式.球体仿真和平板实验验证了所提出的声场重建方法的可行性和稳健性.基于源强密度声辐射模态的声场重建方法简单,利用较少测点数据就可以获得较高的声场重建精度,特别适合于复杂结构的低频声场重建.     

A threshold for the use of Tikhonov regularization in inverse force determination

In the analysis of structure-borne sound from installed machinery, it is important to be able to estimate the operational forces. Assuming that their location is known, indirect approaches based on matrix inversion can be used to reconstruct the operational forces from a set of measured operational responses and corresponding matrix of frequency response functions. In common with many such inverse problems, matrix ill-conditioning can affect the reliability of the results. Methods such as pseudo-inversion of over-determined matrices, singular value rejection, and Tikhonov regularization have been used previously to overcome this and it has been found that Tikhonov regularization generally performs well in reducing the errors in the reconstructed forces. However, full-rank pseudo-inversion (unregularized solution) gives better results than Tikhonov regularization in some cases, particularly with low condition numbers. Since the need for regularization is greatest when the matrix is ill-conditioned, this suggests the introduction of a threshold above which Tikhonov regularization is used and below which pseudo-inversion is used. In this study, the extent to which response errors are amplified in the force estimates is considered and plotted against the matrix condition number. This allows a threshold condition number to be identified above which Tikhonov regularization gives improved results. It is found that the threshold is related not only to the condition number but also to the matrix dimensions including the extent of over-determination. A simple empirical formula is obtained for this threshold that is usable for matrices in a wide range of matrix dimensions.     

The quantification of structure-borne transmission paths by inverse methods. Part 2: Use of regularization techniques

The inversion of an ill-conditioned matrix of measured data lies at the heart of procedures for the quantification of structure-borne sources and transmission paths. In an earlier paper the use of over-determination, singular value decomposition and the rejection of small singular values was discussed. In the present paper alternative techniques for regularizing the matrix inversion are considered. Such techniques have been used in the field of digital image processing and more recently in relation to nearfield acoustic holography. The application to structure-borne sound transmission involves matrices, which vary much more with frequency and from one element to another. In this study Tikhonov regularization is used with the ordinary cross-validation method for selecting the regularization parameter. An iterative inversion technique is also studied. Here a form of cross-validation is developed allowing an optimum value of the iteration parameter to be selected. Simulations are carried out using a rectangular plate structure to assess the relative merits of these techniques. Experiments are also performed to validate the results. Both techniques are found to give considerably improved results compared to singular value rejection.     

联合波叠加法的全息理论与实验研究

当空间声场中同时存在多个相干声源时，运用常规近场声全息方法无法重建每个相干声源表面的声学信息，当然也无法预测每个声源单独产生的空间声场，相干声场的全息重建与预测已成为全息技术推广应用过程中亟待解决的问题.在提出联合波叠加法并将其应用于空间声场变换的基础上，对其进行了实验研究.通过对实际相干声场的全息重建与预测，验证了常规波叠加法在相干声场重建中的局限性、联合波叠加法在相干声场全息重建与预测过程的可行性和准确性，还研究了Tikhonov正则化方法在抑制声学逆问题的非适定性中的有效性和滤波系数的选择原则的可行性，以提高全息重建与预测的精度. 关键词： 近场声全息 联合波叠加 相干声场 Tikhonov正则化     

声学CT复杂温度场重建研究

为提高声学CT复杂温度场重建能力,提出一种利用Markov径向基函数逼近和Tikhonov正则化的温度场重建算法,简称MTR算法。该算法首先用Markov径向基函数的线性组合,逼近介质中的复杂声速场分布,然后利用介质中多路径声波传播时间和Tikhonov正则化法,求解声速场分布,进而利用声速与温度的关系获得温度分布。对单热点、三热点和五热点温度场模型进行了仿真重建,结果表明MTR算法热点定位精度高,重建误差小。开发了声学CT温度场重建实验系统,用电加热器在内装1200 kg大豆的实验粮仓中形成热点,MTR重建结果能正确反映热点位置,热点温度重建误差1.3%。可见,MTR算法复杂温度场重建能力强,可望用于实际储粮温度分布监测。      

The application of regularization technique based on partial optimization in the nearfield acoustic holography

The regularization technique for stabilizing the reconstruction based on the nearfield acoustic holography(NAH) was investigated on the basis of the equivalent source method.In order to obtain higher regularization effect,a regularization method based on the idea of partial optimization was proposed,which inherits the advantages of the Tikhonov and another regularization method—truncated singular value decomposition(TSVD).Through the numerical simulation,it is proved that the proposed method is stabler than the Tikhonov,and more precise than the TSVD.Finally the validity and the feasibility of the proposed method are demonstrated by an experiment carried out in a semi-anechoic room with two speakers.     

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

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

Estimation of electron fluxes in the Earth's geostationary orbit with Tikhonov regularization during a magnetically quiet period

Electron energy fluxes in the Earth's outer radiation belt are estimated using an inverse theory of the Tikhonov regularization based upon observations in geostationary orbits. Particle Detector (PD) experiment aboard a geostationary satellite GEO-KOMPSAT-2A (GK2A) at a geographic longitude of 128.2°E provided observations of electrons within a 150–2,400 keV energy range with an unprecedented energy resolution of ΔE/E in the range of 5–25%. Instrument response functions, calculated with Monte-Carlo simulations, are deconvoluted with electron observations. Using regularization parameters determined from Generalized Cross Validation (GCV), the Tikhonov method was applied to observations made during a geomagnetically quiet period. This Tikhonov regularization method, now possible for observations in the Earth's radiation belt for the first time, allows direct inference of electron fluxes without resorting to predetermined functional forms. Comparisons of our results with those from conventional methods indicate differences among the results as large as ∼200%.     

板结构的声辐射模态分析及其高分辨率法向振动速度重建

提出了一种基于声辐射模态的速度基向量构建方法,该速度基向量不受网格划分的影响,可用于高分辨率的板结构法向振动速度重建。首先对板表面稀疏网格的声辐射模态进行计算,再以声辐射模态和模态系数构建板法向振动速度分布的基向量,然后由声场测量声压求解基向量系数,最后由该系数和加密网格的速度基向量重建高分辨率的板法向振动速度分布。以简支板声源进行仿真计算,当测量声压信噪比为30 dB时,低频的法向振动速度重建误差最低可达3.7%;以固支板声源在消声室中进行实验验证,131.5 Hz振动频率下的重建误差低于7%。该方法实现了只需要少量声压测量点即可精确重建板声源更高分辨率的法向振动速度分布。      

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 vibroacoustic fields in half-space by using hybrid near-field acoustical holography

In this paper we examine the accuracy and efficiency of reconstructing the vibroacoustic quantities generated by a vibrating structure in half-space by using hybrid near-field acoustic holography (NAH) and modified Helmholtz equation least squares (HELS) formulations. In hybrid NAH, we combine modified HELS with an inverse boundary element method (IBEM) to reconstruct a vibroacoustic field. The main advantage of this approach is that the majority of the input data can be regenerated but not measured, thus the efficiency is greatly enhanced. In modified HELS, we expand the field acoustic pressure in terms of outgoing and incoming spherical waves and specify the corresponding expansion coefficients by solving a system of equations obtained by matching the assumed-form solution to the measured acoustic pressure. Here the system of equations is ill conditioned and Tikhonov regularization is implemented through singular value decomposition (SVD) and the generalized cross-validation (GCV) method. Numerical examples of a dilating and oscillating spheres and finite cylinder are demonstrated. Test results show that hybrid NAH can yield a more accurate reconstruction than does a modified HELS, but a modified HELS is more efficient than is hybrid NAH [Work supported by NSF].