关于科学计算中小数点位数处理的问题

lightninng

暂时找不到原因为何，但是感觉是numpy在处理float变量时精度方面的问题,遇到下面的种奇怪的情况

1. import pickle
   
2. with open("num.txt", "rb") as fp:
   
3.     temp = pickle.load(fp)
   
4. print(temp[0], temp[1])
   
5. print(temp[0] == temp[1])

复制代码

输出结果为

1. >>>
   
2. 2.87734318586 2.87734318586
   
3. False
   
4. >>>

复制代码

肉眼观察两个数字完全相等，但是判断两数是否相等时显示False
另外做了一个实验如下

1. >>> a = 2.87734318586
   
2. >>> b = 2.87734318586
   
3. >>> a == b
   
4. True
   
5. >>>

复制代码

PS：说一下出现这个问题的背景，在我的程序中得到了一个实对称的矩阵K(n*n)，自己然后另外一个值全为1/n的矩阵one(n*n)，然后得到新的矩阵K_new = K - K*one - one*K + one*K*one，这里的*表示矩阵乘法，按道理说这里的K_new矩阵应该是实对称矩阵，特征值应该为实数，但是我分析得到的特征值为复数，然后我将K_new中对应位置不相等的两个数组成元组用pickle.dump存放在test.txt中，然后作了上面的实验
附件中是我的num.txt文件 num.rar (18.81 KB, 下载次数: 2)

2015-4-15 13:53 上传

点击文件名下载附件
num.txt

wei_Y

你那个txt里的文件数据本身就是不一样的啦。
(2.877343185855147, 2.8773431858551475)
这个是我下载的数据。
print是四舍五入了。
不过若是这样一个元组打印的话不会四舍五入。但是以文件形式打开就会四舍五入，不知道为啥。。

lightninng

求解答啊~~

lightninng

wei_Y 发表于 2015-4-15 18:18
你那个txt里的文件数据本身就是不一样的啦。
(2.877343185855147, 2.8773431858551475)
这个是我下载的数 ...

谢啦，原来print会自动四舍五入啊，长见识了~~我下午倒腾了一下午这个问题，三点多的时候才想起来，把所有的矩阵元素都转化成str类型发现确实不一样，然后用round函数把所有的矩阵元素四舍五入到小数点后8位，终于判断实对称矩阵成功了~~
费了我一下午的时间

xubo

wei_Y 发表于 2015-4-15 13:53
你那个txt里的文件数据本身就是不一样的啦。
(2.877343185855147, 2.8773431858551475)
这个是我下载的数 ...

奇犽

print四舍五入引起的