已知函数零点在[-1,1]区间,只有一个零点,求函数零点,精度在1e-4.
注意点:
- 使用fabs而不是使用abs,abs返回的是int!!!这是个隐藏的bug
- fabs(mid_val) < thr)即可,这就相当于是精度了,因为不用|mid_val|<thr的话,只能用mid_val == 0了,而要求的精度只到thr
- if (mid_val * left_val < 0) { r = mid; } 只需要计算一边,这里可以找一个示例图画图看一下
#include <vector>
#include <iostream>
#include<cmath>
float func(double x) {
return x * x * x - x + 0.2;
}
float find_zero(float l, float r, float thr)
{
if (func(l) == 0) {
return l;
}
if (func(r)==0) {
return r;
}
while(r - l > thr) { // 这里以thr为判别条件。
float mid = (r+l) / 2;
float mid_val = func(mid);
float left_val = func(l);
// 注意,一定要用fabs!!!因为abs的输入是int,会导致abs(0.2)<1e-4不成立!!!
if (fabs(mid_val) < thr) { //这就相当于是精度了,因为不用|mid_val|<thr的话,只能用mid_val == 0了,而要求的精度只到thr
return mid;
} else if (mid_val * left_val < 0) { //注意,只需要计算一边,这里可以找一个示例图画图看一下
r = mid;
} else {
l = mid;
}
}
return (r + l) / 2;
}
int main(){
float res = find_zero(-1,1,1e-4);
std::cout<< res<<" "<< func(res) <<std::endl;
return 0;
}