为什么这两个矩阵在 R 中不相同？

如何解决为什么这两个矩阵在 R 中不相同？

目前正在学习 R 中的一些基础知识。我遇到了一个问题，即 R 无法将两个（看似）相同的矩阵识别为相同的。

K <- rbind(c(63,32,32),c(69,20,31),c(58,41,30))
L <- rbind(c(36,35,65),c(61,39,36),c(43,67,25))

c3 <- solve(K%*%L)
c4 <- (solve(L)%*%solve(K))

这是c3的回报

> c3
             [,1]         [,2]         [,3]
[1,] -0.008210912  0.002098811  0.006187541
[2,]  0.019277529 -0.007768868 -0.011693360
[3,] -0.010079058  0.005248454  0.005097116

这是c4的回报

> c4
             [,] -0.010079058  0.005248454  0.005097116

c3和c4的typeof和class是一样的。然而这种情况发生了：

> c3 == c4
      [,1]  [,2]  [,] FALSE FALSE FALSE
[2,] FALSE FALSE FALSE
[3,] FALSE FALSE FALSE

解决方法

代码使用的是不准确的浮点运算。如果您查看足够多的小数，不同的方法可能会导致略有不同的值。

c3-c4
##               [,1]          [,2]          [,3]
## [1,] -1.214306e-17  1.734723e-18  1.040834e-17
## [2,]  6.938894e-18  5.204170e-18 -1.214306e-17
## [3,]  5.204170e-18 -7.806256e-18  8.673617e-19

忽略细微的差异，它们是相同的：

all.equal(c3,c4)
## [1] TRUE

见http://www.hep.by/gnu/r-patched/r-faq/R-FAQ_82.html