2022年现代机械优化设计Matlab编程资料 .pdf

《2022年现代机械优化设计Matlab编程资料 .pdf》由会员分享，可在线阅读，更多相关《2022年现代机械优化设计Matlab编程资料 .pdf（9页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、现代机械优化设计方法程序指导老师：学生：学号：机电工程学院2011年 12 月名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 9 页 - - - - - - - - - 1 一、一维优化问题练习题1、2min( )1036f ttt黄金分割法黄金分割法的基本思想：黄金分割法是通过不断单峰区间的长度来搜索极小点的一种有效方法， 它是搜索区间比例缩小， 通过计算比较函数值， 以确定取舍区间。黄金分割法程序框图黄金分割法matlab程序：clcsyms t ;f=t2-10*t

2、+36;a=2;b=10;c=a+0.382*(b-a);d=a+0.618*(b-a);f1=subs(f,c);f2=subs(f,d);while(abs(b-a)=0.01)if(f1f2) a=c;c=d;f1=f2; d=a+0.618*(b-a);f2=subs(f,d);elseif(f1f2) b=d;d=c;f2=f1;名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 9 页 - - - - - - - - - 2 c=a+0.382*(b-a);f1=

3、subs(f,c);else a=c;b=d; c=a+0.382*(b-a); d=a+0.618*(b-a);endendendt=0.5*(b+a)f0=subs(f) 结果分析， 给出的结果为最优解：，由于定义为单精度类型， 故满足了精度要求。Matlab 运行结果如下：二、无约束优化问题1、2212min()4(5)(6)f Xxx初始点：(0)(0)89 ,()45TXf X（1）梯度法梯度法的基本思想： 梯度方向就是函数值变化最大的方向，沿着梯度问题， 函数值上升最快， 而负梯度方向是函数值下降最快的方向，此方法是沿着负梯度方向进行搜索的无约束方法。梯度法程序框图名师资料总结 -

4、 - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 9 页 - - - - - - - - - 3 梯度法在 matlab 程序：(1)function y=fun(x1,x2)y=4*(x1-5)2+(x2-6)2;%此定义子函数，并保存为fun.m文件。(2) clcsyms x1x2d;f=4*(x1-5)2+(x2-6)2;fx1=diff(f,x1);fx2=diff(f,x2);x1=8;x2=9;for n=1:100 f0=subs(f); f1=subs(fx1); f2=

5、subs(fx2);if(double(sqrt(f12+f22)=0.000000001) n x1=vpa(x1) x2=vpa(x2) f0=vpa(f0)break;else名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 9 页 - - - - - - - - - 4 D=fun(x1-d*f1,x2-d*f2); Dd=diff(D,d); dd=solve(Dd); x1=x1-dd*f1; x2=x2-dd*f2;endend%此为梯度法的主函数，并保存为t

6、df.m文件。经matlab 运行结果如下：(2)牛顿法与一维搜索讲的牛顿法公式推导一样，只不过推广到 n 维。牛顿法收敛速度快。牛顿法程序框图名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 9 页 - - - - - - - - - 5 牛顿法在 matlab 中的程序如下：clcsyms x1x2 ;f=4*(x1-5)2+(x2-6)2;fx1=diff(f,x1);fx2=diff(f,x2);fx1x1=diff(fx1,x1);fx1x2=diff(fx1,x

7、2);fx2x1=diff(fx2,x1);fx2x2=diff(fx2,x2);x1=8;x2=9;for n=1:100 f0=subs(f); f1=subs(fx1); f2=subs(fx2);if(double(sqrt(f12+f22)=0.000000001) n x1=vpa(x1) x2=vpa(x2) f0=vpa(f0)break;else名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 9 页 - - - - - - - - - 6 X=x1 x2

8、-inv(fx1x1 fx1x2;fx2x1 fx2x2)*f1 f2; x1=X(1,1); x2=X(2,1);endend经matlab运行结果如下：三、约束优化方法内点法基本思想这种方法是将新目标函数定义于可行域内，迭代过程均在可行域内进行， 逐步逼近最优点。内点法只能用来求解具有不等式约束的优化问题。1）在可行域内选一个初始点内点X0，初始点应选择一个离约束边界较远的可行点。2）选取适当的罚因子初值r0，降低系数 c（0c1)，计算精度。内点法程序框图名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - -

9、 - - - - - 第 7 页，共 9 页 - - - - - - - - - 7 Matlab程序：clc;count=0;n=100;r=1;c=0.5;flag=0;x0=3,3;esp=1.0e-5;while(flag=0)&(countn) syms x1x2 ; f=(x1-2)2+(x2-1)2; g1=x12-x2; g2=x1+x2-2; F=(x1-2)2+(x2-1)2+r*(log(g1)+log(g2); f1=diff(F,x1); f2=diff(F,x2); s=solve(f1,f2,x1, x2); x1=s.x1 s.x2; norm=subs(sqr

10、t(x1(1)-x0(1)2+(x1(2)-x0(2)2);if(normesp) flag=1;break;end名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 8 页，共 9 页 - - - - - - - - - 8 x0=x1; r=c*r; count=count+1;endx1=x0(1)x2=x0(2)f=subs(f)count程序的运行结果：名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 9 页，共 9 页 - - - - - - - - -