1、直接进行关系比较的错误(==)
浮点数可以进行比较,但是由于不同浮点类型中表示精度的差异,所以会引起一些错误。
例1 :
#include <iostream.h>
void main()
{
float f1=7.123456789;
float f2=7.123456787;cout<<(f1!=f2?”not same\n”:”same\n”);
float g=1.0/3.0;
double d=1.0/3.0;
cout<<(g = =d?”same\n”:”not same\n”);
}
运行结果:same //f1与f2相同 not same //g与d不同
f1与f2的前6位有效数位相等而后面的数不同,在计算机中可能被表示为同一个数。因float型精度有限,不能分辨其差异,造成判断有误(认为是相同的数)。解决的方法是使用double型数据。
例2:
#include <iostream.h>
#include <math.h>
#include <iomanip.h>
void main()
{
double d1=123456789.9*9;
double d2=1111111109.1;
cout<<(d1==d2?" same\n":"not same\n");cout<< ((fabs(d1-d2)<1e-6)?"same\n":"not same\n");
cout<<setprecision(9);
cout<<setiosflags(ios::fixed)<<d1<<"\n"<<d2<<"\n";
}
运行结果: not same same1111111109.100000143 1111111109.099999905
2、使用浮点数进行相等(==)和不相等(!=)比较的操作通常是有问题的。浮点数的相等比较,一般用两者相减的值是否落在0的邻域中来判断。
3、测试一个浮点值与零的关系也存在误差。
在林锐的<<高质量C/C++编程>>中提到:不可将浮点变量用“==”或“!=”与任何数字比较。原因是有精度限制。
应该转换成以下形式
if ((x>=-EPSINON) && (x<=EPSINON)) 其中EPSINON是允许的误差(即精度)。 浮点数与0的比较:之所以要这样if(fabs(i)<=1e-6),是因为对于一个数,在数学上,只要这个数小于一个任意小的正数,我们就可以说这个数等于零。这是与极限相关的一个概念。如果想深入了解为什么,可以看看数学分析。
所以,在计算机编程时,我们之所以用这样的方法来比较实型,是因为这样可以任意达到所需的精度,也就是说fabs(i)可以小于更大或者更小的一个数来达到确认是否等于零的目的。
浮点数与零的比较:
等于0关系: fabs(i)<=1e-6 大于0关系: i>1e-6 小于0关系: i<1e-6 例:求一元二次方程的解
a=0 不是二次方程,是一次方程 b2-4ac=0 有两个相等实根 b2-4ac>0 有两个不等实根 b2-4ac<0 有两个共轭复根
//c代码
#include <stdio.h>
#include <math.h>
main( )
{
double a,b,c,disc,x1,x2,realpart,imagpart;
printf("请输入a,b,c:");
scanf("%lf%lf%lf",&a,&b,&c);
if(fabs(a)<=1e-6)
printf("不是2次方程, 是1次方程\n");
else {
disc=b*b-4*a*c;
if(fabs(disc)<=1e-6)
printf("有两个相等的实根:%f",(-b/(2*a)));
else if(disc>1e-6) {
x1=(-b+sqrt(disc))/(2*a);
x2=(-b-sqrt(disc))/(2*a);
printf("有两个不等的实根:\n");
printf("x1=%f\n",x1);
printf("x2=%f\n",x2);
}
else {
realpart=-b/(2*a);
imagpart=sqrt(-disc)/(2*a);
printf("有两个复根:\n");
printf("%f+%fi\n",realpart,imagpart);
printf("%f-%fi\n",realpart,imagpart);
}
}
}
//C++代码#include <iostream.h>
#include <math.h>
#include <iomanip.h>
void main( )
{
double a,b,c,disc,x1,x2,realpart,imagpart;
cout<<"请输入a,b,c:";
cin>>a>>b>>c;
if(fabs(a)<=1e-6)
cout<<"不是2次方程, 是1次方程\n";
else {
disc=b*b-4*a*c;
if(fabs(disc)<=1e-6)
cout<<"有两个相等的实根:"<<(-b/(2*a));
else if(disc>1e-6) {
x1=(-b+sqrt(disc))/(2*a);
x2=(-b-sqrt(disc))/(2*a);
cout<<"有两个不等的实根:\n";
cout<<setprecision(5)<<"x1="<<x1<<endl;
cout<<setprecision(5)<<"x2="<<x2<<endl;
}
else {
realpart=-b/(2*a);
imagpart=sqrt(-disc)/(2*a);
cout<<"有两个复根:\n";
cout<<setprecision(5)<<realpart<<"+"<<imagpart<<"i"<<endl;
cout<<setprecision(5)<<realpart<<"-"<<imagpart<<"i"<<endl;
}
}
}