解决Python中的modf()函数取小数部分不准确问题


Posted in Python onMay 28, 2021

使用math.modf()对一个浮点数进行拆分时经常会遇到如下情况

如下

import math
print(math.modf(2.4)) # 输出 (0.3999999999999999, 2.0)

我们会发现对2.4进行拆分得到的小数并不是0.4,这是因为什么呢?

这是因为计算机采用的是二进制代码,而二进制代码由于计算上的误差无法准确表示某些十进制数的小数部分。

下面我们具体来讲一下。

我们知道一个十进制数转化为二进制数需要分为两部分进行计算:整数部分和小数部分。

整数部分采用“除二取余法”。将这个整数除以2,得到它的余数,然后将商再除以3,直到商为0为止,并将各个得到的余数按照相反的顺序进行排列。

小数部分采用“乘2取整法”,将这个小数乘2,将新得到的数的整数部分取出,再用2乘余下的小数部分,如此往复直到乘积中的小数部分为0或者达到要求的精度为止。并将取出的整数部分按照取出的先后顺序从前到后排列。

123/2=61...1
 61/2=30...1
 30/2=15...0
 15/2=7...1
 7/2=3...1
 3/2=1...1
 1/2=0...1
 (123)10=(1111011)2
0.4*2=0.8...0
0.8*2=1.6...1
0.6*2=1.2...1
0.2*2=0.4...0
0.4*2=0.8...0
0.8*2=1.6...1
0.6*2=1.2...1
(0.4)10=(0.011001100110011001100110...)2
(123.4)2=(1111011.011001100110011...)2

十进制的0.4转化为二进制时会出现重复循环“0110”的情况,但是目前计算机存储浮点数是按照电器和电子工程师协会的标准(IEEE754浮点数存储格式标准)来进行存储的。

IEEE规定单精度浮点最多存储32位(4个字节):

31位是符号位。1表示该数为负,0为正

30~23是指数位(-128-127)

22~0是尾数位,尾数的编码一般是源码和补码

IEEE标准从逻辑上用三元组{S,E,M}表示一个数,如图所示:

解决Python中的modf()函数取小数部分不准确问题

也就是说上面将0.4转换出的二进制代码,我们只能存储23位,即使数据类型为double,也只能存储52位,这样大家便能看出问题出现的原因了。23位的数据显然无法完整表示0.4的二进制数据,于是误差产生了。

那所有的十进制小数都无法完整表示吗?不是的,只要小数部分乘上2最终小数部分能够得到0就不会出现这种问题,比如0.5,0.75。

import math
print(math.modf(1.5))  # 输出 (0.5,1.0)
0.5*2=1...1
(0.5)10=(0.1)2
0.75*2=1.5...1
0.5*2=1...1
(0.75)10=(0.11)2

补充:Python中“取整”的各种问题

一、初衷:

有时候我们分页展示数据的时候,需要计算页数。一般都是向上取整,例如counts=205 pageCouts=20 ,pages= 11 页。

一般的除法只是取整数部分,达不到要求。

二、方法:

1、通用除法:

UP(A/B) = int((A+B-1)/B)

取临界值,计算下A+B-1的范围就OK.

2 、Python除法:

首先要说的是python中的除法运算,

当使用x/y形式进行除法运算时,那么会进行所谓的true除法,比如2.0/3的结果是 0.66666666666666663。

另外一种除法是采用x//y的形式(向下取整),那么这里采用的是所谓floor除法,即得到不大于结果的最大整数值,这个运算时与操作数无关的。比如2//3的结果是0,-2//3的结果是-1,-2.0//3的结果是-1.0。

在python 3.0中,x/y将只执行true除法,而与操作数无关;x//y则执行floor除法。

Python运算向上取整方法:(A+B-1)/B

3、Python match.ceil函数 np.ceil函数

ceil(x)函数是向上取整,即取大于等于x的最接近整数。

import math 
math.ceil(float(205)/20) 
import numpy as np
np.ceil(float(205)/20)

Python 向上取整的算法

1、一般使用floor除法 (np.floor()或者math.floor())

import numpy as np 
n = np.array([-1.7, -2.5, -0.2, 0.6, 1.2, 2.7, 11]) 
floor = np.floor(n)
print(floor)  # [ -2.  -3.  -1.   0.   1.   2.  11.]

2、一般除法/

A=100
B=16
c=100//16
(c=6)

3、round()四舍五入函数。

np.around 返回四舍五入后的值,可指定精度。

around(a, decimals=0, out=None)

a 输入数组

decimals 要舍入的小数位数。 默认值为0。 如果为负,整数将四舍五入到小数点左侧的位置

import numpy as np 
n = np.array([-0.746, 4.6, 9.4, 7.447, 10.455, 11.555])
 
around1 = np.around(n)
print(around1)  # [ -1.   5.   9.   7.  10.  12.]
 
around2 = np.around(n, decimals=1)
print(around2)  # [ -0.7   4.6   9.4   7.4  10.5  11.6]
 
around3 = np.around(n, decimals=-1)
print(around3)  # [ -0.   0.  10.  10.  10.  10.]

Python 分别取整的算法

math模块中的 modf()方法

将整数部分和小数部分分别取出,可以使用math模块中的 modf()方法

例如:

>>> math.modf(4.25)
(0.25, 4.0)
>>> math.modf(4.33)

以上为个人经验,希望能给大家一个参考,也希望大家多多支持三水点靠木。

Python 相关文章推荐
python 递归遍历文件夹,并打印满足条件的文件路径实例
Aug 30 Python
浅谈Pandas中map, applymap and apply的区别
Apr 10 Python
深入分析python数据挖掘 Json结构分析
Apr 21 Python
Django接收自定义http header过程详解
Aug 23 Python
python实现监控阿里云账户余额功能
Dec 16 Python
python实现飞机大战游戏(pygame版)
Oct 26 Python
Python可变对象与不可变对象原理解析
Feb 25 Python
Python作用域与名字空间原理详解
Mar 21 Python
通过实例简单了解python yield使用方法
Aug 06 Python
Python基于Serializer实现字段验证及序列化
Nov 04 Python
只需要这一行代码就能让python计算速度提高十倍
May 24 Python
浅析Django接口版本控制
Jun 26 Python
利用Python+OpenCV三步去除水印
python实现自定义日志的具体方法
May 28 #Python
python 爬取京东指定商品评论并进行情感分析
python b站视频下载的五种版本
May 27 #Python
教你怎么用python selenium实现自动化测试
Python Django框架介绍之模板标签及模板的继承
May 27 #Python
python 算法题——快乐数的多种解法
May 27 #Python
You might like
索尼SONY ICF-7600A(W)电路分析
2021/03/01 无线电
15种PHP Encoder的比较
2007/03/06 PHP
php 设计模式之 工厂模式
2008/12/19 PHP
php实现无限级分类查询(递归、非递归)
2016/03/10 PHP
浅谈php中变量的数据类型判断函数
2017/03/04 PHP
Javascript与flash交互通信基础教程
2008/08/07 Javascript
40个新鲜出炉的jQuery 插件和免费教程[上]
2012/07/24 Javascript
js跨浏览器实现将字符串转化为xml对象的方法
2013/09/25 Javascript
js中split函数的使用方法说明
2013/12/26 Javascript
jQuery中[attribute!=value]选择器用法实例
2014/12/31 Javascript
javascript中的正则表达式使用详解
2015/08/30 Javascript
详解Bootstrap的iCheck插件checkbox和radio
2016/08/24 Javascript
js select实现省市区联动选择
2020/04/17 Javascript
AngularJS的ng-repeat指令与scope继承关系实例详解
2017/01/21 Javascript
jQuery插件FusionCharts实现的MSBar2D图效果示例【附demo源码】
2017/03/24 jQuery
Angular 4 指令快速入门教程
2017/06/07 Javascript
实现jquery放大镜的两种方法
2018/02/22 jQuery
JavaScript实现创建自定义对象的常用方式总结
2018/07/09 Javascript
vuejs点击class变化的实例
2018/09/05 Javascript
详解node字体压缩插件font-spider的用法
2018/09/28 Javascript
JS实现从对象获取对象中单个键值的方法示例
2019/06/05 Javascript
[13:40]TI3青蛙君全程回顾 DOTA2我们为梦想再战
2013/09/13 DOTA
[01:10:58]KG vs TNC 2019国际邀请赛小组赛 BO2 第二场 8.15
2019/08/16 DOTA
python三元运算符实现方法
2013/12/17 Python
用Python制作检测Linux运行信息的工具的教程
2015/04/01 Python
python遍历文件目录、批量处理同类文件
2019/08/31 Python
Python编写打字训练小程序
2019/09/26 Python
Python shutil模块用法实例分析
2019/10/02 Python
使用TensorFlow-Slim进行图像分类的实现
2019/12/31 Python
python 实现图像快速替换某种颜色
2020/06/04 Python
浅谈keras.callbacks设置模型保存策略
2020/06/18 Python
Python unittest装饰器实现原理及代码
2020/09/08 Python
自荐信模板大全
2015/03/27 职场文书
2015团员个人年度总结
2015/11/24 职场文书
团队拓展训练心得体会
2016/01/12 职场文书
nginx proxy_cache 缓存配置详解
2021/03/31 Servers