<p>We propose a way to unify the construction of Newton’s methods for matrix decompositions. This construction is based on a perturbation analysis of a suitable system associated with matrix decompositions. Then it appears that the resolution of a linear Sylvester equation permits defining Newton’s method. We give a general result to analyze the quadratic convergence of this Newton’s method. We apply it to classical matrix decompositions: <i>LU</i> decomposition, <i>QR</i> decomposition, eigenproblem in the diagonalizable case, singular value decomposition, and Schur decomposition. Finally, we propose for future work a generalization of this construction to approximate matrix decomposition by high-order methods.</p>

错误:搜索内容不能为空,请输入英文关键词
错误:关键词超出字数限制,请精简
高级检索

Unified Newton’s method for matrix decompositions

  • Jean-Claude Yakoubsohn

摘要

We propose a way to unify the construction of Newton’s methods for matrix decompositions. This construction is based on a perturbation analysis of a suitable system associated with matrix decompositions. Then it appears that the resolution of a linear Sylvester equation permits defining Newton’s method. We give a general result to analyze the quadratic convergence of this Newton’s method. We apply it to classical matrix decompositions: LU decomposition, QR decomposition, eigenproblem in the diagonalizable case, singular value decomposition, and Schur decomposition. Finally, we propose for future work a generalization of this construction to approximate matrix decomposition by high-order methods.