问题描述
我正在尝试从头开始实现二分根查找算法,我传递的方程是 x^2 - 42 = 0。我的代码不断给出输出“42 的根在 1e 的误差内为 0 -05. ",意思是我的二分算法的递归栈由于某种原因一直返回0。代码如下:
#include <iostream>
#include <cmath>
using namespace std;
class bisection
{
public:
double a;
double b;
double c;
double epsilon;
void set_values(double,double,double);
double f(double x)
{
return x*x - 42;
}
double solve(double a,double b,double epsilon)
{
c = (a+b)/2.0;
if(abs(f(c)) < epsilon)
{
return c;
}
else
{
if (f(a) * f(c) < 0.0)
{
b = c;
solve(a,b,epsilon);
}
else if (f(c) * f(b) < 0.0)
{
a = c;
solve(a,epsilon);
}
}
return 0;
}
};
void bisection::set_values(double left,double right,double error)
{
a = left;
b = right;
epsilon = error;
}
int main()
{
bisection myObj;
myObj.set_values(0.0,10.0,0.00001);
//cout << myObj.f(7);
cout << "The root of 42 is " << myObj.solve(myObj.a,myObj.b,myObj.epsilon) << " within an error of " << myObj.epsilon << ".";
return 0;
}
我必须包含“return 0;”行,因为我的编译器在没有它的情况下抛出“非空函数不会在所有控制路径中返回值”。如何更改我的代码以包含所有递归路径的返回值但仍然产生正确的答案?我对算法的基本设计是否有错误,或者是一个简单的修复?
解决方法
除了代码的其他问题之外,您还需要对 solve
的结果进行实际操作。通常这意味着 return
直接从 solve
获取结果,直到达到基本情况并返回:
double solve(double a,double b,double epsilon)
{
c = (a+b)/2.0;
if(abs(f(c)) < epsilon)
{
return c; //this is your base case
}
else
{
if (f(a) * f(c) < 0.0)
{
b = c;
return solve(a,b,epsilon);
}
else if (f(c) * f(b) < 0.0)
{
a = c;
return solve(a,epsilon);
}
}
}