numpy的问题点积不起作用

Question

我目前正在学习线性代数，我正在尝试尝试我在代码中学到的东西。我在 python 中使用 numpy。

问题在于，正如您所知，任何矩阵乘以它的逆矩阵都会得到单位矩阵。

但在这里它不会发生。

B= np.array([[3.0, 0.0, 2.0], [2.0, 0.0, -2.0], [0.0, 1.0, 1.0]])
Binv = np.linalg.inv(B)
Bbits= np.dot(B,Binv)

我已经手动尝试过了。它正确地得到逆，但点积不是。 查看 Bbits 的输出

[[ 1.00000000e+00  0.00000000e+00 -5.55111512e-17]
 [ 0.00000000e+00  1.00000000e+00  0.00000000e+00]
 [ 0.00000000e+00  0.00000000e+00  1.00000000e+00]]

Answer 1

你的结果实际上是正确的！

运行您提供的完全相同的代码，但添加了np.around()，

B= np.array([[3.0, 0.0, 2.0], [2.0, 0.0, -2.0], [0.0, 1.0, 1.0]])
Binv = np.linalg.inv(B)
Bbits= np.around(np.dot(B,Binv),2)

我们得到以下输出：

array([[ 1.,  0.,  0.],
       [-0.,  1.,  0.],
       [ 0.,  0.,  1.]])

在您的结果中，-5.55111512e-17 只是一个舍入错误！