基本思想
方程的一种求解方法:
把函数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;
}