最优化方法 Armijo非精确搜索+DFP

成 绩： 指导教师签名： 批阅日期： 附录1：源 程 序

clc; clear; syms x1,x2; f=(x1-1)^2+5*(x2-x1^2)^2; v=[x1,x2]; df=jacobian(f,v); df=df'; x0=[2,0]'; g1=subs(df,{x1,x2},{x0(1,1),x0(2,1)}); k=0; H0=[1,0;0,1]; while(norm(g1)>0.001&k<300) if k==0 p=H0*g1; else vk=sk/(sk'*yk)-(H0*yk)/(yk'*H0*yk); w1=(yk'*H0*yk)*vk*vk'; H1=H0-(H0*yk*yk'*H0)/(yk'*H0*yk)+(sk*sk')/(sk'*yk)+w1; p=H1*g1; H0=H1; end x00=x0; result=Usearch1(f,x1,x2,df,x0,p); arf=result(1); x0=x0+arf*p; g0=g1; g1=subs(df,{x1,x2},{x0(1,1),x0(2,1)}); p0=p; yk=g1-g0;sk=x0-x00; k=k+1; end Function result=Usearch1(f,x1,x2,df,x0,p) mu=0.001;sgma=0.99;a=0;b=inf;arf=1; pk=p; x3=x0; x4=x3+arf*pk; 4

f1=subs(f,{x1,x2},{x3(1,1),x3(2,1)}); f2=subs(f,{x1,x2},{x4(1,1),x4(2,1)}); gk1=subs(df,{x1,x2},{x3(1,1),x3(2,1)}); gk1=subs(df,{x1,x2},{x4(1,1),x4(2,1)}); while(f1-f2<=-mu*arf*gk1'*pk) b=arf;arf=(a+arf)/2;x4=x3+arf*pk; f2=subs(f,{x1,x2},{x4(1,1),x4(2,1)}); gk2=subs(df,{x1,x2},{x4(1,1),x4(2,1)}); end while(1>0) if(gk2'*pk

附录2：实验报告填写说明

1．实验项目名称：要求与实验教学大纲一致。

2．实验目的：目的要明确，要抓住重点，符合实验教学大纲要求。 3．实验原理：简要说明本实验项目所涉及的理论知识。 4．实验环境：实验用的软、硬件环境。

5．实验方案（思路、步骤和方法等）：这是实验报告极其重要的内容。概括整个实验过程。

5

对于验证性实验，要写明依据何种原理、操作方法进行实验，要写明需要经过哪几个步骤来实现其操作。对于设计性和综合性实验，在上述内容基础上还应该画出流程图、设计思路和设计方法，再配以相应的文字说明。对于创新性实验，还应注明其创新点、特色。

6．实验过程（实验中涉及的记录、数据、分析）：写明具体实验方案的具体实施步骤，包括实验过程中的记录、数据和相应的分析。

7．实验结论（结果）：根据实验过程中得到的结果，做出结论。 8．实验小结：本次实验心得体会、思考和建议。

9．指导教师评语及成绩：指导教师依据学生的实际报告内容，给出本次实验报告的评价。

6