Analyzing Weighted $\ell_1$ Minimization for Sparse Recovery With Nonuniform Sparse Models

In this paper, we introduce a nonuniform sparsity model and analyze the performance of an optimized weighted ℓ1 minimization over that sparsity model. In particular, we focus on a model where the entries of the unknown vector fall into two sets, with entries of each set having a specific probability of being nonzero. We propose a weighted ℓ1 minimization recovery algorithm and analyze its performance using a Grassmann angle approach. We compute explicitly the relationship between the system parameters-the weights, the number of measurements, the size of the two sets, the probabilities of being nonzero-so that when i.i.d. random Gaussian measurement matrices are used, the weighted ℓ1 minimization recovers a randomly selected signal drawn from the considered sparsity model with overwhelming probability as the problem dimension increases. This allows us to compute the optimal weights. We demonstrate through rigorous analysis and simulations that for the case when the support of the signal can be divided into two different subclasses with unequal sparsity fractions, the weighted ℓ1 minimization outperforms the regular ℓ1 minimization substantially. We also generalize our results to signal vectors with an arbitrary number of subclasses for sparsity.

[1]  Wei Lu,et al.  Modified compressive sensing for real-time dynamic MR imaging , 2009, 2009 16th IEEE International Conference on Image Processing (ICIP).

[2]  Wei Lu,et al.  Modified-CS: Modifying compressive sensing for problems with partially known support , 2009, 2009 IEEE International Symposium on Information Theory.

[3]  B. Hassibi,et al.  Compressed sensing over the Grassmann manifold: A unified analytical framework , 2008, 2008 46th Annual Allerton Conference on Communication, Control, and Computing.

[4]  Rolf Schneider,et al.  Random projections of regular simplices , 1992, Discret. Comput. Geom..

[5]  A. Salman Avestimehr,et al.  Breaking the ℓ1 recovery thresholds with reweighted ℓ1 optimization , 2009, 2009 47th Annual Allerton Conference on Communication, Control, and Computing (Allerton).

[6]  Chiara Sabatti,et al.  Empirical Bayes Estimation of a Sparse Vector of Gene Expression Changes , 2005, Statistical applications in genetics and molecular biology.

[7]  Emmanuel J. Candès,et al.  Decoding by linear programming , 2005, IEEE Transactions on Information Theory.

[8]  B. Grünbaum Grassmann angles of convex polytopes , 1968 .

[9]  Babak Hassibi,et al.  On the Reconstruction of Block-Sparse Signals With an Optimal Number of Measurements , 2008, IEEE Transactions on Signal Processing.

[10]  Babak Hassibi,et al.  Recovering Sparse Signals Using Sparse Measurement Matrices in Compressed DNA Microarrays , 2008, IEEE Journal of Selected Topics in Signal Processing.

[11]  David L. Donoho,et al.  High-Dimensional Centrally Symmetric Polytopes with Neighborliness Proportional to Dimension , 2006, Discret. Comput. Geom..

[12]  Wei Lu,et al.  Regularized Modified-BPDN for compressive sensing with partially known support , 2010, arXiv.org.

[13]  H. Hadwiger Gitterpunktanzahl im Simplex und Wills'sche Vermutung , 1979 .

[14]  P Mcmullen,et al.  Non-linear angle-sum relations for polyhedral cones and polytopes , 1975, Mathematical Proceedings of the Cambridge Philosophical Society.

[15]  L. Santaló Geometría integral en espacios de curvatura constante , 1952 .

[16]  Volkan Cevher,et al.  Model-Based Compressive Sensing , 2008, IEEE Transactions on Information Theory.

[17]  Michael Rabbat,et al.  Compressed Network Monitoring , 2007, 2007 IEEE/SP 14th Workshop on Statistical Signal Processing.

[18]  Weiyu Xu,et al.  Improved sparse recovery thresholds with two-step reweighted ℓ1 minimization , 2010, 2010 IEEE International Symposium on Information Theory.

[19]  Lawrence Carin,et al.  Bayesian Compressive Sensing , 2008, IEEE Transactions on Signal Processing.

[20]  Stephen P. Boyd,et al.  Enhancing Sparsity by Reweighted ℓ1 Minimization , 2007, 0711.1612.

[21]  D. Donoho,et al.  Neighborliness of randomly projected simplices in high dimensions. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[22]  Richard G. Baraniuk,et al.  Detection and estimation with compressive measurements , 2006 .

[23]  D. Donoho For most large underdetermined systems of linear equations the minimal 𝓁1‐norm solution is also the sparsest solution , 2006 .

[24]  Weiyu Xu,et al.  Compressed sensing - probabilistic analysis of a null-space characterization , 2008, 2008 IEEE International Conference on Acoustics, Speech and Signal Processing.

[25]  G. C. Shephard,et al.  Convex Polytopes , 1967 .