二分法求解函数零点

已知函数零点在[-1,1]区间,只有一个零点,求函数零点,精度在1e-4.
注意点:

  1. 使用fabs而不是使用abs,abs返回的是int!!!这是个隐藏的bug
  2. fabs(mid_val) < thr)即可,这就相当于是精度了,因为不用|mid_val|<thr的话,只能用mid_val == 0了,而要求的精度只到thr
  3. 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;
}
  • 9
    点赞
  • 7
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值