PSD矩阵+对角矩阵的LU分解

计算科学 线性代数 密集矩阵
2021-12-21 09:36:28

如果我有一个 psd,对称矩阵A我需要做 LU decomps onBi=A+Di(在哪里Di是一个对角 psd 矩阵,其中Di随机变化)。有没有一种简单的方法可以从 LU 出发A+ 的一些功能Di得到 LUBi?

文献中的任何参考资料将不胜感激。

1个回答

使用 Cholesky 分解或 LDL 分解代替 LU。

从以下答案来看:

没有好的方法可以更新任何这些分解。除非对角矩阵的秩非常低,否则每次都从头开始重新计算分解不会更糟,在这种情况下,也许您可​​以执行类似于 Sherman-Morrison(-Woodbury) 类型的更新。