最优化期末大作业.docxVIP

西安电子科技大学 ——电子工程学院 最优化方法及应用 上 机 报 告 科 目：_______最优化方法及应用__________ 任课老师：_________ _ _______________ 姓 名：__________ ____ ___________ 学 号:___________ ________________ 班 级：_________ ______________ 问题描述(一) 用牛顿法和最速下降法求解以下优化问题： mi QUOTE  最速下降法 1）、算法程序 function min=grads() syms x1 x2 x3 x4 p f=(x1+10*x2)^2+5*(x3-x4)^2+(x2-2*x3)^4+10*(x1-x4)^4; e=0.01; k=1; x=[-0.1 , -2.0 , 0.5 , -2.5]; df=[diff(f,x1) diff(f,x2) diff(f,x3) diff(f,x4)]; dk=-subs(df,[x1,x2,x3,x4],x); while (norm(dk))= e (k150) f=(x1+10*x2)^2+5*(x3-x4)^2+(x2-2*x3)^4+10*(x1-x4)^4; df=[diff(f,x1) diff(f,x2) diff(f,x3) diff(f,x4)]; dk=-subs(df,[x1,x2,x3,x4],x); dk=subs(dk); x=yw(x,dk) ; k=k+1; end f=subs(f,[x1,x2,x3,x4],x); disp( 初始点为： ) disp( x = [-0.1 , -2.0 , 0.5 , -2.5] ) disp( 迭代次数为： ) disp(k) disp( 最小点是： ) disp(x) disp( 最小值是： ) disp(f) return end 一维搜索（0.618法）： function m=yw(x,dk) e=0.01; syms x1 x2 x3 x4 f=(x1+10*x2)^2+5*(x3-x4)^2+(x2-2*x3)^4+10*(x1-x4)^4; a=x; b=x+dk; while norm(b-a)e q1=a+0.328*(b-a); q2=a+0.618*(b-a); f1=subs(f,[x1,x2,x3,x4],q1); f2=subs(f,[x1,x2,x3,x4],q2); if f1f2 a=q1; else b=q2; end f=(x1+10*x2)^2+5*(x3-x4)^2+(x2-2*x3)^4+10*(x1-x4)^4; end m=(a+b)/2; 2）、程序运行结果 3）运行结果分析： 经算得，得到的最小值为0.0094 牛顿法 1）、算法程序 function min=new() e=0.001; k=1; syms x1 x2 x3 x4 ; f=(x1+10*x2)^2+5*(x3-x4)^2+(x2-2*x3)^4+10*(x1-x4)^4; df=[diff(f,x1) diff(f,x2) diff(f,x3) diff(f,x4)]; x=[0.1 , -2 , 1 , -1]; dd=subs(df,[x1,x2,x3,x4],x); while norm(dd)=e (k=1000) dm1=[diff(df(1),x1) ,diff(df(1),x2) ,diff(df(1),x3), diff(df(1),x4)]; dm2=[diff(df(2),x1) ,diff(df(2),x2), diff(df(2),x3), diff(df(2),x4)]; dm3=[diff(df(3),x1), diff(df(3),x2), diff(df(3),x3), diff(df(3),x4)]; dm4=[diff(df(4),x1) ,diff(df(4),x2), diff(df(4),x3), diff(df(4),x4)]; dm=[dm1 ,dm2, dm3, dm4]; dm=subs(dm,[x1,x2,x3,x4],x); dk=-dm\dd; x=yiwei(x,dk); dd=subs(df,[x1,x2,x3,x4],