Partial gradient optimal thresholding algorithms for a class of sparse optimization problems

The optimization problems with a sparsity constraint is a class of important global optimization problems. A typical type of thresholding algorithms for solving such a problem adopts the traditional full steepest descent direction or Newton-like direction as a search direction to generate an iterate on which a certain thresholding is performed. Traditional hard thresholding discards a large part of a vector when the vector is dense. Thus a large part of important information contained in a dense vector has been lost in such a thresholding process. Recent study [Zhao, SIAM J Optim, 30(1), pp. 31-55, 2020] shows that the hard thresholding should be applied to a compressible vector instead of a dense vector to avoid a big loss of information. On the other hand, the optimal k-thresholding as a novel thresholding technique may overcome the intrinsic drawback of hard thresholding, and performs thresholding and objective function minimization simultaneously. This motivates us to propose the so-called partial gradient optimal thresholding method in this paper, which is an integration of the partial gradient and the optimal k-thresholding technique. The solution error bound and convergence for the proposed algorithms have been established in this paper under suitable conditions. Application of our results to the sparse optimization problems arising from signal recovery is also discussed. Experiment results from synthetic data indicate that the proposed algorithm called PGROTP is efficient and comparable to several existing algorithms. This work was founded by the National Natural Science Foundation of China (NSFC) under the grant 12071307. Nan Meng School of Mathematics, University of Birmingham, Edgbaston, Birmingham B15 2TT, United Kingdom E-mail: nxm563@bham.ac.uk Yun-Bin Zhao Shenzhen Research Institute of Big Data, Chinese University of Hong Kong, Shenzhen, China E-mail: yunbinzhao@cuhk.edu.cn Michal Kočvara School of Mathematics, University of Birmingham, Edgbaston, Birmingham B15 2TT, United Kingdom E-mail: m.kocvara@bham.ac.uk 2 Nan Meng et al.

[1]  Zhi-Quan Luo,et al.  Constructing New Weighted ℓ1-Algorithms for the Sparsest Points of Polyhedral Sets , 2017, Math. Oper. Res..

[2]  Yun-Bin Zhao,et al.  Analysis of Optimal Thresholding Algorithms for Compressed Sensing , 2019, ArXiv.

[3]  M. Fornasier,et al.  Iterative thresholding algorithms , 2008 .

[4]  Simon Foucart,et al.  Hard Thresholding Pursuit: An Algorithm for Compressive Sensing , 2011, SIAM J. Numer. Anal..

[5]  Olgica Milenkovic,et al.  Subspace Pursuit for Compressive Sensing: Closing the Gap Between Performance and Complexity , 2008, ArXiv.

[6]  YUN-BIN ZHAO,et al.  Optimal k-thresholding algorithms for sparse optimization problems , 2019, SIAM J. Optim..

[7]  Knud D. Andersen,et al.  The Mosek Interior Point Optimizer for Linear Programming: An Implementation of the Homogeneous Algorithm , 2000 .

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

[9]  S. Foucart,et al.  Hard thresholding pursuit algorithms: Number of iterations ☆ , 2016 .

[10]  Shenglong Zhou,et al.  Subspace Newton method for the $\ell_0$-regularized optimization , 2020 .

[11]  T. Blumensath,et al.  Theory and Applications , 2011 .

[12]  Dong In Kim,et al.  Compressed Sensing for Wireless Communications: Useful Tips and Tricks , 2015, IEEE Communications Surveys & Tutorials.

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

[14]  Duan Li,et al.  Reweighted 1-Minimization for Sparse Solutions to Underdetermined Linear Systems , 2012, SIAM J. Optim..

[15]  Jan Vybíral,et al.  Compressed Sensing and its Applications , 2015 .

[16]  Naihua Xiu,et al.  Global and Quadratic Convergence of Newton Hard-Thresholding Pursuit , 2019, J. Mach. Learn. Res..

[17]  Guoxian Huang,et al.  Soft-thresholding Orthogonal Matching Pursuit for efficient signal reconstruction , 2013, 2013 IEEE International Conference on Acoustics, Speech and Signal Processing.

[18]  Michael A. Saunders,et al.  Atomic Decomposition by Basis Pursuit , 1998, SIAM J. Sci. Comput..

[19]  Mike E. Davies,et al.  Iterative Hard Thresholding for Compressed Sensing , 2008, ArXiv.

[20]  Qingshan Liu,et al.  Newton Greedy Pursuit: A Quadratic Approximation Method for Sparsity-Constrained Optimization , 2014, 2014 IEEE Conference on Computer Vision and Pattern Recognition.

[21]  Rama Chellappa,et al.  Sparse Representations, Compressive Sensing and dictionaries for pattern recognition , 2011, The First Asian Conference on Pattern Recognition.

[22]  Yun-Bin Zhao,et al.  Newton-Step-Based Hard Thresholding Algorithms for Sparse Signal Recovery , 2020, IEEE Transactions on Signal Processing.

[23]  Holger Rauhut,et al.  A Mathematical Introduction to Compressive Sensing , 2013, Applied and Numerical Harmonic Analysis.

[24]  Z. Luo,et al.  Improved RIP-Based Bounds for Guaranteed Performance of Several Compressed Sensing Algorithms , 2020, 2007.01451.

[25]  Joel A. Tropp,et al.  Signal Recovery From Random Measurements Via Orthogonal Matching Pursuit , 2007, IEEE Transactions on Information Theory.

[26]  Di Guo,et al.  Projected Iterative Soft-Thresholding Algorithm for Tight Frames in Compressed Sensing Magnetic Resonance Imaging , 2015, IEEE Transactions on Medical Imaging.

[27]  Olgica Milenkovic,et al.  Subspace Pursuit for Compressive Sensing Signal Reconstruction , 2008, IEEE Transactions on Information Theory.

[28]  Lie Wang,et al.  Orthogonal Matching Pursuit for Sparse Signal Recovery With Noise , 2011, IEEE Transactions on Information Theory.

[29]  Emmanuel J. Candès,et al.  Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information , 2004, IEEE Transactions on Information Theory.

[30]  Thomas Blumensath,et al.  Accelerated iterative hard thresholding , 2012, Signal Process..

[31]  Yun-Bin Zhao,et al.  A New Computational Method for the Sparsest Solutions to Systems of Linear Equations , 2015, SIAM J. Optim..

[32]  Nan Meng,et al.  Newton-Type Optimal Thresholding Algorithms for Sparse Optimization Problems , 2021, Journal of the Operations Research Society of China.

[33]  Yun-Bin Zhao,et al.  Sparse Optimization Theory and Methods , 2018 .

[34]  Rahul Garg,et al.  Gradient descent with sparsification: an iterative algorithm for sparse recovery with restricted isometry property , 2009, ICML '09.

[35]  Deanna Needell,et al.  CoSaMP: Iterative signal recovery from incomplete and inaccurate samples , 2008, ArXiv.

[36]  David L. Donoho,et al.  De-noising by soft-thresholding , 1995, IEEE Trans. Inf. Theory.

[37]  Mike E. Davies,et al.  Normalized Iterative Hard Thresholding: Guaranteed Stability and Performance , 2010, IEEE Journal of Selected Topics in Signal Processing.

[38]  Deanna Needell,et al.  Signal Recovery From Incomplete and Inaccurate Measurements Via Regularized Orthogonal Matching Pursuit , 2007, IEEE Journal of Selected Topics in Signal Processing.