基于Matlab的机械优化设计课后题.docVIP

基于Matlab的机械优化设计课后题 1.一维搜索法 说明：采用0.618法进行编程 代码如下：syms t f=（需要计算的函数） x1=0;h=2; f1=subs(f,x1);x2=x1+h; f2=subs(f,x2); f3=f2-1;t=1; if (f1-f2)0 while f3f2 f3=subs(f,x2+t*h); t=t+1; end x3=x2+t*h; else f3=f2;f2=f1;t=1;f1=f2-1;x3=x2;x2=x1; while f1f2 f1=subs(f,x2-t*h); t=t+1; end x1=x2-t*h; end a=x1;b=x3; e=1e-05;k=0.618; a1=b-k*(b-a); a2=a+k*(b-a); f1=subs(f,a1); f2=subs(f,a2); c=(b-a)/2; while ce if f1f2 a=a1;a1=a2; a2=a+k*(b-a); f1=f2; f2=subs(f,a2); else b=a2;a2=a1; a1=b-k*(b-a); f2=f1; f1=subs(f,a1); end c=(b-a)/2; end 实际运行结果如下： t=a1=5；f=11; t=a1=3.2796;f=22.6590; t=a1=2;f=1.1122e-011; 2.无约束的优化设计 说明：一共采用了三种算法的编程 首先要建立步长函数，步长函数代码如下： function az=buchang(x,d,f) b=symvar(f); syms aa xn=x-d*aa; pa=subs(f,b,xn); pd=diff(pa,aa); aa=solve(pd); aa=real(aa); aa=subs(aa); n=size(aa,1);z=zeros(1,n); for i=1:n xn(:,i)=x-aa(i)*d; end for i=1:n z(i)=subs(f,b,xn(:,i)); end s=z(1);c=1; for i=1:n-1 if sz(i+1) s=z(i+1);c=i+1; end end az=aa(c); end DFP：·¨ syms （函数中的变量） f=（所需要优化的函数） x0=?????;b=symvar(f); n=size(b,2);g=cell(n,1);G=cell(n,n);H0=eye(n,n); for i=1:n g{i,1}=diff(f,b(i)); end for i=1:n %2úéúo￡è????ó for j=1:n G{j,i}=diff(g{i},b(j)); end end g0=subs(g,b,x0); d0=H0*g0; a=buchang(x0,d0,f); x1=x0-a*d0; xcha=x1-x0; dis=mo(xcha);k=1 while dis1e-5 g1=subs(g,b,x1); if g1==0 Gy=subs(G,b,x1) break end y0=g1-g0; s0=x1-x0; H1=H0+(s0*s0)/(s0*y0)-(H0*y0*y0*H0)/(y0*H0*y0); d1=H1*g1; a=buchang(x1,d1,f); x2=x1-a*d1; g0=g1;x0=x1;x1=x2;H0=H1; xcha=x1-x0;dis=mo(xcha); x1 subs(f,b,x1) k=k+1 end 2.坐标轮换法： syms t1 t2 f=(t1^2+t2-11)^2+(t1+t2^2-7)^2; b=symvar(f);n=size(b,2); A=eye(n,n);x0=[1;1]; dis=1;k=0; while dis1e-5 if k==30 break end e=A(:,1); a=buchang(x0,e,f); x1=x0-e*a;i=2; while in+1 e=A(:,i); a=buchang(x1,e,f);