ΰ

来自「计算二维无约束优化的Powell法和坐标轮换法,以及一维搜索的黄金分割法和搜索区」· 代码 · 共 32 行

clear;
h0=0.1;
x0=[0,0];
e=[1,0;0,1];
xp=[0,0];

for k=1:20
    k
    for i=1:2
        i
        x01=x0(1);x02=x0(2);d=e(i,:);
        alpha=goldmethod2(x01,x02,d,h0);  % 确定最优步长 alpha 
        alpha
        x1=x0+alpha*d; 
        x0=x1;
    end
    xn=x1;
    if abs(xn-xp)<1e-3
        break;
    end
    xp=xn;
end