goldstein.m

来自「黄金分割法的matlab程序」· M 代码 · 共 34 行

%goldtein.m,goldstein法的通用子程序,输入的参数为初始步长
function a1=goldstein(a1,p,t)
%a,b为搜索区间的上下界
a=0;
b=10;
m=f(0);
n=ff(0);
%循环的标志变量
flag=0;
%循环求解最佳的步长
while flag==0 
    %搜索步长满足goldstein第一个准则
    if f(a1)<=m+p*a1*n
        %搜索步长满足goldstein第二个准则
        if f(a1)>=m+(1-p)*a1*n
            %输出最佳的步长
            a1=a1;
            flag=1;
            return;
        %搜索步长不满足goldstein准则，继续迭代
        else
            a=a1;
            a1=0.5*(a+b);
            if b<inf
                a1=0.5*(a+b);
            else
                a1=t*a1;
            end
        end
    else
        b=a1;
        a1=0.5*(a+b);
    end
end