ruby - 这是 Matrix 类中的 ruby​​ 错误吗？

Question

我使用 ruby​​ 中的 Matrix 类来计算特征值和特征向量。使用此代码:

m = Matrix[[0.6868,0.6067],[0.6067,0.5978]]
v, d, v_inv = m.eigensystem

正确的结果应该是:

[ 1.25057433  0.03398123]
[[ 0.73251454 -0.68075138]
 [ 0.68075138  0.73251454]]

我使用 Python 通过 numpy 确认了这一点。

但是，我得到的结果如下:

d=[[0.033970204576497576,   0],
 [0,    1.2506297954235022]]

v=[[0.6807528514962294, 0.7325131774785713],
 [-0.7325131774785713,  0.6807528514962294]]

这是 ruby​​ 的错误吗？我的 ruby 版本如下:

ruby 2.1.2p95 (2014-05-08 revision 45877) [x86_64-linux]

Answer 1

这绝不是一个答案，但为了格式化，我会把它放在这里:

说到数学，我会用 r检查结果的正确性:

> entries <- c(0.6868, 0.6067, 0.6067, 0.5978)
> input <- matrix(entries, nrow=2, byrow=TRUE)
> input

#        [,1]   [,2]
# [1,] 0.6868 0.6067
# [2,] 0.6067 0.5978

> input_eigen <- eigen(input)
> input_eigen 

$values
# [1] 1.2506298 0.0339702

$vectors
#            [,1]       [,2]
# [1,] -0.7325132  0.6807529
# [2,] -0.6807529 -0.7325132

我相信上面的结果，这意味着 ruby​​ 可能比 python/numpy 做得更好。