相关疑难解决方法(0)

我想反转一个方形对称正定矩阵.我知道R中有两个功能solve(),chol2inv()但结果不同.我需要知道为什么会这样？

谢谢.

我比较了计算对称矩阵的逆的各种方法：

• 解决（来自 LAPCK 包）
• 解决（但使用更高的机器精度）
• qr.solve（据说更快）
• ginv（MASS 包，Moore-Penrose 算法的实现）
• chol2inv（使用 Cholesky 分解）

通过特征值比较逆矩阵：

R
library(MASS)

## Create the matrix
m = replicate(10, runif(n=10)) 
m[lower.tri(m)] = t(m)[lower.tri(m)]

## Inverse the matrix
inv1 = solve(m)
inv2 = solve(m, tol = .Machine$double.eps)
inv3 = qr.solve(m)
inv4 = ginv(m)
inv5 = chol2inv(m)

## Eigenvalues of the inverse
em1=eigen(inv1)
em2=eigen(inv2)
em3=eigen(inv3)
em4=eigen(inv4)
em5=eigen(inv5)

## Plot the abs of the eigenvalues (may be complex)
myPch=c( 20, 15, 17, 25, 3 )
plot(abs(em1$values),pch=myPch[1],cex=1.5)
points(abs(em2$values),pch=myPch[2], cex=1.5)
points(abs(em3$values),pch=myPch[3], cex=1.5) …