geta.asv

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

function y=geta(x0,e1)

a1=0;%设置k的区间
a2=1;
a=a1+0.382*(a2-a1);
b=a1+0.618*(a2-a1);

while b-a>0.001
    
    xa=x0+k*e1;
    xb=x0+k*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+k*e1;
        xb=x0+k*e1;
        fa=getf(xa);
        fb=getf(xb);
    else
        a1=a1;
        a2=b;
        a=a1+0.382*(a2-a1);
        b=a1+0.618*(a2-a1);

        fa=getf(a);
        fb=getf(b);
    end
end
y=(a+b)/2
end