yijietidu.m

来自「用一阶梯度法」· M 代码 · 共 42 行

%多变量的一阶梯度寻优
%变量数为n，允许迭代的最大次数M,精度要求：0.001

n=4;                  %变量数
x0=[-3;-1;-3;-1];     %初始值
x1=x0;                 %初始点
M=2000;
e=0.001;

df0=dfm(x0);
s=-df0;                   %初始方向

nx1=1;                    %终止判据

tic

for i=1:M
    
    if nx1<e              %终止判断
        break;
    end
    
    for j=1:n             %对各变量方向寻优
       
    a=fpart618(x1,s);     %采用618法进行单变量寻优
    x2=x1+a*s;            %代表k＋1次迭代点
     
    x0=x1;
    x1=x2;
    
    s=-dfm(x1);               %k+1次搜索方向
    
    nx1=norm(x1-x0);
    
    end
    
end

i
x1
fm(x1)
toc