实用最优化方法编程大作业讲述.docxVIP

实用最优化方法编程大作业班 级：姓 名：指导老师：学 号：大连理工大学2015年11月27日版本号：MATLAB 7.11.0 (R2010b)【文件名WP.m】function x = WP(f,x0,var,s0,eps)clcsyms x1 x2 ;c1=0.1;c2=0.5;b=inf;lambda=1;ifnargin==4eps=1.0e-6;endgradf=jacobian(f,var);g0=subs(gradf,var,x0);f0=subs(f,var,x0);gs=g0*s0;a=0;j=0;while j1000new_x=x0+lambda*s0;new_f=subs(f,var,new_x);left=f0-new_f;new_g=subs(gradf,var,new_x);new_gs=new_g*s0;right=-1*c1*lambda*gs; if leftright %不满足第一个条件 j=j+1; b=lambda;lambda=0.5*(lambda+a); else %满足第一个条件 left2=new_gs; right2=c2*g0; if left2right2