The improved separation method of coherent sources with two measurement surfaces based on statistically optimized near-field acoustical holography

Jin Mao1 , Zhongming Xu2 , Zhifei Zhang3 , Yansong He4 , Shu Li5

1, 2, 3, 4, 5State Key Laboratory of Mechanical Transmission, Chongqing University, Chongqing, 400030, China

1Corresponding author

Journal of Vibroengineering, Vol. 17, Issue 2, 2015, p. 674-681.
Received 16 July 2014; received in revised form 13 September 2014; accepted 5 October 2014; published 31 March 2015

Copyright © 2015 JVE International Ltd.

Table of Contents Abstract Full-text  Download PDF References
Views 0
Reads 0
Downloads 21
since May, 2018

A technique for separating coherent sources measured by two parallel arrays is proposed. The two measurement surfaces located in the opposite directions of the coherent sources. Similar to separate the aim source from background noise, this method can separate the single source from coherent sources, which makes the sound field information of single source in complex environment more accurate. Such improved separation method based on statistically optimized near-field acoustical holography, according to the sound pressure relationship between measurement surfaces and reconstruction surfaces to separate the sources, reduces the measurement data and obtains higher precision of reconstruction. The present paper uses the improved separation method to obtain the single source results from numerical simulations, gives the relative reconstruction errors with frequency from 100 Hz to 1400 Hz, and practical measurement.

Keywords: statistically optimized, sound field separation, acoustical holography.

1. Introduction

Near-field acoustical holography (NAH) is a powerful tool that can reconstruct a three-dimensional (3-D) space sound field, i.e., sound pressure, particle velocity, and sound intensity, by propagating sound pressure data measured on a two-dimensional (2-D) measurement surface [1, 2]. Statistically optimized near-field acoustical holography (SONAH) has higher precision of NAH proposed by Hald to reduce the window effects associated with the spatial FFT used in conventional NAH [3-5].

A set of techniques about NAH have been developed to avoid the use of the discrete Fourier transform (DFT) [6], which more or less improve the computational efficiency and reduce the measurement data. There are several methods such as SONAH, inverse boundary element method (IBEM) [7, 8], and equivalent source method (ESM) [9, 10]. Partial field decomposition to identify the coherent sources by means of reference signal was proposed by J. Stuart et al., but it needed multi-reference microphones and to do continuous and more than once measurements [11]. Zhang et al. used orthogonal spherical wave put on maximum value of sound pressure over the sound source plane and double times calendar to confirm the source position, which would lose efficiency under the circumstance that interaction of multiple sources produced maximal values of pseudo sources [12]. Jerome made use of Bayesian theory to reconstruct the coherent sources and then obtained the source position by variable radius of a 2-D Hanning window, which also wasn’t suitable for the maximums of pseudo sources [13]. From the above discussion, those methods could all identify the coherent sources in some cases, but cannot be applied in any case.

Sound field separation technique with two surfaces makes a contribution to separate sound field information of the aim source from background noise [14-17], which applies a way to separate the coherent sources in the same direction of measurement surface. The present paper separates the coherent sources with two parallel measurement surfaces used improved method with SONAH, which would solve sources located at any position, reduce the measurement data and improve the computational time.

2. General underlying formulation

2.1. Theory of SONAH

Basic theory of SONAH: Assume that complex time-harmonic sound pressure pri has been measured in free and homogeneous region Ω by a set of positions ri, i=1, 2,…, N. The amplitude and phase of each elementary wave in the sound field is written as:

p r i = m = 1 M a m Φ m r i ,               i = 1 , , N .

Here, am are the complex expansion coefficients. When NM, the am can be calculated by the least-squares equivalent wave model. There are M wave functions, Φm, measured at N hologram locations.

Generally, the elementary waves are plane and evanescent, and they are considered to be a free-field condition. Sources behind a source plane, radiate into the source free, homogeneous half space. When sources are not restricted to a half space, as illustrated in Fig. 1, where measurements are taken in a free and homogeneous source field between two parallel source planes, z0+zz0-. In this space, the radiated sound field can be represented by plane wave functions of the combined set:

Φ k 1 r = e - j k x x + k y y + k z z - z 0 - ,
Φ k 2 r = e - j k x x + k y y + k z z - z 0 + ,

where z0- and z0+ are source positions in z axis, r=x,y,z is a position vector, and k=kx,ky,kz is a wave number vector with:

k z = k 2 - k x 2 + k y 2 ,                   k x 2 + k y 2 k 2 ,   - j k x 2 + k y 2 - k 2 ,         k x 2 + k y 2 > k 2 .

Matrix A is the formulation of wave function values at the measurement positions, defined as:

A = Φ 1 r 1 Φ 1 r 2         Φ 1 r N Φ 2 r 1 Φ 2 r 2         Φ 2 r N                                                                   Φ M r 1 Φ M r 2         Φ M r N .

Vector p and the matrix αr at reconstruction positions can be written by the following form:

p = p r i = p r 1 p r 2 p r N ,       α r = Φ 1 r s Φ 2 r s   Φ M r s ,

here, rs is a position vector at the reconstruction locations.

To obtain the pressure at the reconstruction location, prs, we use Tikhonov regularized solution to resolve the single value which would come out during the inverse operation of matrices. After a set of calculation, prs is given by:

p r s = p T A H A + ε I - 1 A H α r ,

where AH is the conjugate transpose of A, I is a unite diagonal matrix and ε is positive regularization parameter, provided by Hald [4], is given by:

ε = 1 + 1 2 k d 2 10 - S N R 10 ,

where SNR is the signal to noise ratio including random error and noise, and d is the distance between measurement surface and reconstruction surface.

Fig. 1. Two parallel source planes limiting a homogeneous and source-free region

Fig. 2. The positions of measurement surfaces, reconstruction surfaces and sources

2.2. Sound field separation technique based on SONAH

In order to obtain the information of a single source in an acoustic field where exist other sources, we should separate the target source from the complex environment. The method of the present paper is to measure the sound pressures on two spaced parallel measurement surfaces located at the opposite sides of the sources and then realize the separation by SONAH.

Assume that there are two source points in the acoustic field, since sound pressure is a scalar quantity which is measured in two parallel planes, it can be expressed as in Eq. (9) and Eq. (10) in each plane:

p 1 r 1 = p 11 r 1 + p 21 r 1 ,
p 2 r 2 = p 12 r 2 + p 22 r 2 .

The p1 and p2 refer to the sound pressure on two hologram surfaces, and r1 and r2 refer to the positions of sources on hologram surfaces. The pi1 and pi2 contain respectively the pressure produced from the source i to hologram surface H1 and H2.

Define the transfer vector between measurement surface H1 and reconstruction surface R1 as c11, the transfer vector between measurement surface H2 and reconstruction surface R2 as c22, the transfer vector between measurement surface H1 and reconstruction surface R2 as c12, and the transfer vector between measurement surface H2 and reconstruction surface R1 as c21. Then, the reconstruction proceedings respectively from measurement surface H1 to reconstruction surface R1and from measurement surface H2 to reconstruction surface R2 as in Eq. (11) and Eq. (12):

p 1 c 11 = p 11 c 11 + p 21 c 11 ,
p 2 c 22 = p 12 c 22 + p 22 c 22 .

According to the acoustic field reconstruction principle, it could be obtained that the two sound pressures, locating at the same distances away from two sides of one source point which propagates the spherical wave, are seen as the same size. During the calculation, the distance between H2 and R1 has been converted in terms of the H2 and R1 in the same direction of sound. In addition, the distance between H1 and R2 has the same process. Then, p11c11 could be replaced by p12c21, and also p22c22 could be replaced by p21c12. We can obtain another form of Eq. (11) and Eq. (12), respectively:

p 1 c 11 = p 12 c 21 + p 21 c 11 ,
p 2 c 22 = p 12 c 22 + p 21 c 12 .

This allows Eq. (13) and Eq. (14) to be written as follows:

  p 21 = p 1 c 11 c 0 - p 2 c 22 c 11 c 0 - c 12 ,

where c0=c21-1c22, the c21-1 refers to the generalized inverse. c0 could be worked out by the Tikhonov regularization method according to the references [5] and [18], then the regularized least squares solution as the Eq. (7) is:

c 0 = c 21 + c 21 + ε I - 1 c 21 + c 22 .

Similarly, c21+ is the conjugate transpose of c21 and the value of ε is chosen by the Eq. (8) on the basis of the regularization method of SONAH [4].

As a result, we obtain the p11 expression as follows:

  p 11 = p 1 - p 21 = p 1 - p 1 c 11 c 0 - p 2 c 22 c 11 c 0 - c 12 ,

where  p11 is the pressure of source point S1 on the measurement surface H1, thus the pressure of single source on one measurement could be separated from coherent sources. In addition, the value of S1 on the construction surface prs could be calculated by  p11 according to the Eq. (7), where pT refers to the  p11. Above all, the separation of coherent sources and reconstruction of the single source are realized. In the same way, we can also obtain the relevant acoustic field information of p22.

3. Numerical simulations

According to the set-up illustrated in Fig. 2, a set of measurements were simulated. There, the grid shows an 8×8 element microphone array with 6 cm grid spacing, the microphone being at the corners of the grid. The centers of two coherent in-phase monopole point sources of equal strength are located at (–0.15, 0, 0.11) m and (0.15, 0, 0.06) m, respectively. The distance between measurement surface H1 and source point S1 was 0.06 m, the distance between measurement surface H2 and source point S1 was 0.11 m, and the gap between the two measurement surfaces was 0.17 m. SONAH calculations were executed in the measurement surfaces and construction surfaces. The calculation grid in the construction surfaces had the gird by 13×13 data points with the same grid spacing as measurement grid. The regularization parameter in Eq. (8) was set to 40 dB according to SNR, and the source distance d was set to 0.06 m. In addition, the relative average error level during calculation proceeding was given by the formula:

  L e r r = 10 l o g 10 p i t u r e - p i 2 p i t u r e 2 ,

where  Lerr is the relative average error level and piture, pi are the truth value of pressure and the calculation value of pressure, respectively.

Fig. 3. Pressure amplitude: a) the measurement surface H1; b) theoretical value on reconstruction surface R1; c) calculated value on reconstruction surface R1; d) comparison section view on xoz plane





At first, a simulation at a single frequency of 1000 Hz was investigated, i.e. the aim to separate the source point S1 from the coherent sources.

Figs. 3(a)-3(d) depict the pressure amplitude of measurement surface H1, theoretical value of reconstruction surface R1, calculated value of reconstruction surface R1 and comparison sectional view of theoretical and calculated value in reconstruction surface R1, respectively.

Fig. 4. Relative average error level of SONAH calculation for pressure in the reconstruction surface

After the sound field separation process, Fig. 3(b) shows the good agreement with the theoretical value depicted in Fig. 3(c). Here, the relative average error is –14.98 dB. It is obvious that, by using the improved sound field separation technique with SONAH, the pressure produced by a single source from the coherent sources on the measurement surface H1 can be separated out correctly.

In order to further verify the performance of sound field separation technique with SONAH, the relative average errors at different frequencies (from 100 Hz to 1400 Hz) were studied by simulation. Fig. 4 shows the relationship between frequency and relative average error. Here, the relative average error is all below –10 dB, which is acceptable.

4. Practical measurement

A practical measurement to examine the method proposed in this paper was conducted. The set-up consisted of two loudspeakers radiating steady sound field with 1000 Hz as two coherent sources. Other instruments were set according to simulation conditions. Figs. 5(a) and 5(b) indicate the schematic diagram and photo of the experimental setup.

Fig. 5. Experimental setup: a) the schematic diagram; b) practical photo



Fig. 6. Reconstruction result: a) three-dimensional diagram; b) contour plot



The pressures of time domain on the two measurement surfaces were measured by a scanning microphone composed of 3 microphones a few times and a fixed reference microphone. Every measurement surface puts 81 measuring points, which were composed of 3 microphones with 27 times, respectively. During the whole mensuration, the reference microphone collected the data as the 3 microphones all the way. Then, the complex pressures of measurement surfaces were calculated by the measured data and reference data. The amplitudes of pressures and the phased of pressures were got through FFT transform, further, the phase differences of pressure were calculated. According to the amplitudes and phase differences of pressure, the complex pressure data were worked out.

Figs. 6(a) and 6(b) show that the position of target source S1 on the reconstruction surface R1 is same as the position on the measurement surface H1 and S1 is completely separated from the coherent sources, although there are some errors, which do not disturb the recognition of the target source that is needed to be identified.

5. Conclusions

The improved sound field separation method of coherent sources based on statistically optimized near-field acoustical holography has been proposed and examined in this paper. The method, based on the measurement of pressures in two parallel layers and reconstruction of pressures on the measurement surfaces, respectively, makes use of the relationship between data on the measurement surfaces and on the reconstruction surfaces to separate the single source from coherent sources. This method could not only reduce the measurement data, but also improve the computational accuracy. Their performance in the sound field space has been examined numerically and experimentally.

This separation method of coherent sources would obtain the sound field information of single source, which avoids the interference of other sources with the same frequency. In practical applications, such method could locate the source position more precisely and identify the size of source more accurately.


  1. Maynard J. D., Williams E. G., Lee Y. Nearfield acoustic holography: I. Theory of generalized holography and the development of NAH. Journal of the Acoustical Society of America, Vol. 78, Issue 4, 1985, p. 1395-1413.
  2. Veronesi W. A., Maynard J. D. Nearfield acoustic holography (NAH) II. Holographic reconstruction algorithms and computer implementation. Journal of the Acoustical Society of America, Vol. 81, Issue 5, 1987, p. 1307-1322.
  3. Hald J. Planar Near-field acoustical holography with arrays smaller than the sound source. Proceedings of the International Congresses on Acoustics, Rome, Italy, 2001.
  4. Hald J. Patch near-field acoustical holography using a new statistically optimal method. The 32nd International Congress and Exposition on Noise Engineering Jeju International Convention Center, Seogwipo, Korea, 2003.
  5. Hald J. Basic theory and properties of statistically optimized near-field acoustical holography. Journal of the Acoustical Society of America, Vol. 125, Issue 4, 2009, p. 2105-2120.
  6. Steiner R., Hald J. Near-field acoustical holography without the errors and limitations caused by the use of spatial DFT. International Journal of Acoustics and Vibration, Vol. 6, Issue 2, 2001, p. 83-89.
  7. Oey A., Jang H. W., Ih J. G. Effect of sensor proximity over the non-conformal hologram plane in the near-field acoustical holography based on the inverse boundary element method. Journal of Sound and Vibration, Vol. 329, Issue 11, 2010, p. 2083-2098.
  8. Chen M. Y., Shang D. J., Li Q., et al. Nearfield acoustic holography based on inverse boundary element method for moving sound source identification. Acta Acustica, Vol. 36, Issue 5, 201, p. 489-4951.
  9. Zhang X. Z., Thomas J. H., Bi C. X., et al. Reconstruction of nonstationary sound fields on the time domain plane wave superposition method. Journal of the Acoustical Society of America, Vol. 132, Issue 4, 2012, p. 2427-2436.
  10. Wang Z. T., Wang R. J., Yang D. G. Quantitative measurement method for moving sound sources using acoustical holography based on wave superposition. 20th International Congress on Sound and Vibration, Vol. 3, 2013, p. 2012-2018.
  11. Bolton J. S., Kim Y. J., Niu Y., et al. Multi-reference methods for nearfield acoustical holography. Journal of the Acoustical Society of America, Vol. 129, Issue 4, 2011, p. 2491-2508.
  12. Zhang Y., Bi C.-X., Chen J., et al. Separation method of multiple coherent sources based on hybrid holographic algorithm. Chinese Journal of Mechanical Engineering, Vol. 43, Issue 9, 2007, p. 173-178.
  13. Antoni J. A Bayesian approach to sound source reconstruction: Optimal basis, regularization, and focusing. Journal of the Acoustical Society of America, Vol. 131, Issue 4, 2012, p. 2873-2890.
  14. Finn J., Chen X., Jaud V. A comparison of statistically optimized near field acoustic holography using single layer pressure-velocity measurements and using double layer pressure measurements. Journal of the Acoustical Society of America, Vol. 123, Issue 4, 2008, p. 1842-1845.
  15. Efren F. G., Finn J. Sound field separation with a double layer velocity transducer array. Journal of the Acoustical Society of America, Vol. 130, Issue 1, 2011, p. 5-8.
  16. Efren F. G., Finn J., Quentin L. Sound field separation with sound pressure and particle velocity measurements. Journal of the Acoustical Society of America, Vol. 132, Issue 6, 2012, p. 3818-3828.
  17. Bi C. X., Chen X. Z., Chen J. Sound field separation technique based on equivalent source method and its application in nearfield acoustic holography. Journal of the Acoustical Society of America, Vol. 123, Issue 3, 2008, p. 1472-1478.
  18. Calvetti D., Morigi S., Reichel L., et al. Tikhonov regularization and the L-curve for large discrete ill-posed problems. Journal of Computational and Applied Mathematics, Vol. 123, Issue 1, 2000, p. 423-446.