Abstract:Finding sparse solutions to a system of equations and/or inequalities is an important topic in many application areas such as signal processing, statistical regression and nonparametric modeling. Various continuous relaxation models have been proposed and widely studied to deal with the discrete nature of the underlying problem. In this paper, we propose a quadratically constrained [Formula: see text] (0 < q < 1) minimization model for finding sparse solutions to a quadratic system. We prove that solving the proposed model is strongly NP-hard. To tackle the computation difficulty, a first order necessary condition for local minimizers is derived. Various properties of the proposed model are studied for designing an active-set-based descent algorithm to find candidate solutions satisfying the proposed condition. In addition to providing a theoretical convergence proof, we conduct extensive computational experiments using synthetic and real-life data to validate the effectiveness of the proposed algorithm and to show the superior capability in finding sparse solutions of the proposed model compared with other known models in the literature. We also extend our results to a quadratically constrained [Formula: see text] (0 < q < 1) minimization model with multiple convex quadratic constraints for further potential applications. Summary of Contribution: In this paper, we propose and study a quadratically constrained [Formula: see text] minimization (0 < q < 1) model for finding sparse solutions to a quadratic system which has wide applications in sparse signal recovery, image processing and machine learning. The proposed quadratically constrained [Formula: see text] minimization model extends the linearly constrained [Formula: see text] and unconstrained [Formula: see text]-[Formula: see text] models. We study various properties of the proposed model in aim of designing an efficient algorithm. Especially, we propose an unrelaxed KKT condition for local/global minimizers. Followed by the properties studied, an active-set based descent algorithm is then proposed with its convergence proof being given. Extensive numerical experiments with synthetic and real-life Sparco datasets are conducted to show that the proposed algorithm works very effectively and efficiently. Its sparse recovery capability is superior to that of other known models in the literature.

Algorithms on the sparse solution of under-determined linear systems

CONVERGENCE AND STABILITY OF ITERATIVELY REWEIGHTED LEAST SQUARES FOR LOW-RANK MATRIXRECOVERY

Stability and robustness of the l2/lq-minimization for block sparse recovery.

On the Schatten P-Quasi-norm Minimization for Low-Rank Matrix Recovery

Smoothed Low Rank and Sparse Matrix Recovery by Iteratively Reweighted Least Squares Minimization

Sparse Solutions by a Quadratically Constrained ℓq (0 <<i>q</i>< 1) Minimization Model

Sparse Estimation Via ℓ_q Optimization Method in High-Dimensional Linear Regression

Sparse Solutions by a Quadratically Constrained ℓ q (0 <

Generalized Nonconvex Nonsmooth Low-Rank Minimization

Sparse Estimation Via Lower-Order Penalty Optimization Methods in High-Dimensional Linear Regression.

Subspace Weighted L(2,1) Minimization For Sparse Signal Recovery

Sparse Solutions by a Quadratically Constrained <i>l<sub>q</sub></i> (0 < <i>q</i> < 1) Minimization Model

Nonnegative Iterative Reweighted Method for Sparse Linear Complementarity Problem

Lp quasi-norm minimization: algorithm and applications

Minimization Over the Nonconvex Sparsity Constraint Using A Hybrid First-order method

Reduced-Space Iteratively Reweighted Second-Order Methods for Nonconvex Sparse Regularization

Efficient Low-rank Identification via Accelerated Iteratively Reweighted Nuclear Norm Minimization

l0 Norm Constraint LMS Algorithm for Sparse System Identification.

Subspace weighted ℓ 2,1 minimization for sparse signal recovery

A Fast Algorithm for Sparse Reconstruction Based on Shrinkage, Subspace Optimization, and Continuation

Recovery Conditions of Sparse Signals Using Orthogonal Least Squares-Type Algorithms