导读：浮点数运算是一个非常有技术含量的话题不太容易掌握。许多程序员都不清楚使用==操作符比较float/double类型的话到底出现什么问题 许多人使用float/double进行货币计算时经常会犯錯。这篇文章是这一系列中的精华所有的软件开发人员都应该读一下。

　　随着你经验的增长你肯定 想去深入了解一些常见的东西的細节，浮点数运算就是其中之一

　　在计算机系统的发展过程中，曾经提出过多种方法表达实数

　　【1】典型的比如相对于浮点数的萣点数（Fixed Point Number）。在这种表达方式中小数点固定的位于实数所有数字中间的某个位置。货币的表达就可以使用这种方式比如 99.00 或者 00.99 可以用于表达具有四位精度（Precision），小数点后有两位的货币值由于小数点位置固定，所以可以直接用四位数值来表达相应的数值SQL 中的 NUMBER 数据类型就昰利用定点数来定义的。

　　【2】还有一种提议的表达方式为有理数表达方式即用两个整数的比值来表达实数。

　　定点数表达法的缺點在于其形式过于僵硬固定的小数点位置决定了固定位数的整数部分和小数部分，不利于同时表达特别大的数或者特别小的数最终，絕大多数现代的计算机系统采纳了所谓的浮点数表达方式

　　【3】浮点数表达方式， 这种表达方式利用科学计数法来表达实数即用一個尾数（Mantissa ），一个基数（Base）一个指数（Exponent）以及一个表示正负的符号来表达实数。比如 123.45 用十进制科学计数法可以表达为 1.2345 × 102 其中 1.2345 为尾数，10 為基数2 为指数。浮点数利用指数达到了浮动小数点的效果从而可以灵活地表达更大范围的实数。提示: 尾数有时也称为有效数字（Significand）尾数实际上是有效数字的非正式说法。

　　同样的数值可以有多种浮点数表达方式比如上面例子中的 123.45 可以表达为 12.345 × 101，0.12345 × 103 或者 1.2345 × 102因为这種多样性，有必要对其加以规范化以达到统一表达的目标规范的（Normalized）浮点数表达方式具有如下形式：

　　其中 d.dd...d 即尾数，β 为基数e 为指數。尾数中数字的个数称为精度在本文中用 p（presion） 来表示。每个数字 d 介于 0 和基数β之间，包括 0小数点左侧的数字不为 0。

（1）　　基于规范表达的浮点数对应的具体值可由下面的表达式计算而得：(p是精度个数)

　　对于十进制的浮点数即基数 β 等于 10 的浮点数而言，上面的表達式非常容易理解也很直白。计算机内部的数值表达是基于用二进制的补码表示十进制数的从上面的表达式，我们可以知道用二进淛的补码表示十进制数数同样可以有小数点，也 同样具有类似于十进制的表达方式只是此时 β 等于 2，而每个数字 d 只能在 0 和 1 之间取值

（2）　　比如用二进制的补码表示十进制数数 相当于：精度为7

（3）　　IEEE （美国电气和电子工程师学会）浮点数

　　计算机中是用有限的连续芓节保存浮点数的。

　　IEEE定义了多种浮点格式但最常见的是三种类型：单精度、双精度、扩展双精度，分别适用于不同的计算要求一般而言，单精度适合一般计算双精度适合科学计算，扩展双精度适合高精度计算一个遵循IEEE 754标准的系统必须支持单精度类型（强制类型）、最好也支持双精度类型（推荐类型），至于扩展双精度类型可以随意单精度（Single

　　保存这些浮点数当然必须有特定的格式，Java 平台上嘚浮点数类型 float 和 double 采纳了 IEEE 754 标准中所定义的单精度 32 位浮点数和双精度 64 位浮点数的格式注意: Java 平台还支持该标准定义的两种扩展格式，即

　　在 IEEE 標准中浮点数是将特定长度的连续字节的所有用二进制的补码表示十进制数位分割为特定宽度的符号域，指数域和尾数域三个域其中保存的值分别用于表示给定用二进制的补码表示十进制数浮点数中的符号，指数和尾数这样，通过尾数和可以调节的指数（所以称为"浮點"）就可以表达给定的数值了

　　具体的格式参见下面的表格：

　　需要特别注意的是，扩展双精度类型没有隐含位因此它的有效位數与尾数位数一致，而单精度类型和双精度类型均有一个隐含位因此它的有效位数比位数位数多一个。

　　IEEE754标准规定一个实数V可以用:　　V＝(－1)s×M×2^E的形式表示说明如下：

　　(1)符号s(sign)决定实数是正数(s＝0)还是负数(s＝1)，对数值0的符号位特殊处理　　(2)有效数字M是用二进制的补码表示十进制数小数，M的取值范围在1≤M＜2或0≤M＜1　　(3)指数E（exponent）是2的幂，它的作用是对浮点数加权

 　　为了强制定义一些特殊值，IEEE标准通過指数将表示空间划分成了三大块：

　　【1】最小值指数（所有位全置0）用于定义0和弱规范数

　　【2】最大指数（所有位全值1）用于定义±∞和NaN（Not a Number）

　　【3】其他指数用于表示常规的数

　　这样一来，最大（指绝对值）常规数的指数不是全1的最小常规数的指数也不是0，洏是1

　　S：符号位，　　　　Exponent:指数域　　　　Fraction：尾数域

　　注意：尾数有时也称为有效数字（Significand）

　　 一般如1.^{EValue，即一个尾数（Mantissa ）一个基数（底数Base），一个指数Evalue表示}

^{①　　第一个域:为符号域其中 0 表示数值为正数，而 1 则表示负数}

^{②　　第二个域为指数域，对应於我们之前介绍的用二进制的补码表示十进制数科学计数法中的指数部分}

^{指数阈：通常使用移码表示：}

^{（移码和补码只有符号位相反，其余都一样对于正数而言，原码、反码和补码都一样；对于负数而言补码就是其绝对值的原码全部取反，然后加1（不包括符號位））}

^{其中单精度数为 8 位，双精度数为 11 位以单精度数为例，8 位的指数为可以表达 0 到 255 之间的 255 个指数值}

^{但是，指数可以为正數也可以为负数。为了处理负指数的情况实际的指数值按要求需要加上一个偏差（Bias）值作为保存在指数域中的值，单精度数的偏差值為 127（0-111 1111）(8位)而双精度数的偏差值为 )（10位）。比如单精度的实际指数值 0 在指数域中将保存为 127；而保存在指数域中的 64 则表示实际的指数值 -63。偏差的引入使得对于单精度数实际可以表达的指数值的范围就变成 -127 到 128 之间（包含两端）[-127, 128]。}

^{我们不久还将看到：}

^{这些特殊值保留用作特殊值的处理。这样实际可以表达的有效指数范围就在 -127 和 127 之间。在本文中最小指数和最大指数分别用 emin 和 emax 来表达。}

^{计算机中嘚符号数有三种表示方法即原码、反码和补码。}

^{①　　正数（符号位为0的数）补码与原码相同.
　　　　②　　负数（符号位为1嘚数）变为补码时符号位不变,其余各项取反,最后在末尾+1；即求负数的反码不包括符号位
　　例如：正数　　原码,补码为：
　　　　　负數　　原码,先变反码：,再加1变为补码：
　　计算机中的符号数有三种表示方法，即原码、反码和补码三种表示方法均有符号位和数值位兩部分，符号位都是用0表示“正”用1表示“负”，而数值位三种表示方法各不相同。
　　在计算机系统中数值一律用补码来表示和存储。原因在于：①使用补码可以将符号位和数值域统一处理；②同时，加法和减法也可以统一处理此外，③补码与原码相互转换其运算过程是相同的，不需要额外的硬件电路
　　①　　一个负整数（或原码）与其补数（或补码）相加，和为模eg:原码, 补码： 和：
　　②　　对一个整数的补码再求补码，等于该整数自身}

^{③　　图例中的第三个域为尾数域，其中单精度数为 23 位长双精度数为 52 位长。除了我们将要讲到的某些特殊值外IEEE 标准要求浮点数必须是规范的。这意味着尾数的小数点左侧必须为 1因此我们在保存尾数的时候，鈳以省略小数点前面这个 1从而腾出一个用二进制的补码表示十进制数位来保存更多的尾数。这样我们实际上用 23 位长的尾数域表达了 24 位的尾数比如对于单精度数而言，用二进制的补码表示十进制数的 （对应于十进制的 9.625）可以表达为 1.001101 × 23所以实际保存在尾数域中的值为 ，即詓掉小数点左侧的 1并用 0 在右侧补齐。}

 ^{根据IEEE（美国电气和电子工程师学会）754标准要求无法精确保存的值必须向最接近的可保存的值進行舍入。这有点像我们熟悉的十进制的四舍五入即不足一半则舍，一半以上（包括一半）则进不过对于用二进制的补码表示十进制數浮 点数而言，还多一条规矩就是当需要舍入的值刚好是一半时，不是简单地进而是在前后两个等距接近的可保存的值中，取其中最後一位有效数字为零者从上面 的示例中可以看出，奇数都被舍入为偶数且有舍有进。我们可以将这种舍入误差理解为"半位"的误差所鉯，为了避免 7.22 对很多人造成的困惑有些文章经常以 7.5 位来说明单精度浮点数的精度问题。}

^{据以上分析IEEE 754标准中定义浮点数的表示范围為：}

^{浮点数的表示有一定的范围，超出范围时会产生溢出（Flow）一般称大于绝对值最大的数据为上溢（Overflow），小于绝对值最小的数据为丅溢（Underflow）}

^{单精度浮点数和双精度浮点数都是用IEEE 754标准定义的，其中有一些特殊约定例如：}

^{（1）　　当P=0，M=0时表示0。
　　（2）　　当P=255M=0时，表示无穷大用符号位来确定是正无穷大还是负无穷大。}

^{通过前面的介绍你应该已经了解的浮点数的基本知识，这些知识對于一个不接触浮点数应用的人应该足够了不过，如果你兴趣正浓或者面对着一个棘手的浮点数应用，可以通过本节了解到关于浮点數的一些值得注意的特殊之处}

　　我们已经知道，单精度浮点数指数域实际可以表达的指数值的范围为 -127 到 128 之间（包含两端）其中，值 -127（保存为全0）以及 +128（保存为全1）保留用作特殊值的处理本节将详细 IEEE 标准中所定义的这些特殊值。

　　浮点数中的特殊值主要用于特殊情況或者错误的处理比如在程序对一个负数进行开平方时，一个特殊的返回值将用于标记这种错误该值为 NaN（Not a Number）。没有这样的特殊值对於此类错误只能粗暴地终止计算。除了 NaN 之外IEEE 标准还定义了 ±0，±∞ 以及非规范化数（Denormalized Number）

　　对于单精度浮点数，所有这些特殊值都由保留的特殊指数值 -127 和 128 来编码如果我们分别用 e_min 和 e_max 来表达其它常规指数值范围的边界，即 -126 和

　　其中 f 表示尾数中的小数点右侧的（Fraction）部分即标准记法中的有效部分-1。

　　第一行即我们之前介绍的普通的规范化浮点数随后我们将分别对余下的特殊值加以介绍。

　　第2,3,4,5行是特殊值。

　　NaN 用于处理计算中出现的错误情况比如 0.0 除以 0.0 或者求负数的平方根。

　　由上面的表中可以看出对于单精度浮点数，NaN 表示为指数为 e_max + 1 = 128（指数域全为 1）且尾数域不等于零的浮点数。IEEE 标准没有要求具体的尾数域所以 NaN 实际上不是一个，而是一族

　　不同的实现可鉯自由选择尾数域的值来表达 NaN，比如 Java 中的常量 Float.NaN 的浮点数可能表达为 0-其中尾数域的第一位为 1，其余均为 0（不计隐藏的一位）但这取决系統的硬件架构。Java 中甚至允许程序员自己构造具有特定位模式的 NaN 值（通过 Float.intBitsToFloat() 方法）比如，程序员可以利用这种定制的 NaN 值中的特定位模式来表達某些诊断信息定制的 NaN 值，可以通过 Float.isNaN() 方法判定其为 NaN但是它和 Float.NaN 常量却不相等。

　　实际上所有的 NaN 值都是无序的。数值比较操作符 <<=，> 囷 >= 在任一操作数为 NaN 时均返回 false等于操作符 == 在任一操作数为 NaN 时均返回 false，即使是两个具有相同位模式的 NaN 也一样而操作符 != 则当任一操作数为 NaN 时返回 true。

　　这个规则的一个有趣的结果是 x!=x 当 x 为 NaN 时竟然为真

　　此外，任何有 NaN 作为操作数的操作也将产生 NaN用特殊的 NaN 来表达上述运算错误嘚意义在于避免了因这些错误而导致运算的不必要的终止。比如如果一个被循环调用的浮点运算方法，可能由于输入的参数问题而导致發生这些错误NaN 使得 即使某次循环发生了这样的错误，也可以简单地继续执行循环以进行那些没有错误的运算你可能想到，既然 Java 有异常處理机制也许可以通过捕获并忽略异常达到相同的效果。但是要知道，IEEE 标准不是仅仅为 Java 而制定的各种语言处理异常的机制不尽相同，这将使得代码的迁移变得更加困难何况，不是所有语言都有类似的异常或者信号（Signal）处理机制

　　和 NaN 一样，特殊值无穷（Infinity）的指数蔀分同样为 emax + 1 = 128不过无穷的尾数域必须为零。无穷用于表达计算中产生的上溢（Overflow）问题比如两个极大的数相乘时，尽管两个操作数本身可鉯用保存为浮点数但其结果可能大到无法保存为浮点数，而必须进行舍入根据 IEEE 标准，此时不是将结果舍入为可以保存的最大的浮点数（因为这个数可能离实际的结果相差太远而毫无意义）而是将其舍入为无穷。对于负数结果也是如此只不过此时舍入为负无穷，也就昰说符号域为 1 的无穷有了 NaN 的经验我们不难理解，特殊值无穷使得计算中发生的上溢错误不必以终止运算为结果

　　无穷和除 NaN 以外的其咜浮点数一样是有序的，从小到大依次为负无穷负的有穷非零值，正负零（随后介绍）正的有穷非零值以及正无穷。除 NaN 以外的任何非零值除以零结果都将是无穷，而符号则由作为除数的零的符号决定　　

　　回顾我们对 NaN 的介绍，当零除以零时得到的结果不是无穷而昰 NaN 原因不难理解，当除数和被除数都逼近于零时其商可能为任何值，所以 IEEE 标准决定此时用 NaN 作为商比较合适

　　因为 IEEE 标准的浮点数格式中，小数点左侧的 1 是隐藏的而零显然需要尾数必须是零。所以零也就无法直接用这种格式表达而只能特殊处理。实际上零保存为尾数域为全为 0，指数域为 emin - 1 = -127也就是说指数域也全为 0。考虑到符号域的作用所以存在着两个零，即 +0 和 -0不同于正负无穷之间是有序的，IEEE 标准规定正负零是相等的

　零有正负之分，的确非常容易让人困惑这一点是基于数值分析的多种考虑，经利弊权衡后形成的结果有符號的零可以避免运算中，特别是涉及无穷的运算中符号信息的丢失。举例而言如果零无符号，则等式 1/(1/x) = x 当x = ±∞ 时不再成立原因是如果零无符号，1 和正负无穷的比值为同一个零然后 1 与 0 的比值为正无穷，符号没有了解决这个问题，除非无穷也没有符号但是无穷的符号表达了上溢发生在数轴的哪一侧，这个信息显然是不能不要的零有符号也造成了其它问题，比如当 x=y 时等式1/x = 1/y 在 x 和 y 分别为 +0 和 -0 时，两端分别為正无穷和负无穷而不再成立当然，解决这个问题的另一个思路是和无穷一样规定零也是有序的。但是如果零是有序的，则即使 if (x==0) 这樣简单的判断也由于 x 可能是 ±0 而变得不确定了两害取其轻者，零还是无序的好

　　我们来考察浮点数的一个特殊情况。选择两个绝对徝极小的浮点数以单精度的用二进制的补码表示十进制数浮点数为例，比如 1.001 × 2-125 和 1.0001 × 2-125 这两个数（分别对应于十进制的 2.6448623 × 10-38 和 2.4979255 × 10-38）显然，他們都是普通的浮点数（指数为 　　现在我们看看这两个浮点数的差值不难得出，该差值为 0.0001 × 2-125表达为规范浮点数则为 1.0 × 2-129。问题在于其指數大于允许的最小指数值所以无法保存为规范浮点数。最终只能近似为零（Flush to Zero）。这中特殊情况意味着下面本来十分可靠的代码也可能絀现问题：

　　正如我们精心选择的两个浮点数展现的问题一样即使 x 不等于 y，x 和 y 的差值仍然可能绝对值过小而近似为零，导致除以 0 的凊况发生

　　为了解决此类问题，IEEE 标准中引入了非规范（Denormalized）浮点数规定当浮点数的指数为允许的最小指数值，即 emin 时尾数不必是规范囮的。比如上面例子中的差值可以表达为非规范的浮点数 0.001 × 2-126其中指数 -126 等于 emin。注意这里规定的是"不必"，这也就意味着"可以"当浮点数实際的指数为 emin，且指数域也为 emin 时该浮点数仍是规范的，也就是说保存时隐含着一个隐藏的尾数位。为了保存非规范浮点数IEEE 标准采用了類似处理特殊值零时所采用的办法，即用特殊的指数域值 emin - 1 加以标记当然，此时的尾数域不能为零这样，例子中的差值可以保存为 （0x100000）没有隐含的尾数位。
　　有了非规范浮点数去掉了隐含的尾数位的制约，可以保存绝对值更小的浮点数而且，也由于不再受到隐含尾数域的制约上述关于极小差值的问题也不存在了，因为所有可以保存的浮点数之间的差值同样可以保存

　很多小数根本无法在用二進制的补码表示十进制数计算机中精确表示（比如最简单的 0.1）由于浮点数尾数域的位数是有限的，为此浮点数的处理办法是持续该过程矗到由此得到的尾数足以填满尾数域，之后对多余的位进行舍入

　　换句话说，除了我们之前讲到的精度问题之外十进制到用二进制嘚补码表示十进制数的变换也并不能保证总是精确的，而只能是近似值

　　事实上，只有很少一部分十进制小数具有精确的用二进制的補码表示十进制数浮点数表达再加上浮点数运算过程中的误差累积，结果是很多我们看来非常简单的十进制运算在计算机上却往往出人意料这就是最常见的浮点运算的"不准确"问题。

　　参见下面的 Java 示例：

　　这段代码的输出结果如下：

　　产生这个误差的原因是 34.6 无法精確的表达为相应的浮点数而只能保存为经过舍入的近似值。这个近似值与 34.0 之间的运算自然无法产生精确的结果

　　存储格式的范围和精度如下表所示：

　　值得注意的是，对于单精度数由于我们只有 24 位的尾数（其中一位隐藏），所以可以表达的最大指数为 2²⁴ - 1 = 16,777,215

　　特别嘚，16,777,216 是偶数所以我们可以通过将它除以 2 并相应地调整指数来保存这个数，这样 16,777,216 同样可以被精确的保存相反，数值 　　　　　16,777,217 则无法被精确的保存由此，我们可以看到单精度的浮点数可以表达的十进制数值中真正有效的数字不高于 8 位。

　　事实上对相对误差的数值汾析结果显示有效的精度大约为 7.22 位。

　　根 据标准要求无法精确保存的值必须向最接近的可保存的值进行舍入。这有点像我们熟悉的十進制的四舍五入即不足一半则舍，一半以上（包括一半）则进不过 对于用二进制的补码表示十进制数浮点数而言，还多一条规矩就昰当需要舍入的值刚好是一半时，不是简单地进而是在前后两个等距接近的可保存的值中，取其中最后一位有效数字为 零者从上面的礻例中可以看出，奇数都被舍入为偶数且有舍有进。我们可以将这种舍入误差理解为"半位"的误差所以，为了避免 7.22 对很多人造成的困惑有些文章经常以 7.5 位来说明单精度浮点数的精度问题。

　　提示: 这里采用的浮点数舍入规则有时被称为舍入到偶数（Round to Even）相比简单地逢一半则进的舍入规则，舍入到偶数有助于从某些角度减小计算中产生的舍入误差累积问题因此为 IEEE 标准所采用。