回答

收藏

浮点数被破坏了吗?

技术问答 技术问答 366 人阅读 | 0 人回复 | 2023-09-11

考虑以下代码:
1 @' t# ]& C/ x6 S. f" F, y0.1   0.2 == 0.3  ->  false0.1   0.2                ->  0.30000000000000004为什么会出现这些错误?
: b% b; L: ~6 [- N                                                                / {" S% T6 @; I
    解决方案:                                                               
, p6 Q5 J. R$ q+ w! G                                                                二进制浮点数学就是这样。它基于大多数编程语言IEEE 754 标准。问题的关键在于,数字以这种格式表示为整数乘以 2 的幂。不能准确表示分母不是 2 次方的有理数(如0.1,即)。1/10
) b8 b/ i( d% g因为0.1在标准binary在64格式中,表示可以完全写成
/ k. ^3 ?/ x9 |% ?6 \0.1000000000055151257270218340454015625 十进制* s) S. H+ b- c' J% N/ ?
0x1.999999999999ap-4在C99 hexfloat 在表示法中。
相反,有理数0.1,即1/10,可以完全写成
- Z% l/ \7 L' [. \% o0.或
4 ]- Q8 w6 B0 {0x1.99999999999999...p-4在 C99 hexfloat 表示法的类似物,其中...表示 9 无休止序列。
常数0在程序中.2和0.3也类似于它们的真实值。碰巧最接近double的0.2大于有理数0.2,但最接近double的0.3小于有理数0.3。0.1和的总和0.最终大于理数0.因此,它与代码中的常数不一致。
& {, Y& m% P  {' o5 p1 g9 R8 P( ~浮点算术问题的综合处理是每个计算机科学家都应该知道浮点算术的知识。请参考更容易理解的解释floating-point-gui.de。5 {; a2 d+ I: a7 _0 _
旁注:所有位置(base-N)数字系统都精确地共享这个问题
* ]+ s& ~8 ]* ]老十进制(以 10 为底)的普通数字也有同样的问题,这就是为什么 1/3 这样的数字最终会变成 0.333333333…
. ?: w8 @* B3 Z1 C& l* E你只是偶然发现了一个数字 (3/10),很容易用十进制系统表示,但不适合二进制系统。它也是双向的(在某种程度上):1/16 是十进制的丑数(0.0625),但在二进制中,它看起来像十进制的万分之一(0.0001)** - 如果我们在日常生活中使用 2 的数字系统,你甚至会看到这个数字,本能地理解你可以通过减半一次又一次地到达那里。( w) |4 S6 Z# w: L: g
** 当然,内存中浮点的存储方法并不完全正确(它们使用科学的记数法)。然而,它确实表明,这经常发生在二进制浮点的精度错误中,因为我们通常对现实世界感兴趣的数字通常是10米 -只是因为我们使用了十进制数字系统 -今天。这就是为什么我们会说 71% ,而不是每 7 中有 5 (71% 是近似值,因为 5/7 不能用任何十进制数准确表示)。1 P0 _* P; q  E& L7 H4 |
所以不:二进制浮点数没有被破坏,只是碰巧和其他一切都是基于 N 数字系统同样不完美9 T3 ?" |8 q# J0 @& R7 Z3 e
旁注:编程中使用浮点数: ?; ~' J" @, Z$ D) H9 E/ e3 o. g+ i; C
在实践中,这个精度问题意味着你需要用舍入函数四舍五入你感兴趣的小数字,然后显示它们。, j% E& j3 c+ ~0 L, Y; }
您还需要使用允许一定容差的比较来替换相等测试,这意味着:
  _* [- t% _( Y8 a5 B( j不要做_if (x == y) { ... }+ m+ E' ]- u* f
而是做if (abs(x - y)
! T. T+ Z/ c6 o" v2 b其中abs是绝对值。myToleranceValue您需要为您的特定应用程序选择 - 这与你准备的摆动空间和你想要比较的最大数字(精度损失)有很大关系)。请注意您选择的语言中的epsilon样式常量不能作为公差值。
分享到:
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则