从0.1加0.2不等于0.3谈Python浮点数的前世今生

1. 0.1加0.2不等于0.3？

什么？0.1 加 0.2 不等于 0.3？你确定没有搞错？真的，这是千真万确的事实。不仅 Python 如此，所有浮点数规范遵从IEEEE754二进制浮点数算术标准（ANSI/IEEE Std 754-1985）的编程语言，比如 C，同样如此（如果想在C环境中验证的话，请使用 double 类型）。我们在 Python 的 IDLE 中验证一下：

>>> 0.1 + 0.2 == 0.3
False
>>> 0.1 + 0.2
0.30000000000000004

为什么会这样呢？ 0.1 + 0.2 和 0.3 究竟有什么差别？让我们用 Python 的方式来探究一番吧。我们知道，Python 的世界里，万物皆对象，浮点数也是如此，浮点数对象的方法里面究竟隐藏着怎样的秘密呢？

>>> a = 0.1 + 0.2
>>> b = 0.3
>>> a.hex()
'0x1.3333333333334p-2'
>>> b.hex()
'0x1.3333333333333p-2'

我们用对象 a 表示 0.1 + 0.2，用对象 b 表示 0.3， 浮点数对象的 hex() 方法，返回的是浮点数的十六进制表示。p-2，意思是为小数点左移2位，如果p+3，则表示右移3位。上面的代码，清晰显示了 a 和 b 的十六进制表示确实不一样。将 a 和 b 的 hex() 翻译成二进制是这样的：

a.hex(): 0.01 0011 0011 0011 0011 0011 0011 0011 0011 0011 0011 0011 0011 0100
b.hex(): 0.01 0011 0011 0011 0011 0011 0011 0011 0011 0011 0011 0011 0011 0011

2. 为什么要使用浮点数？

天地茫茫，宇宙无限，星际距离动辄以万亿公里计。然若专注于微观世界，人类认知已经进入了纳米范围之内（氢原子基态的电子轨道半径0.0528纳米）。如果使用定点数表示如此巨大的动态范围，势必耗费大量的存储资源，因此，浮点数就应运而生了。所谓浮点数，简单理解，就是科学计数法。比如：

太阳与地球之间的平均距离：150000000000米，记作 $1.5\times 10^{11}$ 米
氢原子基态的电子轨道半径：0.000000000528米，记作 $5.28\times 10^{-11}$ 米

Python 的浮点数为双精度浮点数，遵从 IEEEE754 标准，使用8个字节共64位表示一个浮点数，其中：

• 符号位：1位
• 指数位：11位
• 尾数位：52位

我们来看看双精度浮点数的有关信息。

>>> import sys
>>> sys.float_info
sys.float_info(
max=1.7976931348623157e+308, 
max_exp=1024, 
max_10_exp=308, 
min=2.2250738585072014e-308, 
min_exp=-1021, 
min_10_exp=-307, 
dig=15, 
mant_dig=53, 
epsilon=2.220446049250313e-16, 
radix=2, 
rounds=1)

Python 的双精度浮点数，可以表示的值域范围：

最小值：$2.2250738585072014\times 10^{-308}$
最大值：$1.7976931348623157\times 10^{308}$

1位符号位表示浮点数的正负，很容易理解；11位指数，对应的底数是2，可以表示的指数的整数范围应该是-1024~1023，为什么 sys.float_info 显示的 min_exp 和 max_exp 却分别是 -1021 和 1024 呢？为什么尾数是看起来这么奇怪呢？

3. 浮点数的二进制和十进制是怎样转换的？

我们知道，二进制整数从低位到高位，权重分别是1、2、4、8等；二进制小数也类似，从小数点后第1位向右，权重分别是0.5、0.25、0.125、0.0625等。据此，可以很容易地将浮点数的二进制转换成十进制。比如，0b0.1，代表十进制的0.5，0b0.01，代表十进制的0.25，0b0.11，代表十进制的0.75。
我们以十进制浮点数 9.25 为例，看看浮点数对象是如何存储的。9.25 写成二进制是：

0b1001.01

小数点左移3位，整数位只保留1位（类似于十进制的科学记数法）：

0b1.00101

现在的尾数是00101，转成十六进制（需要后补3个0变成8位）是0x28。因为小数点左移3位，指数是3。我们来验证一下：

>>> a = 9.25
>>> a.hex()
'0x1.2800000000000p+3'

结果和我们手工转换完全一致。双精度浮点数的尾数，因为默认整数位有一个1，所以双精度浮点数所能表示的最大数的指数，也就从1023变成了1024。当52位尾数全部位1，指数位1024时，代表的最大数是$2^{1024}$，变成以10为底的数，尾数应当是：

>>> import math
>>> math.log10(math.pow(2,1000)) + math.log10(math.pow(2,24)) # log10(2^1024)
308.25471555991675
>>> _ - 308
0.25471555991674677
>>> math.pow(10, _)
1.7976931348623277

这与 sys.float_info 显示的尾数 1.7976931348623157 几乎完全一致。双精度浮点数所能表示的最小数的指数和尾数，和刚才的推导基本相似，我们就不再赘述了。

下面给出浮点数的二进制和十进制互转函数：

# -*- coding: utf-8 -*-

import math

def decimal_bin(number, precision=52):
    """十进制浮点数转二进制字符串，precision为二进制小数位位数"""
    
    num_int = int(number)
    num_dec = number - num_int
    
    s = ''
    while num_dec > 0 and len(s)<precision:
        num_dec *= 2
        if num_dec >= 1:
            s += '1'
            num_dec -= 1
        else:
            s += '0'
 
    return '%s.%s'%(bin(num_int), s)
    
def bin_decimal(bin_str):
    """二进制字符串转十进制浮点数"""
    
    num_int, num_dec = bin_str.split('.')
    
    s = int(num_int, base=2)
    for i, bit in enumerate(num_dec):
        if bit == '1':
            s += math.pow(2, -(i+1))
        
    return s

4. 如何实现0.1加0.2等于0.3？

通过刚才的分析，我们已经认识到了浮点数的本质：在有限的精度内，我们可以精确的设置浮点数对象的值，但一旦参与运算，因为精度误差的原因，将会导致计算结果和预期的数学结果不能完全一致。那么，在浮点运算时，我们有没有方法实现计算结果和预期的数学结果保持一致呢？答案是：有！万能的 Python 已经为我们准备了模块 decimal，专门解决这个问题。请看演示：

>>> from decimal import Decimal
>>> a = 0.1
>>> b = 0.2
>>> c = float(Decimal(str(a)) + Decimal(str(b)))
>>> c
0.3

decimal 模块还有很多功能，有兴趣的同学可以在网上找到更多的学习资料。