Article ID | Journal | Published Year | Pages | File Type |
---|---|---|---|---|
6420606 | Applied Mathematics and Computation | 2015 | 11 Pages |
In this paper, we propose a numerical method for computing the nearest low-rank correlation matrix (LRCM). Motivated by the fact that the nearest LRCM problem can be reformulated as a standard nonlinear equality constrained optimization problem with matrix variables via the Gramian representation, we propose a new algorithm based on the sequential quadratic programming (SQP) method. On each iteration, we do not solve the quadratic program (QP) corresponding to the exact Hessian, but a modified QP with a simpler Hessian. This QP subproblem can be solved efficiently by equivalently transforming it to a sparse linear system. Global convergence is established and preliminary numerical results are presented to demonstrate the proposed method is potentially useful.