A Robust and Efficient Compressed Sensing Algorithm for Wideband Acoustic Imaging

Wideband acoustic imaging, which combines compressed sensing (CS) and microphone arrays, is widely used for locating acoustic sources. However, the location results of this method are unstable, and the computational efficiency is low. In this work, in order to improve the robustness and reduce the computational cost, a DCS-SOMP-SVD compressed sensing method, which combines the distributed compressed sensing using simultaneously orthogonal matching pursuit (DCS-SOMP) and singular value decomposition (SVD) is proposed. The performance of the DCS-SOMP-SVD is studied through both simulation and experiment. In the simulation, the locating results of the DCS-SOMP-SVD method are compared with the wideband BP method and the DCS-SOMP method. In terms of computational efficiency, the proposed method is as efficient as the DCS-SOMP method and more efficient than the wideband BP method. In terms of locating accuracy, the proposed method can still locate all sources when the signal to noise ratio (SNR) is − 20 dB, while the wideband BP method and the DCS-SOMP method can only locate all sources when the SNR is higher than 0 dB. The performance of the proposed method can be improved by expanding the frequency range. Moreover, there is no extra source in the maps of the proposed method, even though the target sparsity is overestimated. Finally, a gas leak experiment is conducted to verify the feasibility of the DCS-SOMP-SVD method in the practical engineering environment. The experimental results show that the proposed method can locate both two leak sources in different frequency ranges. This research proposes a DCS-SOMP-SVD method which has sufficient robustness and low computational cost for wideband acoustic imaging.


Introduction
Acoustic imaging, which uses planar microphone array and beamforming methods [1][2][3][4][5], is widely employed for locating acoustic sources. As for wideband acoustic imaging, one way is to operate the time-domain beamforming technique. The time-domain method has a relatively high computational efficiency and a very good application for non-stationary and strongly transient signals [6]. Zhao using a Tap delay line structure or finite impulse response (FIR) filter to achieve different frequency beamforming [7]. Wilkins presented a true time delay (TTD) beamformer bank in the beamspace and permitted the directions of arrival of broadband sources to be estimated accurately, efficiently, and noniteratively [8]. However, runtime delay quantization effects can cause a need for high sampling and processing of huge amounts of data [9]. Another way to deal with wideband acoustic imaging is to operate in the frequency domain. For wideband signal, we can use the fast Fourier transform (FFT) to divide the array outputs into many narrowband frequency bins and apply some beamforming methods to each narrowband frequency bins, such as delay-and-sum (DAS) beamforming [10], standard capon beamforming(SCB) method [11] and robust capon beamforming method(RCB) [12]. However, the problem of varying mainlobe width as a function of frequency exists on the results when

Open Access
Chinese Journal of Mechanical Engineering those methods are applied to wideband acoustic imaging. Wang et al. [13] have proposed the shaded robust capon beamformer method (SRCB), which can obtain approximately constant mainlobe width for wideband acoustic imaging. In Ref. [14], robust Capon beamforming with pre-steering produces can locate acoustic emissions with significant accuracy and less ghost than ordinary beamforming. The wideband acoustic imaging also can be performed by jointing different narrowband frequency bins. Kassis et al. [15] have proposed the wideband zero-forcing MUSIC (ZF-MUSIC) method for locating aeroacoustic sources. The wideband ZF-MUSIC criterion was proposed to avoid the maximization of the criterion at each frequency bin. The ZF-MUSIC method increases the ability to separate sources having different powers than the MUSIC method, but the width of the frequency band cannot be too large in case the low frequencies degrade the resolution. Guo developed a robust nearfield wideband beamformer design approach based on adaptive-weighted convex optimization. The method employs the adaptive array signal processing theory, adjusts weights flexibly, and improves the beamforming performance [16]. He proposed a new direction of arrival (DOA) estimation method of wideband source, which is based on iterative adaptive spectral reconstruction, which can be applied to coherent sources and improve the accuracy of DOA estimation [17]. But they all increase computational complexity and time costs, cannot meet the real-time requirements. Donoho [18], Candès, Romberg, and Tao [19][20][21][22] proposed the theory of compressed sensing (CS). CS shows that a signal having a sparse representation can be recovered exactly from a small set of linear, nonadaptive measurements. If most entries of a signal are zeros, the signal is sparse. As for the acoustic imaging problem, the number of sources that are usually assumed to be point sources is much less than the node number behind the grid. Therefore, compressed sensing can be easily used for acoustic imaging. Chu et al. [23] have applied the Bayesian CS method to nearfield wideband aeroacoustic imaging. The method has good robustness in poor signal to noise ratio (SNR) cases and can obtain a wide dynamic range. However, it has more computational costs than beamforming methods. Chaturvedi used CS to reconstruct cross-correlation of wideband signals from the cross-correlation of sub-Nyquist samples to estimate the DOA [24], but cross-correlation provides inferior accuracy in the experiment [25]. Boufounos et al. combined joint sparsity models and CoSaMP algorithm to wideband array processing. However, it is known that CoSaMP algorithm fails to provide satisfactory performance in the source location and spectral estimation applications, especially in the presence of closely spaced sources [26].
In this paper, we propose a new CS method called the DCS-SOMP-SVD method for wideband acoustic imaging, which combines the distributed compressed sensing using simultaneously orthogonal matching pursuit (DCS-SOMP) method and singular value decomposition (SVD). In this study, we will study the performance of the DCS-SOMP-SVD method for wideband acoustic imaging by comparing that with the wideband basis pursuit (BP) method and the DCS-SOMP method.
This paper is organized as follows: Section 2 describes the observation model of acoustic signal propagation and the math model of the wideband acoustic imaging. Then our proposed method is presented in Section 3. Subsequently, the performance of the proposed method is compared with the other two methods by simulations in Section 4. The analysis of the DCS-SOMP-SVD method is also shown in this section. Section 5 provides a gas leakage experiment to verify the feasibility of the DCS-SOMP-SVD method in actual application. Finally, we conclude this paper in Section 6.

Observation Model and Math Model
for Wideband Acoustic Imaging [·] T denotes the transpose operator. The source plane z 0 is discretized into N = u × u equidistant grids at known (2020) 33:95 discrete positions P = [P 1 , · · ·, P N ] . Let Y be the vector of wavefield measurement at M microphones of the array in the frequency domain, and the unknown vector X comprise source strengths at all N grid nodes. With the help of the microphone array, we can get the pressure fields of the microphones and obtain Y = [Y 1 , · · ·, Y M ] T by operating FFT. The n-th element of X equals zero if there is no source at the nth gird node. Otherwise it is nonzero. The sources in our model are supposed to be uncorrelated monopoles in order to simplify the physical process and build up the acoustic propagation model explicitly [27].

Observation Model for Acoustic Imaging
The pressure field at the m-th microphone is given by: where r mn = � − P m − P n � denotes the distance between the mth microphone and the nth grid node, ω = 2πf with f being the frequency, the wave number k = ω/c with c being the sound speed, and X n is the amplitude of the nth grid node.
The model can be compactly expressed in matrix form: where A is a M × N matrix and defined as the measurement matrix.
In the practical engineering environment, there are often measurement errors and background noise. In this paper, we assume errors and background noise as additive Gaussian white noise (AGWN), which is mutually independent and identically distributed and independent to sources. Thus, a more realistic propagation model can be depicted as: where e denotes background noise and errors.
By solving the linear system Eq.
(2), we can recover the signal x ∈ C N . Mathematically speaking, if we want to solve the linear system Eq. (2) without any distortion, the number of measurements M, i.e., the number of microphones in acoustic imaging, should be at least as large as the signal length N . Otherwise, the system will be severely underdetermined and have no unique solution.
Thanks to sparsity, one can perfectly recover X by solving the following optimization problem: Unfortunately, the l 0 -minimization problem is an NPhard problem [28][29][30] and thus computationally intractable. Candès and Tao [31,32] have proved that under certain condition, Eq. (4) is equivalent to the following l 1 -optimization problem: where �X� 1 = n i=1 |x i |. As to Eq. (3) which takes noise into account, we can solve it by solving the following second order cone programming (SOCP): where ε is a specified tolerance for noise e.
Eqs. (5) and (6) are the convex relaxation of their according to original NP-Hard problems and can be solved by basis pursuit (BP) (also known as l 1 -minimization method) with polynomial computational time [31,[33][34][35]. The BP method has both merits and drawbacks. BP provides theoretical performance guarantees, but its computational costs may be a limitation.
Apart from convex relaxation, several greedy methods, which compute the support set of the signal iteratively and approximate the sparse signal of Eq. (3) until a preset stopping condition [36][37][38][39][40][41][42], are also widely used. Greedy methods have the advantages of easy implementation, fast convergence, and low complexity.

Joint Sparsity Model for Wideband Acoustic Imaging
For the narrow acoustic imaging, we have developed the CS algorithm based on the greedy algorithm [43][44][45]. For a wideband signal, we choose its characteristic frequency band, and its lower and upper bound frequency are f min and f max . A simple method is to divide the chosen wideband into several narrowband, then for each frequency f j we have: where J is the number of frequencies within [ f min , f max ].
We can directly solve Eq. (7) by the CS algorithms, but the computational efficiency can be a problem as the bandwidth grows.
However, the positions of the sound sources do not change with frequency. In other words, the signals at different frequencies satisfy the simultaneous sparse approximation. So, we combine the signal vectors at each frequency into a signal matrix and solve it jointly [46]. So, the problem can be transformed into an optimization problem with the help of joint sparsity: T denotes the energy vector, which is the mean square of source power on each node, i.e., X (l 2 ) n = �X n f 1 , X n f 2 , . . . , X n (f J )� 2 . − X and − Y can be obtained by stacking the data Y (f j ) and signal vec- A is a block-diagonal matrix with each measurement matrix A(f j ) as its element, and − ε is a specified tolerance for noise − e.
Then the original problem has been simplified to a SOCP problem. Thus the BP can be used. Furthermore, the greedy algorithm can also be applied. Simultaneous orthogonal matching pursuit (SOMP) is a kind of greedy pursuit algorithm proposed by Tropp et al. [47] on the basis of OMP. It can compute provably good solutions to several simultaneous sparse approximation problems. Sarvotham et al. [48,49] expanded the theory for distributed compressed sensing (DCS) and achieved joint recovery of multiple signals from incoherent projections through SOMP.
DCS-SOMP has the same advantages as OMP and has higher computational efficiency than BP. However, we found that DCS-SOMP could not locate sources in low SNR. So, we propose a robust method combining DCS-SOMP and SVD, which can be used for wideband acoustic imaging even though the SNR is low.

DCS-SOMP-SVD Method for Wideband Acoustic Imaging
In this section, we present a practical approach for wideband acoustic imaging in low SNR environments, which combines the DCS-SOMP algorithm and the SVD. The measurements of sound pressure received by the microphone array in time-domain are divided into B blocks, where each block contains L data points and has 50% overlap. We perform L-point discrete Fourier transform (DFT) and choose the data between the lower and the upper bound frequency f min and f max . We can obtain a M × B data matrix y(f j ) under each frequency f j from the B blocks of data: Similarly, the source strength can also be divided as a N × B matrix x f j , which consists of source strengths under the frequency each frequency f j . We employ SVD for y f j : We define the reduced M × K dimensional matrix y f j SV , which involves most of the signal power as Here K is the source sparsity, which is the actual number of sources. I K is a K × K identity matrix and 0 is a K × (B − K ) zero matrix. In addition, we transform the signal matrix x f j at each frequency f j as x f j SV = x f j V D K , and let e f j SV = e f j V D K , to obtain the system: After SVD, the signal subspace is reserved, and the noise subspace is abandoned. Then we can apply the DCS-SOMP algorithm to solve the problem in Eq. (11). Same as the OMP algorithm, the DCS-SOMP algorithm also has two ways to terminate iterations. We have introduced the difference between the two ways in our previous work [45]. Unfortunately, the source sparsity K is unknown a priori in many cases. In this work, we use the target sparsity K T , which is larger than K , as the stopping condition, same as the approach of the OMP-SVD [45,50].
Donoho et al. [51] gave a phase diagram to depict the performance of CS. The diagram shows that high-accuracy reconstruction can be obtained for small ρ(= K /M) and large δ(= M/N ) , while for large ρ and small δ the reconstruction fails. The phase transition analysis [52] shows that the maximum source sparsity K can be accurately reconstructed with an empirical formula M ≈ 2K log(N) [42]. To obtain high-accuracy reconstruction (small ρ and large δ ) of the DCS-SOMP-SVD, the maximum possible number of sources K is obtained by By applying DCS-SOMP, we can obtain an approximate solution X f j at each frequency f j and transform X f j into the matrix X f j SV .Then we can get the strengths and locations of sources by averaging the source strengths of all columns of signal subspace X f j SV . The (discrete) source strength X * j can be obtained by denotes the source power covariance matrix and X f j SV (k) denotes the kth column of the matrix X f j SV . As the location of unknown sources is assumed to be the same, we can depict the maps with location and overall (11)  where x * jn is the power of the nth node for the jth frequency within the frequency range. Now we are ready to present our DCS-SOMP-SVD method for wideband acoustic imaging. The sequence of steps is as follows. [48,49] Step 1: Obtain the measurement of sound pressure in time-domain from the microphone array.
Step 2: Give the lower and upper frequency bound f min and f max . For each frequency f ∈ [f min , f max ] , construct the measurement matrix A according to the frequency f , the node position, and the distance z 0 with Eq. (2).
Step 3: Divide the measurement of sound pressure into B blocks at each frequency f j and each block contains data block lengths L , where each block has 50% overlap.
Step 4: Obtain the M × B data matrix y(f j ) at each frequency f j by performing k-point DFT on the data in step 1.
Step 5: Compute the source sparsity K with Eq. (12) as the maximum number of sources can be located. On this basis, we use the target sparsity K T as the stopping condition.
Step 6: Discretize the source plane, which is z 0 away from the array plane, with u × v = N nodes.
Step 7: Perform SVD for the data matrix y(f j ) with Eq. (10).
Step 8: Transform the data matrix y as y(f j ) SV = U ΛD K = y f j ΛD K .
Step 9: Repeat from step 3 to step 8 for signal data and error data to obtain x(f j ) SV and e(f j ) SV , respectively.
Step 11: Set the iteration counter l = 1 . For each signal index j ∈ {1, · · ·, J } , initialize the orthogonalized coefficient vectors β j = 0 , also initialize the set of selected indices Ω = ∅ . Let r j,l denote the residual of the measurement y(f j ) remaining after the first l iterations, and initialize r j,0 = y f j .
Step 12: Select the dictionary vector that maximizes the value of the sum of the magnitudes of the projections of the residual at each narrow band, and add its index to the set of selected indices: where A j,n is the nth column of the measurement matrix A f j . (14) OASPL n = 10log 10 10 x * 1n 10 + 10 x * 2n 10 + · · · + 10 x * Jn 10 (15) n l = argmax n=1,2,...,N J j=1 �r j,l−1 , A j,n � �A j,n � 2 (16) Ω = Ω n l Step 13: Operate Schmidt regularization and orthogonalize the selected basis vector against the orthogonalized set of previously selected dictionary vectors: where γ k,l is the regularization result of selected column A j,n j , which equals to the multiplication of the amplitude of A j,n j and Q after QR factorization.
Step 14: Update the estimate of the coefficients β j for the selected vector and residuals r j,l : where β j = R j,l ∼ X j,l with ∼ X j,l being the least square solution of linear system r j,l−1 = A j,n l X j,l and R j,l being the division of R after QR factorization of the selected column A j,n l and the amplitude of A j,n l .
Step 15: l = l + 1 . Return to step 12 until l = K T .
Step 16: Apply QR factorization on the mutilated basis A j, = Q j R j = Γ j R j . Since in each narrow-band y f j = Γ j R j = A j, X j, = Γ j R j, , where X j, is the mutilated coefficient vector, we can compute the signal estimates ∼ X j as: Step 17: End iteration and let X f j =X j .
Step 18: Reshape X f j as the N × K matrix X f j SV .
Step 20: Find the source positions using the indices of the nonzero elements in OASPL and the observation model in Section 2.

Simulation Results and Analysis
In this section, the wideband BP method, DCS-SOMP method, and DCS-SOMP-SVD method are compared in terms of computational effort, reconstruction precision of OASPL, and the robustness to signal to noise ratio (SNR). More, we investigate the effect of frequency range and iteration number on our method. The source maps for the wideband BP method are obtained by employing convex optimization (CVX) toolbox, a package for (17)  specifying and solving convex programs [53], to solve the Eq. (8). The value of ε in Eq. (8) is set to be the norm of background noise roughly, where the background noise is the difference between the measurement data with the AGWN and those without the AGWN.

Simulation Configuration
In the simulation, an optimized random plane microphone array with a circular aperture of 1.0 m is used to obtain the measurement data. The plane microphone array, which is shown in Figure 2, consists of 60 microphone sensors. The interested sources plane is a 0.8×0.8 m 2 rectangular zone, which is 0.8 m away from the array plane. A piece of 14 s music is regarded as the source signal added additional Gaussian white noise (AGWN), and its spectrum, which is shown in Figure 3, is a typical wideband spectrum.

Method Comparison
In  Figure 4 shows the source maps obtained by the wideband BP method, the DCS-SOMP method and the DCS-SOMP-SVD method. The squares in the source maps indicate the actual source positions, while the small colored circular dots indicate the reconstructed source positions. Figure 4(b) and (c) indicate that the DCS-SOMP and the DCS-SOMP-SVD method can obtain the super-resolution source maps like the wideband BP method, and the whole sources can be located by the three methods exactly. The CPU-time of methods is obtained using a Core i7 multi-core PC. The CPU-time of the DCS-SOMP-SVD method is 2.463 s, which is a little slower than the DCS-SOMP method, but is far less than the CPU-time 126.5 s of the wideband BP method. More, the DCS-SOMP-SVD method has a smaller OASPL error ratio among each source in comparison with the wideband BP method, which is shown in Table 1. We conclude that the DCS-SOMP-SVD method is as efficient as the DCS-SOMP method and much more efficient than the wideband BP method.

Effect of SNR on Source Maps
In this section, we investigate the effect of SNR on the source maps of the wideband BP method, the DCS-SOMP method, and the DCS-SOMP-SVD method. We choose 2500-3000 Hz as the frequency range, which has good performance in Section 4.2. We choose the same sources, as shown in Section 4.2. The performances of the wideband BP method, DCS-SOMP method, and our method are studied as SNR equals 0 dB, −10 dB, and −20 dB.  (2020) 33:95 The source maps obtained by the wideband BP method, the DCS-SOMP method, and the DCS-SOMP-SVD method at different SNRs are shown in Figure 5, Figure 6, and Figure 7. When the SNR equals 0 dB, the DCS-SOMP-SVD method can locate the whole sources exactly, while the wideband BP method can only locate three sources, and the DCS-SOMP method has 1 grid error in locating sources. Besides, the max OASPL error of DCS-SOMP-SVD method among the whole sources being 0.0041% is less than the max OASPL error of the wideband BP method being 10.6171% and the max OASPL error of DCS-SOMP method being 1.4726%. As the SNR reduces to −10 dB, the DCS-SOMP-SVD method can still locate whole sources, while the DCS-SOMP method and the wideband BP method can only locate part of the sources. Besides, some reconstructed sources outside of the expected source positions exist on the source maps of wideband BP method, and one of reconstructed sources outside of the expected source positions has larger OASPL than the reconstructed source located at the actual source position, as shown in Figure 6(a). The DCS-SOMP-SVD method has better OASPL reconstruction precision than the DCS-SOMP method and the wideband BP method. When the SNR equals −20 dB, the DCS-SOMP-SVD method can still locate all sources accurately, while the DCS-SOMP method and the wideband BP method have failed to locate any sources. Compared with the wideband BP method and the DCS-SOMP method, our method has better robustness in low SNRs.

Effect of Frequency Range on Sources Maps
In this section, we investigate the effect of the frequency range on the source maps. The wideband BP method  has a relatively large computational effort in Section 4.2, which limits the use of the wideband BP method at a wider frequency range. Therefore, we focus on studying the effect of frequency range on the DCS-SOMP method and the DCS-SOMP-SVD method. Here, the frequency ranges are chosen as 1000-1500 Hz, 1000-2000 Hz, and 1000-3000 Hz, and the SNR is set to 0 dB. Figure 8 shows the source maps obtained by the DCS-SOMP method and the DCS-SOMP-SVD method at different frequency ranges. When the frequency range equals 1000-1500 Hz, Figure 8(b) shows there are only three minor errors on the source maps. However, Figure 8(a) shows that the DCS-SOMP method can only locate one source. There are over ten grid errors in locating other sources, which is considered that the source having the max minor error cannot be located by the DCS-SOMP method. When the frequency range broadens to 1000-2000 Hz, the position errors become smaller. The whole sources are located by the DCS-SOMP-SVD method, in which the max error is merely one grid as shown in Figure 8(d). Figure 8(c) shows the DCS-SOMP method can also locate all sources while the max errors are four grids. Furthermore, as the frequency range equals 1000-3000 Hz, Figure 8(f ) shows that the whole sources can be located by the DCS-SOMP-SVD method correctly. The DCS-SOMP method can also have a satisfying result when the frequency range equals 1000-3000 Hz. There is only one grid error in the map obtained by the DCS-SOMP method shown in Figure 8(e). From the above comparison, as the frequency range broadens, the location results gradually become better. Table 2 shows the error of the locating results of the DCS-SOMP method and the DCS-SOMP-SVD method (2020) 33:95 under different frequency ranges. We can see that the results obtained by the DCS-SOMP-SVD are better than that of the DCS-SOMP method at the same frequency ranges.
In order to identify the performance of the DCS-SOMP-SVD method after broadening frequency range, the max OASPL error being merely 0.0008% among whole sources is computed as the frequency range equals 1000-3000 Hz. It can conclude that the DCS-SOMP-SVD method can obtain better performance by broadening the frequency range.

Effect of the Target Sparsity on Sources Maps
In the above study, the iterative termination conditions of both the DCS-SOMP method and the DCS-SOMP-SVD method are set to be the target sparsity, which are consistent with the actual source sparsity. However, in real cases, the source sparsity is not known in advance, and we use the overestimated target sparsity, which is larger than the actual source sparsity. Therefore, it is necessary to study the effect of the target sparsity of the DCS-SOMP and the DCS-SOMP-SVD method on source maps. In this section, the performance of the DCS-SOMP and the DCS-SOMP-SVD method are investigated as the target sparsity equals 10. The frequency range and SNR are set to 2500-3000 Hz and 0 dB, respectively. Figure 9 shows the sources maps obtained by the DCS-SOMP method and the DCS-SOMP-SVD method as the target sparsity equals 10. When the target sparsity equals 10, there are six extra sources on the map obtained by the DCS-SOMP method shown in Figure 9(a), which has a bad effect on determining the location of the actual (2020) 33:95 source. Figure 9(b) shows that the whole sources can be located exactly by the DCS-SOMP-SVD method when the target sparsity equals 10. Besides, there are no extra sources on the map.

Experimental Results and Analysis
In order to verify the feasibility of the DCS-SOMP-SVD method in the real application, we conducted an experiment of gas leakage for wideband acoustic imaging at the Northwestern Polytechnical University.

Experiment Configurations
In the experiment, a planar microphone array composed of 24 microphone sensors shown in Figure 10 was used to measure the audio data. Figure 11 shows the experimental configuration taken by a camera installed in the center of the array. The microphone array is 0.7 m away from the source plane. A 10 s recordings were stored for each measurement with sampling frequency being 44.1 kHz. Thus we can obtain many time-domain sampling data for each microphone. The sampling data were divided into 50% overlapping data block, where every data block contains 1024 sampling points. We performed FFT after applying the Hanning window to each data block. Then we can obtain the measurement data of microphone array in frequency domain. The observation zone, which is 0.5708 × 0.4280 m 2 , is discretized into N = 61 × 46 grids.
The spectrum of measurement data is shown in Figure 12, where we can find that the frequency of sources locates at low and middle frequency range mainly. Due to (2020) 33:95 the long calculation time of the wideband BP method, the experimental results of the method will not be shown.
The experimental source maps of the wideband DCS-SOMP method and DCS-SOMP-SVD method were obtained as the frequency ranges equal 3500-4000 Hz, 6500-7000 Hz, and 11000-11500 Hz, which were shown in Figure 13, Figure 14, and Figure 15.   Figure 13 shows the source maps obtained by the DCS-SOMP method and the DCS-SOMP-SVD method when the frequency range equals 3500-4000 Hz. It can be seen from Figure 12 that there is a peak near 4000 Hz, so the SNR of the experiment is relatively high. In this frequency range, both the DCS-SOMP method and the DCS-SOMP-SVD method can locate two sources. Figure 14 shows the source maps obtained by the DCS-SOMP method and the DCS-SOMP-SVD method when the frequency range equals 6500-7000 Hz. When the frequency range is 6500-7000 Hz, the amplitude of the signal is further reduced from Figure 15. In this case, the DCS-SOMP method can only locate the left source and fails to locate the right one,  (2020) 33:95 while our DCS-SOMP-SVD method can still locate two sources. Figure 15 shows the source maps obtained by the DCS-SOMP method and the DCS-SOMP-SVD method when the frequency range equals 11000-11500 Hz. As the frequency increases, the amplitude of the signal gradually decreases. When the frequency range is 11000-11500 Hz, the DCS-SOMP method can only locate the right source, and the DCS-SOMP-SVD method can locate two sources. Therefore, we conclude that our method has better localization performance under different frequency ranges than the DCS-SOMP method.  In this paper, we have proposed a DCS-SOMP-SVD method for wideband acoustic imaging. The performance of the proposed method has been studied through both simulations and experiments. We have also compared the source maps obtained by the DCS-SOMP-SVD method with those obtained by the wideband BP method and the DCS-SOMP method. The main conclusions are as follows.

Source Maps Using Experiment Results of the Gas Leakage
(1) The simulation results show that the DCS-SOMP-SVD method, as well as the wideband BP method and the DCS-SOMP method, can locate all sources as the SNR equals 30 dB. The CPU-time of the DCS-SOMP-SVD method is 2.463 s, which is a little slower than the DCS-SOMP method, but is far less than the CPU-time 126.5 s of the wideband BP method. (2) However, when the SNR decreases to −20 dB, the DCS-SOMP-SVD method can still locate all sources accurately, while the DCS-SOMP method and the wideband BP method have failed to locate any sources. (5) In the gas leak experiment, the results show that the DCS-SOMP-SVD method can locate all leak sources in 3500-4000 Hz, 6500-7000 Hz, and 11000-11500 Hz, while the DCS-SOMP method can only locate the leak source when the SNR is high.