二进制浮点数表示

分享这篇笔记的目的是想驱动自己多动手,多动脑,开放自己,认识别人。如有问题,请指正。

从整数到小数

先回忆下如何用二进制表示十进制整数。例如,5 = 1×2^2 +0×2^1 +1×2^0 = (101)B
以此类推,
5.3 = 1×2^2+0×2^1+1×2^0+1×2 ^{-1}+1×2^{-2} = (101.11)B
这样的表示方法,会被我们理所应当地想到。但是这样表示小数存在不足

存在的问题

  1. 能精确表示的小数有限,只能表示 被x×2^y表示的数字。如1/3、 1/5等无法被上述方式表示。
  2. 在定长位数下(float是32位),不能同时表示极大数和及其接近0的数。(ps,如果小数点不定位置,那么要描述小数点位置又需要额外的位数,更加浪费空间)。

IEEE标准表示浮点数

要解决的问题:在给定区间内,整数有限个,但实数无限个
基本矛盾:有限个二进制位去表示无限多的实数
表示:使用类似科学计数法方式表示。
本文只谈float32类型情况下。

表示公式

V=(-1)^S+F×2^E
这样的表示小数点是浮动的,所以称为浮点数。(相对的,小数点确定的叫做定点数,整数也是一种定点数)
S:标志位(1位),大家都懂
E:阶码(8位),作用是对浮点数加权
F:尾数(23位),这个数有些特别,它的取值范围[1,2)或[0,1),俩个取值范围跟E的取值相关,后文详细解释。

来源 深入理解计算机系统v3

阶码E 和 尾数F 的取值

阶码状态 尾数状态 值的表示 备注
E∈(0,255) 任意 1+\frac{F}{2^{-23}}×2^{E-127} 规格化
E=0 任意 \frac{F}{2^{-23}}×2^{E-126} 非规格化
E=255 0 无穷大、小
E=255 ≠0 NaN

1.规格化情况

阶码E的二进制位既不全是0,也不全是1。
因为E∈(0,255)的,E-127 刚好平分E的值域。这样,E既可以表示正指数,也可表示副指数。
尾数F的值域[1,2)

2. 非规格化情况

阶码E的二进制位全0
尾数F的值域[0,1)

阶码E和尾数F的含义

这样的表示公式是非常反直觉的,为什么要将非规格化情况单独提出,这驱使我探究其中含义。

阶码E

先定义阶数组。根据E取值不同,2^{E-127}的值可以组成的数组是阶数组,阶数组可见下表

阶数组元素 元素值
A[0] 2^{-126}
A[1] 2^{-126}
A[2] 2^{-125}
A[3] 2^{-124}
…… ……
A[127] 2^{0}
A[128] 2^{1}
…… ……
A[253] 2^{126}
A[254] 2^{127}
A[255]

阶码E就是阶数组的序号(数组下标)
阶数组元素算出后,每个阶的值域也我们自然能算出。

阶数组元素 元素值 表示范围 精度
A[0] 2^{-126} [0,2^{-126}) 2^{-149}
A[1] 2^{-126} [2^{-126},2^{-125}) 2^{-149}
A[2] 2^{-125} [2^{-125},2^{-124}) 2^{-148}
A[3] 2^{-124} [2^{-124},2^{-123}) 2^{-147}
…… …… …… ……
A[127] 2^{0} [2^{0},2^{1}) 2^{-23}
A[128] 2^{1} [2^{1},2^{2}) 2^{-22}
…… …… …… ……
A[150] 2^{23} [2^{23},2^{24}) 2^{0}=1
…… …… …… ……
A[254] 2^{127} [2^{127},2^{128}) 2^{104}
A[255] [0,2^{-126}) -

这样的表示,就将整个实数域都囊括了。非规格化情况(A[0]阶),巧妙地填充了0附近的空白。

尾数F

尾数F理解成一个分数的分子,分母是2^23,分数的“1”是E阶的值域。打个比方,一条线段的长度为A[E],平均分成2^23段,取其前F段。

计算出浮点数的阶号E和尾数F

计算F时,采用二分法。
下面以 N=25f 进行举例

1.找到N的阶码

2^k <= N < 2^{k+1}
2^4 <= N < 2^{5}
阶码E= k+127 = 131=1000 0011

2.二分法找尾数F

判断 Y or N F从右向左的二进制占位
IF 25 >=24 Y 1
IF 25 >=28 N 10
IF 25 >=26 N 100
IF 25 >=25 Y 1001

解释一下,阶区间[24,25),中位数的十进制表示是24,25>24,判断为真,F的最高位填1。
剩余区间中位数28,25<28,判断为假,F的下一位填0。依次类推,最终F为1001时。
所以最终25f的二进制表示为:
高16位:0 1000 0011 1001 000
低16位:0000 0000 0000 0000
十六进制表示:41 C8 00 00

浮点数表示带来的精度问题

  1. python中: 0.3+0.6=0.899999

  2. 32位浮点数,超过2的24次方之后,精度大于1,运算 准确性不如整数。

    2,000,000,000f+1f = 2,000,000,000f

解决方法:Kahan Summation算法

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

推荐阅读更多精彩内容