最优化方法中关于黄金分割课题设计.docVIP

阐述黄金分割的基本思路及原理 基本思路：黄金分割法适用于[a,b]区间上的任何单股函数求极小值问题,对函数除要求”单峰”外不做其他要求,甚至可以不连续.因此,这种方法的适应面非常广,黄金分割法也是建立在区间消去法原理基础上的试探方法,即在搜索区间[a,b内适当插入两点a1,a2,并计算其函数值。a1,a2将区间分成三段，应用函数的单峰性质，通过函数值大小的比较，删除其中一段，是搜索区间得以缩小。然后再在保留下来的区间上作同样的处理，如此迭代下去，是搜索区间无限缩小，从而得到极小值的数值近似解。 基本原理：在单谷区间内适当插入两点，由此把区间分为三段，然后再通过比较这两点函数值的大小，就可以确定是删去最左端还是最右端，或者同时删去左右两端保留中间段.如此继续下去可将单谷区间无缩小. 基本原理：所谓黄金分割就是将一线段分为两段时，要求整段长与较长段的比值正好等于较长段与较短段的比值（如图所示），即 于是有，解出其正根 . 由此可见长段的长度应为全长的倍，而短段的长度应为全长的倍.因为古代的人们认为按的比率来分割线段时最协调，胜似黄金，故称之为黄金分割. 黄金分割的算法步骤. （1）给定初始区间，精度要求。 令，，并计算与。令。 若，停止，且。否则，当时，转；当时，转。 令，，计算，令，转。 令，计算，令，转。 黄金分割即0.618的算法流程图 黄金分割即0.618的算法的MATLAB程序 function [x,fval,iter]=f618(a,b) iter=0; while abs(b-a)1e-5 iter=iter+1; lambda=a+0.382*(b-a); miu=a+0.618*(b-a); [dy,f1]=fun1(lambda); [dy,f2]=fun1(miu); if f1f2 a=lambda; disp([thenum2str(iter)iteration search area is :[num2str(a),num2str(b)]]) else b=miu; disp([thenum2str(iter)iteration search area is:[num2str(a),num2str(b)]]) end end x=(a+b)/2; [dy,fval]=fun1(x); 根据0.618算法编写程序,求函在计算机上实现黄金分割 数值实验（5、6个表达式也写出） 例1求数在区间[0,1]上的极大值。 解：令 Fun1程序为： function [dy,val]=fun1(x) val=-((sin(x))^6*tan(1-x)*exp(30*x)); dy=-(6*(sin(x))^5*cos(x)*tan(1-x)*exp(30*x)+(sin(x)^6*(1/cos(1-x))^2*(-1)*exp(30*x)+(sin(x))^6*tan(1-x)*exp(30*x)*30); end 运行结果： [x,fval,iter]=f618(0,1) The1iteration search area is:[0.382,1] The2iteration search area is:[0.61808,1] The3iteration search area is:[0.76397,1] The4iteration search area is:[0.85413,1] The5iteration search area is:[0.90985,1] The6iteration search area is:[0.94429,1] The7iteration search area is:[0.94429,0.97872] The8iteration search area is:[0.95744,0.97872] The9iteration search area is:[0.96557,0.97872] The10iteration search area is:[0.96557,0.9737] The11iteration search area is:[0.96867,0.9737] The12iteration search area is:[0.96867,0.9737] The13iteration search area is:[0.96986,0.97178] The14iteration search area is:[0.96986,0.97104] The15iteration search area is:[0.97031,0.971