Robust principal component analysis using facial reduction

被引:3
作者
Ma, Shiqian [1 ]
Wang, Fei [2 ]
Wei, Linchuan [3 ]
Wolkowicz, Henry [3 ]
机构
[1] Univ Calif Davis, Dept Math, Davis, CA 95616 USA
[2] Royal Inst Technol, Dept Math, Div Optimizat & Syst Theory, Stockholm, Sweden
[3] Univ Waterloo, Fac Math, Dept Combinator & Optimizat, Waterloo, ON N2L 3G1, Canada
关键词
Robust principal component analysis; Semidefinite cone; Facial reduction; SEMIDEFINITE; ALGORITHMS;
D O I
10.1007/s11081-019-09476-9
中图分类号
T [工业技术];
学科分类号
08 ;
摘要
We introduce a novel approach for robust principal component analysis (RPCA) for a partially observed data matrix. The aim is to recover the data matrix as a sum of a low-rank matrix and a sparse matrix so as to eliminate erratic noise (outliers). This problem is known to be NP-hard in general. A classical approach to solving RPCA is to consider convex relaxations. One such heuristic involves the minimization of the (weighted) sum of a nuclear norm part, that promotes a low-rank component, with an l(1) norm part, to promote a sparse component. This results in a well-structured convex problem that can be efficiently solved by modern first-order methods. However, first-order methods often yield low accuracy solutions. Moreover, the heuristic of using a norm consisting of a weighted sum of norms may lose some of the advantages that each norm had when used separately. In this paper, we propose a novel nonconvex and nonsmooth reformulation of the original NP-hard RPCA model. The new model adds a redundant semidefinite cone constraint and solves small subproblems using a PALM algorithm. Each subproblem results in an exposing vector for a facial reduction technique that is able to reduce the size significantly. This makes the problem amenable to efficient algorithms in order to obtain high-level accuracy. We include numerical results that confirm the efficacy of our approach.
引用
收藏
页码:1195 / 1219
页数:25
相关论文
共 26 条
[1]  
Aybat NS, 2016, ALGORITHMS STABLE PC
[2]  
Blair J.R., 1993, Graph theory and sparse matrix computation, P1
[3]   Proximal alternating linearized minimization for nonconvex and nonsmooth problems [J].
Bolte, Jerome ;
Sabach, Shoham ;
Teboulle, Marc .
MATHEMATICAL PROGRAMMING, 2014, 146 (1-2) :459-494
[4]   Robust PCA via Principal Component Pursuit: A review for a comparative evaluation in video surveillance [J].
Bouwmans, Thierry ;
Zahzah, El Hadi .
COMPUTER VISION AND IMAGE UNDERSTANDING, 2014, 122 :22-34
[5]   Robust Principal Component Analysis? [J].
Candes, Emmanuel J. ;
Li, Xiaodong ;
Ma, Yi ;
Wright, John .
JOURNAL OF THE ACM, 2011, 58 (03)
[6]  
CHANDRASEKARAN V, 2009, AARXIV09062220 MIT
[7]  
Cheng J., 2015, P SEG TECH PROGR AUG, P4646
[8]  
Drusvyatskiy Dmitriy, 2017, Foundations and Trends in Optimization, V3, P77, DOI 10.1561/2400000011
[9]   NOISY EUCLIDEAN DISTANCE REALIZATION: ROBUST FACIAL REDUCTION AND THE PARETO FRONTIER [J].
Drusvyatskiy, D. ;
Krislock, N. ;
Voronin, Y. -L. ;
Wolkowicz, H. .
SIAM JOURNAL ON OPTIMIZATION, 2017, 27 (04) :2301-2331
[10]   COORDINATE SHADOWS OF SEMIDEFINITE AND EUCLIDEAN DISTANCE MATRICES [J].
Drusvyatskiy, Dmitriy ;
Pataki, Gabor ;
Wolkowicz, Henry .
SIAM JOURNAL ON OPTIMIZATION, 2015, 25 (02) :1160-1178