最优化方法及其matlab程序设计 马昌凤版 课后答案 杭电课件

ρ=0.55;σ=0.4;µk=norm(feval(Fk,x0));k=0;epsilon=1e-6;n=length(x0);while(k<maxk)

fk=feval(Fk,x0);%计算函数值jfk=feval(JFk,x0);%计算Jacobi阵gk=jfk’*fk;

dk= (jfk′ jfk+µk eye(n))\gk;%解方程组Gk*dk=-gk,计算搜索方向if(norm(gk)¡epsilon),break;end%检验终止准则m=0;mk=0;

while(m<20)%用Armijo搜索求步长

newf=0.5 norm(feval(Fk,x0+ρm dk))2;oldf=0.5 norm(feval(Fk,x0))2;if(newf<oldf+sigma ρm gk′ dk)mk=m;break;endm=m+1;end

x0=x0+ρmk dk;muk=norm(feval(Fk,x0));k=k+1;endx=x0;val=0.5*µ2k;

20