Cholesky 因子的计算

计算科学 线性代数 算法
2021-12-18 02:16:36

因此,Cholesky 分解定理表明,任何实对称正定矩阵M有 Cholesky 分解M=LL在哪里L是一个下三角矩阵。

给定M,我们已经知道有快速算法来计算其 Cholesky 因子L.

现在,假设给了我一个矩形m×n矩阵A,我知道AA是肯定的。有没有办法计算 Cholesky 因子LAA无需计算AA显式然后应用 Cholesky 分解算法?

如果A是一个非常大的矩形矩阵AA明确地似乎非常昂贵,因此是这个问题。

2个回答

是的,您可以使用 QR 分解获得因子(直到条目的符号);看到这个答案请注意,如果您感兴趣的只是解决导致涉及的正规方程的最小二乘问题ATA, 可以直接使用 QR 分解。

是的。计算QR因式分解和取L=RT; 重新调整行R如有必要(通过更改它们的某些符号)使对角线的符号为非负(因为 Cholesky 因子被定义为具有非负对角线)。

对于稀疏 QR 分解,请参见 http://dl.acm.org/citation.cfm?id=174408