工程优化及应用第三章.ppt

求出： 极值的条件： 极值充分条件为： 将极值点方程带入上式 仅取正号 二点三次插值法一般流程： 编写程序应用时建议结合教材框图编写，其更具普适性、鲁棒性。 例　用三次插值法求解下列问题： 教材P56-58的 D.S.C.法、Powell法及其组合法是区间搜索一二次插值法的结合！ P59-P64介绍了有理插值、连分式方法属特殊方法，含教材作者的一些研究成果，大家参阅教材，注意其适应条件，必要时可选用. 方法综述 （1）如目标函数能求二阶导数：用Newton法 收敛快 （2）如目标函数能求一阶导数 ：首先考虑用三次插值法，收敛 较快；二分法、收敛速度慢， 但可靠；二次插值法也可选择. （3）只需计算函数值的方法 ：首先考虑用二次插值法, 收敛快 黄金分割法收敛速度较慢，但可靠 作业：P65 3.1, 3.2, 3.3, 3.7， 各种方法的例题。 定理：设 f：R→R 在[a, b ]上是单峰函数， a≤ x1 x2 ≤b 。 那么 1°若 f (x1)≥ f (x2)，则 x* ∈[x1 , b] ,如左下图 2°若 f (x1) f (x2) ，则 x* ∈[a, x2 ], 如右下图 0.618法（黄金分割法） f(x) x a b (I) 消去[a, x1 ] x* x1 x2 f(x) x a b (II) 消去[x2, b] x* x2 x1 (II) 若f(x1) f(x2), x*∈[a,x2] （I） 若f(x1)≥f(x2)，x*∈[x1,b] 通过上述定理，选二点 x1 x2 , 比较 f (x1 ) 与 f (x2 ) ,可去掉 [a , x1 ] 或者[x2 , b]. 考虑条件： 1°对称原则： x1 – a = b- x2 ……(1) （使“坏”的情况去掉，区间长度不小于“好”的情况） 2°保持缩减比原则 t =(保留的区间长度／原区间长度) 不变。 （使每次保留下来的节点， x1或 x2 ，在下一次的比较中成 为一个相应比例位置的节点 ）。 推导缩减比 t : 如图设第一次保留[a, x2 ] (去掉[x2 , b]), 第二次保留的长度为[a, x1 ], 则 α x1 x2 b 0.618法（黄金分割法） 整理(2) ： x2 = a + t(b-a) x1 = a + t ( x2 -a)= a + (1-t) ( b -a ) 结合(1)式： 故 t≈0.618 注意: 上式有 t 2 = 1-t , 故有 x1 = a+(1-t)(b-a)=a+0.328(b-a) x2 = a+t(b-a)= a+0.618(b-a) 0.618法（黄金分割法） t 2 + t – 1 = 0 当区间的长度充分小时，可将区间中点取做极小点的近似点。 0.618 法----算法 初始[a, b], ε0 x1 = α + （1- t）(b -α ) x2 =α +t (b -α ) b -α ε? STOP; x* =(α+b)/2 yes f（ x1 ）-f（ x2 ）0? No α= x1 , x1 = x2 x2 =α +t ( b –α) yes b= x2 , x2 = x1 x1 = α + （1-t）( b -α ) No α x1 x2 b x2 b α x1 x2 b α x1 优点：不要求函数可微，且每次迭代只需计算一 个函数值，计算量小，程序简单 缺点：收敛速度慢。 黄金分割法（0.618 法）的优缺点 0.618法（黄金分割法） 例 ：试用0.618法求目标函数