$\ell _{0}$ Factor Analysis: A P-Stationary Point Theory

Factor analysis is a widely used modeling technique for stationary time series, which achieves dimensionality reduction by revealing a hidden low-rank plus sparse structure of the covariance matrix. Such an idea of parsimonious modeling has also been important in the field of systems and control. In this article, a nonconvex nonsmooth optimization problem involving the <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"><tex-math notation="LaTeX">$\ell _{0}$</tex-math></inline-formula> norm is constructed in order to achieve the low-rank and sparse additive decomposition of the sample covariance matrix. We establish the existence of an optimal solution and characterize these solutions via the concept of proximal stationary points. Furthermore, an ADMM algorithm is designed to solve the <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"><tex-math notation="LaTeX">$\ell _{0}$</tex-math></inline-formula> optimization problem, and a subsequence convergence result is proved under reasonable assumptions. Finally, numerical experiments demonstrate the effectiveness of our method in comparison with some alternatives in the literature.