解决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使用MONGODB入门实例
May 11 Python
python下MySQLdb用法实例分析
Jun 08 Python
用python找出那些被“标记”的照片
Apr 20 Python
Python实现PS图像明亮度调整效果示例
Jan 23 Python
Python实现抢购IPhone手机
Feb 07 Python
python 列表,数组,矩阵两两转换tolist()的实例
Apr 04 Python
Python适配器模式代码实现解析
Aug 02 Python
python 定时器每天就执行一次的实现代码
Aug 14 Python
Python读取配置文件(config.ini)以及写入配置文件
Apr 08 Python
python函数中将变量名转换成字符串实例
May 11 Python
使用tensorflow根据输入更改tensor shape
Jun 23 Python
Django中日期时间型字段进行年月日时分秒分组统计
Nov 27 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
基于mysql的bbs设计(二)
2006/10/09 PHP
简单介绍下 PHP5 中引入的 MYSQLI的用途
2007/03/19 PHP
php页面缓存方法小结
2015/01/10 PHP
php页面跳转session cookie丢失导致不能登录等问题的解决方法
2016/12/12 PHP
PHP基于rabbitmq操作类的生产者和消费者功能示例
2018/06/16 PHP
浅谈Laravel中的三种中间件的作用
2019/10/13 PHP
js获取单选按钮的数据
2006/11/27 Javascript
js 实现打印网页中定义的部分内容的代码
2010/04/01 Javascript
关于js注册事件的常用方法
2013/04/03 Javascript
jquery js 重置表单 reset()具体实现代码
2013/08/05 Javascript
JavaScript代码简单实现求杨辉三角给定行的最大值
2013/10/29 Javascript
倒记时60刷新网页的js代码
2014/02/18 Javascript
使用不同的方法结合/合并两个JS数组
2014/09/18 Javascript
浅谈下拉菜单中的Option对象
2015/05/10 Javascript
jQuery 3.0 的 setter和getter 模式详解
2016/07/11 Javascript
js实现精确到毫秒的倒计时效果
2016/08/05 Javascript
AngularJS实现使用路由切换视图的方法
2017/01/24 Javascript
react脚手架如何配置less和ant按需加载的方法步骤
2018/11/28 Javascript
在Python的Django框架中编写编译函数
2015/07/20 Python
python之matplotlib学习绘制动态更新图实例代码
2018/01/23 Python
Python基于Flask框架配置依赖包信息的项目迁移部署
2018/03/02 Python
python画柱状图--不同颜色并显示数值的方法
2018/12/13 Python
Python操作Word批量生成合同的实现示例
2020/08/28 Python
canvas烟花特效锦集
2018/01/17 HTML / CSS
HTML5 form标签之解放表单验证、增加文件上传、集成拖放的使用方法
2013/04/24 HTML / CSS
C语言笔试题回忆
2015/04/02 面试题
Unix里面如何在后台运行程序
2016/10/14 面试题
历史专业个人求职信范文
2013/12/07 职场文书
建筑工地文明标语
2014/10/09 职场文书
2014年十八届四中全会思想汇报范文
2014/10/17 职场文书
高中校园广播稿
2014/10/21 职场文书
简历自我评价模板
2015/03/11 职场文书
圣贤教育改变命运观后感
2015/06/16 职场文书
丧事主持词
2015/07/02 职场文书
工作报告范文
2019/06/20 职场文书
如何设计高效合理的MySQL查询语句
2021/05/26 MySQL