最优化方法4第四章剖析.ppt

搜索区间及其确定 对分法 Newton切线法 黄金分割法 抛物线法 三， 即使曲线比较正常，在[a,b]中或者凹或者凸，初始点的选取必须适当，在图4.10(a)的情况，曲线凹，应选b作为初始点；而在图4.10(b)的情况，曲线凸，应选a作为初始点； 单变量函数求最小值的标准形式为 s.t x = fminbnd(fun,x1,x2,options) % options为指定优化参数选项 [x,fval] = fminbnd(…) % fval为目标函数的最小值 [x,fval,exitflag] = fminbnd(…) %xitflag为终止迭代的条件 [x,fval,exitflag,output] = fminbnd(…) % output为优化信息 说明 若参数exitflag0，表示函数收敛于x，若exitflag=0，表示超过函数估计值或迭代的最大数字，exitflag0表示函数不收敛于x；若参数output=iterations表示迭代次数，output=funccount表示函数赋值次数，output=algorithm表示所使用的算法。 例1 计算下面函数在区间(0，1)内的最小值。 §4.4 黄金分割法 一、黄金分割法基本原理 要介绍黄金分割法有必要回顾一下古老的黄金分割问题．所谓黄金分割就是将一线段分为二段的方法．这样分后，要求整段长L与较长段x的比值正好等于较长段x与较短段 的比值（如图4.11所示），即 ． 于是有 ，解出其正根 ． 由此可见长段的长度应为全长的0.618倍，而短段的长度应为全长的0.382倍．因为古代的人们认为按0.618的比率来分割线段是最协调，胜似黄金，故称之为黄金分割． 图4.11 用黄金分割法进行一维搜索，其基本思想是在单谷区间内适当插入两点，由此把区间分为三段，然后再通过比较这两点函数值大小，就可以确定是删去最左段还是最右段，或者同时删去左右两段保留中间段．如此继续下去可将单谷区间无限缩小． 二、黄金分割法迭代步骤 现在提出一个问题，就在 上如何选取二点 使得迭代次数最小而区间缩短最快？要解决这个问题，人们想到对区间 选二点等价于将区间长度 进行黄金分割， 也就是将第一个搜索点 取在 的0.618处，第二个搜索点 取成 的对称点即 的0.382处（如图4.12所示），亦即 ， ． 计算 与 的值，并根据 与 的值的大小关系分情况讨论： （1） 若 ，说明 是好点，于是把区间 划掉，保留 ，则 内有一保留点 ，置新的区间 ； （2）若 ，说明 是好点，于是应将 划掉，保留 ，则 内有保留点 ，置新的区间 ． （3）若 ，则应具体分析，看极小点可能在哪一边再决定取舍，在一般情况下，可同时划掉 ，和 ，仅保留中间的 置新的区间 ． 接下来是在留下的区间 内找好点．重复上面的步骤，直到搜索区间 小于给定的允许误差 为止．这样就得到黄金分割法迭代算法： 图4.12 已知 ，常数 0.382，终止限 ． (1)确定 的初始搜索区间 ． (2)计算 ． (3)计算 ． (4) 若