R中的solve(a)和ginv(a)有什么区别?

时间:2014-09-08 10:28:47

标签: r rstudio

我正在研究R中的矩阵运算。我遇到了两种不同的方法来找到矩阵的逆矩阵:

  • MASS::ginv(a)
  • solve(a)

它们返回相同的结果或看起来似乎但当我尝试使用ginv(a) == solve(a)等同它们时,只有前导对角线的元素似乎相等。它为这些元素返回True,否则返回零。

我对这种行为感到困惑。即使我可以看到它们具有完全相同的数值,也会发生这种情况。

我正在使用R Studio和Windows。

 library(MASS)
 a1 <- rbind(c(1,-0.25),c(-0.25,1))
 solve(a1)
 ginv(a1)
 solve(a1) == ginv(a1)

请有人向我解释发生了什么。是因为伪逆?

0 个答案:

没有答案