上次在校园招聘的时候,问了几个学生一个关于浮点数计算的问题,就是下面的代码为什么第一行返回false,而第二行和第三行都返回true。
Console.WriteLine("1.123f + 1.345f == 2.468f ? {0}",
1.123f + 1.345f == 2.468f);// False
Console.WriteLine("1.123f + 1.344f == 2.467f ? {0}",
1.123f + 1.344f == 2.467f);// True
Console.WriteLine("1.123 + 1.345 == 2.468 ? {0}",
1.123 + 1.345 == 2.468);// True
我们知道,integer类型占用的是4个字节,可以表示的数字范围在2**-32到2**32之间.而且每一个整数在计算机里面都有相应的内存表示方法,例如数字12345的表示方法就是:
00000000 00000000 00110000 00111001
即根据二进制到十进制的转换方法: 12345 = 213+ 212+ 25+ 24+ 23+ 20
而在计算机float型,这种表示方法最大的问题就是如何在二进制里面保存小数点的位置--因为小数点的位置不固定.于是计算机专家们就想到了用科学计数法的方式来表示浮点数,因为在科学计数法里面,小数点的位置总是固定,就是在第一个数字的后面,例如12345的科学计数法的表示为:1.2345 * 105。而0.012345的科学计数法表示方式为:1.2345 * 10-2.
使用公式来表示的话,科学计数法的公式是。
因为我们的计算机比较笨,只能处理0和1,所以在计算机里面表示浮点数的时候,上面的公式中的基数b是2,而不是10.在计算机内存当中,保存的实际是浮点数的计算公式,而不是确切的值,所以说计算机里面浮点数都是近似值,而不是确切的值.在计算机中,以float类型为例,内存中32个位所代表的内容分别是:
SEEEEEEE EMMMMMMM MMMMMMMM MMMMMMMM
其中S代表符号位,1代表这个浮点数是一个负数,而0则表示是正数。而E就是科学计数法里面的指数e,由于e既有可能是正数,也有可能是负数,因为1234.5的科学计数法是1.2345 * 104,而0.12345的科学计数法是1.2345 * 10-1。所以e的计算规则是EEEEEEE E表示的数减去27- 1(7是因为我们有8个位来表示指数),这样8个表示指数的位就可以用来表示负数和正数了。
而M就是用来计算科学计数法里面的m,计算规则是,从左往右开始,第一个M代表2-1,第二个M表示2-2,依次类推。由于m要么就是大于1的小数,要么就是小于1的小数(记住,我们计数法里面的基数是2,而不是10),如果e的值大于0,那么我们就可以加上这个1,如果e的值小于0,那么我们可以省略这个1,这是因为我们总是可以通过调整e的值来做到这一点。因此m里面小数点前面的1被省略掉了,这个省略掉的1可以根据e的值来推算出m的小数点前面到底有没有这个1。
而对于double类型来说,在计算机中各个位表示的信息如下所示:
SEEEEEEE EEEEMMMM MMMMMMMM MMMMMMMM MMMMMMMM MMMMMMMM MMMMMMMM MMMMMMMM
现在我们再回过头来看看上面的代码:
doublea = 1.345f;
doubleb = 1.123f;
doubleresult = a + b;
doubleexpected = 2.468f;
Console.WriteLine("result == expected ? {0}", result == expected);
调试这个程序,在内存中看一看a的值是3ff5851ec0000000(为了说明方便,我将文章开头的程序改头换面了一下),在计算器里面将这个十六进制的值换算成二进制的是:
11111111110101100001010001111011000000000000000000000000000000
因为计算器会将前缀零省略掉,因此上面的值实际上是:
00111111 11110101 10000101 00011110 11000000 00000000 00000000 00000000
由于e部分的值大于等于0,因为127 – (27– 1) = 1,所以我们在计算m的时候需要加上隐含的1,也就是说上面黄色背景的值实际上表示的是:
(十进制的)1.0 + 0101 10000101 00011110 11000000 00000000 00000000 00000000所表示的小数
= 1.0 + 2**-2 + 2**-4 + 2**-5 + 2**-10 + 2**-12 + 2**-16 + 2**-17 + 2**-18 + 2**-19 + 2**-21 + 2**-22
= 1.0 + 0.25
+ 0.0625
+ 0.03125
+ 0.0009765625
+ 0.000244140625
+ 0.0000152587890625
+ 0.00000762939453125
+ 0.000003814697265625
+ 0.0000019073486328125
+ 0.000000476837158203125
+ 0.0000002384185791015625
= 1.0 + 0.3450000286102294921875
= 1.3450000286102294921875
因此按照上面的算法,b的值实际上应该是1.1230000257492065,因此a + b实际上是2.468000054359436,而上面的程序expected的实际值应该是2.4679999351501465,这就是为什么上面1.345f + 1.123f != 2.468f的原因。
而1.344f + 1.123f == 2.467f,纯粹是巧合,因为前两者相加的值恰好等于后者在计算机里面的表现形式。
最后洋洋洒洒一大篇的结论就是,浮点类型不能用==号来判断,因为浮点数是一个近似值,只能通过两者相减小于一个可以接受的误差来判断。