Augmented _1 and nuclear-norm models with a globally linearly convergent algorithm

From MaRDI portal
Publication:2873229

DOI10.1137/120863290zbMATH Open1279.68329arXiv1201.4615OpenAlexW2075826622MaRDI QIDQ2873229FDOQ2873229


Authors: Ming-Jun Lai, Wotao Yin Edit this on Wikidata


Publication date: 23 January 2014

Published in: SIAM Journal on Imaging Sciences (Search for Journal in Brave)

Abstract: This paper studies the long-existing idea of adding a nice smooth function to "smooth" a non-differentiable objective function in the context of sparse optimization, in particular, the minimization of ||x||1+1/(2alpha)||x||22, where x is a vector, as well as the minimization of ||X||+1/(2alpha)||X||F2, where X is a matrix and ||X||* and ||X||F are the nuclear and Frobenius norms of X, respectively. We show that they can efficiently recover sparse vectors and low-rank matrices. In particular, they enjoy exact and stable recovery guarantees similar to those known for minimizing ||x||1 and ||X||* under the conditions on the sensing operator such as its null-space property, restricted isometry property, spherical section property, or RIPless property. To recover a (nearly) sparse vector x0, minimizing ||x||1+1/(2alpha)||x||22 returns (nearly) the same solution as minimizing ||x||1 almost whenever alphage10||x0||infty. The same relation also holds between minimizing ||X||+1/(2alpha)||X||F2 and minimizing ||X||* for recovering a (nearly) low-rank matrix X0, if alphage10||X0||2. Furthermore, we show that the linearized Bregman algorithm for minimizing ||x||1+1/(2alpha)||x||22 subject to Ax=b enjoys global linear convergence as long as a nonzero solution exists, and we give an explicit rate of convergence. The convergence property does not require a solution solution or any properties on A. To our knowledge, this is the best known global convergence result for first-order sparse optimization algorithms.


Full work available at URL: https://arxiv.org/abs/1201.4615




Recommendations





Cited In (36)





This page was built for publication: Augmented \(\ell_1\) and nuclear-norm models with a globally linearly convergent algorithm

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2873229)