golden.m

黄金分割法 最优化 求极值 ddddddddddddddddddddddddddd

function [x,f]=golden(fun,eps,a,b)
x1=a+0.382*(b-a);
x2=a+0.618*(b-a);
f1=feval(fun,x1);
f2=feval(fun,x2);
k=1;
while b-a>=eps
    if(f1>f2)
        a=x1;
        x1=x2;
        x2=a+0.618*(b-a);
        f1=f2;
        f2=feval(fun,x2);
    else
        b=x2;
        x2=x1;
        x1=a+0.382*(b-a);
        f2=f1;
        f1=feval(fun,x1);
    end
    k=k+1;
end

x=(b+a)/2;
f=feval(fun,x);