关于浮点数的比较

我们都知道,浮点数比较,应该是比较一个范围,如果在这个范围内,我们认为相等,为什么不能用"=="比较呢?下面给出几段代码,欢迎大家讨论.

#include <stdio.h>

int main()
{
    float f = 1/3.0;
    printf("f = %f\n", f);
    if (1/3.0 == f)
    {
        printf("Hello World!\n");
    }

    return 0;
}

这段代码,编译会有警告:warning C4305: 'initializing' : truncation from 'const double ' to 'float '.也就是说,初始化变量的时候,将double型截断为float型,有精度损失.

同理,则if语句里面的1/3.0也是double型的了,那么,后面的f则会由float型提升为double型,这里,精度不一样的问题就出来了,此代码不会输出"Hello World!\n".

那么我们再看一段代码:

#include <stdio.h>

int main()
{
    float f = 1/3.0f;
    printf("f = %f\n", f);
    if (1/3.0f == f)
    {
        printf("Hello World!\n");
    }

    return 0;
}

这段代码,编译没有警告,而且,能够输出"Hello World!\n".

#include <stdio.h>

int main()
{
    double f = 1/3.0;
    printf("f = %f\n", f);
    if (1/3.0 == f)
    {
        printf("Hello World!\n");
    }
    return 0;
}

这段代码,编译没有警告,而且,也能够输出"Hello World!\n".

但是,请注意,这里都是没有对数据进行运算过,下面的一个例子

#include <stdio.h>

int main()
{
    float f = 1.111f+2.222f;
    printf("f = %f\n", f);
    if (3.333f == f)
    {
        printf("Hello World!\n");
    }
    return 0;
}

就已经 不能够输出"Hello World!\n"了!

当然,如果编译器将前面的加法给优化了的话,可能会输出 "Hello World!\n". 

因计算机里面存放数据都是二进制的(具体存储方式请查询别人的其它文章),所以表示的精度和十进制所表示的精度不一样,所以,有些我们十进制能够很方便表示的数据,使用二进制却不能够完整的表示出来,只能够表示个近似值,就如同使用十进制数表示无理数一样,我们也需要取个合理的近似值(如π我们常用3.14).

当对数据进行运算,因有些数据不能使用二进制表示,则存在误差,不能表示的数据越多,则误差越大,以致不能使用"=="来比较浮点数,就是存在这方面的误差.

最后比较浮点数,还是使用范围比较,范围的精度看实际情况.

(文章中如有什么不正确的地方,欢饮指正,不甚感激!)