Python如何执行精确的浮点数运算


Posted in Python onJuly 31, 2020

问题

你需要对浮点数执行精确的计算操作,并且不希望有任何小误差的出现。

解决方案

浮点数的一个普遍问题是它们并不能精确的表示十进制数。 并且,即使是最简单的数学运算也会产生小的误差,比如:

>>> a = 4.2
>>> b = 2.1
>>> a + b
6.300000000000001
>>> (a + b) == 6.3
False
>>>

这些错误是由底层CPU和IEEE 754标准通过自己的浮点单位去执行算术时的特征。 由于Python的浮点数据类型使用底层表示存储数据,因此你没办法去避免这样的误差。

如果你想更加精确(并能容忍一定的性能损耗),你可以使用 decimal 模块:

>>> from decimal import Decimal
>>> a = Decimal('4.2')
>>> b = Decimal('2.1')
>>> a + b
Decimal('6.3')
>>> print(a + b)
6.3
>>> (a + b) == Decimal('6.3')
True

初看起来,上面的代码好像有点奇怪,比如我们用字符串来表示数字。 然而, Decimal 对象会像普通浮点数一样的工作(支持所有的常用数学运算)。 如果你打印它们或者在字符串格式化函数中使用它们,看起来跟普通数字没什么两样。

decimal 模块的一个主要特征是允许你控制计算的每一方面,包括数字位数和四舍五入运算。 为了这样做,你先得创建一个本地上下文并更改它的设置,比如:

>>> from decimal import localcontext
>>> a = Decimal('1.3')
>>> b = Decimal('1.7')
>>> print(a / b)
0.7647058823529411764705882353
>>> with localcontext() as ctx:
...   ctx.prec = 3
...   print(a / b)
...
0.765
>>> with localcontext() as ctx:
...   ctx.prec = 50
...   print(a / b)
...
0.76470588235294117647058823529411764705882352941176
>>>

讨论

decimal 模块实现了IBM的”通用小数运算规范”。不用说,有很多的配置选项这本书没有提到。

Python新手会倾向于使用 decimal 模块来处理浮点数的精确运算。 然而,先理解你的应用程序目的是非常重要的。 如果你是在做科学计算或工程领域的计算、电脑绘图,或者是科学领域的大多数运算, 那么使用普通的浮点类型是比较普遍的做法。 其中一个原因是,在真实世界中很少会要求精确到普通浮点数能提供的17位精度。 因此,计算过程中的那么一点点的误差是被允许的。 第二点就是,原生的浮点数计算要快的多-有时候你在执行大量运算的时候速度也是非常重要的。

即便如此,你却不能完全忽略误差。数学家花了大量时间去研究各类算法,有些处理误差会比其他方法更好。 你也得注意下减法删除以及大数和小数的加分运算所带来的影响。比如:

>>> nums = [1.23e+18, 1, -1.23e+18]
>>> sum(nums) # Notice how 1 disappears
0.0
>>>

上面的错误可以利用 math.fsum() 所提供的更精确计算能力来解决:

>>> import math
>>> math.fsum(nums)
1.0
>>>

然而,对于其他的算法,你应该仔细研究它并理解它的误差产生来源。

总的来说, decimal 模块主要用在涉及到金融的领域。 在这类程序中,哪怕是一点小小的误差在计算过程中蔓延都是不允许的。 因此, decimal 模块为解决这类问题提供了方法。 当Python和数据库打交道的时候也通常会遇到 Decimal 对象,并且,通常也是在处理金融数据的时候。

以上就是Python如何执行精确的浮点数运算的详细内容,更多关于Python执行精确的浮点数运算的资料请关注三水点靠木其它相关文章!

Python 相关文章推荐
用Python编写一个简单的俄罗斯方块游戏的教程
Apr 03 Python
使用Python编写一个简单的tic-tac-toe游戏的教程
Apr 16 Python
Python学习笔记之解析json的方法分析
Apr 21 Python
Python中GIL的使用详解
Oct 03 Python
Python Flask 搭建微信小程序后台详解
May 06 Python
教你一步步利用python实现贪吃蛇游戏
Jun 27 Python
kali中python版本的切换方法
Jul 11 Python
python Gunicorn服务器使用方法详解
Jul 22 Python
分享PyCharm的几个使用技巧
Nov 10 Python
Python input函数使用实例解析
Nov 22 Python
python3.8与pyinstaller冲突问题的快速解决方法
Jan 16 Python
浅谈keras中的Merge层(实现层的相加、相减、相乘实例)
May 23 Python
Python使用shutil模块实现文件拷贝
Jul 31 #Python
Python基于pyjnius库实现访问java类
Jul 31 #Python
Python如何将字符串转换为日期
Jul 31 #Python
Python在字符串中处理html和xml的方法
Jul 31 #Python
python中selenium库的基本使用详解
Jul 31 #Python
Python过滤序列元素的方法
Jul 31 #Python
python中的django是做什么的
Jul 31 #Python
You might like
ThinkPHP验证码和分页实例教程
2014/08/22 PHP
php进行ip地址掩码运算处理的方法
2016/07/11 PHP
Thinkphp5框架实现获取数据库数据到视图的方法
2019/08/14 PHP
屏蔽网页右键复制和ctrl+c复制的js代码
2013/01/04 Javascript
给js文件传参数(详解)
2014/07/13 Javascript
Nodejs学习笔记之Global Objects全局对象
2015/01/13 NodeJs
PHP守护进程实例
2015/03/06 Javascript
jquery 中ajax执行的优先级
2015/06/22 Javascript
js鼠标按键事件和键盘按键事件用法实例汇总
2016/10/03 Javascript
vue实现可增删查改的成绩单
2016/10/27 Javascript
BootStrap 实现各种样式的进度条效果
2016/12/07 Javascript
Vue下的国际化处理方法
2017/12/18 Javascript
node.js基础知识小结
2018/02/26 Javascript
Vue自定义指令实现checkbox全选功能的方法
2018/02/28 Javascript
vue router的基本使用和配置教程
2018/11/05 Javascript
JS栈stack类的实现与使用方法示例
2019/01/31 Javascript
Vue表情输入组件 微信face表情组件
2019/02/11 Javascript
详解vue中axios请求的封装
2019/04/08 Javascript
JavaScript使用表单元素验证表单的示例代码
2019/08/20 Javascript
js 判断当前时间是否处于某个一个时间段内
2019/09/19 Javascript
基于javascript的无缝滚动动画实现2
2020/08/07 Javascript
解决vue 退出动画无效的问题
2020/08/09 Javascript
Python Web框架Flask下网站开发入门实例
2015/02/08 Python
Python3.5基础之NumPy模块的使用图文与实例详解
2019/04/24 Python
pytorch索引查找 index_select的例子
2019/08/18 Python
python_array[0][0]与array[0,0]的区别详解
2020/02/18 Python
python实现翻译word表格小程序
2020/02/27 Python
纽约复古灵感的现代珠宝品牌:Lulu Frost
2018/03/03 全球购物
墨尔本最受欢迎的复古风格品牌:Princess Highway
2018/12/21 全球购物
C#中有没有静态构造函数,如果有是做什么用的?
2016/06/04 面试题
支部鉴定材料
2014/06/02 职场文书
贷款工资证明范本
2015/06/12 职场文书
Ajax实现局部刷新的方法实例
2021/03/31 Javascript
Python OpenCV 图像平移的实现示例
2021/06/04 Python
Oracle 死锁的检测查询及处理
2021/09/25 Oracle
科学家研发出新型速效酶,可在 24 小时内降解塑料制品
2022/04/29 数码科技