浮点数为什么不精确?

最近在项目中遇到了浮点数相加结果显示异常的问题,如下图所示


demo测试结果

理论上来说,d应该是等于0的,为什么结果和预期的不一样呢,所以抽空研究了一下float

原理

浮点数的表示有一个 IEEE 的标准,它定义了两个基本的格式。

  • 用 32 比特表示单精度的浮点数,即 float 或者 float32 类型。
  • 用 64 比特表示双精度的浮点数,即 double 或者 float64 类型。
    这两种原理差不多,所以只要看单精度类型的就可以了
浮点数的标准定义

第一部分是一个符号位,用来表示是正数还是负数,用 s 来表示,所有的浮点数都是有符号的。

接下来是一个 8 个比特组成的指数位。一般用 e 来表示。8 个比特能够表示的整数空间,就是 0~255。我们在这里用 1~254 映射到 -126~127 这 254 个有正有负的数上。浮点数希望能够表示很小的数,所以指数位也会有负数。

最后,是一个 23 个比特组成的有效数位。我们用 f 来表示。综合科学计数法,我们的浮点数就可以表示成下面这样: (−1)s × 1.f × 2e

以 0.5 为例子。0.5 的符号为 s 应该是 0,f 应该是 0,而 e 应该是 -1,也就是0.5 = (−1)0 × 1.0 × 2−1 = 0.5,对应的浮点数表示,就是 32 个比特。如下图所示
0.5

s=0,e=2−1,需要注意,e 表示从 -126 到 127 个,-1 是其中的第 126 个数,这里的 e 如果用整数表示,就是 26 + 25 + 24 + 23 + 22 + 21 = 126,1.f=1.0。

为什么会不精确

10进制小数转换为2进制的方法:乘2,然后取整(整数部分最大为1,最小为0),小数部分继续乘2,取整,直到小数部分0为止。
以0.1为例,转化如下

乘2 整数部分 二进制表示
1 0.1*2=0.2 0 0
2 0.2*2=0.4 0 0
3 0.4*2=0.8 0 0
4 0.8*2=1.6 1 1
5 0.6*2=1.2 1 1
6 0.2*2=0.4 0 0

从第6行开始就会重复2~5行的运算,0.1就会变成一个无限循环的小数:0.000110011……,这里的0011无限循环。浮点数其实是用二进制的科学计数法来表示的,所以我们可以把小数点左移三位,这个数就变成了:1.10011001100110011…× 2-4

这种科学计数法对应打上面的格式为:符号位 s = 0,对应的有效位 f=1001100110011…。因为 f 最长只有 23 位,那这里“0011”无限循环,最多到 23 位就截止了。于是,f=10011001100110011001100。对应的指数为 e,代表的应该是 -4。其对应的是123,转化成二进制就是 01111011。那么最终的二进制表示为:001111011 10011001100110011001100

现在再将二进制浮点数转化为十进制,结果为0.0999999940395

知道这些有什么用?

在讨论这个之前,先看一下浮点数如何相加。首先要把两个的指数位对齐,也就是把指数位都统一成两个其中较大的那个,那么指数位较小的数,需要进行有效位右移,在右移的过程中,最右侧的有效位就被丢弃掉了。这会导致对应的指数位较小的数,在加法发生之前,就丢失精度。如果两个数的指数位差出 23 位,较小的数右移 24 位之后,所有的有效位就都丢失了。这也就意味着,虽然浮点数可以表示上到 3.40×1038,下到 1.17×10−38 这样的数值范围。但是在实际计算的时候,只要两个数,差出 224,也就是 16777216倍,那这两个数相加之后,结果完全不会变化。如下图所示

浮点数相加

这就是浮点数运算大数吃小数的问题,那么怎么解决呢,方法其实也很简单,就是记录一下误差用于下次计算
gap

举个例子,比如我们需要循环计算2000w个1.0,根据上面的可知,当计算到16777216就会出现大数吃小数的情况,那么就会造成严重的误差
2000w个1.0相加

如果我们记住每次的误差,就能解决这个问题

参考文章:深入浅出计算机组成原理

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 215,384评论 6 497
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 91,845评论 3 391
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 161,148评论 0 351
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 57,640评论 1 290
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 66,731评论 6 388
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,712评论 1 294
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,703评论 3 415
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,473评论 0 270
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,915评论 1 307
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,227评论 2 331
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,384评论 1 345
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,063评论 5 340
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,706评论 3 324
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,302评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,531评论 1 268
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,321评论 2 368
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,248评论 2 352

推荐阅读更多精彩内容