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简单判断序列是否为空的方法
Jun 30 Python
Golang与python线程详解及简单实例
Apr 27 Python
python构建自定义回调函数详解
Jun 20 Python
Python人脸识别初探
Dec 21 Python
如何安装并在pycharm使用selenium的方法
Apr 30 Python
Python3 socket即时通讯脚本实现代码实例(threading多线程)
Jun 01 Python
python编写一个会算账的脚本的示例代码
Jun 02 Python
Python类super()及私有属性原理解析
Jun 15 Python
Python爬虫基于lxml解决数据编码乱码问题
Jul 31 Python
Python使用urlretrieve实现直接远程下载图片的示例代码
Aug 17 Python
python解决OpenCV在读取显示图片的时候闪退的问题
Feb 23 Python
Python还能这么玩之只用30行代码从excel提取个人值班表
Jun 05 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
一周学会PHP(视频)Http下载
2006/12/12 PHP
AJAX的使用方法详解
2017/04/29 PHP
js自定义事件及事件交互原理概述(二)
2013/02/01 Javascript
JavaScript中几个重要的属性(this、constructor、prototype)介绍
2013/05/19 Javascript
JavaScript中的Math.sin()方法使用详解
2015/06/15 Javascript
Jquery1.9.1源码分析系列(六)延时对象应用之jQuery.ready
2015/11/24 Javascript
基于JavaScript实现简单的随机抽奖小程序
2016/01/05 Javascript
JS实现关闭当前页而不弹出提示框的方法
2016/06/22 Javascript
Web打印解决方案之证件套打的实现思路
2016/08/29 Javascript
jquery判断类型是不是number类型的实例代码
2016/10/07 Javascript
Bootstrap 网格系统布局详解
2017/03/19 Javascript
Vue响应式添加、修改数组和对象的值
2017/03/20 Javascript
vue+mockjs模拟数据实现前后端分离开发的实例代码
2017/08/08 Javascript
Bootstrap treeview实现动态加载数据并添加快捷搜索功能
2018/01/07 Javascript
vue 动态绑定背景图片的方法
2018/08/10 Javascript
vue-cli4.x创建企业级项目的方法步骤
2020/06/18 Javascript
vue Treeselect 树形下拉框:获取选中节点的ids和lables操作
2020/08/15 Javascript
[49:13]DOTA2上海特级锦标赛C组资格赛#1 OG VS LGD第一局
2016/02/27 DOTA
跟老齐学Python之Import 模块
2014/10/13 Python
Python常用算法学习基础教程
2017/04/13 Python
Python+Turtle动态绘制一棵树实例分享
2018/01/16 Python
python编写暴力破解zip文档程序的实例讲解
2018/04/24 Python
用python-webdriver实现自动填表的示例代码
2021/01/13 Python
canvas实现扭蛋机动画效果的示例代码
2018/10/17 HTML / CSS
您在慕尼黑的跑步商店:Lauf-bar
2019/10/11 全球购物
初中毕业典礼演讲稿
2014/09/09 职场文书
查摆问题自查报告范文
2014/10/13 职场文书
校本课程教学计划
2015/01/19 职场文书
大学生毕业个人总结
2015/02/15 职场文书
计算机专业自荐信
2015/03/05 职场文书
2015年事业单位工作总结
2015/04/27 职场文书
2016秋季运动会前导词
2015/11/25 职场文书
导游词之安徽醉翁亭
2020/01/10 职场文书
总结Pyinstaller打包的高级用法
2021/06/28 Python
Spring Security中用JWT退出登录时遇到的坑
2021/10/16 Java/Android
SQL Server使用CROSS APPLY与OUTER APPLY实现连接查询
2022/05/25 SQL Server