现代机械优化设计Mtlab编程.docVIP

《现代机械优化设计方法》 程 序 指导老师： 学 生： 学 号： 机电工程学院 2011年12月 一、一维优化问题练习题 1、 黄金分割法 黄金分割法的基本思想：黄金分割法是通过不断单峰区间的长度来搜索极小点的一种有效方法，它是搜索区间比例缩小，通过计算比较函数值，以确定取舍区间。 黄金分割法程序框图 黄金分割法 matlab程序： clc syms t; f=t^2-10*t+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); else if(f1f2) b=d;d=c;f2=f1; c=a+0.382*(b-a);f1=subs(f,c); else a=c;b=d; c=a+0.382*(b-a); d=a+0.618*(b-a); end end end t=0.5*(b+a) f0=subs(f) 结果分析，给出的结果为最优解：，由于定义为单精度类型，故满足了精度要求。Matlab运行结果如下： 二、无约束优化问题 1、 初始点： （1）梯度法 梯度法的基本思想：梯度方向就是函数值变化最大的方向，沿着梯度问题，函数值上升最快，而负梯度方向是函数值下降最快的方向，此方法是沿着负梯度方向进行搜索的无约束方法。 梯度法程序框图 梯度法在matlab程序： (1)function y=fun(x1,x2) y=4*(x1-5)^2+(x2-6)^2;%此定义子函数，并保存为fun.m文件。 (2) clc syms x1 x2 d; 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=subs(fx2); if(double(sqrt(f1^2+f2^2))=0.000000001) n x1=vpa(x1) x2=vpa(x2) f0=vpa(f0) break; else D=fun(x1-d*f1,x2-d*f2); Dd=diff(D,d); dd=solve(Dd); x1=x1-dd*f1; x2=x2-dd*f2; end end %此为梯度法的主函数，并保存为tdf.m文件。 经matlab运行结果如下： (2)牛顿法 与一维搜索讲的牛顿法公式推导一样，只不过推广到n维。牛顿法收敛速度快。 牛顿法程序框图 牛顿法在matlab中的程序如下： clc syms x1 x2; f=4*(x1-5)^2+(x2-6)^2; fx1=diff(f,x1); fx2=diff(f,x2); fx1x1=diff(fx1,x1); fx1x2=diff(fx1,x2); 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(f1^2+f2^2))=0.000000001) n x1=vpa(x1) x2=vpa(x2) f0=vpa(f0) break; else X=[x1 x2]-inv([fx1x1 fx1x2;fx2x1 fx2x2])*[f1 f2]; x1=X(1,1); x2=X(2,1); end end 经matlab运行结果如下： 三、约束优化方法 内点法基本思想 这种方法是将新目标函数定义于可