优化方法 上机大作业优化方法 上机大作业机械工程与材料能源学部能源与动力学院能源与环境工程 联系方式：15326178310x0=0;1T;%初始值s0=-1;1T;%初始搜索方向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*deps if n3g=gfun(x0+d1*s0);d= double(solve(s0*g);x1=x0+d*s0;g1=gfun(x1);if norm(g1)=epss0=-g0;g=gfun3_1(x0+d1*s0);d= double(solve(s0*g);x1=x0+d*s0;g1=gfun3_1(x1);if( norm(g1)=eps d=-g20g0; x1=x0+d;g1=gfun3_1(x1);if( norm(g1)=epss0=-h0*g0;g=gfun3_1(x0+d1*s0);d= double(solve(s0*g);x1=x0+d*s0;g1=gfun3_1(x1);if( norm(g1)bi(i)+epsilon), index(i)=0; endendwhile (k0), Aee=Ae; end for(j=1:ni) if(index(j)0), Aee=Aee; Ai(j,:); end end gk=H*x+c; m1,n1 = size(Aee); dk,lamk=qsubp(H,gk,Aee,zeros(m1,1); if(norm(dk)ne) y,jk=min(lamk(ne+1:length(lamk); end if(y=0) exitflag=0; else exitflag=1; for(i=1:ni) if(index(i) & (ne+sum(index(1:i)=jk) index(i)=0; break; end end end k=k+1; else exitflag=1; alpha=1.0; tm=1.0; for(i=1:ni) if(index(i)=0)&(Ai(i,:)*dk0) tm1=(bi(i)-Ai(i,:)*x)/(Ai(i,:)*dk);