[已解决]为啥我使用numpy计算行列式时，为0的行列式它不等于0而是等于一个很小的数啊

ALNY

import numpy as np
import numpy.linalg as nl
x1=np.array([[1.,2.,3],[4.,5.,6.],[7.,8.,9,]])
print(x1)
print(nl.det(x1))

[[1. 2. 3.]
[4. 5. 6.]
[7. 8. 9.]]
-9.51619735392994e-16

最佳答案

月排行榜 / 总排行榜

Twilight6

2021-9-16 12:31:22

ALNY 发表于 2021-9-16 12:20
那请问有什么处理方法吗

呃，这误差没办法处理了，这是计算机存储数据时候产生的误差，但是对结果进行处理还是可以的：

1. import numpy as np
   
2. import numpy.linalg as nl
   
3. x1=np.array([[1.,2.,3],[4.,5.,6.],[7.,8.,9,]])
   
4. print(x1)
   
5. print(round(nl.det(x1),10))

复制代码

楼上的结果为 0.0 ，可能是因为计算机存储方式不同的原因吧，就比如：

1. >>> 0.1+0.1-0.2
   
2. 0.0
   
3. 
   
4. >>> 0.1+0.7-0.8
   
5. -1.1102230246251565e-16

复制代码

跳转到最佳答案楼层

Twilight6

因为 Python 浮点数是有误差的

导致最终结果为一个非常小的数，这么小的数视为零即可

ALNY

Twilight6 发表于 2021-9-16 11:59
因为 Python 浮点数是有误差的

导致最终结果为一个非常小的数，这么小的数视为零即可

那请问有什么处理方法吗

傻眼貓咪

為什麼我用你的代碼執行結果為0.0？

1. import numpy as np
   
2. import numpy.linalg as nl
   
3. x1=np.array([[1.,2.,3],[4.,5.,6.],[7.,8.,9,]])
   
4. print(x1)
   
5. print(nl.det(x1))

复制代码

1. [[1. 2. 3.]
   
2. [4. 5. 6.]
   
3. [7. 8. 9.]]
   
4. 0.0

复制代码

ALNY

傻眼貓咪 发表于 2021-9-16 12:21
為什麼我用你的代碼執行結果為0.0？

我也很奇怪 是numpy版本问题么

Twilight6

ALNY 发表于 2021-9-16 12:20
那请问有什么处理方法吗

呃，这误差没办法处理了，这是计算机存储数据时候产生的误差，但是对结果进行处理还是可以的：

1. import numpy as np
   
2. import numpy.linalg as nl
   
3. x1=np.array([[1.,2.,3],[4.,5.,6.],[7.,8.,9,]])
   
4. print(x1)
   
5. print(round(nl.det(x1),10))

复制代码

楼上的结果为 0.0 ，可能是因为计算机存储方式不同的原因吧，就比如：

1. >>> 0.1+0.1-0.2
   
2. 0.0
   
3. 
   
4. >>> 0.1+0.7-0.8
   
5. -1.1102230246251565e-16

复制代码

ALNY

Twilight6 发表于 2021-9-16 12:31
呃，这误差没办法处理了，这是计算机存储数据时候产生的误差，但是对结果进行处理还是可以的：

确实只能处理结果了....
刚刚试了下  这个好神奇啊

傻眼貓咪

学到了

1. print(0.1+0.7-0.8)

复制代码

1. -1.1102230246251565e-16

复制代码

好像只有十进制 0.1 才会这样