二分法

基本思想

方程的一种求解方法：
把函数f(x)的零点所在的区间[a，b](满足f(a）●f(b）<0)“一分为二”，得到[a，m]和[m，b]。根据“f(a)●f(m)<0”是否成立，取出零点所在的区间[a，m]或[m，b]，仍记为[a，b]。所对得的区间[a，b]重复上述步骤，直到包含零点的区间[a，b]“足够小”，则[a，b]内的数可以作为方程的近似解。

举例

用二分法求方程2x^3+4x*x+3x-6=0的一个根，并要求绝对误差不超过0.001.

算法

给定精确度ξ,用二分法求函数f(x)零点近似值的步骤如下:
1 确定区间[a,b],验证f(a)•f(b)<0,给定精确度ξ.
2 求区间(a,b)的中点d.
3 计算f(d).
(1) 若f(d)=0,则c就是函数的零点;
(2) 若f(a)•f(d)<0,则令b=d;
(3) 若f(d)•f(b)<0,则令a=d.
(4) 判断是否达到精确度ξ:即若|a-b|<ξ,则得到零点近似值a(或b),否则重复2-4.

程序

double funx(double  x)
{   return(2*x*x*x-4*x*x+3*x-6);  }
double fun( double  a, double  b)
{
    double d;
    d=(a+b)/2;
    while(fabs(a-b)>0.001)
    {   if(funx(a)*funx(d)<0) b=d;
		if(funx(d)*funx(b)<0) a=d;
	d=(a+b)/2;
    }
    return  d;
}