try1.m

来自「用随即搜索方法编写的m文件」· M 代码 · 共 36 行

clear;
x0=[3,2]';
e1=[0.5,0.5]';
a1=0;%设置k的区间
a2=1;
a=a1+0.382*(a2-a1);
b=a1+0.618*(a2-a1);

while b-a>0.001
    
    xa=x0+a*e1;
    xb=x0+b*e1;
    fa=getf(xa);
    fb=getf(xb);
    
    if fa>fb
        a1=a;
        a2=a2;
        a=a1+0.382*(a2-a1);
        b=a1+0.618*(a2-a1);
        xa=x0+a*e1;
        xb=x0+b*e1;
        fa=getf(xa);
        fb=getf(xb);
    else
        a1=a1;
        a2=b;
        a=a1+0.382*(a2-a1);
        b=a1+0.618*(a2-a1);
        xa=x0+a*e1;
        xb=x0+b*e1;
        fa=getf(xa);
        fb=getf(xb);
    end
end
y=(a+b)/2