jintuifa.asv

黄金分割法

%宋斌
%进退法求函数极小值
a0=5;                               %初始点的位置
h0=1.0e-5;                          %初识步长
c_jingdu=1.0e-6;                    %精度
x=a0;
h=h0;
f1=x^4-4*x^3-6*x^2-16*x+4;          %计算f(a0)
x1=x+h;
f2=x1^4-4*x1^3-6*x1^2-16*x1+4;      %计算f(a0+h0)
while(1)
    if f2<=f1
        f1=f2;
        x=x+h;
        f2=x^4-4*x^3-6*x^2-16*x+4;
    else
        if f2>=f1
            if abs(h)<=c_jingdu
                break
            else
            f1=f1;
            h=-h/4;
            x=x+h;
            f2=x^4-4*x^3-6*x^2-16*x+4;
            end
        end
    end
end
x