The Noise-Sensitivity Phase Transition in Compressed Sensing

David L. Donoho, Arian Maleki, Andrea Montanari
2011 IEEE Transactions on Information Theory  
Consider the noisy underdetermined system of linear equations: y = Ax 0 + z, with A an n × N measurement matrix, n < N , and z ∼ N(0, σ 2 I) a Gaussian white noise. Both y and A are known, both x 0 and z are unknown, and we seek an approximation to x 0 . When x 0 has few nonzeros, useful approximations are often obtained by 1 -penalized 2 minimization, in which the reconstructionx 1,λ solves min{ y − Ax 2 2 /2 + λ x 1 }. Consider the reconstruction mean-squared error MSE = E x 1,λ − x 0 2 2 /N
more » ... and define the ratio MSE/σ 2 as the noise sensitivity. Consider matrices A with iid Gaussian entries and a largesystem limit in which n, N → ∞ with n/N → δ and k/n → ρ. We develop exact expressions for the asymptotic MSE ofx 1,λ , and evaluate its worst-case noise sensitivity over all types of k-sparse signals. The phase space 0 ≤ δ, ρ ≤ 1 is partitioned by the curve ρ = ρ MSE (δ) into two regions. Formal noise sensitivity is bounded throughout the region ρ < ρ MSE (δ) and is unbounded throughout the region ρ > ρ MSE (δ). The phase boundary ρ = ρ MSE (δ) is identical to the previously-known phase transition curve for equivalence of 1 − 0 minimization in the k-sparse noiseless case. Hence a single phase boundary describes the fundamental phase transitions both for the noiseless and noisy cases. Extensive computational experiments validate the predictions of these predictions, including the existence of game theoretical structures underlying it (saddlepoints in the payoff, least-favorable signals and maximin penalization). Underlying our formalism is an approximate message passing soft thresholding algorithm (AMP) introduced earlier by the authors. Other papers by the authors detail expressions for the formal MSE of AMP and its close connection to 1 -penalized reconstruction. The focus of the present paper is on computing the minimax formal MSE within the class of sparse signals x 0 . 1 1 Here and below we write a ∼ b if a/b → 1 as both quantities tend to infinity.
doi:10.1109/tit.2011.2165823 fatcat:cheheznymvd75memxzrlgodzpq