【Matlab学习手记】非线性方程求解：二分法

功能：二分法求解非线性方程的一个解，采用递归的方式。

• 源码

function root = HalfInterval_Search(fun, a, b, eps)
% 二分法求函数fun在区间[a b]上的一个零点
% 输入参数说明：
%  fun --- 目标函数
%  a --- 搜索区间的下界
%  b --- 搜索区间的上界
%  eps --- 区间的最小阀值长度
%  输出参数说明：
%  root --- 函数的一个根

root = [];   %  返回空值
if fun(a) == 0
   root = a;
end;
if fun(b) == 0
   root = b;
end
if fun(a)*fun(b) > 0
   disp('两端点函数值乘积大于0，无解！');
   return;
else 
   root = FindRoot(fun, a, b, eps);
end
end

function r = FindRoot(fun, a, b, eps)
fa = fun(a);
x = 0.5*(a + b);
fx = fun(x);
if fa*fx > 0                % a和中点之间没有根
   t = 0.5*(a + b);
   r = FindRoot(fun, t, b, eps);
else
   if fa*fx == 0
      r = 0.5*(a + b);   % 输出根
   else
      if abs(b - a) <= eps
         r = 0.25*(b + 3*a);
      else                     % a和中点之间有根  
         s = 0.5*(a + b);
         r = FindRoot(fun, a, s, eps);
      end
   end
end  
end