考虑以下代码:
0.1 + 0.2 == 0.3 -> false 0.1 + 0.2 -> 0.30000000000000004
为什么会出现这些错误?
二进制浮点数学就是这样。在大多数编程语言中,它基于IEEE754标准。问题的症结在于数字以这种格式表示为整数乘以2的幂。分母不是2的幂的有理数(例如0.1,是1/10)无法精确表示。
0.1
1/10
对于0.1标准binary64格式,表示形式可以完全按照
binary64
0.1000000000000000055511151231257827021181583404541015625 以十进制表示,或 0x1.999999999999ap-4以C99十六进制表示法表示。 相比之下,合理数量0.1,这是1/10可以完全按照书面
0.1000000000000000055511151231257827021181583404541015625
0x1.999999999999ap-4
0.1 以十进制表示,或 0x1.99999999999999...p-4以C99十六进制表示法的类似形式表示,其中…表示9的无休止序列。 常量0.2和0.3程序中的常量也将接近其真实值。碰巧的是,最近double以0.2低于合理数量较大0.2但最近double以0.3低于合理数量较小0.3。的总和0.1和0.2卷起比有理数较大0.3,并因此与在代码中不同意恒定。
0x1.99999999999999...p-4
double
0.3
0.2
每个计算机科学家都应该了解浮点算术,这是对浮点算术问题的相当全面的处理。有关更易于理解的说明,请参见float-point-gui.de。
float-point-gui.de
旁注:所有位置(以N为底的)数字系统均会精确地共享此问题
普通的旧十进制数(以10为底)有相同的问题,这就是为什么像1/3这样的数字最终会变成0.333333333 ...
0.333333333 ...
您偶然发现了一个数字(3/10),该数字很容易用十进制表示,但不适合二进制。它也是双向的(在某种程度上):1/16是一个丑陋的数字,十进制(0.0625),但是在二进制中,它看起来像十进制的10,000十分整洁(0.0001)****-如果我们在习惯在我们的日常生活中使用以2为底的数字系统,您甚至会看到该数字,并且本能地了解到您可以通过将某物减半,一次又一次减半而到达那里。
**当然,这并不完全是将浮点数存储在内存中的方式(它们使用科学计数法的形式)。但是,它确实说明了二进制浮点精度误差趋于增加的观点,因为我们通常感兴趣的“真实世界”数通常是10的幂-但这仅仅是因为我们使用了十进制数天-今天。这也是为什么我们说类似71%而不是“每7个中有5个”(71%是近似值,因为5/7不能用任何十进制数字精确表示)的原因。
否:二进制浮点数没有坏,它们恰好与其他所有base-N数系统一样不完美:)
侧面说明:在编程中使用浮点数
实际上,这种精度问题意味着您需要使用舍入函数将浮点数四舍五入为您想要的任意小数位,然后再显示它们。
您还需要用允许一定程度的容忍的比较替换相等性测试,这意味着:
千万不能做if (x == y) { ... }
if (x == y) { ... }
反而做if (abs(x - y) < myToleranceValue) { ... }。
if (abs(x - y) < myToleranceValue) { ... }
abs绝对值在哪里。myToleranceValue需要为您的特定应用选择-这与您准备允许多少“摆动空间”以及要比较的最大数字有很大关系(由于精度问题) )。当心所选语言中的“ epsilon”样式常量。这些不得用作公差值。
myToleranceValue