机械优化设计黄金分割法实验报告资料.docVIP

实验1 一维搜索(黄金分割)法实验报告 实验日期 201X 年 12 月 11 日 报告日期 201X 年 12 月 12 日 班级 XXXXX级机自XXXX班 姓名 XXXXXX 学号 XXXXXXXXXXXXXXX 1、实验目的 了解黄金分割法的基本原理； 熟悉matlab程序使用方法； 学习上机调试、运行所编写的程序。 2、黄金分割法原理 该法适用于[a,b]区间上单谷函数极小值问题。在搜索区间[a,b]内按照0.618比例加入两点α1,α2，并计算其函数值。α1,α2将区间分成三段，然后利用区间消去法，通过比较函数值大小，删除其中一段，使搜索区间缩短，在保留区间进行同样处理，直到搜索区间缩小到指定精度为止。 3、编制MATLAB优化程序 编写函数文件，并命名为fx.m保存，程序代码如下： function f=fx(w) %f=w^2-10*w+36; %f=w^4-5*w^3+4*w^2-6*w+60; %f=((w+1)^4)*((w-2)^2); 注：上述“%”后面分别为要求解的三个方程，求解该方程式把相应方程式前面的“%”删除，点击保存，并运行下面的hjf.m文件，输入相应的初始步长h0、初始点x0、收敛法则epsilan的值 编写进退法程序文件，命名为ab1.m保存，程序代码如下： function [a,b]=ab1(h0,x0) h=h0;x1=x0;f1=fx(x1); x2=x1+h;f2=fx(x2); if f2f1 h=-h;x3=x1;f3=f1;x1=x2;f1=f2;x2=x3;f2=f3; end x3=x2+h;f3=fx(x3); while f2=f3 x1=x2;f1=f2;x2=x3;f2=f3; x3=x2+2*h;f3=fx(x3); end if h0 a=x3;b=x1; else a=x1;b=x3; end 编写黄金分割法程序文件，命名为hjfgf.m保存，程序代码如下： function hjf clear h1 = input(h0=?);x1=input(x0=?); epsilan=input(epsilan=?); [a,b]=ab1(h1,x1); x1=a+0.382*(b-a);f1=fx(x1); x2=a+0.618*(b-a);f2=fx(x2); while abs(b-a)epsilan if f1f2 a=x1;x1=x2;f1=f2; x2=a+0.618*(b-a);f2=fx(x2); else b=x2;x2=x1;f2=f1; x1=a+0.382*(b-a);f1=fx(x1); end end xm=(a+b)/2; [Optimal result:,blanks(3),xm=[,... num2str(xm),],blanks(6),fm=,num2str(fx(xm))] 4、实验结果 结果: = 5.0001 11 （h0=0.3、x0=0、epsilan=0.001） 结果: = 3.2795 22.659 （h0=0.3、x0=0、epsilan=0.001） 结果: =__-0.99989__ __1.4253e-015__ （h0=0.3、x0=0、epsilan=0.001） 励志美文美句摘抄 　　1、不要放弃自己就是真正的坚强，虚心就是坚强，努力就是坚强，从头再来就是坚强，正直就是坚强，学会坚强之前要学会如何爱惜自己。 　　2、人生，就没有，永远的悲痛；也没有，永远的欢欣。能使我们坚强的，往往不是顺境，而是逆境；能让我们醒悟的，往往不是高兴，而是伤心。学会忍受，懂得艰辛，于曲折中前进。 　　3、人都说比天空和大地更远的距离是人与人的距离，因为人心里都会藏匿太多的猜忌和戒备，想要快乐就甩开生命中这些过于沉重，却又不必要的行李吧，生命中有爱就足够了。试着给周围你所熟识的还有你还陌生的人一个真诚无惕的微笑吧，它可以触摸到他人的心灵，微笑是有感染力有连带性的，它会无声的渗透进每个易感的心灵，更会让更多的心灵为之感动，心中有爱就会快乐，就会让微笑发自心底，灿烂在脸上。 　　4、千万不要因为自己已经到了结婚年龄而草率结婚。想结婚，就要找一个能和你心心相印相辅相携的伴侣。不要因为放纵和游戏而恋爱，不要因为恋爱而影响工作和事业，更不要因一桩草率而失败的婚姻而使人生受阻。感情用事往往会因小失大。 　　5、你要从现在开始，微笑着面对生活，不要抱怨生活给了你太多的磨难，不要抱怨生活中有太