gongetidu.m

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

%多变量寻优的共轭梯度
%允许迭代的最大次数M,变量数N，

n=4;        %变量数
x0=[-3;-1;-3;-1];
x1=x0;      %代表k次迭代点

M=2000;
e=0.0001;

df0=dfm(x0);

s0=-df0;
s=s0;

ndf0=norm(df0);
ndf1=ndf0;

tic

for i=1:M

    if ndf1<e
        break;
    end

    for j=1:n

    a=fpart618(x1,s);     %采用618法进行单变量寻优
    x2=x1+a*s;            %代表k＋1次迭代点

    x0=x1;
    x1=x2;

    df0=dfm(x0);
    df1=dfm(x1);

    ndf0=ndf1;
    ndf1=norm(df1);

    s0=s;
    s=-df1+s0*ndf1^2/ndf0^2;            %搜索方向

    end

    s=-df1;

end

i
x1
fm(x1)
toc