hooke.txt

来自「hooke算法的自我编写的小小的程序能够用来计算问题的最优值」· 文本 代码 · 共 27 行

function hooke(x1)
u=0.001;t=1;b=0.5;a=1.5;
y1=x1;
e=[1,0;0,1];
while t>u
    for j=1:2
        if fun(y1+t*e(j,:))<fun(y1)
            y2=y1+t*e(j,:);
            y1=y2;
        elseif fun(y1-t*e(j,:))<fun(y1)
            y2=y1-t*e(j,:);
            y1=y2;
        else
            y2=y1;
        end
        if fun(y2)<fun(x1)
            x2=y2;
            y1=x2+a*(x2-x1);
            x1=x2;
        elseif t<u
            return
        else
            t=b*t;y1=x1;x2=x1;
        end
    end
end
x1