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 01 Python
使用Python解析JSON数据的基本方法
Oct 15 Python
简单实现Python爬取网络图片
Apr 01 Python
python smtplib模块自动收发邮件功能(二)
May 22 Python
win7下python3.6安装配置方法图文教程
Jul 31 Python
python学生信息管理系统(完整版)
Apr 05 Python
分享Python切分字符串的一个不错方法
Dec 14 Python
python爬虫开发之使用python爬虫库requests,urllib与今日头条搜索功能爬取搜索内容实例
Mar 10 Python
Python PyQt5整理介绍
Apr 01 Python
Pytest allure 命令行参数的使用
Apr 18 Python
Python实现为PDF去除水印的示例代码
Apr 03 Python
Python Pandas解析读写 CSV 文件
Apr 11 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邮件发送,php发送邮件的类
2011/03/24 PHP
PHP教程之PHP中shell脚本的使用方法分享
2012/02/23 PHP
解析php 版获取重定向后的地址(代码)
2013/06/26 PHP
jquery不支持toggle()高(新)版本的问题解决
2016/09/24 PHP
PHP脚本自动识别验证码查询汽车违章
2016/12/20 PHP
php模拟post提交请求调用接口示例解析
2020/08/07 PHP
用jQuery模拟select下拉框的简单示例代码
2014/01/26 Javascript
jquery处理json对象
2014/11/03 Javascript
jQuery实现表格行上下移动和置顶效果
2015/06/05 Javascript
jQuery代码实现发展历程时间轴特效
2015/07/30 Javascript
jQuery添加和删除指定标签的方法
2015/12/16 Javascript
js精准的倒计时函数分享
2016/06/29 Javascript
Nodejs下DNS缓存问题浅析
2016/11/16 NodeJs
ES6概念 ymbol.for()方法
2016/12/25 Javascript
使用Browserify来实现CommonJS的浏览器加载方法
2017/05/14 Javascript
深入探究angular2 UI组件之primeNG用法
2017/07/26 Javascript
轻松理解vue的双向数据绑定问题
2017/10/30 Javascript
在vue中通过axios异步使用echarts的方法
2018/01/13 Javascript
JS定时器如何实现提交成功提示功能
2020/06/12 Javascript
原生JS实现微信通讯录
2020/06/18 Javascript
Vue解决移动端弹窗滚动穿透问题
2020/12/15 Vue.js
[02:21]十步杀一人,千里不留行——DOTA2全新英雄天涯墨客展示
2018/08/29 DOTA
python操作日期和时间的方法
2014/03/11 Python
Python程序中的观察者模式结构编写示例
2016/05/27 Python
Python的Flask框架应用程序实现使用QQ账号登录的方法
2016/06/07 Python
老生常谈Python之装饰器、迭代器和生成器
2017/07/26 Python
解决Python下imread,imwrite不支持中文的问题
2018/12/05 Python
在python中实现同行输入/接收多个数据的示例
2019/07/20 Python
python和php学习哪个更有发展
2020/06/17 Python
python利用线程实现多任务
2020/09/18 Python
搭建pypi私有仓库实现过程详解
2020/11/25 Python
同学会邀请书大全
2014/01/12 职场文书
网上开店必备创业计划书
2014/01/26 职场文书
四风存在的原因分析
2014/02/11 职场文书
中学生秋季运动会广播稿
2014/09/21 职场文书
【海涛解说】暗牧也疯狂,牛蛙成配角
2022/04/01 DOTA