A wave field synthesis approach to reproduction of spatially correlated sound fields

This article discusses an open-loop wave field synthesis (WFS) approach for the reproduction of spatially correlated sound fields. The main application concerns laboratory reproduction of turbulent boundary layer wall pressure on aircraft fuselages and measurement of their sound transmission loss. The problem configuration involves reconstruction of random sound pressure distributions on a planar reproduction surface using a planar array of reproduction monopoles parallel to the reproduction plane. In this paper, the WFS formulation is extended to sound fields with imposed time and spatial correlation properties (or equivalently imposed cross-spectral density in the frequency and wave number domains). Numerical examples are presented for the reproduction of a propagating plane wave, diffuse acoustic field and wall pressure in subsonic or supersonic turbulent boundary layers. The reproduction accuracy is examined in terms of the size of the source plane and reproduction plane, their separation, and the number of reproduction sources required per acoustic wavelength. While the reproduction approach cannot reconstruct sub-wavelength correlation scales of subsonic turbulent boundary layers, it effectively reconstructs correlation scales larger than the acoustic wavelength, making it appropriate for diffuse acoustic field and supersonic turbulent layers.     

Wave field synthesis of moving virtual sound sources with complex radiation properties

An approach to the synthesis of moving virtual sound sources with complex radiation properties in wave field synthesis is presented. The approach exploits the fact that any stationary sound source of finite spatial extent radiates spherical waves at sufficient distance. The angular dependency of the radiation properties of the source under consideration is reflected by the amplitude and phase distribution on the spherical wave fronts. The sound field emitted by a uniformly moving monopole source is derived and the far-field radiation properties of the complex virtual source under consideration are incorporated in order to derive a closed-form expression for the loudspeaker driving signal. The results are illustrated via numerical simulations of the synthesis of the sound field of a sample moving complex virtual source.     

Adaptive wave field synthesis for active sound field reproduction: experimental results

Sound field reproduction has applications in music reproduction, spatial audio, sound environment reproduction, and experimental acoustics. Sound field reproduction can be used to artificially reproduce the spatial character of natural hearing. The objective is then to reproduce a sound field in a real reproduction environment. Wave field synthesis (WFS) is a known open-loop technology which assumes that the reproduction environment is anechoic. The room response thus reduces the quality of the physical sound field reproduction by WFS. In recent research papers, adaptive wave field synthesis (AWFS) was defined as a potential solution to compensate for these quality reductions from which WFS objective performance suffers. In this paper, AWFS is experimentally investigated as an active sound field reproduction system with a limited number of reproduction error sensors to compensate for the response of the listening environment. Two digital signal processing algorithms for AWFS are used for comparison purposes, one of which is based on independent radiation mode control. AWFS performed propagating sound field reproduction better than WFS in three tested reproduction spaces (hemianechoic chamber, standard laboratory space, and reverberation chamber).     

Adaptive wave field synthesis for broadband active sound field reproduction: signal processing

Sound field reproduction is a physical approach to the reproduction of the natural spatial character of hearing. It is also useful in experimental acoustics and psychoacoustics. Wave field synthesis (WFS) is a known open-loop technology which assumes that the reproduction environment is anechoic. A real reflective reproduction space thus reduces the objective accuracy of WFS. Recently, adaptive wave field synthesis (AWFS) was defined as a combination of WFS and active compensation. AWFS is based on the minimization of reproduction errors and on the penalization of departure from the WFS solution. This paper focuses on signal processing for AWFS. A classical adaptive algorithm is modified for AWFS: filtered-reference least-mean-square. This modified algorithm and the classical equivalent leaky algorithm have similar convergence properties except that the WFS solution influences the adaptation rule of the modified algorithm. The paper also introduces signal processing for independent radiation mode control of AWFS on the basis of plant decoupling. Simulation results for AWFS are introduced for free-field and reflective spaces. The two algorithms effectively reproduce the sound field and compensate for the reproduction errors at the error sensors. The independent radiation mode control allows a more flexible tuning of the algorithm.     

水下目标对幅度高斯相关海面环境噪声场扰动特性的研究

理论探讨了水下目标受海面环境噪声场激发产生的散射场对海洋环境噪声场的扰动特性,基于幅度高斯型相关海面噪声源模型,利用任意声源分布声场的积分表达式和单极子源半无限空间格林函数的球面波展开式,获得了刚性球目标在海面单极子和偶极子源海洋环境噪声“照射”下的直达噪声场、散射噪声场、总噪声场,以及二接收点直达、散射和总噪声协方差的理论表达式,表明总噪声场除了和目标阻抗特性、接收点的方向有关外,还受到海面噪声源相关特性,以及直达与散射噪声场之间干涉的较大影响。数值计算结果给出较近距离范围内,刚球目标声学可见度约为4~5dB,并指出在海洋环境噪声场中,刚球目标散射的方向特性明显不同于平面波入射。      

Horizontal refraction of sound rays due to short-period internal waves in the ocean

The effect of a short-period internal wave measured in field conditions on the horizontal (side) refraction of sound rays is estimated. The angle of horizontal refraction, i.e., the angle between the direction of the signal arrival in the horizontal plane and the true direction to the sound source, is determined. The influence of various factors, such as the position of the receiving system in depth with respect to the layer of high sound velocity gradients, the rotation of the transmitter-receiver track with respect to the internal wave front, etc., on the horizontal refraction is estimated. Numerical calculations are carried out. Conclusions about the possible errors that arise in determining the azimuth direction to the sound source because of the effect of short-period internal waves are derived.     

Vertical amplitude phase structure of a low-frequency acoustic field in shallow water

We obtain in integral and analytic form the relations for calculating the amplitude and phase characteristics of an interference structure of orthogonal projections of the oscillation velocity vector in shallow water. For different frequencies and receiver depths, we numerically study the source depth dependences of the effective phase velocities of an equivalent plane wave, the orthogonal projections of the sound pressure phase gradient, and the projections of the oscillation velocity vector. We establish that at low frequencies in zones of interference maxima, independently of source depth, weakly varying effective phase velocity values are observed, which exceed the sound velocity in water by 5–12%. We show that the angles of arrival of the equivalent plane wave and the oscillation velocity vector in the general case differ; however, they virtually coincide in the zone of the interference maximum of the sound pressure under the condition that the horizontal projections of the oscillation velocity appreciably exceed the value of the vertical projection. We give recommendations on using the sound field characteristics in zones with maximum values for solving rangefinding and signal-detection problems.     

三维浅海环境下孤立子内波对低频声能流的传播影响

针对三维浅海环境下孤立子内波对低频声信号传播特性的影响问题,基于Oxyz坐标系下的三维浅海低频声场有限元计算方法,以声能流为研究对象,仿真分析了内波存在对低频声信号传播特性的影响规律.研究结果表明:受内波影响,在xOz平面,声能流垂直分量的传播偏转角度呈现周期性的起伏规律;随着声源深度的增加,内波对声能流偏转角度的影响...     

波场合成与波场分析的有源房间补偿方法

波场合成是一种空间声重放技术,利用扬声器阵列在宽阔的听音区域内重建声场。为消除或者消减听音房间的反射对重建波场的影响,利用测量环绕听音区域的闭合曲线上声压和声压梯度来分析波场,推导了基于圆形阵列进行波域分解的公式,利用多通道逆滤波进行了平面波域的房间补偿,实验结果显示该算法在整个听音区域内都是有效的。与传统补偿方法相比,边界元法所需的测量传声器数目少计算复杂度低,而波域分解具有更充分的波场分析能力,因此是一种更有效的有源房间补偿方法.      

Ultrasound fields in attenuating media

For medical ultrasonic imaging and for nondestructive testing, the attenuation of pressure waves and the resulting shift in wave velocity are important features in commonly used transmission media such as biological tissue. An algorithm for the numerical evaluation of pressure field distributions generated by ultrasonic transducers is presented. The attenuation and dispersion of the sound transmission medium are taken into consideration. The sound fields are computed numerically for continuous wave as well as pulse excitation. The transducer has plane or gently curved geometry and is embedded in a plane rigid baffle. The numerically determined pressure fields are presented as 3D plots, as gray-scale images for a fixed time stamp (like a snapshot), or as isobars regarding the maximum values over time for each local point in the area under investigation. The algorithm described here can be utilized as a tool for design of ultrasound transducers, especially array antennas.     

Low-frequency horizontal acoustic refraction caused by internal wave solitons in a shallow sea

The effect of internal wave solitons on the sound field generated by a point source in a shallow sea is considered. In the framework of the theory of “horizontal rays and vertical modes,” the sound field pattern governed by the aforementioned hydrodynamic effect is investigated. It is shown that solitons can induce time-periodic focusing and defocusing of horizontal rays propagating at shallow angles to the internal wave front. This may result in the formation of “dynamical” horizontal sound channels, which, in its turn, results in considerable temporal fluctuations of the field along the acoustic track oriented along the internal wave front. For the sound field calculations, an approach is developed on the basis of the parabolic approximation in the horizontal plane and the mode representation in the vertical direction. The results obtained can be used for remote monitoring of internal wave packets in a shallow sea.     

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

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

Simulation of the vector wave field of a low-frequency sound source in the ocean: Some important features

An approach to the simulation of low frequency vector wave fields in stratified media (mainly in the ocean) is considered. The approach is characterized by an improved stability with respect to dividing the medium into many layers of arbitrary thickness. The model for the sound field of a point source is based on an integral representation of two-dimensional, cylindrically symmetric vector wave fields in inhomogeneous media, so that the contributions of all types of waves are included automatically. The model medium is subdivided into N horizontally homogeneous layers for which 4(N?1) equations are formulated to satisfy the boundary conditions between adjacent layers. The method of the generalized Schmidt matrix is used to obtain the coefficients of the equations; these coefficients are substituted into the expressions (of the Fourier-Bessel integral type) for the local parameters of the field. The latter are calculated according to the numerical procedure, and the results are used to model the distributions of the acoustic pressure and the horizontal and vertical components of the particle velocity in liquid and elastic media. The instability of the calculation procedure may result in a disagreement between the model and the exact solution. However, the disagreement is shown to occur mainly in models containing excessively thick layers. A way for improving the stability of the numerical model is suggested. The simulation results are compared with the exact analytical solution for the simplest example and with the results obtained according to the commonly used generalized matrix procedure (the benchmark problem). The examples of the practical application of the model for investigating more complex seismoacoustic wave fields in the ocean are presented.     

采用振速测量的改进统计最优近场声全息方法

针对常规统计最优近场声全息在空间多源声场重建过程中所需波函数项数多、重建精度不理想的问题,本文提出了一种基于振速测量的改进统计最优近场声全息算法。与常规算法不同,改进算法主要根据声源特点选取合适的波函数组合来计算声场传递矩阵。通过数值仿真初步验证了该方法的准确性和有效性,并与常规算法进行了详细的对比分析。仿真结果表明,改进算法重建精度高,随频率的变化相对误差波动较小,且随着频率的增大相对误差有逐渐减小的趋势;此外,不同的波函数组合,重建效果差异很大,当选取的波函数与声源共形且数量一致时重建效果最好。     

采用球面波叠加法还原自由声场的方法

为解决非自由声场中近场声全息重建时,干扰声在目标声源表面产生的散射影响,提出一种基于球面波叠加法的自由场还原技术。该技术首先采用基于球面波叠加法的声场分离技术获得向内和向外传播的声场,然后以目标声源的表面导纳作为边界条件,实现目标声源辐射声和散射声的分离,从而获得等效于自由声场的测量条件。该技术为准确实现非自由声场中的噪声源识别创造了条件。文中首先详细描述了该技术的基本原理,并提出一种最优球面波展开项数选取方法,最后通过数值仿真说明了该技术的有效性。结果表明:在频率较低时,散射声影响较小,采用声场分离技术和自由场还原技术效果相当;但随着频率升高,散射声影响逐步增强,必须采用自由场还原技术才能准确获得目标声源辐射声。      

Recovery of the free field in noisy environment by using the spherical wave superposition method

To remove the scattering effect of the disturbing sound on the target source when implementing nearfield acoustic holography in a non-free field, a free field recovery technique based on the spherical wave superposition method is proposed. In the method, the sound field separation technique based on the spherical wave superposition method is first used to separate the incoming and outgoing fields, and a further step for separating the radiated and scattered fields is performed by utilizing the surface admittance of the target source as the boundary condition. The technique makes it possible to correctly identify noise sources in a non-free sound field. The basic principle of the technique is described firstly, a method for choosing the optimal number of spherical wave expansion terms is given, and two numerical simulations are used to demonstrate the validity of this technique. It is shown that, for the lower frequency, the scattering effect can be neglected, and the radiated field of the target source can be obtained by the sound field separation technique, however, as the increasing of the frequency, the scattering effect cannot be neglected, and the free field recovery technique has to be used to obtain the radiated field of the target source.     

迭代重加权时域平面波叠加法的非稳态声场重建

提出了迭代重加权时域平面波叠加法,并应用于非稳态声场重建。该方法首先通过时域离散卷积得到全息面时域声压与虚源面时域波数谱离散状态公式;然后根据稀疏表示和迭代重加权,得到残差范数与加权解范数最小化问题表达式;最后将每一对波数下虚源面时域波数谱与时域传播核卷积叠加得到重建面上的时域声压。为了验证该方法的有效性,给出了激励板声源的仿真验证,并在半消声室内进行了实验验证。讨论了一些参数对该方法的影响,并通过与基于Tikhonov正则化的时域平面波叠加法比较,突出该方法的重建优势。仿真和实验结果表明,该方法对非稳态声场重建具有很好的效果,并且重建精度高。     

一种考虑声线弯曲的温度场重建算法

声波在非均匀温度场中传播时会由于声波的折射而产生声线弯曲现象。为提高非均匀温度场声学CT (Computer Tomography)重建精度,提出一种考虑声线弯曲的重建算法。首先用最小二乘法获得一个不考虑声线弯曲的重建温度场,然后用打靶-插值法确定本征声线出射角,用三角形前向展开法追踪声线,获得此温度场中声波发射/接收器间的本征声线,建立本征声线上声波传播时间与温度分布间的关系,进而实现考虑声线弯曲的最小二乘法温度场重建。单峰、双峰温度场仿真重建结果表明:所采用的本征声线追踪法计算简单运行速度快;考虑声线弯曲后,可明显提高温度场重建精度。因此所提重建算法能更好地适应实际温度场重建需求。      

Spatial-temporal correlation of a diffuse sound field

The spatial correlation function of the sound in a diffuse field is a quantity widely used in many reverbrant room acoustic applications. Although results for the spatial and temporal correlation for pure-tone and band-limited diffuse fields have already been developed, these have not been generalized for other signal types. This work presents a generalized derivation of the diffuse field spatial-temporal correlation which can be used for stationary random signals with given power spectral density. It is shown that the spatial-temporal correlation depends entirely on the temporal correlation of the signal exiting the diffuse field, or alternatively on its power spectral density. A simulation using the plane wave model is presented for tonal and broadband diffuse sound fields.