优化方法 上机大作业.docVIP

优化方法 上机大作业 机械工程与材料能源学部 能源与动力学院 能源与环境工程 联系方式 x0=[0;1]T;%初始值 s0=[-1;1]T;%初始搜索方向 c1=0.1;c2=0.5;a=0;b=inf;d=1;n=0; x1=x0+d*s0; g0=[(x0(2)-x0(1)^2)*x0(1)-2*(1-x0(1));(x0(2)-x0(1)^2)]; g1=[(x1(2)-x1(1)^2)*x1(1)-2*(1-x1(1));(x1(2)-x1(1)^2)]; f1=(x1(2)-x1(1)^2)^2+(1-x1(1))^2; f0=(x0(2)-x0(1)^2)^2+(1-x0(1))^2; while((f0-f1-c1*d*g0*s0)||(g1*s0c2*g0*s0)) if ((f0-f1)(-c1*d*g0*s0)) b=d;d=(d+a)/2; x1=x0+d*s0; g0=[(x0(2)-x0(1)^2)*x0(1)-2*(1-x0(1));(x0(2)-x0(1)^2)]; g1=[(x1(2)-x1(1)^2)*x1(1)-2*(1-x1(1));(x1(2)-x1(1)^2)]; f1=(x1(2)-x1(1)^2)^2+(1-x1(1))^2; f0=(x0(2)-x0(1)^2)^2+(1-x0(1))^2; elseif (((g1)*s0)(c2*(g0)*s0)) a=d; if(2*d=(d+b)/2) d=2*d; else d=(d+b)/2; end x1=x0+d*s0; g0=[(x0(2)-x0(1)^2)*x0(1)-2*(1-x0(1));(x0(2)-x0(1)^2)]; g1=[(x1(2)-x1(1)^2)*x1(1)-2*(1-x1(1));(x1(2)-x1(1)^2)]; f1=(x1(2)-x1(1)^2)^2+(1-x1(1))^2; f0=(x0(2)-x0(1)^2)^2+(1-x0(1))^2; end end x1 f1=(x1(2)-x1(1)^2)^2+(1-x1(1))^2 x1 = -0.0000 1.0000 d = 1.1102e-016 f1 = 2 function f = fun( x ) %UNTITLED3 Summary of this function goes here % Detailed explanation goes here f=x(1)^2-2*x(1)*x(2)+2*x(2)^2+x(3)^2+x(4)^2-x(2)*x(3)+2*x(1)+3*x(2)-x(3); End function g = fun( x ) %UNTITLED4 Summary of this function goes here % Detailed explanation goes here g=[2 -2 0 0;-2 4 -1 0;0 -1 2 0;0 0 0 2]*x+[2;3;-1;0]; end x0=[0;0;0;0]; %初始值 eps=1.0e-4; %精度 g0=gfun(x0); s0=-g0; n=0; syms d1; while norm(g0)eps if n3 g=gfun(x0+d1*s0); d= double(solve(s0*g)); x1=x0+d*s0; g1=gfun(x1); if norm(g1)eps n=n+1; x0=x1; break else s0=-g1+(norm(g1)^2/norm(g0)^2)*s0; x0=x1; g0=g1; end elseif n=3 x0=x1; g0=gfun(x0); s0=-g0; n=0; end n=n+1; end x0 n fun(x0) x0 = -4 -3 -1 0 n = 3 ans = -8 function f= fun3_1(x ) %FUN3 Summary of this function goes here % Detailed explanation goes here f=x(1)+2*x(2)^2+exp(x(1)^2+x(2)^2); end function g= gfun3_1(x) %GFUN3_1 Summary of this function goes here % Detailed explan