机械优化设计的matlab程序

《机械优化设计的matlab程序》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、黄金分割法%黄金分割法f=inline('x^2-7*x+10','x');%目标函数%搜索区间和收敛精度a=2;b=8;epsilon=0.35;x1=b-0.618*(b-a);f1=f(x1);%计算左试点x2=a+0.618*(b-a);f2=f(x2);%计算右试点fork=1:7fprintf(1,'迭代次数k=%3.0f',k)iff1<=f2%舍弃右区间b=x2;x2=x1;f2=f1;x1=b-0.618*(b-a);f1=f(x1);else%舍弃左区间a=x1;x1=x2;f1=f2;x2=a+0.618*(b-a);f2=

2、f(x2);endx=0.5*(b+a);fprintf(1,'迭代区间-左端a=%3.4f',a)fprintf(1,'试点1坐标值x1=%3.4f',x1)fprintf(1,'函数值f1=%3.4f',f(x1))fprintf(1,'迭代区间-右端b=%3.4f',b)fprintf(1,'试点2坐标值x2=%3.4f',x2)fprintf(1,'函数值f2=%3.4f',f(x2))fprintf(1,'区间中点x=%3.4f',x)disp''endfprintf(1,'区间中点函数值fx=%3.4f',

3、f(x))鲍威尔法%鲍威尔法-计算第2环起始点和搜索方向symsx1x2f=60-10*x1-4*x2+x1^2+x2^2-x1*x2;disp'目标函数：'pretty(f);%第1环沿e2方向搜索x01=0;x02=0;X0=[x01x02];%第1环初始点x1=5.0000;x2=0;X=[x1x2];%第1环沿e2方向初始点fy1=60-10*x1-4*x2+x1^2+x2^2-x1*x2;%第1环e2方向初始点函数值e1=[1,0];e2=[0,1];%坐标轴单位方向s1=e2(1);s2=e2(2);S=[s1s2];m=0;[xy1,x

4、y2,fyx]=powell(m,x1,x2,X,s1,s2,S);%计算影射点及其函数值xs1=2*xy1-x01;xs2=2*xy2-x02;Xs=[xs1xs2];fxs=60-10*xs1-4*xs2+xs1^2+xs2^2-xs1*xs2;fx0=60-10*x01-4*x02+x01^2+x02^2-x01*x02;%第1环初始点函数值disp'第1环影射点坐标'disp(Xs)%判断第1环搜索函数值下降最大方向df01=fx0-fy1;%第1环沿e1方向初始点函数值-第1环沿e1方向终点函数值df02=fy1-fyx;%第1环沿e2方向

5、起始点函数值-第1环沿e2方向始点函数值ifdf01>df02dfm=df01;m=1;elsedfm=df02;m=2;end%计算POWELL判别式f1=fx0;f2=fyx;f3=fxs;f123=[f1f2f3];disp'第1环初始点、终点、影射点函数值'disp(f123)fP1=(f1-2*f2+f3)*(f1-f2-dfm)^2;fP2=0.5*dfm*(f1-f3)^2;fP=[fP1fP2];disp'两个POWELL判别式的值'fprintf(1,'POWELL判别式1的值fP1=%10.4f',fP1)fprintf(1,

6、'POWELL判别式1的值fP2=%10.4f',fP2)fprintf(1,'函数值下降最大方向m=%2.0f',m)iff1>f3&fP2>fP1sx1=xy1-X0(1);%同时不满足两个判别式时，产生新方向sx2=xy2-X0(2);sx=[sx1sx2];ifm==1S=[e2,sx];%新方向取代上环搜索中函数值下降最大方向m=1elseS=[e1,sx];%新方向取代上环搜索中函数值下降最大方向m=2endelseS=[e1,e2];%满足某个判别式时，保留上环搜索方向enddisp'@@@@第2环搜索方向@@@@'disp(S

7、)%第1环沿新方向搜索x1=xy1;x2=xy2;X=[x1x2];s1=xy1-x01;s2=xy2-x02;S=[s1s2];[xy1,xy2,fyx]=powell(m,x1,x2,X,s1,s2,S);disp'@@@@第2环搜索起始点@@@@'iffyx

8、次项和二次项系数fxs0=60-10*x1-4*x2+x1^2+x2^2-x1*x2;fxs1=-10*s1