我想对一个方形对称正定矩阵求逆。我知道R中有两个函数solve()
和chol2inv()
,但是它们的结果不同。我需要知道为什么会这样吗?
谢谢。
最佳答案
对于solve
,您需要给出原始矩阵,但是对于chol2inv
,则使用预先计算的cholesky分解:
set.seed(1)
a<-crossprod(matrix(rnorm(9),3,3))
a_chol<-chol(a)
solve(a)
[,1] [,2] [,3]
[1,] 1.34638151 -0.02957435 0.8010735
[2,] -0.02957435 0.32780020 -0.1786295
[3,] 0.80107345 -0.17862950 1.4533671
chol2inv(a_chol)
[,1] [,2] [,3]
[1,] 1.34638151 -0.02957435 0.8010735
[2,] -0.02957435 0.32780020 -0.1786295
[3,] 0.80107345 -0.17862950 1.4533671
关于r - 矩阵求逆,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/15336123/