fsxgold.m

来自「非线性规划中的最速下降法,黄金分割法,阻尼牛顿法,牛顿切线法.」· M 代码 · 共 11 行

function x=fsxgold(f,e,a,b)
% fsxgold函数 黄金分割算法 
% x=fsxgold(f,e,a,b)为输入函数 f为函数 e为允许误差 (a,b)为搜索空间;
L=(5^(1/2)-1)/2;
u=a+L*(b-a);f2=subs(f,u);t=a+b-u;f1=subs(f,t);
while (abs(u-t)>=e)
    if (f1<=f2)  b=u;u=t;f2=f1;t=a+b-u;f1=subs(f,t);
    else  a=t;t=u;f1=f2;u=a+L*(b-a);f2=subs(f,u); 
end;
end;
x=0.5*(t+u);