优化方法-上机大作业

《优化方法-上机大作业》由会员分享，可在线阅读，更多相关《优化方法-上机大作业（20页珍藏版）》请在金锄头文库上搜索。

1、优化方法 上机大作业优化方法 上机大作业机械工程与材料能源学部能源与动力学院能源与环境工程 联系方式：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-

2、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

3、= 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)+epsi

4、lon), 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);