Chinaunix首页 | 论坛 | 博客
  • 博客访问: 617407
  • 博文数量: 172
  • 博客积分: 10010
  • 博客等级: 上将
  • 技术积分: 1252
  • 用 户 组: 普通用户
  • 注册时间: 2009-06-29 22:26
文章分类

全部博文(172)

文章存档

2011年(6)

2010年(7)

2009年(159)

我的朋友

分类: C/C++

2010-06-16 22:24:48

一、浮点数

1、定义

  浮点数是属于有理数中某特定子集的数的数字表示,在计算机中用以近似表示任意某个实数。具体的说,这个实数由一个整数或定点数(即尾数)乘以某个基数(计算机中通常是2)的整数次幂得到,这种表示方法类似于基数为10的科学记数法。

2、浮点计算

  是指浮点数参与的运算,这种运算通常伴随着因为无法精确表示而进行的近似或舍入。
  一个浮点数a由两个数m和e来表示:a = m × b^e。在任意一个这样的系统中,我们选择一个基数b(记数系统的基)和精度p(即使用多少位来存储)。m(即尾数)是形如±d.ddd...ddd的p位数(每一位是一个介于0到b-1之间的整数,包括0和b-1)。如果m的第一位是非0整数,m称作规格化的。有一些描述使用一个单独的符号位(s 代表+或者-)来表示正负,这样m必须是正的。e是指数。

3、结构

  由此可以看出,在计算机中表示一个浮点数,其结构如下:
  尾数部分(定点小数 阶码部分(定点整数)
  
数符± 尾数m 阶符± 阶码e
  这种设计可以在某个固定长度的存储空间内表示定点数无法表示的更大范围的数。

二、实例

1、题目

  例如,一个指数范围为±4的4位十进制浮点数可以用来表示43210,4.321或0.0004321,但是没有足够的精度来表示432.123和43212.3(必须近似为432.1和43210)。当然,实际使用的位数通常远大于4。

2、特别数值

  此外,浮点数表示法通常还包括一些特别的数值:+∞和−∞(正负无穷大)以及NaN('Not a Number')。无穷大用于数太大而无法表示的时候,NaN则指示非法操作或者无法定义的结果。

3、二进制表示

  众所周知,计算机中的所有数据都是以二进制表示的,浮点数也不例外。然而浮点数的二进制表示法却不像定点数那么简单了。

4、浮点数的概念

  先澄清一个概念,浮点数并不一定等于小数,定点数也并不一定就是整数。所谓浮点数就是小数点在逻辑上是不固定的,而定点数只能表示小数点固定的数值,具用浮点数或定点数表示某哪一种数要看用户赋予了这个数的意义是什么。
  C++中的浮点数有6种,分别是:
  float:单精度,32位
  unsigned float:单精度无符号,32位
  double:双精度,64位
  unsigned double:双精度无符号,64位
  long double:高双精度,80位
  unsigned long double:高双精度无符号,80位(嚯,应该是C++中最长的内置类型了吧!)
  然而不同的编译器对它们的支持也略有不同,据我所知,很多编译器都没有按照IEEE规定的标准80位支持后两种浮点数的,大多数编译器将它们视为double,或许还有极个别的编译器将它们视为128位?!对于128位的long double我也仅是听说过,没有求证,哪位高人知道这一细节烦劳告知。
  下面我仅以float(带符号,单精度,32位)类型的浮点数说明C++中的浮点数是如何在内存中表示的。先讲一下基础知识,纯小数的二进制表示。(纯小数就是没有整数部分的小数,讲给小学没好好学的人)
  纯小数要想用二进制表示,必须先进行规格化,即化为 1.xxxxx * ( 2 ^ n ) 的形式(“^”代表乘方,2 ^ n表示2的n次方)。对于一个纯小数D,求n的公式如下:
  n = 1 + log2(D); // 纯小数求得的n必为负数
  再用 D / ( 2 ^ n ) 就可以得到规格化后的小数了。接下来就是十进制到二进制的转化问题,为了更好的理解,先来看一下10进制的纯小数是怎么表示的,假设有纯小数D,它小数点后的每一位数字按顺序形成一个数列:
  {k1, k2, k3, ... , kn}
  那么D又可以这样表示:
  D = k1 / (10 ^ 1 ) + k2 / (10 ^ 2 ) + k3 / (10 ^ 3 ) + ... + kn / (10 ^ n )
  推广到二进制中,纯小数的表示法即为:
  D = b1 / (2 ^ 1 ) + b2 / (2 ^ 2 ) + b3 / (2 ^ 3 ) + ... + bn / (2 ^ n )
  现在问题就是怎样求得b1, b2, b3,……,bn。算法描述起来比较复杂,还是用数字来说话吧。声明一下,1 / ( 2 ^ n )这个数比较特殊,我称之为位阶值。

5、例二

  例如0.456,第1位,0.456小于位阶值0.5故为0;第2位,0.456大于位阶值0.25,该位为1,并将0.456减去0.25得0.206进下一位;第3位,0.206大于位阶值0.125,该位为1,并将0.206减去0.125得0.081进下一位;第4位,0.081大于0.0625,为1,并将0.081减去0.0625得0.0185进下一位;第5位0.0185小于0.03125……
  最后把计算得到的足够多的1和0按位顺序组合起来,就得到了一个比较精确的用二进制表示的纯小数了,同时精度问题也就由此产生,许多数都是无法在有限的n内完全精确的表示出来的,我们只能利用更大的n值来更精确的表示这个数,这就是为什么在许多领域,程序员都更喜欢用double而不是float。
  float的内存结构,我用一个带位域的结构体描述如下:
  struct MYFLOAT
  {
  bool bSign : 1; // 符号,表示正负,1位
  char cExponent : 8; // 指数,8位
  unsigned long ulMantissa : 23; // 尾数,23位
  };
  符号就不用多说了,1表示负,0表示正
  指数是以2为底的,范围是 -128 到 127,实际数据中的指数是原始指数加上127得到的,如果超过了127,则从-128开始计,其行为和X86架构的CPU处理加减法的溢出是一样的。
  比如:127 + 2 = -127;-127 - 2 = 127
  尾数都省去了第1位的1,所以在还原时要先在第一位加上1。它可能包含整数和纯小数两部分,也可能只包含其中一部分,视数字大小而定。对于带有整数部分的浮点数,其整数的表示法有两种,当整数大于十进制的16777215时使用的是科学计数法,如果小于或等于则直接采用一般的二进制表示法。科学计数法和小数的表示法是一样的。
  小数部分则是直接使用科学计数法,但形式不是X * ( 10 ^ n ),而是X * ( 2 ^ n )。拆开来看。
  0 00000000 0000000000000000000000
  符号位 指数位 尾数位
  ----------------------------------------------------------------------------------

三、浮点前导数字的分布

1、简介

  作者:concreteHAM
  
 什么是浮点数,不用我多说,这里我们要讨论的是规格化的任意进制浮点数的前导数字的概率分布。
  在《计算机程序设计艺术》第二卷中做了非常深入的讨论,这里我从中精炼出要点。

2、实例

  例如:
  2.345 E 67
  这是一个十进制规格化浮点数,前导数字就是2 。
  就只有一个“随机”的浮点数而言,讨论其分布式没有意义的,我们要讨论的是充分多个“随机”数进行的一系列运算后产生的浮点结果的前导数字分布。
  假设现在有一巨大的浮点数集,依此对数集中每个浮点数都乘以2,其中有一个十进制浮点数F,它的前导数字是1,那么它底数可能的值范围就是1.000…~1.999…,乘以一个数字2,那么它的底数就变成2.000…~3.999…,很明显乘以2以前前导数字是1的浮点个数与现在前导数字是2、3的浮点个数相同。以此我们接下来分析。
  对于一个b进制的浮点数,它的前导数字x范围就是0 < x < b,设f(x)是上述数集的前导数字的概率密度函数(注:是密度函数),那么它在前导数字u和v之间(0
  ∫[u,v]f(x)dx (1)
  由前面所述的,对于一个充分小增量Δx,f(x)必须满足这样一个公式:
  f(1)Δx = x*f(x)Δx (2)
  因为:
  f(1)Δx是f(1)微分段内的概率,根据前面所述,f(1)Δx概率等于f(1*x)*(x*Δx)
  很明显:
  f(x) = f(1)/x (3)
  两边在[1,b]之间进行积分,等号左边必定为1,右边等于f(1)ln(b):
  1 = f(1)ln(b) (4)
  得:f(1) = 1/ln(b) 带入(3)中:
  f(x) = 1/(x*ln(b))
  那么利用(1)式得:
  ∫[u,v]1/(x*ln(b))dx
  = ln(v/u) / ln(b) (5)
  这就是求前导数字的概率分布函数。
  例如b = 10进制时,前导数字为1的概率就是:
  = ln((1+1)/1) / ln(10)
  ≈ 0.301
  前导数字为9的概率就是:
  = ln((9+1)/9) / ln(10)
  ≈0.0458
  以下是一个测试程序(Mathematica软件):
  T[n_,b_]:=Block[{res={},ran,i,a},
  For[i=1,i
  res=Append[res,0]
  ];
  For[i=0,i
  ran=Random[]*Random[]*Random[]; 充分打乱浮点数
  ran=Log[b,ran];
  a=Floor[b^(ran-Floor[ran])]; 取出前导数字
  res[[a]]++ 对前导数字个数统计
  ];
  Return[res]
  ]
  执行T[100000,10],以10进制测试100000个浮点数,得到一个分布:
  {30149, 18821, 13317, 9674, 7688, 6256, 5306, 4655, 4134}
  和理论值相当接近。
阅读(3099) | 评论(0) | 转发(0) |
给主人留下些什么吧!~~