Passive Source Localization Using Compressive Sensing

This paper presents an underwater passive source localization method by forming an underdetermined linear inversion problem. The signal strength on a specified grid is evaluated using sparse reconstruction algorithms by exploiting the spatial sparsity of the source signals. Our strategy leads to a high ratio of measurements to sparsity (RMS), an increase in the peak sharpness with a low side lobe level, and minimization of the dimensionality of the problem due to the formulation of the system equation of the multiple snapshots based on the data correlation matrix. Furthermore, to reduce the computational burden, pre-locating with Bartlett is presented. Our proposed technique can perform close to Bartlet and white noise gain constraint processes in the single-source scenario, but it can give slightly better results while localizing multiple sources. It exhibits the respective characteristics of traditionally used Bartlett and white noise gain constraint methods, such as robustness to environmental/system mismatch and high resolution. Both the simulated and experimental data are processed to demonstrate the effectiveness of the method for underwater source localization.


Introduction
Matched-field processing (MFP) has been extensively studied for passive source localization [1]. It works by scanning through spatial grids of potential targets while looking for the best matching of the observed data and the sound propagation physical model in an ocean waveguide associated with the scanning target position. In regards to the performance metrics, such as resolution and robustness, existing MFP methods have their own pros and cons [1][2][3]. For example, conventional (Bartlett) MFP is relatively robust, however suffering from some fundamental limitations such as low resolution and a high side lobe level. On the other hand, methods such as Minimum Variance Distortionless Response (MVDR) offer a high resolution and low side lobe level [4,5], however often showing significant sensitivity to environmental/system mismatch; hence, when the assumed and true models differ [6,7], the methods frequently give biased estimates [8].
In many practical situations, the number of source signals is rather small at any time instant compared with the number of scanning grids. On that basis, we can exploit the spatial sparsity feature of the underwater source signals and develop a high-resolution MFP method based on sparse reconstruction algorithms, which have been attracting great interest due to their relationship with Compressive Sensing (CS) [9]. According to the CS theory, the spatial sparsity of a signal 1. The sources are stationary random processes and uncorrelated with each other. A number of independent snapshots are available. 2. Both the signal and noise processes follow a zero-mean complex Gaussian distribution, and the noise is independent of the signal processes. 3. The source motion effect is ignored, and random environmental variation is stationary.
We assume a network of M sensors and Q sources located at r q , q = 1, 2, · · · , Q, and these sources emit a narrow-band signal at frequency f . Suppose we have L independent measurements; the complex envelope of the received signal by M number of sensors at the l th snapshot, y l ( f ), can be expressed as [41,42]: where g f , r q is an M × 1 vector of Green's function from the source located at r q to the receiver array; b l,q ( f ) is a random process of the q th source incorporating amplitude and phase variability; n l ( f ) is a stationary noise vector. For each sample, the received signal can be written in a matrix equation: where G ( f ) is composed of Green's function corresponding to each source location as its columns: , g ( f , r 2 ) , · · · , g f , r Q ; the source signals group into a vector: Collection of all the snapshots results in: where the columns of Y indicate different time samples: Y = [y 1 ( f ) , y 2 ( f ) , · · · , y L ( f )], and B and N are defined similarly. Thus, the goal of underwater source localization is to determine the columns of matrix G, which is related to the source locations, based on the known measurements Y and unknown signals B. Several MFP algorithms either explicitly or implicitly involve the formation of the sample covariance matrixK Y , given by: Some adaptive MFP algorithms need enough snapshots to makeK Y full rank. However, in the time-varying channel and high-speed source cases, the number of snapshots available for accumulatinĝ K Y is limited [6,7]. Thus, diagonal loading can be applied, and the sample covariance matrix is replaced by: where µ is a small constant [1].

The Traditional MFP Methods
The Bartlett MFP correlates the data vector with the replica vector and scans all the search region to find the maximum output [1]. The output of Bartlett MFP can be written as: where the replica vector w c (r) is the normalized Green's function vector, It is well known that Bartlett MFP is relatively easy to implement, ant-interference, and not sensitive to environmental mismatch, but it suffers from a broad main lobe and high sidelobe level, which results in low resolution. The narrow main lobe and sidelobe suppression are important if more than one source is present. Later, we use Bartlett MFP as the pre-treatment of other methods.
The WNC method is a high-resolution method, as it suppresses sidelobes while offering a degree of robustness in frequently-encountered system/environmental mismatch scenarios. The WNC is a potential standard as, unlike MVDR, it is more robust to system/environmental mismatch [1,4]. Its output is given by: where: The adaptive weights w n correspond to diagonally-loaded MVDR weights and are obtained by solving: To satisfy a white noise gain constraint G W NG , this iterative algorithm finds the diagonal loading γ for each replica vector such that: In practice, we normalize it and express it as 10 log 10 (κ 2 /M) ≤ 0 dB [43]. The WNC process is a better choice for baseline comparison because it has been shown in many articles that the MVDR has a better resolution than Bartlett MFP, but it is sensitive to environmental/system mismatch and needs enough snapshots to makeK Y full rank [2,3,5].

Sparse Reconstruction-Based MFP with a Single Snapshot
In this section, we start to formulate the MFP problem as a sparse representation problem. For simplicity, we describe the sparse reconstruction-based MFP technique for a single time sample in this section, then the approach will be extended to the multiple snapshot case in the next section.
Instead of estimating the source locations directly, we discretize the region containing the sources with an N-point grid, {r 1 , r 2 , · · · , r N }, and assume that each source is located on one of the grid points. The signal strength on each grid point forms the unknown vector b CS to be solved, where the n th elementb n is nonzero if some source comes from r n and zero otherwise. In practice, there are usually very few sources within the interested area, with the number much smaller than the grid points. Thus, we have reached a sparse representation of the spatial distribution of sources b CS . We construct a replica matrix based on the N-point grid, where g ( f , r n ) is Green's function from a test source located at r n from the receiver array. Since the number of grid points is always larger than the measurements (that is N M), G CS is overcomplete. In this framework, the measurement matrix is known and does not depend on the actual source locations.
For a single time sample, we reach the measuring equation: Different from (2), the columns of G CS and the rows of b CS are indexed by grid positions {r 1 , r 2 , · · · , r N }. Therefore, the MFP problem is translated into a sparse signal reconstruction problem, which is to determine the source strength on each grid point b CS based on the measurements y and known measurement matrix G CS .
In this case, Equation (12) is underdetermined and does not have a unique solution. Recent developments in the CS theory have shown that a Q-sparse N-dimensional signal (only has Q non-zeros entries) can be uniquely reconstructed by O (Q log (N/Q )) incoherent linear measurements [44]. Therefore, Equation (12) can be solved through sparse reconstruction algorithms, such as Basis Pursuit (BP) [45][46][47], Matching Pursuit (MP) [45,48], and Compressive Sampling Matching Pursuit (CoSaMP) [49]. In this paper, we adopt the BP algorithm, which has been shown to be near optimal and robust to noise and to coherent dictionary entries as well [43]. Thus, in the perspective of spatial sparse constraint, the MFP problem is formulated as: The regularization parameter ε is introduced to tolerate or penalize the noise in the measurements.
ε should be set appropriately: if ε is too small, it will not be able to tolerate noise, causing several peaks in the output; if ε is too large, there will be no penalty for noise, always resulting in a zero solution. In (13), the l 1 -norm is imposed on b CS to emphasize sparsity in space; thus, the sparse-reconstruction-based MFP can achieve high resolution and suppress the side lobes effectively. RMS can be set as a performance indicator of sparse reconstruction-based MFP methods. The higher the RMS, the better the performance of sparse reconstruction [50][51][52]. The RMS of (12) is M/Q (M measurements for Q sparsity).

Sparse Reconstruction-Based MFP with Multiple Snapshots
For multiple snapshots, we can reach: where The matrix B CS is parameterized temporally and spatially, but it is only sparse in space, with QL non-zeros. The numerical solution of (14) is a bit more involved than that of a single sample case, since we should solve the problem in the N × L-dimensional complex space.

Related Work
Malioutov et al. provided some ideas on solving (14) and reducing the problem's dimensionality [11]. The first thought is separate treatment of each snapshot, thus reducing the problem dimensionality to N, and finally averaging the results across all the snapshots. This approach is simple, but it has low performance and is sensitive to noise. The second approach is to treat different time samples in synergy. Since matrix B CS is only sparse in space, one can compute the l 2 -norm of each row of B CS , given by b (l 2 ) , which is a Q-sparse N-dimensional signal. Penalize the l 1 -norm of b (l 2 ) , and the cost function becomes: This treatment can achieve better performance than treating each snapshot separately, but the computational cost for BP is O (N × L) 3 , which is greatly increased by the number of snapshots. Malioutov et al. proposed a Singular-Value Decomposition (SVD) method to reduce the dimensionality. They took the SVD of the measurement data matrix Y = UΣV H and kept the components in signal subspace, which resulted in an M × Q matrix Y SV = YV Q , where V Q is the first Q columns of V corresponding to the Q largest singular values. Thus, the measurement equation is transformed to: where B SV = BV Q and N SV = NV Q . The problem dimensionality is reduced to N × Q; however, the RMS is also retained as M/Q (M × Q measurements for Q 2 sparsity), and it requires sources to be sufficiently separated. We consider this method to be closely related to our work. It is also a good candidate for the comparison with our proposed method. We call this the 1-SVD method in the rest of this paper. Mantzel et al. also proposed CS imaging for MFP named as cMFP, where it focuses more on the computational complexity reduction rather than the imaging model and recovery performance [31]. The method proposed in our work mainly focuses on performance and partially on the reduction of the computational burden. We draw a comparison of this method with our proposed method in terms of performance, and the results show that our method has clear advantages in terms of the localization of the source and main lobe to sidelobe ratio.
where Φ is an M × N encoding matrix. The question remains as to how to choose the encoding matrix Φ so that the solution to this method is similar to the standard MFP, e.g., Bartlett. To reduce the computational burden, the construction of Φ is proposed to be a random linear mapping. The performance of the cMFP will match that of the standard MFP when Φ maintains the energy of the differences between the observations Y and all scalar multiples of Green's function at different points. This method is used for comparison, and it is named cMFP in this paper.

CS-R Method
In this section, we present a new approach, which can extend the use of sparse signal representation ideas to practical underwater source localization problems. It is convenient to write the data correlation matrix of (14) as: In the case of white noise, there isK N =K N − σ 2 n I (where σ 2 n is the noise variance), and hence, all the entries inK N ,K GBN andK NGB are zero mean. Since the sources are uncorrelated with each other,K B only has large entries along its diagonal. Therefore, Equation (18) can be approximated as: where: where vec (·) denotes the vectorization operator (stacking the columns of a matrix on top of each other). It is obvious that d is an N-dimensional real vector with Q sparsity. Since the number of measurements is increased, while the sparsity is retained at the minimum Q, the RMS of (19) is increased.
Since the rows of the measurement matrix A may be linearly dependent and the dimensionality of A can be quite high, SVD is used to reduce the measurement space and increase the SNR. Taking the SVD of A, A = UΣV H , reserving theM largest singular values λ m , m = 1, 2, · · · ,M, and projecting the measurements onto the corresponding output basis vectors, we derive the following equation: whereΣ is the diagonal matrix consisting of theM largest singular values andŪ,V denote the correspondingM columns of U and V, respectively. Generally, after the processing of (21), almost all the signal power is retained, whereas there is only a part of the noise power left. Thus, SNR increases, which will benefit the latter sparse reconstruction process. Furthermore, the problem dimensionality is decreased toM × Q. Equation (21) can be solved through: The schematic diagram of implementing the CS-R method is given in Figure 1.

Restricted Isometry Property for Sparse Recovery
The Restricted Isometry Property (RIP) for sparse recovery is related to the orthonormality of the columns of the A matrix. The RIP condition is a necessary, though not sufficient, condition for signal recovery [46]. Let the sparsity Q of a vector be the number of nonzero elements, and define the RIP constant δ Q < 1 to be the smallest positive number such that the inequality: for every vector d which is Q-sparse. Here, we can define another quantity, the Q,Q-restricted orthogonality constant θ Q,Q , as the smallest number that satisfies: for all d andd such that d andd are Q-sparse andQ-sparse, respectively, and have disjoint supports. The constants δ Q and θ Q,Q are related by the following inequalities, A sensing matrix is said to satisfy the 2Q-RIP if δ 2Q < √ 2 − 1 [52]. In this work, the 2Q-RIP condition is always satisfied for a noiseless environment, where the sparse target is Q = 1. The condition in the presence of bounded error is also fulfilled; the value of 0.98 is incurred, where that condition is given as δ 2Q + θ Q,2Q < 1 [53]. A direct consequence of (23) is that δ 2Q < √ 2 − 1 is in fact a strictly stronger condition than δ 2Q + θ Q,2Q < 1 [54].

Comparison between the CS-R Method and the Bartlett Method
The output of Bartlett is B Bart For the system equation of CS-R (18), ignoring the noiseK N +K GBN +K NGB , we can simplify the system equation aŝ Then,K Y can be expressed as: where Q = diag ( g (r 1 ) 2 , g (r 2 ) 2 , · · · , g (r N ) 2 ). G CS Q −1 denotes normalizing the Green's function of each grid point; the n th column of G CS Q −1 is . Through Least-Squares (LS), we derive the following equations: where F = G CS Q −1 H G CS Q −1 −1 . Furthermore, we obtain: and:K where G CS Q −1 HK Y G CS Q −1 is the objective equation of Bartlett MFP. CS-R translates the system equation into solving (22), wherein the l 1 -norm is used for the penalty of spatial sparsity, ensuring that CS-R can achieve high resolution and suppress side lobes effectively; the system equation of the l 2 -norm can be seen as equivalent to the objective equation of Bartlett MFP, meaning that CS-R can provide robust estimation similar to that of Bartlett.
Compared to the other sparse-reconstruction-based MFP methods, CS-R has three added advantages: the problem dimensionality is reduced to the minimum, and the unknown signal is real, so the computational cost is lower; there are high RMS gains, with better reconstruction performance; robustness, which is rather important for underwater source localization.

Pre-Locating
Since the MFP problem in our framework is performed on a grid, the accuracy depends on the grid size. However, it is computationally expensive to make the grid uniformly fine on the whole region. Malioutov et al. introduced the idea of adaptive grid refinement, in which the grid is refined around the source locations derived by the previous estimation [11]. Inspired by this, if we have obtained approximate knowledge of source locations using low resolution, but simple and robust methods, such as Bartlett MFP, the grid for sparse reconstruction-based MFP can be specified to cover only the regions where sources are located with a high probability, as shown in Figure 2 [10]. Multiple sub-regions may be selected for multiple targets' localization. These regions can be treated as a whole, reconstructing the multi-targets at the same time; whereas the isolated regions can also be treated separately, largely reducing the computational load of sparse reconstruction algorithms. Pre-locating has some limitations while locating multiple sources, especially closely located sources or detecting a weak submerged target in the presence of a loud surface interferer. We observe that pre-locating works fine in the single-source scenario; hence, on a prior knowledge basis, we can skip pre-locating in the multi-source scenario.

Pre-locating with Bartlett MFP
Grid for sparse-reconstruction-based MFP

Simulation Results
Similar to Bartlett, all those sparse reconstruction-based MFP methods do not require the sample covariance matrix to be full rank; thus, they also do not need a large number of snapshots. Our simulation and experimental data processing will focus on the snapshot deficient case, that is L < M. In this section, we will use numerical simulation to verify that the CS-R method can give results similar to other existing adaptive processes, namely Bartlett, WNC, and 1-SVD. Firstly, we compare the ambiguity surface of CS-R to Bartlett MFP. Next, the source location estimation is performed for CS-R in comparison with Bartlett, WNC, and 1-SVD, followed by the performance analysis with environmental uncertainty.
The simulation scenario is shown in Figure 3. The water depth is 102 m, and the acoustic field is received at an array of M = 40 sensors spanning the water column. The source, located at the depth and range of 20 m and 12 km, respectively, emits a narrow-band signal at a frequency of 630 Hz. We further assumed that twenty independent snapshots are available for each estimate (L = 20), and SNR = 5 dB. The location of the peak on each ambiguity surface is shown by a circle in the simulation results.

CS-R vs. Bartlett
The results of Bartlett MFP and CS-R are shown in Figure 4, wherein CS-R can trade off between the l 1 -norm and l 2 -norm through the regularization parameter ε. The top row in Figure 4 shows that the Bartlett MFP has low resolution, and its side lobe level is rather high. As shown in Figure 4e,f, a good choice of ε is two; the ambiguity surface of CS-R has a sharp peak, and the peak corresponds to the source location. When ε is 0.3 times 2, many peaks appear, shown in the Figure 4c,d; when ε is set at a higher value, 3 times 2 for example, CS-R degrades to the Bartlett MFP, and its side lobe level grows and main lobe broadens as shown in Figure 4g,h. Besides, the ambiguity surfaces of Bartlett and CS-R have the same multi-modal structure. The example illustrates two points. First, the choice of the parameter ε is very important; when the parameter ε is too large, CS-R degrades into Bartlett MFP, which is because the l 2 -norm of CS-R can be equivalent to Bartlett. In this case, an optimum value of the parameter ε is estimated to achieve a trade-off between high resolution and robustness.
The choice of the regularization parameter ε plays an important role in our source localization framework, which balances the fit of the solution to the data versus the sparsity prior. This question arises in many practical inverse problems, and the parameter selection can be made if the noise statistics are known or can be estimated. The idea is to select ε for the residuals of the solutiond(ε) for some known statistics of the noise when such statistics are available. If the distribution of the noise is known or can be estimated [55], then the regularization parameter can be calculated such

Source Location Estimation Comparison in a Complex Environment
This simulation is performed using the environment described in Figure 3 for both depth and range estimation of the underwater source location estimation. The results presented in Figure 5 indicate that the Bartlett surface does not have a prominent peak, but provides the right peak position. WNC and 1-SVD provide the right peak position with a sharper surface as compared to Bartlett. cMFP performs closely to Bartlett using random received vector from the source, but a slightly broader main lobe than Bartlett. On the other hand, CS-R achieves a higher resolution than the WNC and 1-SVD with the correct peak position and low side lobe level, which indicates its robustness and high-resolution capabilities. Using the same simulation environment, two sources are localized in Figure 6. In Figure 6a, the locations of two sources, which have the same depth and are 0.4 km apart in range, are estimated; while Source 1 has a higher SNR (30 dB) and Source 2 exhibits weaker SNR (8 dB). It is observed that Bartlett and cMFP show many peaks higher than the Source 2 level. On the other hand, all other processes find the right peak location, and CS-R has lower sidelobes. Note that pre-locating is not implemented in the multi-source scenario for CS-R. It can be seen in Figure 6b that all processes suffer from a high side-lobe level when there is a 1% mismatch. Bartlett and cMFP again show many false peak locations other than sources. The WNC and CS-R processes perform better than 1-SVD.

Performance Analysis with Environmental Mismatch
Performance analysis using a simulated environmental model can provide important guidance for practical applications. We focus on performance analysis with environmental uncertainty here. While the respective true sound speeds in the sediment are 1579 m/s on the top and 1629 m/s on the bottom, we assume they are 1595 m/s on the top and 1645 m/s on the bottom instead to create a case of some environmental mismatch deliberately. Performance analysis in source range estimation is evaluated to compare Bartlett, WNC, 1-SVD, CS-R, and cMFP. Results are shown in Figure 7, which are derived based on 5000 Monte Carlo simulations with 20 snapshots for every run.
1. In Figure 7a, we present plots of the the Mean Square Error (MSE) of each method versus SNR.
The search for the source range location is between 5 and 15 km, and the depth is between 1 and 102 m. 2. Figure 7b gives the evaluation results of the probability of localization error (Pe) versus SNR.
The search for range is between 5 and 15 km, and the depth is between 1 and 102 m.
It is found in Figure 7 that the performance of CS-R and cMFP is much closer to Bartlett than either WNC and 1-SVD, following closely Bartlett in the whole SNR the environmental mismatch scenario for the single source, showing robustness similar to the Bartlett process in the single-source scenario. The performance comparison for source localization is investigated with environmental mismatch while having two sources, and the results are shown in Figure 8. Results of Bartlett, cMFP, and 1-SVD of this investigation were relatively poor; Bartlett and cMFP were unable to estimate the correct location of the second source due to its broad main lobe and high sidelobe level, and 1-SVD's performance was similar to Bartlett due to its high sensitivity to mismatch. Figure 8b also compares the probability of error for the multi-source scenario versus SNR. The WNC and CS-R processes performed much better, e.g., at SNR = 5 dB, the probability of the error of the localization of second source was very low for WNC and CS-R. The performance of CS-R was slightly better in the two-source scenario than WNC, as shown in Figure 8.

Peak to Background Ratio Analysis
To demonstrate the high-resolution capabilities of CS-R in underwater source location estimation, which include a sharp main lobe level and a low side lobe level, the Peak to Background Ratio (PBR) analysis was performed.
PBR of CS-R was compared with WNC, Bartlett, cMFP, and 1-SVD while changing different parameters, as shown in Figure 9. The obtained results showed the high resolution advantage of CS-R over other conventional methods in different challenging scenarios, e.g., snapshot deficient and low SNR. Although all processes can perform in the snapshot-deficient scenario, the purpose here is to show that, unlike some other methods, e.g., MVDR, our proposed method can also perform in those kinds of scenarios. We used 20 snapshots in Figure 9a and 20 dB in Figure 9b. We observe in Figure 9 that CS-R achieves PBR better than Bartlett and WNC and is close to 1-SVD. CS-R only struggled when the SNR was low, as shown in Figure 9a, but still exhibited slightly better results than Bartlett and WNC. Equation (30) below is used to calculate the PBR of different source location estimation methods [56].
where P A is the total integrated area within peak boundaries minus the background area and B A is the rectangular area of background height times the length between peak boundaries. Note that the background height is defined as the minimum intensity at peak boundaries, and peak height is defined as the maximum intensity between peak boundaries minus background height.

Experimental Results
The Asian Seas International Acoustics Experiments (ASIAEX) were carried out in the spring of 2001 in the East China Sea (ECS) [57]. Here, we used part of the data from those experiments for our analysis. The site location is shown in Figure 10. The scenario for the shallow water experiment is shown in Figure 11. The water depth was about 102 m, and the acoustic field was received at a vertical line array consisting of 60 sensors. The source signal frequency was 630 Hz. Conductance Temperature Depth (CTD) sensors were used to calibrate the depth of the acoustic sensors and measure the Sound Speed Profile (SSP) of the water column. The SNR was about 15 dB. The range (determined by GPS coordinates) and depth (prior knowledge) of the source at the time of experiment (see Figure 12 at about 15 hours of time) was about 11.6 km and 40 m, respectively.
After removing noisy channels, we selected 42 sensors and 30 independent snapshots in utilizing experimental data; thus, the dimensionality of the received signal Y was 42 × 30, which is a snapshot-deficient scenario in the statisticalsense. The search interval for the source location was between 11 and 13 km and 1 and 102 m for range and depth, respectively. Pre-locating was implemented for CS-R, and the discretization of the search grid for CS-R remained the same for all other processes. Figure 13 summarizes the ambiguity surfaces of the Bartlett, WNC, 1-SVD, cMFP, and CS-R processes, respectively.    We observe in Figure 13 that the main lobe of the Bartlett process is broad, and the peak to side lobe ratio is rather low. However, it is most likely to provide an optimum peak location. The cMFP's results show that it focuses more on the computational complexity reduction rather than the imaging model and recovery performance; therefore, it suffers from a broad main lobe and high side lobe level. The WNC provides a much sharper peak closer to the true source location, but 1-SVD fails to provide the right peak location due to its sensitivity to environmental/system mismatch. It can be seen that the CS-R achieves a slightly sharper spectrum than WNC, but it also has the same peak corresponding to the source location as Bartlett and WNC.
To demonstrate the properties of our proposed method, we constructed the multi-source scenario using real data. A time-delayed copy of the ASIAEX was added to the original ASIAEX to create a distance of about 700 m between two sources. In Figure 14, the results of two-source localization using experimental data are shown. The first source was located at a range of about 11.6 km and a depth of about 40 m, while the second source was located at a range of about 12.3 km and a depth of also about 40 m. Thirty time snapshots from the vertical line array of 42 sensors were used. The SNR for the first and second sources were 15 dB and 8 dB, respectively. The two source results demonstrate that WNC and CS-R performed closely to each other, and CS-R can achieved slightly better results, as shown in the Figure 14. The Bartlett and cMFP processes were unable to show the weaker source's true location due to the high sidelobe level. On the other hand, 1-SVD remained unable to find the true locations of both sources. In future, the reported technique can be used for different scenarios, i.e., passive detection and tracking of an underwater acoustic target [58].

Conclusions
In this paper, we considered the underwater source localization problems and explored a high-resolution and robust CS-R method based on sparse reconstruction algorithms for the matched field processing application. In our framework, the source localization problem was formulated to solve an underdetermined equation by discretizing the source space into grids. Then, in order to handle the multiple snapshots, to reduce the signal space dimensionality, and to increase RMS, we established the measurement equation according to the data correlation matrix. CS-R uses the l 1 -norm to penalize the spatial sparsity and the l 2 -norm to penalize the error of the measurement equation. Pre-locating with the Bartlett processes was also proposed to reduce the computational load for the single-source scenario. The results of the numerical simulations suggested that the CS-R method allowed achieving a trade-off between high resolution and robustness through the choice of the regularization parameter. The experimental data were used to show that the performance characteristics of our proposed technique were similar to the adaptive process. CS-R exhibited high resolution capabilities and robustness close to both WNC and Bartlett and showed better performance than 1-SVD and cMFP. Both the single source and multi-sources scenario were used to investigate the performance of different processes. In the single-source scenario, the Bartlett process was more robust, and CS-R could achieve results close to it, but at the cost of increased computational load. In the two-source scenario, the CS-R's performance was slightly better than WNC. As future work, the performance of CS-R in data replica mismatch can be improved.