最速下降法牛顿法.docVIP

最速下降法 最速下降法迭代公式是 计算步骤如下： （1）给定初点，允许误差＞0，令k=0。 （2）计算搜索方向 （3）若，则 ，停止；否则令 ，由一维搜索步长 ，使得 （4）令 ，k=k+1，转步骤（2）。 用最速下降法求min f=(x1-4)^4+(x2+2)^2+1,初始点是（1，-3） %允许误差e %zk是给定初点 %z是找到的最优点 %k是迭代次数 %选用步长公式; syms x1 x2 f=(x1-4)^4+(x2+2)^2+1;%目标函数 e=1.0e-4; df_dx1=diff(f,x1); df_dx2=diff(f,x2); f1=[df_dx1;df_dx2]; flag=1; k=0; zk=[1;-3]; while flag g=subs(f1,[x1;x2],zk); %梯度g(zk) norm_g=norm(g); if norm_g=e%如果g的模小于允许误差，则已找到最优点 flag=0; z=zk; zuiyou_f=subs(f,[x1;x2],z); disp(极小点是) z disp(最小值是) zuiyou_f disp(迭代次数) k break else %g的模不小于允许误差的情况 k=k+1; p=-g; %负梯度 h=jacobian(f1); H=subs(h,[x1;x2],zk); %Hesse矩阵 H（xk） t=(g*g)/(g*H*g); %t是搜索步长 zk=zk+t*p; end end 牛顿法 已知目标函数及其梯度，Hesse矩阵，终止限 。 （1）选定初始点；计算 ；置k=0。 （2）计算 。 （3）由方程 解出。 （ 4 ）计算 。 判别终止准则是否满足：如满足，则打印最优解（ ）结束；否则，置k=k+1，转（2）。 用牛顿法求min f=(x1-4)^4+(x2+2)^2+1,初始点是（1，-3） %允许误差e %zk是给定初点 %z是找到的最优点 %k是迭代次数 clear all close all clc syms x1 x2 f=(x1-4)^4+(x2+2)^2+1;%目标函数 e=1.0e-10; df_dx1=diff(f,x1); df_dx2=diff(f,x2); f1=[df_dx1;df_dx2]; flag=1; k=0; zk=[1;-3]; while flag fk=subs(f,[x1;x2],zk); gk=subs(f1,[x1;x2],zk); fk=subs(f,[x1;x2],zk); %f(xk) gk=subs(f1,[x1;x2],zk); %梯度 mo_gk=norm(gk); if mo_gk=e %如果g的模小于允许误差则已找到最优点 z=zk; zuiyou_f=subs(f,[x1;x2],z); disp(极小点是) z disp(最小值是) zuiyou_f disp(迭代次数) k break else %g的模不小于允许误差的情况 k=k+1; Gk=subs(jacobian(f1),[x1;x2],zk); %Hesse矩阵 pk=-inv(Gk)*gk;%搜索步长 zk=zk+pk; end end