浮点数的使用方法图解(正确认识浮点数)
大家在日常开发中,必然使用过浮点数,也会发现浮点数不是精确的,那究竟是什么原因造成的呢?
奇怪的结果
var_dump((1-0.9) == 0.1); //输出:bool(false)
很奇怪吧!1-0.9怎么能不等于0.1呢?这是为什么呢?这要从浮点数的储存标准开始说。
IEEE 754浮点数在计算机中是根据IEEE 754(二进制浮点数算数标准)储存的。计算公式为: (-1)^S x M x 2^E
32位单精度储存结构(对应占位)
符号(S)
阶码(E)
尾数(M)
1
8
23
64位双精度储存结构(对应占位)
符号(S)
阶码(E)
尾数(M)
1
11
52
解释:
- S: 符号(0正,1负)
- E: 阶码(指数)
- M: 尾数(二进制小数,数字的实体部分)
M(尾数)和E(阶码)不同情况需要分别对待
E(阶码)的三种状态及对应的M表示
从图中(截图于深入理解计算机系统)我们可以分为三种情况(第三种又分为两种特殊情况)
规格化E既不等于0也不等于255(将S按十进制计算),这个时候的E=E-127,M的二进制小数默认省略了1.,也就是M=1.M(二进制小数)我们做一个简单的测试看一下二进制00111110001000000000000000000000(32位)表示的对应的浮点数为多少?
- 首先拆分二进制: 0 01111100 01000000000000000000000
- E = 124 = 124 - 127 = -3
- M = 1.01000000000000000000000
- 套公式: 1 x 1.01000000000000000000000 x 2^-3 = 0.00101000000000000000000000 = 2^-3 2^-5 = 0.15625
使用PHP验证一下结果:
var_dump(unpack('f', pack('l', bindec('00111110001000000000000000000000')))[1]); //输出: float(0.15625)
上面的例子没有丢失精度,下面看一个丢失精度的例子:
printf('2s', decbin(unpack('l', pack('f', 1/3))[1])); //输出: 00111110101010101010101010101011 var_dump(unpack('f', pack('l', bindec('00111110101010101010101010101011')))[1]); float(0.33333334326744)
丢失精度最主要原因就在于M(二进制小数),我们只能精确的表示2^n倍数的数(2^-1(0.5),2^-2(0.25),2^-3(0.125)...),丢了在所难免。
非规格化E等于0,这个时候E=-126,M的二进制小数前缀为0.,也就是M=0.M(二进制小数),具体过程就不写了,和上面类似
特殊情况E等于255(全部位都为1),如果M全部为0,那么表示为无穷大,否则表示为NAN(不是一个数)
不要比较浮点数
var_dump(unpack('f', pack('l', bindec('01111111100000000000000000000000')))[1]); //输出: float(INF) var_dump(unpack('f', pack('l', bindec('01111111100000000000000000000110')))[1]); //输出: float(NAN)
总之,浮点数是不准确的。尤其在我们日常工作中,不要比较浮点数的大小,如果需要精确的比较计算,请使用bc*系列函数。还有一点,浮点数不准确和PHP没有任何关系,PHP不背这个锅。
,
免责声明:本文仅代表文章作者的个人观点,与本站无关。其原创性、真实性以及文中陈述文字和内容未经本站证实,对本文以及其中全部或者部分内容文字的真实性、完整性和原创性本站不作任何保证或承诺,请读者仅作参考,并自行核实相关内容。文章投诉邮箱:anhduc.ph@yahoo.com