浮点数指数域(阶码)的存储格式

浮点数公式:V=(-1)M2E 

分2种情况,非规格(-1<V<1这个区间时)和规格(V<=-1或V>=1的区间时)。

1、当表示非规格化小数时,阶码全为0。

2、当为规格化小数时,根据阶码e的计算公式:

E=e - Bias
Bias=2^k-1 - 1
e=E + Bias
e=E + 2^k-1 - 1

举例:(5.25)10 =(101.01)=1.0101 * 22

根据规格数的定义:M=1+f,f=M-1,故f=0.0101

E=2,e=2+28-1 -1 = (129)10 =(10000001)2 

当然,以上是为了准确说明和引导理解。实际中,因为单精度浮点数的价码是约定好的8位,共bias也就是127,换成2进制是01111111,我们可以直接拿2的指数的二进制数和0111111相加就可以得出价码的二进制形式了。

重点:有的人说:价码就是2的指数的移码,经过验证这是错误的。移码的定义:

按照上述定义:e=27+E,而根据上面的公式可知:e=2+ E - 1

所以:正确的求法是:知道原码后,将原码从右往左数,第一个1为分界线,右边(含第一个1)求反,左边最高位求反,其余位不变。

这个方法属于二进制的速算技巧了,可以证明出来的,有机会补充上。

举例:(5.25)2 = (101.01)2 = (1.0101 * 22 ) ,由公式得知:E=2=e-127,e=2+127,2的二进制是10,根据上述速算法:e=(10000001),f=01010000000000000000000

故5.25的浮点存储为:01000000101010000000000000000000

编程检查一下:

#include <cstdio>
int main(){
    float a=5.25;
    unsigned int b=*(int *)&a;  //把b转为无符号整型,取二进制位
    for(int i=0;i<32;i++){
        printf("%d",(b>>(31-i))&1?1:0);
    }
    printf("\nb=%d\n",b);  //看看把浮点格式转为int格式后,数字变成啥,for fun
    return 0;
}

输出:

vagrant@ubuntu-bionic:~$ ./t
01000000101010000000000000000000   //哈哈,激动,和我刚才推测的一样
b=1084751872   //看到没,float的5.25和unsigned int的1084751872一样的格式。

例2:(0.2)10别小看这个数,在二进制中它是一个无限循环小数(0.001100110011...)2,0011是循环节。这样的数如何解析呢?

它仍然是规格数,不要以为0和1之间,就是非规格数,非规格数用于表示0和极小的数,0.2显然不是,所以它还是规格数。

首先求f和e:把(0.00110011...)2左移3位,让它的小数点左边为1,这是规格数的要求,这样变成这样(1.1001100110011...)2*2-3,变大了23,自然要在缩小23,才保证原数不变。

根据这个:(1.1001100110011...)2*2-3  括号内对应M,后面的指数-3对应E,则f=M-1=(0.100110011...),e=E+127=124=(01111100)2,这里在重点说一下,e的快捷求法,-3的二进制为:3的补码(按位求反+1)得(11111101)2,然后,从右边第一个1为分界线,左边最高位求反,其余不变,右边(含第一个1)全部求反,得(01111100)2 ,这个数就是指数域了。

然后小数域:有f=(0.100110011...)2得出24位(多算1位,考虑舍入得问题),是:(0.100110011001100110011001)2,看到最后第24位是1,故向前进1位,得到(0.10011001100110011001101)2,这个数就是小数域的值了。

总上所述:(0.2)10 的二进制存储为:(00111110010011001100110011001101)2

把上面的验证程序改个数,验证一下输出:

vagrant@ubuntu-bionic:~$ ./t
00111110010011001100110011001101   //哈哈,激动,又对了
b=1045220557

PS,还有一个验证程序的写法,就是用printf函数也可以:不过这个输出16进制,需要自己转一下二进制:

#include <cstdio>
using namespace std;
int main(void)
{
float f = 0.2;
int *df = (int *)&f;
printf("%#x\n",*df);
return 0;
}

输出:

vagrant@ubuntu-bionic:~$ ./t1
0x3e4ccccd

 

有兴趣可以把a变成任一浮点数,先自己推一下,看答案。

另外,欢迎有人来查错,但不要再无脑的抄别人了,这样,才能共同进步。

posted @ 2019-03-07 09:10  立体风  阅读(1555)  评论(0编辑  收藏  举报