二分法求函数零点

一、零点及零点存在性定理

零点定义，对于函数$y=f(x)$，使得$f(x)=0$的实数$x$叫做函数$f(x)$的零点。换句话说，函数$y=f(x)$的零点就是方程$f(x)=0$的实数根，也就是函数$y=f(x)$的图像与$x$轴的交点。

零点存在性定理，如果函数$f=f(x)$在区间$[a,b]$上的图像是连续的曲线，并且有$f(a)\cdot f(b)<0$，我们就说函数$y=f(x)$在开区间$(a,b)$内有零点，即存在$c\in (a,b)$使得$f(c)=0$。

注意：满足该定理是函数存在零点的充分不必要条件。如果该函数是一个单调函数，那么零点有且仅有一个。

二、二分法求函数零点

利用零点存在定理，可以用来求取函数的零点的近似值。二分法的基本思想是通过不断地将零点所在的区间一分为二，使得两个端点逐步逼近零点，从而得到零点近似值的方法叫做二分法。

按照上面的流程图实现一下即可：

#include <iostream>

double f(double x)
{
    return 2 * x * x * x - 4 * x * x + 3 * x - 6;
}


double findZero(double a, double b,double err)
{
    double m = 0.0;
    if (fabs(f(a)) < err)
        return a;
    if (fabs(f(b)) < err)
        return b;

    while (fabs(a - b) > err)
    {
        m = (a + b) / 2;
        if (f(a) * f(m) < 0)
        {
            b = m;
        }
        else
        {
            a = m;
        }
        if(fabs(f(m))<err)
		{
			return m;
		}
    }
    return m;
}


int main()
{
    std::cout << findZero(-10, 10, 0.001) << std::endl;
}