黄金分割法求解函数fx-湖南大学.ppt

luluoxx@hnu.cn 第3章常用的一维搜索最优化方法 湖南大学车辆工程 怀旧型缝纫机 FSAE悬架的空间四连杆机构 双横臂独立悬架结构 运动轨迹控制导向作用 弹簧阻尼器：减震吸能作用 悬架的上下轮跳： 控制接地性保证驱动 最优化设计目标：轮距变化最小 减少磨胎 Lunar Exploration Landing 硬着陆 软着陆：CE-II project 转移机构（空间四连杆机构） 月球车 着陆器 曲柄摇杆机构优化设计 讲解内容 导论 搜索区间的确定 黄金分割法的基本原理和算法 插值法的基本原理和算法 算例 总结与讨论 导论 一维问题的作用? 工程中的问题多为复杂问题: 多目标-多维-约束 解决上述问题的根本前提和基础 数值算法的比较与选择(Algorithm) 计算的速度和精度(Convergence) 稳定性(Robust) 适用性(Applicable) 一维搜索方法 一维搜索方法(续) 最优化数学模型 最优步长的计算 在设计空间中选定一个初始设计点X(0) 然后从X(0)这一点出发，按照某一优化方法所规定的原则，确定初始搜索方向S(0)，沿这个方向寻求最优步长α(0)，获得一个目标函数值有所改进的设计点X(1) 然后以X(1)点作为新的始点，再构造此点的新的搜索方向S(1)，求新的最优步长α(1)，求得改进的设计点X(2) 以此类推，重复这种过程，获得目标函数值不断改进的点列：X(0), X(1),…X(k+1)等点，最后可以得到满足所规定的收敛准则或终止准则要求的理论最优点的近似最优点X* 这种寻找最优点的反复过程称为数值迭代方法(定向移动) 一维搜索的基本步骤 探索区间的确定 进退法 外推法 最优步长的计算 解析方法（微分法或变分法） 求函数的一阶导数，甚至二阶导数 只适用于小型（维数n）的最优化问题 数值迭代法(直接法：黄金分割法和平分法，和间接法：二次插值法和三次插值法 ） 不用求导数，只要计算函数值 有收敛速度和收敛精度的问题（求收敛方向S和收敛步长a） 适用于大型的最优化问题（维数n=100) 计算量和存储量问题 初始单峰区间的确定及算法-1 进退法是一种通过比较函数值大小来确定单峰区间的方法 任意给定初始点x1和步长h,算出f(x1) 和x2 =x1 + h点的f(x2)函数值进行比较 初始单峰区间的确定及算法-2 图(a)，f(x1) f(x2),则说明x*x1，将步长增加一倍，取x3 =x2 + 2h 图(b)，f(x1) f(x2),则说明x * x1,则需改变步长符号，取x3 =x2 -h 以此类推，即每跨一步为前一次步长的2倍，直至函数值增加为止 图3-3 (教材pp42) 用进退法 求探索区间[a,b]的 计算程序框图 黄金分割法 基本思想是序列消去原理，即在探索空间内，选取计算点计算函数值并进行比较，消去部分区间，达到缩短探索区间的目的 黄金分割法也称0.618法，是通过对黄金分割点函数值的计算和比较，将初始区间逐次进行缩小，直到满足给定的精度要求，即求得一维极小点的近似解X* 动画 取点规则 黄金分割法的关键是如何不断找出区间内的2个对称点，保证极小点不会丢掉，且收敛快 黄金分割法的基本步骤 1. 确定搜索区间 [a, b], 并计算函数值f (a) and f (b) 2. 令x1=b-0.618(b-a), x2=a+0.618(b-a), 并分别计算函数值f(x1) , f (x2) 3. 比较f(x1) 和f(x2), 直到满足收敛准则 ε≦|b-a| 时终止,则最优化的解为 x*=1/2(b-a) 和f(x*) 4. 如果 f (x1) ＝ f(x2)则舍去右边, 令 b= x2, x2 = x1, f (x2) =f (x1)取点x1 =b-0.618(b-a) 并计算函数值 f (x1) ; 如果 f (x1) f(x2)则舍去左边, 令 a= x1, x1 = x2, f (x1) =f (x2)取点x2 =a+0.618(b-a) 并计算函数值 f (x2) ;然后返回步骤3进行迭代 图3-13 (教材pp52) 不要事先估计n的 黄金分割法 计算程序框图 黄金分割法的特点-1 收敛速度快 区间缩短率为等速均匀的0.618*(b-a) 黄金分割法的关键是如何不断找出区间内的2个对称点，保证最优点不会丢掉，且收敛快 给定收敛精度ε，则计算次数（即区间的缩短次数）N 黄金分割法的特点-2 不必要求f(x)可微，只要利用函数值大小的比较，即可很快地找到最优解 除了第一次缩小区间要计算两个点及其函数值以外，其余每次只要计算一个点及其函数值 可靠性好 算例与结果分析 Maintenance time interval, tp* = 3.7 (mont