We apply Dykstra's alternating projection algorithm to the constrained least-squares matrix problem that arises naturally in statistics and mathematical economics. In particular, we are concerned with the problem of finding the closest symmetric positive definite bounded and patterned matrix, in the
Dykstra's algorithm for constrained least-squares rectangular matrix problems
β Scribed by R. Escalante; M. Raydan
- Publisher
- Elsevier Science
- Year
- 1998
- Tongue
- English
- Weight
- 404 KB
- Volume
- 35
- Category
- Article
- ISSN
- 0898-1221
No coin nor oath required. For personal study only.
β¦ Synopsis
ln a recent paper, the authors applied Dykstra's alternating projection algorithm to solve constrained least-squares n x n matrix problems. We extend these results in two different directions. First, we make use of the singular value decomposition to solve now constrained leastsquares rectangular m x n matrix problems that arise in several applications. Second, we propose a new and improved implementation of the projection algorithm onto the e-positive definite set of matrices. This implementation does not require the computation of all elgenvalues and eigenvectors of a matrix per iteration, and still guarantees convergence. Finally, encouraging preliminary numerical results are discussed.
π SIMILAR VOLUMES
For linear least squares problems min x Ax -b 2 , where A is sparse except for a few dense rows, a straightforward application of Cholesky or QR factorization will lead to catastrophic fill in the factor R. We consider handling such problems by a matrix stretching technique, where the dense rows ar
We present some perturbation results for least squares problems with equality constraints. Relative errors are obtained on perturbed solutions and Lagrange multipliers of the problem, based on the equivalence of the problem to a consistent system of linear equations.