我有n个矩阵,我想找到每个平方根。但是我的算法需要cholesky。我得到的错误矩阵不是正定的。我将对角线元素转换成真正的一个。但我得到了同样的错误。有没有其他方法可以找到矩阵的cholesky?
答案 0 :(得分:0)
如果你的矩阵距离正定很远,那么你无能为力 - Cholesky因子分解是基于它是正定的假设。
通常情况下,矩阵基本上是正定的,但由于一些小的数值问题非常略微不对称。如果这是您遇到的问题,您可以通过说x
强制它(让我们称之为x = (x+x')/2
)对称。
希望有所帮助!