浮点数的精度问题

对浮点数的精度产生了疑问:

Q: 基于IEEE 754标准的double转化为int类型时,数据是如何处理的?请以一种编译器举例说明 

  

从网上找了些资料，又参考了韩卫平的专栏http://blog.csdn.net/akirya/archive/2009/07/19/4360923.aspx,终于理解了一些东西,感谢akirya的帮助!

 

/***********************************************************************
    Copyright (c) 2009,liushac
    All rights no reserved.

    ID:   newfoti.c

    功能: newftoi函数(把单精度浮点数转换成相应的整数)

***********************************************************************/
#include <stdio.h>

void newftoi(float s,int d[],int m);
int power(int base,int n);
int main()
{
    int i=0,arrayindex=0,line=0,number=0;

    int destination[2];
    destination[0]=0;
    destination[1]=0;
    float source;;

    /* ----------test1-------------*/

    printf("------test1:-------/n");
    i++;
    line++;                                             /* 行号从1开始,这是显示格式 */
    while (i<32)
    {   /* 测试 */
        (i%2)?(number=(power(-2,i))+1):(number=(power(-2,i))-1); /* 2的正整数次幂 */
        (float) source=(number);

        newftoi(source,destination,arrayindex);
   
        printf("%d/tn=%d/t%.6f=%c%d/n",line,number,source,destination[arrayindex],destination[arrayindex+1]);
        i++;
        line++;
    }

    /*-----------test2---------------*/

    printf("/n------test2:-------/n");
    i=0;
    line=0;
    line+=1;                                      /* 行从1开始 */
    number=16777212;
    while (i<10)
    {
        (i%2)? number:number*=-1;
        (float) source=(number);
        newftoi(source,destination,arrayindex);
   
        printf("%d/tn=%d/t%.6f=%c%d/n",line,number,source,destination[arrayindex],destination[arrayindex+1]);
        i++;
        (i%2)?(number-=1):(number+=1);;
        line++;
    }

    return 0;
}

/* newftoi函数: 32位float转化为int,IEEE 754标准;
转化后的整数放到入口参数中的数组d[]中,符号放入d[m],整数放入d[m+1];
注意float=0.0时,转化为+2;其余整数部分能精确转换[-16777216 ~ +16777216] */
void newftoi(float s,int d[],int m)
{
   int n=0;
   int fs=*(int*)(&s);
   if (fs>>31)
   {   /* float符号位为负 */
       d[m++]='-';
   }
   else
   {   /* float符号位为正 */
       d[m++]='+';
   }

   n=((fs&0x7fffffff)>>23)-127;                  /* 得到移位的个数 */
   if (n==0)
   {
       d[m]=0;
   }
   if (n<=23)
   {   /* 科学计数转化为正常计数,再去掉小数部分 */
       d[m]=( (fs&(~(~0<<23)))|(1<<23) )>>(23-n);
   }
   else
   {   /* 科学计数转化为正常计数,再加上丢掉的部分 */
       d[m]=( (fs&(~(~0<<23)))|(1<<23) )<<(n-23);
   }
}
/* power: 求底的n次幂: n>=0 这个函数不是一个实用的乘幂函数,它只能处理比较小的整数的正整数次幂 */
int power(int base,int n)
{
    int i,p;
    p=1;
    for (i=1;i<=n;++i)
        p=p * base;
    return p;
}