# The Sparsity Adaptive Reconstruction Algorithm Based on Simulated Annealing for Compressed Sensing

Authors: Yangyang Li; Jianping Zhang; Guiling Sun; Dongxue Lu
Journal: Journal of Electrical and Computer Engineering
 (2019)

---

## Abstract

This paper proposes a novel sparsity adaptive simulated annealing algorithm to solve the issue of sparse recovery. This algorithm combines the advantage of the sparsity adaptive matching pursuit (SAMP) algorithm and the simulated annealing method in global searching for the recovery of the sparse signal. First, we calculate the sparsity and the initial support collection as the initial search points of the proposed optimization algorithm by using the idea of SAMP. Then, we design a two-cycle reconstruction method to find the support sets efficiently and accurately by updating the optimization direction. Finally, we take advantage of the sparsity adaptive simulated annealing algorithm in global optimization to guide the sparse reconstruction. The proposed sparsity adaptive greedy pursuit model has a simple geometric structure, it can get the global optimal solution, and it is better than the greedy algorithm in terms of recovery quality. Our experimental results validate that the proposed algorithm outperforms existing state-of-the-art sparse reconstruction algorithms.


---

## Body

## 1. Introduction

In recent years, the research of compressed sensing (CS) [1–3] has received more attention in many fields. CS is a novel signal processing theory, which can efficiently reconstruct a signal by finding solutions to underdetermined linear equations [4]. It is very magical that CS theory can exactly restore the sparse signal or compressible signals by a sampling rate which does not satisfy the Nyquist–Shannon sampling theorem. And many papers have demonstrated that CS can effectively get the key information from a few noncorrelative measurements.The core idea of CS is that sampling and compression are done simultaneously. CS technology can reduce the hardware requirements, further reduce the sampling rate, improve the signal restoration quality, and save the cost of signal processing and transmission. Currently, CS has been widely used in wireless sensor networks [5, 6], information theory [7], image processing [8–10], earth science, optical/microwave imaging, pattern recognition [11], wireless communications [12, 13], atmosphere, geology, and other fields.CS theory is mainly divided into three aspects: (1) sparse representation; (2) uncorrelated sampling; (3) sparse reconstruction. However, how to reconstruct sparse signals is crucial. It is a huge challenge for the researchers to propose a reconstruction algorithm with reliable accuracy. Since reconstruction methods need to recover the original signal from the low-dimensional measurements, the signal reconstruction requires solving an underdetermined equation. That may include infinite solutions. Mathematically, the issue can be regarded as thel0-norm minimization. However, it is an NP-hard problem [14]. It requires exhaustively listing all possible solutions of the original signal. That is unrealistic for large-scale data processing problems.At present, some scholars have proposed many reconstruction algorithms. Mostly, the nonconvexl0-norm minimization will be replaced with the convex l1-norm minimization to design the reconstruction model. There are two important kinds of reconstruction algorithms. One is greedy pursuit algorithms, and the other is the convex optimization algorithm. The convex optimization algorithm solves a much easier l1-norm minimization problem (e.g., basis pursuit algorithm). But it requires high computational complexity to achieve exact restoration. However, a series of iterative greedy algorithms receive great attention due to their low computational complexity and simple geometric interpretation. The greedy algorithms include the orthogonal matching pursuit (OMP) [15], the stagewise OMP (StOMP) [4], the regularized OMP (ROMP) [16], the subspace pursuit (SP) [17], the compressive sampling OMP (CoSaMP) [18], and so on. The key step of these methods is to estimate the support sets with high probability. At each iteration, one or several indexes of the sparse signal are added into the support sets for subsequent processing. If those indexes are reliable, they will be added to the final support sets to recover the original signal.In this paper, we propose a novel sparsity adaptive simulated annealing method (SASA) to reconstruct the sparse signal efficiently. This new method combines the advantage of the SAMP [19] algorithm and the simulated annealing (SA) algorithm to solve the CS reconstruction problem. The simulated annealing algorithm is known for its efficient global optimization strategies and superior results in solving nonconvex problems. The proposed algorithm can estimate the sparsity level, and it achieves superior performance in finding the optimization solution. The experimental results have shown that the proposed SASA method outperforms existing state-of-the-art sparse reconstruction algorithms.The main contributions of this paper are listed as follows:(1)
We take advantage of the sparsity adaptive simulated annealing algorithm in global searching to guide the sparse reconstruction(2)
The proposed algorithm can reconstruct a sparse signal accurately, and it does not require sparsity as a priori information(3)
The proposed algorithm has both the simple geometric structure of the greedy algorithm and the global optimization ability of the SA algorithmThe remainder of this paper is organized as follows. Section2 introduces the background of CS and some approximation methods that may be applied in the later section. In Section 3, we propose the SASA algorithm and give the algorithm’s frameworks. We discuss the experimental results in Section 4. Finally, Section 5 draws the conclusion.


## 2. Background and Related Work

In this section, we introduce the compressive sensing model and the sparse signal approximation method.


### 2.1. Compressive Sensing Model

Sparse reconstruction is to recover the high dimensional original signalxx∈RN×1. yy∈RM×1 is the measurement vector. The compressive sensing model can be defined as(1)y=Φx,where Φ is the sensing matrix Φ∈RM×N, M≪N. Mathematically, this problem can be solved by the l0-norm minimization. It can be written as follows:(2)minxx0,s.t.y=Φx,where ∗0 is the zero norm. However, the l0-norm minimization is the NP-hard problem. It needs to exhaustively list all CNK (K is the sparsity) possible solutions for determining the locations of the nonzero entries in x. Therefore, it is usually converted to the convex l1-norm minimization problem as follows:(3)minxx1,s.t.y=Φx.The above model is a convex optimization problem. It can resort to linear programming methods or Bregman split algorithm. And linear programming methods have shown to be effective in solving such problems with high accuracy. Theoretically, the CS has explained that the sparse signalx can be exactly reconstructed from only a few observation y when the sensing matrix Φ satisfies the restricted isometry property (RIP). This is an important property for us to analyze and design the reconstruction algorithm. The RIP is given as follows:(4)1−δsx22≤Φx22≤1+δsx22,where δs is a constant 0≤δs<1. Denote x22 as the l2 norm of a vector x.Generally speaking, ifδS is very close to one, then it is possible that y cannot preserve any information on x (i.e., Φx22≈0, x is the null space of the sensing matrix Φ). In this case, it is nearly impossible to reconstruct the sparse signal x by using any greedy algorithms.There are many kinds of sensing matrixΦ. It is mainly divided into deterministic matrices and learning matrices. Based on signal characteristics, the learning matrix can be continuously updated by using a machine learning method. However, the commonly used sensing matrices include the Gaussian random matrix, Fourier random matrix, and Toeplitz matrix.


### 2.2. Sparse Signal Approximation Method

In this section, we give some notations and matrix operators to facilitate the summary of the proposed algorithm.Definition 1 (projection and residue).
Lety∈RM×1 and ΦS∈RM×S. If Φ is the full column rank, then Φ′Φ−1Φ′ is the Moore–Penrose pseudoinverse of Φ, and Φ′ represents the matrix transpose. The projection of the measurements y onto span ΦS is written as(5)yp=ΦSΦS†y,where ΦS† is the pseudoinverse of the matrix ΦS and is defined as(6)ΦS†=ΦS′ΦS−1ΦS′.
It can be seen from equation (5) that we can get the restoration of the sparse signal x by using the following formulation:(7)x^=ΦS†y=ΦS′ΦS−1ΦS′y.
The residue of the projection vectory is expressed as(8)yr=E−ΦSΦS†y.
Suppose thatΦS′Φs is invertible. The projection of y onto span (ΦS) is defined as(9)yp=ΦSΦS†y.
E is the unit matrix. To illustrate the relationship between projection and residue in iteration, we draw Figure 1 [17]. Refer to Figure 1 to help better understand the definition.
As can be seen from Figure1, the residue is perpendicular to the hyperplane which can be ensured by using equation (7). It means that the residue will be minimized in each iteration. And, that can speed up the convergence of the algorithm.Figure 1
In each iteration, aT-dimensional hyperplane closer to y is obtained.


## 2.1. Compressive Sensing Model

Sparse reconstruction is to recover the high dimensional original signalxx∈RN×1. yy∈RM×1 is the measurement vector. The compressive sensing model can be defined as(1)y=Φx,where Φ is the sensing matrix Φ∈RM×N, M≪N. Mathematically, this problem can be solved by the l0-norm minimization. It can be written as follows:(2)minxx0,s.t.y=Φx,where ∗0 is the zero norm. However, the l0-norm minimization is the NP-hard problem. It needs to exhaustively list all CNK (K is the sparsity) possible solutions for determining the locations of the nonzero entries in x. Therefore, it is usually converted to the convex l1-norm minimization problem as follows:(3)minxx1,s.t.y=Φx.The above model is a convex optimization problem. It can resort to linear programming methods or Bregman split algorithm. And linear programming methods have shown to be effective in solving such problems with high accuracy. Theoretically, the CS has explained that the sparse signalx can be exactly reconstructed from only a few observation y when the sensing matrix Φ satisfies the restricted isometry property (RIP). This is an important property for us to analyze and design the reconstruction algorithm. The RIP is given as follows:(4)1−δsx22≤Φx22≤1+δsx22,where δs is a constant 0≤δs<1. Denote x22 as the l2 norm of a vector x.Generally speaking, ifδS is very close to one, then it is possible that y cannot preserve any information on x (i.e., Φx22≈0, x is the null space of the sensing matrix Φ). In this case, it is nearly impossible to reconstruct the sparse signal x by using any greedy algorithms.There are many kinds of sensing matrixΦ. It is mainly divided into deterministic matrices and learning matrices. Based on signal characteristics, the learning matrix can be continuously updated by using a machine learning method. However, the commonly used sensing matrices include the Gaussian random matrix, Fourier random matrix, and Toeplitz matrix.


## 2.2. Sparse Signal Approximation Method

In this section, we give some notations and matrix operators to facilitate the summary of the proposed algorithm.Definition 1 (projection and residue).
Lety∈RM×1 and ΦS∈RM×S. If Φ is the full column rank, then Φ′Φ−1Φ′ is the Moore–Penrose pseudoinverse of Φ, and Φ′ represents the matrix transpose. The projection of the measurements y onto span ΦS is written as(5)yp=ΦSΦS†y,where ΦS† is the pseudoinverse of the matrix ΦS and is defined as(6)ΦS†=ΦS′ΦS−1ΦS′.
It can be seen from equation (5) that we can get the restoration of the sparse signal x by using the following formulation:(7)x^=ΦS†y=ΦS′ΦS−1ΦS′y.
The residue of the projection vectory is expressed as(8)yr=E−ΦSΦS†y.
Suppose thatΦS′Φs is invertible. The projection of y onto span (ΦS) is defined as(9)yp=ΦSΦS†y.
E is the unit matrix. To illustrate the relationship between projection and residue in iteration, we draw Figure 1 [17]. Refer to Figure 1 to help better understand the definition.
As can be seen from Figure1, the residue is perpendicular to the hyperplane which can be ensured by using equation (7). It means that the residue will be minimized in each iteration. And, that can speed up the convergence of the algorithm.Figure 1
In each iteration, aT-dimensional hyperplane closer to y is obtained.


## 3. SASA: Sparsity Adaptive Simulated Annealing Matching Pursuit Algorithm

The greedy pursuit algorithms for sparse recovery have low computational complexity and perfect reconstruction quality. Thus, how to design a stable and efficient CS reconstruction algorithm is a meaningful research issue. But some existing greedy algorithms often get a local optimization solution to the reconstruction problem. To better solve this issue, we introduce an intelligent optimization algorithm to get better reconstruction result.Simulated annealing (SA) is a famous heuristic optimization algorithm, and it has been empirically indicated that SA has superior performance in finding global optimization solutions. As a special optimization reconstruction problem, we can directly solve the problem (2) by using the SA algorithm. However, due to the high computational complexity, it is not practical to directly use SA to solve the sparse recovery problem.As we all know, the greedy algorithms have low complexity and simple geometric interpretation. But greedy algorithms do not have global optimization capabilities, and it is easy to get a suboptimal solution. This means that the greedy strategy and SA can complement each other. We combine the advantage of the SAMP algorithm and the SA algorithm to recover a sparse signal from only a few measurementsy.Inspired by this, we propose a new method to solve the sparse reconstruction issue. The proposed algorithm not only does not require sparsity as a necessary prior but also has superior optimization performance. This paper proposes a novel sparse signal recovery algorithm called sparsity adaptive simulated annealing matching pursuit algorithm (SASA).


### 3.1. Brief of the Simulated Annealing Algorithm

Kirkpatrick proposed the SA algorithm in 1983. SA was firstly proposed to solve thermodynamic problems. The SA algorithm can escape from a local optimization solution. Specially, since it can add a few interference items with a certain probability, SA is more likely to find global optimums or approximate optimization solutions. So far, the SA algorithm is still receiving great attention due to its special strategy.The detail procedure of the standard SA is summarized in Algorithm1 [20]. In Algorithm 1, f⋅ is the objective function and NH⋅ is a function that returns a random neighbor solution to neighborhood structure H.Algorithm 1:Simulated annealing algorithm.
(1)
Input:
(2)
Initial solution S;
(3)
Initial the temperaturet0;
(4)
Cooling rateα;
(5)
The outer-loop iterationoutmax;
(6)
The inner-loop iterationinnermax.
(7)
S∗⟵S, t⟵t0
(8)
Iteration:
(9)
Whilei<outmax do
(10)
forj<innermax do
(11)
S′⟵ random neighbor of a random neighborhood k, S′∈NHS
(12)
Δ⟵fS′−fS
(13)
if Δ≤0 then
(14)
S⟵S′
(15)
iffS′<fS∗ then
(16)
S∗⟵S′
(17)
else
(18)
Take a randomx∈0,1
(19)
ifx<e−Δ/t then
(20)
S⟵S′
(21)
t⟵α×t: adjusting temperature
(22)
if t<ε then t⟵t0
(23)
Output: S∗


### 3.2. Proposed SASA

In this part, we introduce the processing of the SASA algorithm in detail. In real signal processing, the sparsity is usually unavailable for restoration task. We cannot know the signal sparsity in advance. If the sparsity levelK can be accurately obtained, we can solve the problem (2) easily.Theorem 1.
Suppose the support setI^ is the unique solution for the ground truth x^ and K<sparkΦ/2. The sparkΦ is the minimum number of columns from Φ, and those columns are linearly dependent. Then, we can get the estimation solution x˜=x^ if and only if ΦI˜ΦI˜†y=y. I˜ is the estimated solution. And, the x˜ can be obtained by the least squares.From theorem 1, we can define the cost function as(10)fI=ΦIΦI†y−y2,where I denotes the support set in each iteration. In the following, we give the pseudocode of the SASA algorithm, and the main steps of signal reconstruction are summarized in Algorithm 2.Algorithm 2:SASA algorithm.
(1)
Input:
(2)
Measurement signaly;
(3)
Sensing matrixΦ;
(4)
Initializeoutmax, innermaxsq, β.
(5)
Iteration:
(6)
Estimated sparsity levelK.
(7)
Initialize theK support set I0.
(8)
Itemp=I0;
(9)
InitializeT0=10∗fI0;
(10)
While i<outmax do
(11)
forj<innermax do
(12)
Iij=Itemp;
(13)
Chooseq elements from C−Iij and combined with Iij as E;
(14)
Calculatexij=ΦE†y;
(15)
Choose theK largest elements from xij as Itemp;
(16)
Calculateδ=fItemp−fIij;
(17)
if δ<0 then
(18)
Itemp=Itemp
(19)
else
(20)
Calculateβ=min1,exp−absδ/T0
(21)
Take a randomra∈0,1
(22)
if β>ra then Itemp=Itemp;
(22)
elseItemp=Iij
(23)
Calculateerror=fItemp
(24)
if error<β then break
(25)
end for
(26)
UpdateT0=T0/logi+1
(27)
End while
(28)
Output:
x
˜
=
Φ
I
temp
†
y;


## 3.1. Brief of the Simulated Annealing Algorithm

Kirkpatrick proposed the SA algorithm in 1983. SA was firstly proposed to solve thermodynamic problems. The SA algorithm can escape from a local optimization solution. Specially, since it can add a few interference items with a certain probability, SA is more likely to find global optimums or approximate optimization solutions. So far, the SA algorithm is still receiving great attention due to its special strategy.The detail procedure of the standard SA is summarized in Algorithm1 [20]. In Algorithm 1, f⋅ is the objective function and NH⋅ is a function that returns a random neighbor solution to neighborhood structure H.Algorithm 1:Simulated annealing algorithm.
(1)
Input:
(2)
Initial solution S;
(3)
Initial the temperaturet0;
(4)
Cooling rateα;
(5)
The outer-loop iterationoutmax;
(6)
The inner-loop iterationinnermax.
(7)
S∗⟵S, t⟵t0
(8)
Iteration:
(9)
Whilei<outmax do
(10)
forj<innermax do
(11)
S′⟵ random neighbor of a random neighborhood k, S′∈NHS
(12)
Δ⟵fS′−fS
(13)
if Δ≤0 then
(14)
S⟵S′
(15)
iffS′<fS∗ then
(16)
S∗⟵S′
(17)
else
(18)
Take a randomx∈0,1
(19)
ifx<e−Δ/t then
(20)
S⟵S′
(21)
t⟵α×t: adjusting temperature
(22)
if t<ε then t⟵t0
(23)
Output: S∗


## 3.2. Proposed SASA

In this part, we introduce the processing of the SASA algorithm in detail. In real signal processing, the sparsity is usually unavailable for restoration task. We cannot know the signal sparsity in advance. If the sparsity levelK can be accurately obtained, we can solve the problem (2) easily.Theorem 1.
Suppose the support setI^ is the unique solution for the ground truth x^ and K<sparkΦ/2. The sparkΦ is the minimum number of columns from Φ, and those columns are linearly dependent. Then, we can get the estimation solution x˜=x^ if and only if ΦI˜ΦI˜†y=y. I˜ is the estimated solution. And, the x˜ can be obtained by the least squares.From theorem 1, we can define the cost function as(10)fI=ΦIΦI†y−y2,where I denotes the support set in each iteration. In the following, we give the pseudocode of the SASA algorithm, and the main steps of signal reconstruction are summarized in Algorithm 2.Algorithm 2:SASA algorithm.
(1)
Input:
(2)
Measurement signaly;
(3)
Sensing matrixΦ;
(4)
Initializeoutmax, innermaxsq, β.
(5)
Iteration:
(6)
Estimated sparsity levelK.
(7)
Initialize theK support set I0.
(8)
Itemp=I0;
(9)
InitializeT0=10∗fI0;
(10)
While i<outmax do
(11)
forj<innermax do
(12)
Iij=Itemp;
(13)
Chooseq elements from C−Iij and combined with Iij as E;
(14)
Calculatexij=ΦE†y;
(15)
Choose theK largest elements from xij as Itemp;
(16)
Calculateδ=fItemp−fIij;
(17)
if δ<0 then
(18)
Itemp=Itemp
(19)
else
(20)
Calculateβ=min1,exp−absδ/T0
(21)
Take a randomra∈0,1
(22)
if β>ra then Itemp=Itemp;
(22)
elseItemp=Iij
(23)
Calculateerror=fItemp
(24)
if error<β then break
(25)
end for
(26)
UpdateT0=T0/logi+1
(27)
End while
(28)
Output:
x
˜
=
Φ
I
temp
†
y;


## 4. Simulation Results and Analysis

In this section, a series of simulations are explained to evaluate the superior performance of the proposed SASA algorithm. To demonstrate the efficiency of the proposed algorithm, the SASA is compared with those existing state-of-the-art reconstruction methods. Those methods include SP, StOMP, CoSaMP, SWOMP [21], and ROMP. We performed simulation experiments on one-dimensional data and two-dimensional data, respectively. We adopted the Gaussian random matrix as the sampling operator for all methods.All the simulations were performed on the MATLAB (R2017b) on Windows 10 system with a 3.20 GHz Intel Core i7 processor and 8 GB memory.


### 4.1. One-Dimensional Signal

In the experiment, the Gaussian random sparse signal length is set toN = 256. The sensing matrix uses the random matrix to get the observation y. To evaluate the reconstruction performance of all algorithms, we conducted a variety of simulation experiments. Simultaneously, we calculated the average of 500 independent experiments as the final experimental result to avoid randomness of the experiments.In Figure2, we give the curve of the exact reconstruction rate under different measurement number M. The exact reconstruction rate can be interpreted as the number of correctly reconstructed signal in 500 independent experiments. If a signal is correctly reconstructed, we can think x−x^<10−6. x^ is the reconstruction signal. The sparsity K is set to 20. We can observe that the proposed algorithm SASA has a higher exact reconstruction rate for the different M. Especially, even when the measurement number M is less than 55, the proposed SASA is significantly superior to other algorithms.Figure 2
Curve of the exact reconstruction rate under different measurement numbersM.In Figure3, we provide the curve of the exact reconstruction rate under different sparsity level K. The measurement number is set to 60. K varies between 5 and 35. From Figure 3, we find that SASA is far superior to other algorithms under different sparsity levels. Even if K > 25, the proposed SASA is better than other classical algorithms.Figure 3
Curve of the exact reconstruction rate under different sparsity levelK. K=5,10,…,25,30.In Figure4, we provide the curve of the exact reconstruction rate under different signal length N. K = 10 and M = 25 are fixed, and signal length varies between 30 and 120. We observe that the proposed SASA performs better than the other five methods in the exact recovery rate. In general, the proposed SASA achieves higher exact reconstruction rate under different signal lengths.Figure 4
Curve of the exact reconstruction rate under different signal lengthK=10,M=25.In Figure5, we tested the exact reconstruction performance of the SASA algorithm under different sparsity levels. We give the exact reconstruction rate as a function of measurement number M. K varies between 5 and 25, and M varies between 10 and 200. We find that the SASA can exactly reconstruct the signal of different sparsity levels.Figure 5
Curve of the exact reconstruction rate under different measurement numbers.Through the above various simulation experiments, we can demonstrate that the proposed SASA algorithm has superior reconstruction performance under different reconstruction conditions.


### 4.2. Two-Dimensional Image

To validate the performance of the proposed SASA algorithm for real data, we use standard test images with a size of256×256. In addition, we need to note that the input image will be normalized, which facilitates sparse processing by using the discrete wavelet transform (DWT). We adopted a Gaussian random matrix as the compression sampling operator Φ. M = 190 and N = 256 are fixed. Figure 6 shows the image reconstructed by CoSaMP, SP, SWOMP, OMP, and our SASA algorithms. To assess the superior performances of all algorithms, the peak signal-to-noise ratio (PSNR) is employed as the reconstruction image quality indices. The PSNR is defined as follows:(11)PSNR=10log10M×N×2552x−x^2.Figure 6
Images reconstructed by different algorithms.All methods are evaluated on three different images, and those images include the Lena, Barche, and Camer. As shown in Figure6, we can find that the proposed SASA method gets better image reconstruction quality than other algorithms. The proposed SASA achieves the best PSNR. The PSNR also demonstrates that the SASA algorithm can make the reconstructed image smoother and achieve more details.Conclusively, various experiments show that the proposed SASA algorithm achieves superior performance for different types of data. Compared with other algorithms, the proposed SASA got higher exact reconstruction rate and PSNR.


## 4.1. One-Dimensional Signal

In the experiment, the Gaussian random sparse signal length is set toN = 256. The sensing matrix uses the random matrix to get the observation y. To evaluate the reconstruction performance of all algorithms, we conducted a variety of simulation experiments. Simultaneously, we calculated the average of 500 independent experiments as the final experimental result to avoid randomness of the experiments.In Figure2, we give the curve of the exact reconstruction rate under different measurement number M. The exact reconstruction rate can be interpreted as the number of correctly reconstructed signal in 500 independent experiments. If a signal is correctly reconstructed, we can think x−x^<10−6. x^ is the reconstruction signal. The sparsity K is set to 20. We can observe that the proposed algorithm SASA has a higher exact reconstruction rate for the different M. Especially, even when the measurement number M is less than 55, the proposed SASA is significantly superior to other algorithms.Figure 2
Curve of the exact reconstruction rate under different measurement numbersM.In Figure3, we provide the curve of the exact reconstruction rate under different sparsity level K. The measurement number is set to 60. K varies between 5 and 35. From Figure 3, we find that SASA is far superior to other algorithms under different sparsity levels. Even if K > 25, the proposed SASA is better than other classical algorithms.Figure 3
Curve of the exact reconstruction rate under different sparsity levelK. K=5,10,…,25,30.In Figure4, we provide the curve of the exact reconstruction rate under different signal length N. K = 10 and M = 25 are fixed, and signal length varies between 30 and 120. We observe that the proposed SASA performs better than the other five methods in the exact recovery rate. In general, the proposed SASA achieves higher exact reconstruction rate under different signal lengths.Figure 4
Curve of the exact reconstruction rate under different signal lengthK=10,M=25.In Figure5, we tested the exact reconstruction performance of the SASA algorithm under different sparsity levels. We give the exact reconstruction rate as a function of measurement number M. K varies between 5 and 25, and M varies between 10 and 200. We find that the SASA can exactly reconstruct the signal of different sparsity levels.Figure 5
Curve of the exact reconstruction rate under different measurement numbers.Through the above various simulation experiments, we can demonstrate that the proposed SASA algorithm has superior reconstruction performance under different reconstruction conditions.


## 4.2. Two-Dimensional Image

To validate the performance of the proposed SASA algorithm for real data, we use standard test images with a size of256×256. In addition, we need to note that the input image will be normalized, which facilitates sparse processing by using the discrete wavelet transform (DWT). We adopted a Gaussian random matrix as the compression sampling operator Φ. M = 190 and N = 256 are fixed. Figure 6 shows the image reconstructed by CoSaMP, SP, SWOMP, OMP, and our SASA algorithms. To assess the superior performances of all algorithms, the peak signal-to-noise ratio (PSNR) is employed as the reconstruction image quality indices. The PSNR is defined as follows:(11)PSNR=10log10M×N×2552x−x^2.Figure 6
Images reconstructed by different algorithms.All methods are evaluated on three different images, and those images include the Lena, Barche, and Camer. As shown in Figure6, we can find that the proposed SASA method gets better image reconstruction quality than other algorithms. The proposed SASA achieves the best PSNR. The PSNR also demonstrates that the SASA algorithm can make the reconstructed image smoother and achieve more details.Conclusively, various experiments show that the proposed SASA algorithm achieves superior performance for different types of data. Compared with other algorithms, the proposed SASA got higher exact reconstruction rate and PSNR.


## 5. Conclusions

In this paper, we propose a new sparsity adaptive simulated annealing algorithm. It can solve the sparse reconstruction issue efficiently. The proposed algorithm considers both adaptive sparsity and global optimization for sparse signal reconstruction. A series of experimental results show that the proposed SASA algorithm achieves the best reconstruction performance and is superior to the existing state-of-the-art methods in terms of reconstruction quality. For those advantages, the proposed SASA algorithm has broad application prospects and a higher guiding significance for the sparse signal reconstruction.



---

*Source: 6950819-2019-07-14.xml*