Abstract
The compressiveequivalent source method nearfield acoustic holography (CESM) is disturbed by reverberation in the enclosed space such as room and cabin, which leads to large reconstruction error and disturbs the judgment of sound source position. In order to solve this problem, an improved CESM algorithm based on room impulse response (RIR) is proposed to filter out reverberation interference in this paper. Different from the original equivalent source method, the improved algorithm constructs the transfer function through the room impulse response to establish the relationship between the equivalent source and the sound pressure on any plane in space, and the sparse signal reconstruction method of the compressive sensing technology is used to obtain the strength of the equivalent source. Then the transfer function to any plane of space is established according to the free field Green’s function, to eliminate the interference of reverberation and improve the effect of sound source location. The accuracy and effectiveness of the improved method are verified by preliminary numerical simulation. And the results show that compared with the original algorithm, this method has obvious advantages in sound source localization in a reverberant field.
Highlights
 An acoustic holographic algorithm suitable for use in a reverberant environment is proposed
 The mirror source method is used to establish the sensing matrix in the sound field reconstruction formula
 The reverberation environment interference is compared by simulation modeling, and the sound source recognition effect of the improved algorithm is verified
1. Introduction
With the rapid development of computer and other technologies, Monitor Technology of Mechanical Equipment has made remarkable progress in all aspects. Especially in recent years, the advance of acoustic monitoring technology, which has the advantage of noncontact intuitive, attracting a large number of domestic and foreign scholars to pay attention to the research [1]. An important research direction of acoustic monitoring is the location of sound source.
At present, sound source location technology mainly includes: traditional identification method, timedomain analysis method, acoustic holography and beamforming [2], among which acoustic holography method has many advantages such as fast measurement speed, high imaging efficiency and good spatial resolution [3]. However, most of the researches are carried out in free and semifree acoustic fields such as laboratories with silencer materials or open outdoor areas. These research methods perform better in the free field, but in the reverberation environment, the performance decreases due to the distortion caused by the reflected sound waves from the surrounding walls [46]. Most of the mechanical equipment are working in the indoor workshop, engine room and other enclosed space, so the sound source location of this equipment needs to consider environmental noise and reverberation interference factors.
In reverberation environment, the room impulse response represents the acoustic channel from the sound source to the sensor and describes the characteristics of the acoustic environment. Therefore, research on eliminating reverberation interference is generally realized by reducing or directly removing the influence of room impulse response on the received signal. The dereverberation method can be classified into reverberation suppression and reverberation cancellation according to whether or not the impulse response is estimated [7]. Among them, the reverberation cancellation method is based on reverberation estimation and inversion, and the dereverberation problem is regarded as a class of inverse problems. The room impulse response is calculated by analyzing the environment, and then solving inverse of the estimated impulse response to obtain the sound source signal [8].
The principle of reverberation cancellation is the same as the equivalent source method of nearfield acoustic holography. In the equivalent source method, the radiated sound field can be expressed as the superposition of sound fields by several equivalent sources (point sources), and the sound source intensity is obtained by inverting the transfer matrix formed by Green’s function of free field, so as to reconstruct the entire sound field [9, 10]. The solution of equivalent source intensity is the key to sound field reconstruction because it is severely illconditioned. In order to increase the accuracy and robustness of traditional methods, regularization methods [1113] were proposed. In addition, [14] introduced the compressive sensing principle and proposed the CompressiveEquivalent Source Method (CESM) which effectively improved the positioning accuracy of sound source location.
Zhao estimated the position of the sparse sound source by using the compressed sensing technology, that the dictionary is constructed based on the room impulse response of different positions in space, improving antireverberation capability of the positioning algorithm [15]. However, the positioning accuracy of this method is not high, which is related to the spatial grid division.
Considering the problem of reverberation interference in the environment, in the study we take the indoor impulse response as the transfer function to establish the sound source transfer model, estimate the equivalent source intensity through the compressive sensing method, and then reconstruct the sound field through the free field Green’s function. This method can be regarded as a kind of filtering processing, which can eliminate the reflected sound wave interference and improve the sound source location effect of the algorithm under the condition of reverberation.
The rest of this article is organized as follows. Section 2 introduces the basic principle of equivalent source nearfield acoustical holography, and the principle of CESM estimates the equivalent source strength by compressive sensing method. Section 3 briefly reviews three measurement methods of room impulse response and introduces Image Source Method (ISM) [16] in detail. Section 4 proposes an equivalent source algorithm for reverberation environment based on room impulse response. Section 5 compared and verified the correctness and feasibility of the improved algorithm through COMSOL Multiphysics software simulation modeling, and the final conclusion was given in Section 6.
2. Compressiveequivalent source method (CESM)
The equivalent source method (ESM) is a near field sound holographic algorithm with superior performance. The basic idea is to replace the actual sound field by the superposition of the sound field generated by several simple sources. It is critical to obtain the appropriate equivalent source strength by solving the equation based on the sound pressure value of the measuring surface. This is usually an underdetermined problem, which can be solved by Tikhonov regularization. However, this method is easy to be disturbed under low SNR conditions and affects the positioning effect.
The equivalent source method is to use a series of simple point sources as the equivalent source, which is evenly distributed behind the source, as shown in Fig. 1.
In a free field, the sound pressure on the measuring surface can be regarded as the superposition of the sound pressure propagated by the continuum of sources on the sound source surface:
where $p\left({\mathbf{r}}_{h}\right)$ denotes the measured surface sound pressure, $q\left({\mathbf{r}}_{v}\right)$ denotes the equivalent source intensity, $g({\mathbf{r}}_{h},{\mathbf{r}}_{v})$ denotes the free space Green’s function of the equivalent source to the holographic surface, $\rho $ is the medium density, and $c$ is the speed of sound in the medium. where $k=\omega /c$ and $\omega $ is the angular frequency.
Fig. 1Schematic diagram of equivalent source method principle, where Sv is the equivalent sound source surface, Ss is the real sound source surface, and Sh is the holographic measurement surface
After discretization, it can be expressed as:
where ${\mathbf{Q}}_{v}=[{q}_{1},\cdots ,{q}_{N}{]}^{T}\in {C}^{N}$ is the strength vector of $N$ equivalent sources are uniformly distributed in the sound source, ${\mathbf{P}}_{h}=\left[p\right({r}_{1},\omega ),\cdots ,p({r}_{M},\omega ){]}^{T}\in {C}^{M}$ is thess measured pressure vector by $M$ sensor arrays, and ${\mathbf{r}}_{M}$ is the position vector of the sensor. ${\mathbf{G}}_{hv}$ is the freefield transfer function from the equivalent source to the sensor array, and the element ${G}_{hv}(m,n)$ in the matrix represents transfer function from the $n$th equivalent source to the $m$th sensors:
The solution of Eq. (2) is the estimation of equivalent sound source strength, which is the basis of deriving the sound pressure distribution of the whole sound field. However, the number of equivalent sources is generally more than the number of array elements, $N>M$. In this case, Eq. (2) is an underdetermined equation with a nonunique solution.
The solution of Eq. (2) can constitute an optimization problem:
The common method is to solve the ${\Vert \mathbf{Q}\Vert}_{2}$ minimization by Tikhonov regularization and obtain the equivalent source strength. It is simplified as:
However, the selection of regularization parameters $\lambda $ in Ep. (6) needs to know the relative noise level and is greatly affected by the relative noise level.
Compressive sensing technology has great advantages in dealing with the reconstruction of sparse signals. Considering the fact that the sound source is not densely spaced in space and is sparsely spaced. The compressive sensing theory can be applied to the solution of equivalent source strength [17]. In this method, the source surface is equivalent to a number of uniformly distributed point sources, in which the strength of most equivalent sound sources is zero and sparsity exists.
This method considers the sparsity of sound sources in space to solve Eq. (2) under the ${\Vert \mathbf{Q}\Vert}_{0}$ minimization constraint, that is the solution with the minimum number of sound sources in space. It is very difficult to directly solve this problem, but when the problem of research is sparse, the minimization of ${\Vert \mathbf{Q}\Vert}_{0}$ can be relaxed to the ${\Vert \mathbf{Q}\Vert}_{1}$ minimization problem, and the original problem is transformed into:
The above problems can be solved by the convex optimization algorithm, such as Base Pursuit algorithm or Matching Pursuit algorithm. The partial sound radiating field can be simulated accurately by solving the equivalent source strength by compressive sensing. Compared with the traditional equivalent source method, this method has a smaller location error in a larger reconstruction volume, and the broadband sound field can be effectively reconstructed.
3. Establishment of the room impulse response model
In the actual engine room or workshop, the radiated sound wave of a sound source will produce reverberation because of the reflection on the walls, which will bring interference to the location of the sound source. Reverberation can be seen as the superposition of sound waves after multiple reflections in space. Therefore, the signal received by the sensor in the closed environment can be regarded as the convolution of a sound source signal and room impulse response, and plus environmental noise. The relationship between the source signal and the sensor in closed environment can be established by the impulse response.
For the establishment of impulse response function in reverberation environment, several methods are proposed, such as wavebased modeling, raybased modeling, statistical modeling and so on [18].
The wavebased modeling methods are mainly applied to the simple case such as rigid wall room. This method can obtain an accurate impulse response function to solve the wave equation of the room sound field by the finite element method. But it is difficult to determine the boundary conditions, and the number of computations is large.
The statistical modeling methods, such as the Statistical Energy Analysis, are widely used in highfrequency noise analysis, but this method cannot simulate the change of sound field state with time, and that is not suitable for the auralization of a sound field.
The raybased modeling methods were developed from geometrical room acoustics. The methods require simulating all possible sound reflection paths from a source to a receiver for an ideal impulse response. Such as the image source method. By setting the mirror source to indicate the source direction of the reflected sound wave, the image source method can calculate the path length of the reflected waves of each level and simulate the sound field environment to obtain the room impulse response [19]. As shown in Fig. 2.
In a cuboid closed environment, the sound wave reaches the sensor position through multiple reflections from the walls, and the discrete signal received by the sensor can be expressed as:
And $t=n\cdot T$, where $T$ is the sampling period, which is the time between the samples, and $n$ is an integer that establishes the time position of each sample. $h\left(n\right)$ is the impulse response function from the sound source $s\left(n\right)$ to the sensor, $\mathrm{*}$ is the convolution operation, and $w\left(n\right)$ is noise.
Fig. 2Reverberation model of enclosed space
After Fourier transform:
Generally, the path length of the reflection can be obtained by calculating the distance between the image source and the sensor. The room impulse response is obtained by superimposing the propagation functions of each direct wave and all reflected wave paths:
where $R=\Vert {\mathbf{R}}_{\theta}+{\mathbf{R}}_{f}\Vert $, ${\mathbf{R}}_{p}$ is the distance from the sound source or image source to the sensor, and $\theta $ is the eight sound source locations (include image source), where $q$, $j$, $\upsilon $ takes on values 0 or 1, and $\gamma $ is the order of reflection:
The position coordinates of the sound source $S$ is ($x$, $y$, $z$), the sensor position coordinates ($x\mathrm{\text{'}}$, $y\mathrm{\text{'}}$, $z\mathrm{\text{'}}$); ${\mathbf{R}}_{f}$ is the virtual room size of the highorder reflection:
where (${L}_{x}$, ${L}_{y}$, ${L}_{z}$) is the size of the room, the range of each element of (${m}_{x}$, ${m}_{y}$, ${m}_{z}$) is a positive integer set.
In fact, the reflection coefficient of each wall should also be considered. Eq. (10) is amended as:
where ${\beta}_{x,\xi}^{}$ represents the reflection coefficient of two walls in the $x$axis, $\xi =$ 1 represents the wall close to the origin of the coordinate, and $\xi =$ 2 represents the wall away from the origin of the coordinate, the coefficient being related to the sound absorption coefficient $\alpha $ of the wall material, $\beta =\pm \sqrt{1\alpha}$.
According to the specific environment, the impulse response from the sound source to any position in the closed space can be obtained through Eq. (13).
4. Equivalent source method based on room impulse response
In the reverberation environment, the acoustic holography algorithm based on compressive sensing is seriously affected by the reflected sound wave interference.
By analyzing the principle of ESM, the improved algorithm replaces the transfer function in the free field with the room impulse response function to describe the signal received by the sensor more accurately, and then reconstructs the sound field through the freefield transfer function to achieve the filtering of the reflected sound wave.
4.1. Constructing the transfer function based on the room impulse response
As shown in Fig. 3, it is assumed that $N$ simple sound sources of different intensities are uniformly distributed in the sound source, and the position of the point source is marked as ${q}_{i}$. The point source of position ${q}_{i}$ to the $m$th array element on the sensor array has a room impulse response of ${H}_{m,{q}_{i}}\left(k\right)$, which together form a reverberation transfer function ${\mathbf{G}}_{h}\left(k\right)$.
Fig. 3Measurement schematic diagram of CS equivalent source method different grayscale of point source indicates different strength
${H}_{m,{q}_{i}}\left(k\right)$ is obtained by image sound method, and the specific process is shown in Fig. 4.
Fill the calculated ${H}_{m,{q}_{i}}\left(k\right)$ into Eq. (14), and get ${\mathbf{G}}_{h}\left(k\right)$:
4.2. Calculate equivalent source strength
The received signal of $M$ sensors can form a vector $\mathbf{X}=[{x}_{1},{x}_{2},\dots ,{x}_{M}]$, which is expressed in the wavenumber domain as follows:
where $\mathbf{W}\left(k\right)$ represents environmental noise. The Eq. (15) is consistent with the Eq. (2), and the number of equivalent sources is sparse in the spatial domain. The equivalent source strength vector $\mathbf{S}\left(k\right)$ can be reconstructed by the method of compressive sensing, and the elements of the vector $\mathbf{S}\left(k\right)$ represent the strength of the equivalent sound source at different locations.
Fig. 4A flow diagram for constructing the impulse response function from a sound source point to the sensor
The equivalent source strength can be solved by solving the following equation:
This paper uses orthogonal matching pursuit algorithm to solve [20], and the general steps are as follows:
1. Initialization: ${\mathbf{r}}_{0}=\mathbf{X}$, ${\mathrm{\Lambda}}_{0}=\mathrm{\xd8}$, ${G}_{h}=\xd8$, $t=$ 1;
2. Search index ${\lambda}_{t}$: ${\lambda}_{t}=\mathrm{a}\mathrm{r}\mathrm{g}\underset{j=\mathrm{1,2},\cdots ,N}{\mathrm{m}\mathrm{a}\mathrm{x}}\left\u27e8{\mathbf{r}}_{t1},{\mathbf{\alpha}}_{j}\u27e9\right$;
3. Settings ${\mathrm{\Lambda}}_{t}={\mathrm{\Lambda}}_{t1}\cup \left\{{\lambda}_{t}\right\}$, ${\mathbf{G}}_{h,t}={\mathbf{G}}_{h,t1}\cup {\mathbf{\alpha}}_{{\lambda}_{t}}$;
4. Find the least squares solution of $\mathbf{X}={\mathbf{G}}_{h,t}{\mathbf{\theta}}_{t}$:${\widehat{\mathbf{\theta}}}_{t}=\mathrm{a}\mathrm{r}\mathrm{g}\underset{{\theta}_{t}}{\mathrm{m}\mathrm{i}\mathrm{n}}\Vert \mathbf{X}{\mathbf{G}}_{h,t}{\mathbf{\theta}}_{t}\Vert $;
5. Update residual: ${\mathbf{r}}_{t}=\mathbf{X}{\mathbf{G}}_{h,t}{\widehat{\mathbf{\theta}}}_{t}$;
6. Judge $\Vert {\mathbf{r}}_{t}\Vert \le \epsilon $, if true, outputs the estimated value ${\widehat{\mathbf{\theta}}}_{t}$ of $\mathbf{S}$, otherwise it returns to step 2.
4.3. Reconstructing the sound pressure distribution of any plan
According to the case of the equivalent source in the free sound field, the transfer matrix ${G}_{D}$ is constructed with the green’s function of the equivalent source:
where ${g}_{m,{q}_{n}}\left(k\right)$ is the Green’s function in the free field, and is calculated according to Eq. (4).
After obtaining the equivalent source strength, the sound pressure of any plane in the closed space can be reconstructed according to Eq. (15):
5. Simulation experiment
The experiment simulates the nearfield acoustic radiation environment generated by metal disk under excitation force. The reverberant field model (Fig. 5) in the closed rectangular space and the free sound field model (Fig. 6) were established respectively. In the model of free sound field, all boundaries absorb sound waves completely without any reflection. In the reverberation model, the sound reflection coefficient of boundaries is set in accordance with Table 1. The setting parameters of the two models are identical except the boundary condition.
Table 1The reflection coefficient of the inner wall of the model room
The surface of the space  Sound reflection coefficient 
The surrounding walls  0.9 
Floor  0.89 
Ceiling  0.99 
Fig. 5The reverberation sound field model
Fig. 6The free sound field model
Set a rectangular parallelepiped space with a length, width and height of 3×4×3 m, with the corner of the cuboid as the origin, and the edge line as the coordinate system to establish the coordinate system, and fill air in the space. Place a disc with a radius of 0.05 m and a thickness of 0.01 m at the (1.5, 1, 1.5) point of the room. The material of the disc is structural steel. The edge of the disc is fixed, and the sound field is generated by applying an excitation force at its center. The excitation force is comprising a series of harmonic signals over a frequency range of 502000 Hz and an interval of 50 Hz.
The pressure acoustic simulation environment was selected for simulation calculation in COMSOL Multiphysics software. The simulation model is divided with the free tetrahedral mesh, the frequency domain compilation equation is selected for compilation, and the steadystate solver is used to solve the equation.
Fig. 7Data collection array
Fig. 8The sound pressure distribution at different frequencies in the free field and reverberation field
a)
b)
c)
d)
e)
f)
g)
h)
i)
The left column represents the sound pressure distribution at $y=$ 1.1 m in the free field, and the other two columns represent the sound pressure distribution at $y=$ 1.1 m and $y=$ 1.2 m in the reverberation field.
At a position 0.2 m away from the disk, a data collection array is placed to collect the sound pressure data in the sound field, as shown in Fig. 7. The spacing of array elements is 0.05 m, and the number of array elements is 21×21. The sound pressure value in the sound field is collected by the array.
From the simulation data, the sound pressure distribution of the plane at $y=$ 1.1 m is obtained in the free field environment and the reverberation environment at 400 Hz, 600 Hz, and 1400 Hz, as shown in Fig. 8.
It can be seen from the distribution of sound pressure that the reverberation seriously interferes with the judgment of the position of the sound source. As the frequency increases, the sound pressure distribution becomes more complicated, the sound wave has stronger directivity, and the reflection efficiency increases in the closed space, which makes the distribution of sound pressure in the space more confusing, and it is difficult to judge the sound source position.
Fig. 9 shows sound pressure distribution images at 400 Hz, 600 Hz, and 1400 Hz, which are obtained by processing the sound pressure data collected in the reverberant field by the CESM method and the improved method, respectively, and thereby determining the sound source position. The equivalent source number of both methods is set to 31×31.
Fig. 9Sound pressure image at y= 1.1 m obtained by the two algorithms
a)
b)
c)
d)
e)
f)
It can be seen from Fig. 9 that the sound pressure distribution image obtained by the CESM algorithm determines that the sound source position will cause a serious positional deviation, and the positioning method fails. The sound pressure image obtained by the improved algorithm more accurately indicates the sound source position, that is, the hot spot position of the sound pressure image can be discriminated as the sound source position, indicating that the algorithm effectively removes the environmental reverberation interference. The improved algorithm is suitable for sound source localization in a reverberant environment.
In order to compare the reconstructed sound pressure with the value of the sound pressure in the free field, Eq. (19) is used to describe the reconstruction error quantitatively, which reflects the approximate degree of the sound pressure value on the whole:
where ${p}_{t}$ represents the reconstructed sound pressure value matrix, and ${p}_{r}$ represents the sound pressure value matrix measured in the free field environment. According to Eq. (19), the reconstruction errors of the two methods at 502000 Hz frequency are obtained, as shown in Fig. 10. If the reconstruction error is less than 30 %, the sound source location can be clearly identified, indicating that the reconstruction is more accurate, and the misjudgment location will occur when the error is more than 40 %.
It can be seen that in the reverberation environment, the CESM method still has limited localization effect at low frequencies. However, as the frequency increases, the algorithm cannot accurately determine the location of the sound source. Especially when the frequency is higher than 350 Hz, the original algorithm completely fails.
In a wide frequency band, the reconstruction error of the improved algorithm is about 30 %, and the sound source location can be determined stably and accurately. This algorithm has obvious advantages in the location of sound source in reverberation environment, but with the increase of frequency, the location effect decreases, and the location of sound source will be misjudged.
Fig. 10Reconstruction error at different frequencies
6. Conclusions
The localization of sound source in a reverberant environment is studied，based on the equivalent source method of nearfield sound holography. By modifying the construction of the transfer matrix, the accuracy of solving the equivalent source strength is improved. Then, the sound pressure field distribution without reverberation interference is deduced by using the transfer function constructed by the freefield transfer function, so as to direct the location of the sound source. The whole process of the algorithm can be seen as filtering the reverberation interference, leaving only the direct wave signal and reconstructing the free environment sound field.
Theoretical analysis and simulation results show that, compared with the CESM method, in the reverberation environment, the improved algorithm can accurately guide the location judgment of sound source in a wide frequency range, broaden the application scope of the CESM method, and is an exploration of the application in the actual engineering environment. However, in general, the reconstruction error will increase with the increase of frequency, and the situation of the highfrequency band needs further improvement.
References

Mba D., Rao Raj B. K. N. Development of acoustic emission technology for condition monitoring and diagnosis of rotating machines: bearings, pumps, gearboxes, engines, and rotating structures. The Shock and Vibration Digest, Vol. 38, Issue 1, 2006, p. 316.

Chen Xinzhao Progress of noise source identification technology. Journal of Hefei University of Technology, Vol. 32, Issue 5, 2009, p. 609614, (in Chinese).

Jia Wenqiang, Chen Jin, Li Jiaqing, et al. Local sound field reconstruction technique based on combined wave superposition and beamforming method. Journal of Vibration and Shock, Vol. 29, Issue 1, 2010, p. 125127, (in Chinese).

Park S. H., Kim Y. H. Visualization of passby noise by means of moving frame acoustic holography. The Journal of the Acoustical Society of America, Vol. 110, Issue 5, 2001, p. 2326.

Xi Jiaxin, Feng Zengming, Wang Guoqiang, et al. Vibration and noise source identification methods for a diesel engine. Journal of Mechanical Science and Technology, Vol. 29, Issue 1, 2015, p. 181189.

XinZhao C., ChuanXing B., WeiBing L. NAH technology and its application in the noise source. Modern Technology of Vibration and Noise, Vol. 4, 2005, p. 311.

Habets E. E. Single and multimicrophone speech dereverberation using spectral enhancement. Eindhoven, Eindhoven University of Technology, Netherlands, 2007.

Huang Y., Benesty J., Chen J. Acoustic MIMO signal processing. Signals and Communication Technology, Vol. 65, Issue 47, 2006, p. 233236.

Valdivia N. P., Williams E. G., Herdic P. C. Equivalent sources method for supersonic intensity of arbitrarily shaped geometries. Journal of the Acoustical Society of America, Vol. 135, Issue 4, 2014, p. 4662.

Koopmann G. H., Song L., Fahnline J. A method for computing acoustic fields based on the principle of wave superposition. Journal of the Acoustical Society of America, Vol. 86, Issue 5, 1989, p. 24332438.

Leclere Q. Acoustic imaging using underdetermined inverse approaches: frequency limitations and optimal regularization. Journal of Sound and Vibration, Vol. 321, Issues 35, 2009, p. 605619.

Kim Y., Nelson P. A. Optimal regularisation for acoustic source reconstruction by inverse methods. Journal of Sound and Vibration, Vol. 275, Issues 35, 2004, p. 463487.

Tikhonov A. N. Solution of incorrectly formulated problems and the regularization method. Soviet Mathematics Doklady, Vol. 4, 1963, p. 10351038.

FernandezGrande E., Xenaki A., Gerstoft P. A sparse equivalent source method for nearfield acoustic holography. The Journal of the Acoustical Society of America, Vol. 141, Issue 1, 2017, p. 532542.

XiaoYan Z., Lin Z., ZhenYang W. Sound source localization algorithm for microphone array based on compressed sensing. Journal of Southeast University (Natural Science), Vol. 45, Issue 2, 2015, p. 203207, (in Chinese).

Habets E. A. P. Room Impulse Response Generator. Technical Report, Technische Universiteit Eindhoven, 2006.

Hu Dingyu, Li Hebing, Hu Yu, et al. Sound field reconstruction with sparse sampling and the equivalent source method. Mechanical Systems and Signal Processing, Vol. 108, 2018, p. 317325.

Savioja L., Huopaniemi J., Lokki T., Väänänen R. Creating interactive virtual acoustic environments. Journal of the Audio Engineering Society, Vol. 47, Issue 9, 1999, p. 675705.

Allen J., Berkley D. Image method for efficiently simulating small room acoustics. Journal of the Acoustical Society of America, Vol. 65, Issue 4, 1979, p. 943950.

Tropp J. A., Gilbert A. C. Signal recovery from random measurements via orthogonal matching pursuit. IEEE Transactions on Information Theory, Vol. 53, Issue 12, 2007, p. 46554666.