C语言浮点数运算Word文档格式.docx
- 文档编号:17371076
- 上传时间:2022-12-01
- 格式:DOCX
- 页数:11
- 大小:27.90KB
C语言浮点数运算Word文档格式.docx
《C语言浮点数运算Word文档格式.docx》由会员分享,可在线阅读,更多相关《C语言浮点数运算Word文档格式.docx(11页珍藏版)》请在冰豆网上搜索。
IEEE754是最广泛使用的二进制浮点数算术标准,被许多CPU与浮点运算器所采用。
IEEE754规定了多种表示浮点数值的方式,在本文档里只介绍32bits的float浮点类型。
它被分为3个部分,分别是符号位S(signbit)、指数偏差E(exponentbias)和小数部分F(fraction)。
2012-05-2219:
21上传
下载附件(7.41KB)
其中S位占1bit,为bit31。
S位为0代表浮点数是正数,S位为1代表浮点数是负数,比如说0x449A522C的S位为0,表示这是一个正数,0x849A522C的S位为1,表示这是一个负数。
E位占8bits,为bit23~bit30。
E位代表2的N次方,但需要减去127,比如说E位为87,那么E位的值为2(87-127)=9.094947017729282379150390625e-13。
F位占23bits,为bit0~bit22。
F位是小数点后面的位数,其中bit22是2-1=0.5,bit21是2-2=0.25,以此类推,bit0为2-23=0.00000011920928955078125。
但F位里隐藏了一个1,也就是说F位所表示的值是1+(F位bit22~bit0所表示的数值),比如说F位是0b10100000000000000000001,只有bit22、bit20和bit0为1,那么F位的值为1+(2-1+2-3+2-23),为1.62500011920928955078125。
综上所述,从二进制数换算到浮点数的公式为:
(-1)S×
2E-127×
(1+F)。
但还有几个特殊的情形:
若E位为0并且F位也为0时表示浮点数0,此时浮点数受S位影响,表现出+0和-0两种0,但数值是相等的。
比如二进制数0x00000000表示+0,二进制数0x80000000表示-0。
若E位为0并且F位不为0时浮点数为(-1)S×
2-126×
F,注意,E位的指数是-126,而不是0-127=-127,而且F位是0.xx格式而不是1.xx格式,比如0x00000001的浮点数为2-126×
2-23=1.4012984643248170709237295832899e-45,而不是20-121×
(1+2-23)。
一旦E为不为0,从0变为1,不是增加2倍的关系,因为公式改变了。
若E位为255并且F位不为0时表示非数值,也就是说是非法数,例如0x7F800001。
若E位为255并且F位为0时表示无穷大的数,此时浮点数受S位影响,例如0x7F800000表示正无穷大,0xFF800000表示负无穷大。
当我们使用1个数除以0时,结果将被记作0x7F800000。
浮点型在多个处理器间通信时,传递的数值是它的二进制数,比如说1234.5678这个浮点数的二进制数是0x449A522B,如果使用串口发送的话,就会发现串口里发送的是0x44、0x9A、0x52和0x2B这4个数(发送的顺序也可能是逆序,这与约定的字节序有关,与浮点格式无关),接收端接收到这4个数字后再组合成0x449A522B,按照IEEE754的定义被解析成1234.5678,这样就实现浮点数通信了。
如果两个处理器所使用的浮点数规则不同,则无法传递浮点数。
浮点数的换算
下面来看看浮点数与二进制数如何转换。
例1,二进制数换算成浮点数:
假如在内存中有一个二进制数为0x449A522C,先将十六进制转换成二进制,如下:
0100
1001
1010
0101
0010
1100
按照SEF的格式分段,如下:
0
10001001
00110100101001000101100
这个数值不是特殊的情形,可以按照公式(-1)S×
(1+F)转换。
S位的值为(-1)0=1,E位的值为2137-127=1024。
F位的值为1+2-3+2-4+2-6+2-9+2-11+2-14+2-18+2-20+2-21=1.205632686614990234375。
最终结果为1×
1024×
1.205632686614990234375=1234.56787109375。
其中F位比较长,使用二进制方式转换比较麻烦,也可以先转换成十六进制再计算,转换为十六进制如下:
0011
1000
0x3
0x4
0xA
0x5
0x8
F位为23bits,需要在最后补一个0凑成24bits,共6个十六进制数。
F位的值为1+3×
16-1+4×
16-2+10×
16-3+4×
16-4+5×
16-5+8×
16-6=1.205632686614990234375,与上面使用二进制方法得到的结果相同。
例2,浮点数换算成二进制数:
下面我们将-987.654e30换算成二进制数。
我们先不看符号位,将987.654e30归一化为整数部分为1的形式,也就是写作987.654e30=2E-127×
(1+F)的标准形式,其中E=log(987.654e30)/log2+127=109.6+127,取E位的整数值为109+127=236,再求F=987.654e30/2236-127-1=0.52172193,这个小数位数保留8位就够了,已经超出了7位的精度。
然后我们求小数部分的二进制数,这个转换就没啥好说的了,依次减去2的幂,从2-1一直到2-23,够减的位置1,不够减的位置0,例如,2-1为0.5,0.52172193-0.5=0.02172193,F位的bit22置1,2-2为0.25,0.02172193不够减,F位的bit21置0,2-3为0.125,0.02172193不够减,F位的bit20置0,2-4为0.0625,0.02172193不够减,F位的bit19置0……,一直算到F位的bit0,这样就得到F位的数值。
如果觉得使用二进制方式转换太麻烦的话也可以使用十六进制进行转换。
16-1为0.0625,0.52172193/0.0625=8.3,说明够减8个,记做0x8,0.52172193-0.0625×
8=0.02172193,16-2为0.00390625,0.02172193/0.00390625=5.6,说明够减5个,加上刚才的0x8记做0x85,以此类推:
16的-N次幂
被减数
十六进制数
减后的数
1
0.0625
0.52172193
0x8
0.02172193
2
0.00390625
0.02172193
0x85
0.00219068
3
0.000244140625
0.00219068
0x858
0.000237555
4
0.0000152587890625
0.000237555
0x858F
0.0000086731640625
5
0.00000095367431640625
0.0000086731640625
0x858F9
0.00000009009521484375
6
0.000000059604644775390625
0.00000009009521484375
0x858F91
一直凑够23bits,也就是6个十六进制,得到0x858F91,换算成二进制如下所示:
1000
1111
0001
由于只有23bits有效,因此需要去掉最后一个bit,二进制本着0舍1入的原则,变成
001
最后需要再补上前面的S位和E位。
由于是负数,S位为1。
E位为236,二进制形式为11101100,将S、E、F位组合在一起就形成了:
11101100
从左边最高位开始,4个一组合并成十六进制:
0110
1100
0111
1001
换算成十六进制为:
0xF
0x6
0x2
0xC
0x7
0x9
综上所述,-987.654e30换算成二进制数为0xF642C7C9。
浮点数的精度
在前面的讲解中可以看到1.xx这个数量级的最小数是2-23,对应的十进制数值为1.00000011920928955078125,可以精确表示到小数点后23位,但有些C语言书上却说float型的有效位只有6~7位,这是为什么?
这是因为二进制小数与十进制小数没有完全一一对应的关系,二进制小数对于十进制小数来说相当于是离散的而不是连续的,我们来看看下面这些数字:
二进制小数
十进制小数
2-23
1.00000011920928955078125
2-22
1.0000002384185791015625
2-21
1.000000476837158203125
2-20
1.00000095367431640625
2-19
1.0000019073486328125
2-18
1.000003814697265625
不看S位和E位,只看F位,上表列出了1.xx这个数量级的6个最小幂的二进制小数,对应的十进制在上表的右边,可以看到使用二进制所能表示的最小小数是1.00000011920928955078125,接下来是1.0000002384185791015625,这两个数之间是有间隔的,如果想用二进制小数来表示8位有效数(只算小数部分,小数点前面的1是隐藏的默认值)1.00000002、1.00000003、1.00000004...这些数是无法办到的,而7位有效数1.0000001可以用2-23来表示,1.0000002可以用2-22来表示,1.0000003可以用2-23+2-22来表示。
从这个角度来看,float型所能精确表示的位数只有7位,7位之后的数虽然也是精确表示的,但却无法表示任意一个想表示的数值。
但还是有一些例外的,比如说7位有效数1.0000006这个数就无法使用F位表示,二进制小数对于十进制小数来说相当于是离散的,刚好凑不出1.0000006这个数,从这点来看float型所能精确表示的位数只有6位。
至于5位有效值的任何数都是可以使用F位相加组合出来的,即便是乘以E位的指数后也是可以准确表示出来的。
因此float型的有效位数是6~7位,但这个说法应该不是非常准确,准确来说应该是6位,C语言的头文件中规定也是6位。
对于一个很大的数,比如说1234567890,它是F位乘上E位的系数被放大了的,但它的有效位仍然是F位所能表示的6位有效数字。
1234567890对应的二进制数是0x4E932C06,其中F位的数值为1.1497809886932373046875,E位的数值为230=1073741824,1073741824×
1.1497809886932373046875=1234567936,对比1234567890,也只有高7位是有效位,后3位是无效的。
int型定点数可以准确的表示1234567890,而float浮点数则只能近似的表示1234567890,精度问题决定了float型根本无法取代int型。
浮点数的比较
从上面的讨论可以看出,float型的有效位数是6位,那么我们在用float型运算时就要注意了,来看下面这段程序:
#include<
stdio.h>
intmain(void)
{
floata=9.87654321;
floatb=9.87654322;
if(a>
b)
{
printf("
a>
b\n"
);
}
elseif(a==b)
a==b\n"
else
a<
return0;
}
按照我们平时的经验来说这段程序应该走a<
b的分支,但程序运行的结果却走了a==b的分支,原因就是float型的精度问题,float型无法区分出小数点后的第8位数,在内存中,a和b的二进制数都是0x411E0652,因此就走了a==b的分支。
某些编译器在编译时会发现a和b的值超出了浮点数的精度,会产生一个告警,提示数据超过精度,但有些编译器则不会做任何提示。
最可怕的是有一类编译器,调试窗口里显示的长度超出float型的精度,比如说a的值显示为9.87654321,b的值显示为9.87654322,但在运行时,硬件可不管这套,硬件认为这2个数都是0x411E0652,因此实际运行结果是a==b的分支。
以前就遇到过一个同学在QQ群里问一个类似的问题,在调试窗口里明明写着a是9.87654321,小于b的9.87654322,但运行结果却是a==b。
当时我给他说了半天也没能让他明白这个问题,希望他能有机会看到这个文档,希望他这次能够明白^_^。
由于精度这个问题的限制,我们在浮点数比较时就需要加一个可接受的精度条件来做判决,比如说上面的这个问题,如果我们认为精度在0.00001就足够了,那么a-b之差的绝对值只要小于0.00001,我们就认为a和b的值是相等的,大于0.00001则认为不等,还要考虑到a-b正负等情况,因此可以将上面的程序改写为:
if(a-b<
-0.00001)
elseif(a-b>
0.00001)
例子中a和b之差的绝对值小于0.00001,因此认为a和b是相等的,运行程序,也正确的打印了a==b。
也许你会觉得费了这么大的劲,最后2个程序运行的结果还是一样的,这不是画蛇添足么?
硬件已经自动考虑到精度问题了,为什么我们还要去再做一个精度的限定?
这是因为我们在应用中的精度往往要低于硬件的6位精度。
比如说我们使用2个AD采集2V的电压,并对这2个电压值做比较,一般要求精确到0.1V即可。
实际情况下AD采集出来的数值都会围绕真实值有上下的波动,比如说AD的精度是0.001V,我们采集出的电压数值就可能是2.003V、2.001V、1.999V等围绕2V波动的数值,如果我们在程序里不对精度加以限制就对这些数值做比较,就会发现几乎每次采集到的电压值都是不同的。
在这种情况下我们就需要将精度设定为0.1V,将上面例子中的0.00001改为0.1,就可以正确的判断出每次采集的电压值都是相同的了。
在实际使用AD采样时可能并不需要使用浮点数,我一般都是使用定点数来代替浮点数进行处理的,请参考另一篇案例《C语言使用定点数代替浮点数计算》。
下面我们再看一个例子:
floata=987654321;
floatb=987654322;
这个例子中的两个数都是很大的数,已经远远超过了0.00001的精度,运行结果是不是应该是a<
b?
但程序运行的结果依然是a==b。
这是因为这个例子里的a和b并不是1.xx的数量级,我们将a和b进行归一化,都除以900000000就会发现a=1.09739369,b=1.097393691,只是在第9位才出现不同,因此在0.00001这个精度下,这2个数还是相等的。
换个角度来看,a和b虽然是很大的数了,但F位仅能表示23bits,有效值仅有6位,a和b的大是因为E位的指数放大F位表现出来的,但有效值依然是6位。
在内存中a和b的二进制数都是0x4E6B79A3。
其中E位为156,2156-127=536870912,F位为0xD6F346,F位1.xx数量级的1.83964955806732177734375被E位放大了536870912倍,E位如果算7位有效精度的话能精确到0.0000001,乘以536870912已经被放大了53倍多,这就说明a和b的个位与十位已经不是有效位数了,所以程序运行的结果表现出a==b也是正常的。
由此可见,设定一个合理的精度是需要结合浮点数的数量级的,这看起来似乎比较难,毕竟在程序运行时十分精确的跟踪每个浮点数的数量级并不容易实现,但实际应用中需要比较的浮点数往往都会有其物理含义,例如上面电压的例子,因此,根据浮点数的物理含义,浮点数的精度还是比较好确定的。
当然在一些复杂的数值运算过程中可能会存在非常复杂的情况,这时浮点数的精度问题就比较棘手了,所幸我所做的都是比较简单的东西,那些复杂的情况就不讨论了,我也没能力讨论^_^。
上面所说的都是同等数量级下的浮点数进行比较的情况,不同数量级的浮点数比较则没有这个限制,比如说1.23456789与12.23456789的比较,在E位已经足以区分大小了,因此F位的精度就没有必要再比较了。
浮点数的加减
二进制小数与十进制小数之间不存在一一对应的关系,因此某些十进制很整的加减法小数运算由二进制小数来实现就表现出了不整的情况,来看下面的例子:
floata=10.2;
floatb=9;
floatc;
c=a-b;
printf("
%f\n"
c);
如果用十进制计算的话变量c应该为1.2,在VisualC++2010环境下实验输出为1.200000,但实际上c变量的值是1.1999998,只不过是在输出时被四舍五入为1.200000罢了。
在内存中c变量的二进制数是0x3F999998,它对应的浮点数是0.19999980926513671875。
如果我们将printf函数%f的格式改为%.7f格式,就会看到c变量输出的值是1.1999998。
两个数量级相差很大的数做加减运算时,数值小的浮点数会受精度限制而被忽略,看下面的例子:
intmain(void
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 语言 浮点 运算