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
python如何查看系统网络流量的信息
Sep 12 Python
利用 Monkey 命令操作屏幕快速滑动
Dec 07 Python
Python实现PS滤镜碎片特效功能示例
Jan 24 Python
将TensorFlow的模型网络导出为单个文件的方法
Apr 23 Python
python issubclass 和 isinstance函数
Jul 25 Python
安装python及pycharm的教程图解
Oct 10 Python
关于tf.TFRecordReader()函数的用法解析
Feb 17 Python
Django单元测试中Fixtures的使用方法
Feb 26 Python
Python动态导入模块:__import__、importlib、动态导入的使用场景实例分析
Mar 30 Python
PyInstaller的安装和使用的详细步骤
Jun 02 Python
python中__slots__节约内存的具体做法
Jul 04 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 header()函数使用说明
2008/07/10 PHP
Yii实现多数据库主从读写分离的方法
2014/12/29 PHP
PHP类相关知识点实例总结
2016/09/28 PHP
PHP实现的常规正则验证helper公共类完整实例
2017/04/27 PHP
PHP堆栈调试操作简单示例
2018/06/15 PHP
php微信公众号开发之音乐信息
2018/10/20 PHP
laravel框架中控制器的创建和使用方法分析
2019/11/23 PHP
PHP执行系统命令函数实例讲解
2021/03/03 PHP
js no-repeat写法 背景不重复
2009/03/18 Javascript
围观tangram js库
2010/12/28 Javascript
jQuery实现公告文字左右滚动的实例代码
2013/10/29 Javascript
seaJs的模块定义和模块加载浅析
2014/06/06 Javascript
JQuery做的一个简单的点灯游戏分享
2014/07/16 Javascript
动态加载jQuery的两种方法实例分析
2015/07/17 Javascript
浅谈JavaScript的Polymer框架中的事件绑定
2015/07/29 Javascript
AngularJS入门心得之directive和controller通信过程
2016/01/25 Javascript
Jquery操作cookie记住用户名
2016/03/29 Javascript
原生js仿淘宝网商品放大镜效果
2017/02/28 Javascript
详解AngularJs ui-router 路由的简单介绍
2017/04/26 Javascript
对于js垃圾回收机制的理解
2017/09/14 Javascript
Vue 动态路由的实现及 Springsecurity 按钮级别的权限控制
2019/09/05 Javascript
JS如何调用WebAssembly编译出来的.wasm文件
2020/11/05 Javascript
[42:20]Secret vs Liquid 2019国际邀请赛小组赛 BO2 第二场 8.15
2019/08/17 DOTA
python scipy求解非线性方程的方法(fsolve/root)
2018/11/12 Python
django认证系统实现自定义权限管理的方法
2019/08/28 Python
python实现多进程按序号批量修改文件名的方法示例
2019/12/30 Python
python小白切忌乱用表达式
2020/05/29 Python
Django filter动态过滤与排序实现过程解析
2020/11/26 Python
HTML5拖拽功能实现的拼图游戏
2018/07/31 HTML / CSS
Python面试题:Python里面如何生成随机数
2015/03/12 面试题
运动会通讯稿50字
2014/01/30 职场文书
网络工程专业自荐信范文
2014/03/16 职场文书
三分钟演讲稿范文
2014/04/24 职场文书
国际商务英语专业求职信
2014/07/08 职场文书
学习退步检讨书
2014/09/28 职场文书
javascript数组includes、reduce的基本使用
2021/07/02 Javascript