The Augmented Lagrange Multiplier Method for Exact Recovery of Corrupted Low-Rank Matrices
Zhouchen Lin, Minming Chen, Yi Ma
TL;DR
This work introduces scalable augmented Lagrange multiplier methods (Exact ALM and Inexact ALM) for robust principal component analysis and extends the approach to matrix completion. By formulating the problems with nuclear-norm and L1-regularization, the ALM framework achieves exact convergence with practical efficiency, often outperforming state-of-the-art APG and SVT-based methods in both speed and accuracy. The authors provide detailed implementation guidelines, convergence proofs, and public Matlab code, making the methods readily usable for large-scale low-rank recovery tasks. The results suggest significant practical impact for applications in image processing, web data analysis, and bioinformatics where large corrupted datasets must be recovered efficiently.
Abstract
This paper proposes scalable and fast algorithms for solving the Robust PCA problem, namely recovering a low-rank matrix with an unknown fraction of its entries being arbitrarily corrupted. This problem arises in many applications, such as image processing, web data ranking, and bioinformatic data analysis. It was recently shown that under surprisingly broad conditions, the Robust PCA problem can be exactly solved via convex optimization that minimizes a combination of the nuclear norm and the $\ell^1$-norm . In this paper, we apply the method of augmented Lagrange multipliers (ALM) to solve this convex program. As the objective function is non-smooth, we show how to extend the classical analysis of ALM to such new objective functions and prove the optimality of the proposed algorithms and characterize their convergence rate. Empirically, the proposed new algorithms can be more than five times faster than the previous state-of-the-art algorithms for Robust PCA, such as the accelerated proximal gradient (APG) algorithm. Moreover, the new algorithms achieve higher precision, yet being less storage/memory demanding. We also show that the ALM technique can be used to solve the (related but somewhat simpler) matrix completion problem and obtain rather promising results too. We further prove the necessary and sufficient condition for the inexact ALM to converge globally. Matlab code of all algorithms discussed are available at http://perception.csl.illinois.edu/matrix-rank/home.html
