%多變量的一階梯度尋優(yōu)
%變量數(shù)為n，允許迭代的最大次數(shù)M,精度要求：0.001

n=4;                  %變量數(shù)
x0=[-3;-1;-3;-1];     %初始值
x1=x0;                 %初始點
M=2000;
e=0.001;

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

nx1=1;                    %終止判據(jù)

tic

for i=1:M
    
    if nx1<e              %終止判斷
        break;
    end
    
    for j=1:n             %對各變量方向尋優(yōu)
       
    a=fpart618(x1,s);     %采用618法進行單變量尋優(yōu)
    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