非线性方程的数值解法——二分法

非线性方程的数值解法——二分法

考虑单变量线性方程：
$$f(x)=0$$
第一步：先通过隔离法选取初始的含根区间。若已知f(x)在(a,∞)上有根，取步长h>0，计算f(a+ih)*f(a+(i+1)h)，i=0,1… 直到出现j，使得f(a+jh)*f(a+(j+1)h)<0，含根区间为(a+jh,a+(j+1)h)，下文记为(a,b)。
第二步：二分法计算根。

二分法总是收敛的。实际计算中，对于给定精度epsilon>0，当区间长度|bk-ak|<2epsilon时终止运算，取x*≈x_k+1=1/2(ak+bk)，此时误差小于半个区间长度：
$$|x_{k+1}-x\ast |<=\frac{1}{2}(b_k-a_k)=\frac{1}{2^{k+1}}(b-a)<ϵ$$
所以可以用上式求出二分的次数k，使算法在有限步终止。

matlab实现：
与上面不同的是，给定已知的求根区间(a,b)；结束计算的条件是函数值的差小于误差而非根的误差。

funciton[r,k] = half(fun,a,b,e)
%Inputs  -fun 需要求根的函数
         -a,b 含根区间的左右端点
         -e   容忍误差
%outputs -r   根
         -k   二分次数
while abs(fun(b)-fun(a))>e    %函数值的差小于等于容忍误差时停止计算
    c = 0.5*(a+b);
    if fun(c)==0
    	k = k+1;
        break;
    end
    if fun(a)*fun(c)<0
        b = c;
    elseif fun(b)*fun(c)<0
        a = c;
    end
    k = k+1; % 二分次数
end
r = c;  % 根

%函数
fun = inline('x^3-2*x-5','x');
%调用half
[r,k] = half(fun,2,3,10^-3);

二分法的优点是算法简单，总是收敛。
缺点是收敛速度慢，一般用于求近似根。

参考资料：《数值计算方法》；倪勤等；高等教育出版社