我正在C ++中应用Regula Falsi方法,但是问题在于==运算符,当F(x3)变为0时,如果(fabs(f(x3))== 0应该停止并从循环中出来,但它并没有停止,为什么为什么为什么....
像第12次迭代后的以下输出f(x3)= 0,但是if(fabs(f(x3)== 0))无法运行。循环不会停止,因此不应进行第13次迭代
float f(float x)
{
float f_x;
f_x= pow(x,3)+(3*x)-5;
return f_x;
}
int main(int argc, char** argv)
{
float a,b,tol,x3;
int itr,n;
cout << "enter the iterations";
cin >> itr;
cout << "enter the interval a";
cin >> a;
cout <<"enter the interval b";
cin >> b;
cout << "enter the toleration";
cin >> tol;
cout.setf(std::ios_base::fixed, std::ios_base::floatfield);
cout.precision(5);
//cout<<"fa="<<f(a)<<"fb"<<f(b);
cout<<"n\t\ta\t\tb\t\tx3\t\tf(a)\t\tf(b)\t\tf(x3)" <<endl;
if (f(a)*f(b)<0 && a<b)
{
for (n=0;n<itr;n++)
{
x3=a-((b-a)*f(a))/(f(b)-f(a));
cout << "xx"<<fabs(f(x3));
if (fabs( f(x3))==0)
{
cout << "Solution"<<fabs(f(x3));
break;
}
else
{
cout<<n+1 <<"\t\t"<<a <<"\t\t"<<b <<"\t\t"<<x3<<"\t\t"<<f(a)
<<"\t"<<f(b)<<"\t\t"<<f(x3) <<endl;
if(f(x3)*f(a)<0)
b=x3;
else
if(f(x3)*f(b)<0)
a=x3;
}
}
}
else
cout<< "No Solution Exist";
return 0;
}
输出值
输入迭代13
输入间隔a1
输入间隔b2
输入公差1
**n a b x3 f(a) f(b) f(x3)**
1 1.00000 2.00000 1.10000 -1.00000 9.00000 -0.36900
2 1.10000 2.00000 1.13545 -0.36900 9.00000 -0.12980
3 1.13545 2.00000 1.14774 -0.12980 9.00000 -0.04487
4 1.14774 2.00000 1.15197 -0.04487 9.00000 -0.01542
5 1.15197 2.00000 1.15342 -0.01542 9.00000 -0.00529
6 1.15342 2.00000 1.15391 -0.00529 9.00000 -0.00181
7 1.15391 2.00000 1.15408 -0.00181 9.00000 -0.00062
8 1.15408 2.00000 1.15414 -0.00062 9.00000 -0.00021
9 1.15414 2.00000 1.15416 -0.00021 9.00000 -0.00007
10 1.15416 2.00000 1.15417 -0.00007 9.00000 -0.00003
11 1.15417 2.00000 1.15417 -0.00003 9.00000 -0.00001
12 1.15417 2.00000 1.15417 -0.00001 9.00000 0.00000
13 1.15417 2.00000 1.15417 -0.00000 9.00000 0.00000
最佳答案
这里的问题不是浮点精度。这就是您愿意接受结果的容忍度。在大多数情况下,随着您进行更多的迭代,Regula falsi将使您越来越接近正确的结果,但它无法提供确切的答案。因此,您要做的决定是,您希望结果接近多少?这是现实世界中对准确性的要求与获得结果所需的时间之间的权衡;更高的精度需要更多的计算时间。因此,请为您的任务选择一个可接受的公差,然后重复循环直到结果在该公差之内。如果结果太慢,则必须增加公差。
关于c++ - C++ ==运算符不起作用,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/18932336/