r - 矩阵求逆

Question

我想对一个方形对称正定矩阵求逆。我知道R中有两个函数solve()和chol2inv()，但是它们的结果不同。我需要知道为什么会这样吗？

谢谢。

Answer 1

对于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