geta.m

用随即搜索方法编写的m文件

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+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;
end