PHP浮点数的一个常见问题


Posted in PHP onMarch 10, 2016

PHP是一种弱类型语言, 这样的特性, 必然要求有无缝透明的隐式类型转换, PHP内部使用zval来保存任意类型的数值, zval的结构如下(5.2为例):

struct _zval_struct {
 /* Variable information */
 zvalue_value value;  /* value */
 zend_uint refcount;
 zend_uchar type; /* active type */
 zend_uchar is_ref;
};

上面的结构中, 实际保存数值本身的是zvalue_value联合体:

typedef union _zvalue_value {
 long lval;     /* long value */
 double dval;    /* double value */
 struct {
  char *val;
  int len;
 } str;
 HashTable *ht;    /* hash table value */
 zend_object_value obj;
} zvalue_value;

今天的话题, 我们只关注其中的俩个成员, lval和dval, 我们要意识到, long lval是随着编译器, OS的字长不同而不定长的, 它有可能是32bits或者64bits, 而double dval(双精度)由IEEE 754规定, 是定长的, 一定是64bits.

请记住这一点, 造就了PHP的一些代码的”非平台无关性”. 我们接下来的讨论, 除了特别指明, 都是假设long为64bits

IEEE 754的浮点计数法, 我这里就不引用了, 大家有兴趣的可以自己查看, 关键的一点是, double的尾数采用52位bit来保存, 算上隐藏的1位有效位, 一共是53bits.

在这里, 引出一个很有意思的问题, 我们用c代码举例(假设long为64bits):

long a = x;
 assert(a == (long)(double)a);

请问, a的取值在什么范围内的时候, 上面的代码可以断言成功?(留在文章最后解答)

现在我们回归正题, PHP在执行一个脚本之前, 首先需要读入脚本, 分析脚本, 这个过程中也包含着, 对脚本中的字面量进行zval化, 比如对于如下脚本:

<?php
$a = 9223372036854775807; //64位有符号数最大值
$b = 9223372036854775808; //最大值+1
var_dump($a);
var_dump($b);

输出:

int(9223372036854775807)
float(9.22337203685E+18)

也就说, PHP在词法分析阶段, 对于一个字面量的数值, 会去判断, 是否超出了当前系统的long的表值范围, 如果不是, 则用lval来保存, zval为IS_LONG, 否则就用dval表示, zval IS_FLOAT.

凡是大于最大的整数值的数值, 我们都要小心, 因为它可能会有精度损失:

<?php
$a = 9223372036854775807;
$b = 9223372036854775808;
 
var_dump($a === ($b - 1));

输出是false.

现在接上开头的讨论, 之前说过, PHP的整数, 可能是32位, 也可能是64位, 那么就决定了, 一些在64位上可以运行正常的代码, 可能会因为隐形的类型转换, 发生精度丢失, 从而造成代码不能正常的运行在32位系统上.

所以, 我们一定要警惕这个临界值, 好在PHP中已经定义了这个临界值:

<?php
 echo PHP_INT_MAX;
 ?>

当然, 为了保险起见, 我们应该使用字符串来保存大整数, 并且采用比如bcmath这样的数学函数库来进行计算.

另外, 还有一个关键的配置, 会让我们产生迷惑, 这个配置就是php.precision, 这配置决定了PHP再输出一个float值的时候, 输出多少有效位.

最后, 我们再来回头看上面提出的问题, 也就是一个long的整数, 最大的值是多少, 才能保证转到float以后再转回long不会发生精度丢失?

比如, 对于整数, 我们知道它的二进制表示是, 101, 现在, 让我们右移俩位, 变成1.01, 舍去高位的隐含有效位1, 我们得到在double中存储5的二进制数值为:

0/*符号位*/ 10000000001/*指数位*/ 0100000000000000000000000000000000000000000000000000
5的二进制表示, 丝毫未损的保存在了尾数部分, 这个情况下, 从double转会回long, 不会发生精度丢失.

我们知道double用52位表示尾数, 算上隐含的首位1, 一共是53位精度.. 那么也就可以得出, 如果一个long的整数, 值小于:

2^53 - 1 == 9007199254740991; //牢记, 我们现在假设是64bits的long
那么, 这个整数, 在发生long->double->long的数值转换时, 不会发生精度丢失.

关于浮点数,还有一点,就是对于如下的这个常见问题的回答:

<?php
 $f = 0.58;
 var_dump(intval($f * 100)); //为啥输出57
?>

为啥输出是57啊? PHP的bug么?

我相信有很多的同学有过这样的疑问, 因为光问我类似问题的人就很多, 更不用说bugs.php.net上经常有人问…

要搞明白这个原因, 首先我们要知道浮点数的表示(IEEE 754):

浮点数, 以64位的长度(双精度)为例, 会采用1位符号位(E), 11指数位(Q), 52位尾数(M)表示(一共64位).

符号位:最高位表示数据的正负,0表示正数,1表示负数。

指数位:表示数据以2为底的幂,指数采用偏移码表示

尾数:表示数据小数点后的有效数字.

这里的关键点就在于, 小数在二进制的表示, 关于小数如何用二进制表示, 大家可以百度一下, 我这里就不再赘述, 我们关键的要了解, 0.58 对于二进制表示来说, 是无限长的值(下面的数字省掉了隐含的1)..

0.58的二进制表示基本上(52位)是: 0010100011110101110000101000111101011100001010001111
0.57的二进制表示基本上(52位)是: 0010001111010111000010100011110101110000101000111101

而两者的二进制, 如果只是通过这52位计算的话,分别是:

0.58 -> 0.57999999999999996
0.57 -> 0.56999999999999995

至于0.58 * 100的具体浮点数乘法, 我们不考虑那么细, 有兴趣的可以看(Floating point), 我们就模糊的以心算来看… 0.58 * 100 = 57.999999999

那你intval一下, 自然就是57了….

可见, 这个问题的关键点就是: “你看似有穷的小数, 在计算机的二进制表示里却是无穷的”

so, 不要再以为这是PHP的bug了, 这就是这样的…..

PHP 相关文章推荐
PHP在字符断点处截断文字的实现代码
Apr 21 PHP
有关php运算符的知识大全
Nov 03 PHP
浅析php中常量,变量的作用域和生存周期
Aug 10 PHP
ThinkPHP CURD方法之limit方法详解
Jun 18 PHP
phpword插件导出word文件时中文乱码问题处理方案
Aug 19 PHP
php简单分页类实现方法
Feb 26 PHP
简单谈谈favicon
Jun 10 PHP
最新最全PHP生成制作验证码代码详解(推荐)
Jun 12 PHP
老生常谈php中传统验证与thinkphp框架(必看篇)
Jun 10 PHP
PHP开发之归档格式phar文件概念与用法详解【创建,使用,解包还原提取】
Nov 17 PHP
php从数据库中获取数据用ajax传送到前台的方法
Aug 20 PHP
PHP压缩图片功能的介绍
Mar 21 PHP
简单谈谈php浮点数精确运算
Mar 10 #PHP
PHP实现仿百度文库,豆丁在线文档效果(word,excel,ppt转flash)
Mar 10 #PHP
Zend Framework教程之Loader以及PluginLoader用法详解
Mar 09 #PHP
php注册登录系统简化版
Dec 28 #PHP
详解WordPress中用于更新和获取用户选项数据的PHP函数
Mar 08 #PHP
Zend Framework教程之Autoloading用法详解
Mar 08 #PHP
Zend Framework教程之Resource Autoloading用法实例
Mar 08 #PHP
You might like
php学习之简单计算器实现代码
2011/06/09 PHP
判断、添加和删除WordPress置顶文章的相关PHP函数小结
2015/12/10 PHP
php设计模式之策略模式应用案例详解
2019/06/17 PHP
解决windows上php xdebug 无法调试的问题
2020/02/19 PHP
用php实现分页效果的示例代码
2020/12/10 PHP
javascript脚本调试方法小结
2008/11/24 Javascript
ASP.NET中使用后端代码注册脚本 生成JQUERY-EASYUI的界面错位的解决方法
2010/06/12 Javascript
JS焦点图切换,上下翻转
2011/05/12 Javascript
jQuery插件-jRating评分插件源码分析及使用方法
2012/12/28 Javascript
鼠标放在图片上显示大图的JS代码
2013/03/26 Javascript
JS获取当前日期时间并定时刷新示例
2021/03/04 Javascript
抛弃Nginx使用nodejs做反向代理服务器
2014/07/17 NodeJs
简介JavaScript中用于处理正切的Math.tan()方法
2015/06/15 Javascript
javascript中判断json的方法总结
2015/08/27 Javascript
关于事件mouseover ,mouseout ,mouseenter,mouseleave的区别
2015/10/12 Javascript
jquery实现简单的全选和反选功能
2016/01/02 Javascript
Three.js利用性能插件stats实现性能监听的方法
2017/09/25 Javascript
JS实现的3des+base64加密解密算法完整示例
2018/05/18 Javascript
Python3操作SQL Server数据库(实例讲解)
2017/10/21 Python
Python实现的生成格雷码功能示例
2018/01/24 Python
python实现两个dict合并与计算操作示例
2019/07/01 Python
python实现的发邮件功能示例
2019/09/11 Python
matplotlib 对坐标的控制,加图例注释的操作
2020/04/17 Python
聊聊python中的循环遍历
2020/09/07 Python
使用python画出逻辑斯蒂映射(logistic map)中的分叉图案例
2020/12/11 Python
HTML5中FileReader接口使用方法实例详解
2017/08/26 HTML / CSS
HTML5 3D衣服摇摆动画特效
2016/03/17 HTML / CSS
Myprotein法国官网:欧洲第一运动营养品牌
2019/03/26 全球购物
GC是什么?为什么要有GC?
2013/12/08 面试题
妇产医师自荐信
2014/01/29 职场文书
2014年采购员工作总结
2014/11/18 职场文书
研究生个人学年总结
2015/02/14 职场文书
教师工作能力自我评价
2015/03/04 职场文书
golang操作redis的客户端包有多个比如redigo、go-redis
2022/04/14 Golang
Python+Pillow+Pytesseract实现验证码识别
2022/05/11 Python
CSS子盒子水平和垂直居中的五种方法
2022/07/23 HTML / CSS