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

功能：黄金分割法求解非线性方程的一个根。

•  源码

function root = Gold_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
   x1 = a + 0.382*(b - a);
   x2 = a + 0.618*(b - a);
   f1 = fun(x1);
   f2 = fun(x2);
   while abs(b - a) > eps      %  控制精度
      if f1*f2 < 0                    %  区间两端点都调整
         a = x1;
         b = x2;
      else
         fa = fun(a);
         if f1*fa > 0                 %  左端点a和x1之间没有根，调整左端点的值
            a = x2;
         else
            b = x1;                    %  左端点a和x1之间有根，调整右端点的值
         end
      end
      x1 = a + 0.382*(b - a);
      x2 = a + 0.618*(b - a);
      f1 = fun(x1);
      f2 = fun(x2);  
   end
   root = (b + a)/2;       % 输出根   
end