第三讲一维优化方法.ppt

黄金分割法可使相邻两次搜索区间都具有相同的缩短率0.618。因而内分点的取点规则为： x1=a+ 0.382(b-a) x2=a+ 0.618(b-a) 终止原则： 最优解： k 为缩短区间的次数 第二十九页，共六十二页，2022年，8月28日 黄金分割法的搜索过程 1、给出初始搜索区间[a, b]及收敛精度 ?； 2、按坐标点计算公式计算，并计算相应的函数值； 3、缩短搜索区间； 4、检查是否满足收敛条件； 5、若满足收敛条件，则取最后两点的平均值作为极小点的近似解。 第三十页，共六十二页，2022年，8月28日 黄金分割法的流程图 第三十一页，共六十二页，2022年，8月28日 例题 例题3.3 试用黄金分割法求目标函数 f(x)=x2-6x+9 的最优解。给定初始区间[1，7]，收敛精度ε=0.4。 解：第一次区间缩短计算过程： 计算两内点及对应函数值： x1=a+0.382(b-a)=3.292 y1=f(x1)=0.085264 x2=a+0.618(b-a)=4.708 y2=f(x2)=2.917264 作数值比较，可见y1y2,再做置换： 用终止准则判断: a(1) ← a=1 b(1) ← x2=4.708 第三十二页，共六十二页，2022年，8月28日 为第二次区间缩短做准备，做置换： x2←x1=3.292 ， y2←y1=0.085264 各次缩短区间的计算数据见下页表。第六次区间缩短的端点 满足精度要求，终止计算。取最优解为 第三十三页，共六十二页，2022年，8月28日 区间缩短次数 a b x1 x2 y1 y2 原区间 1 2 3 4 5 6 1 1 2.416456 2.416456 2.416456 2.750917 2.750917 7 4.708 4.708 3.832630 3.292 3.292 3.085305 3.292 2.416456 3.292 2.957434 2.750917 2.957434 2.878651 4.708 3.292 3.832630 3.292 2.957434 3.085305 2.957434 0.085264 0.340524 0085264 0.001812 0.062044 0.001812 0.014725 2.917264 0.085264 0.693273 0.085264 0.001812 0.007277 0.001812 黄金分割法例题计算数据 第三十四页，共六十二页，2022年，8月28日 格点法 则最优解为：x* ? xm ， y* ? ym 若不能满足精度要求，把当前区间作为初始搜索区间，重复上述步骤直至满足精度为止 设一维函数为 f (x)，搜索区间为[a, b]，一维收敛精度为?。在区间[a, b]的内部取 n 个等分点： x1 , x2 , … , xn。区间被分为 n+1 等分，各分点坐标为: 对应各点的函数值为 y1 , y2 , … , yn。比较其大小，取最小者： 格点法的原理 ym=min { yk , k=1,2,…,n} 则在区间[xm-1 , xm+1]内必包含极小点，取[xm-1 , xm+1]为缩短后新区间，若新区间满足收敛条件： xm+1- x m-1? ? 第三十五页，共六十二页，2022年，8月28日 格点法的区间缩短 新区间 y1 ym-1 ym ym+1 yn x1 a xm-1 xm xm+1 xn b 第三十六页，共六十二页，2022年，8月28日 格点法流程图 第三十七页，共六十二页，2022年，8月28日 例题 例题3.2：用格点法求一维目标函数 f(x) = 4x2-12x+10 的最优解。给定搜索区间[a，b]为[1，2.2]，迭代精度ε=0.2，内分点数 n = 4。 解：计算区间端点的函数值 f(a)=2 f(b)=2.96 由上式确定四个内分点的位置，并计算其函数值，计算结果如表所示。其中最小的函数值为： 对应的点： 第三十八页，共六十二页，2022年，8月28日 区间缩短次数 xk yk xm a b b-a 第一次 1.24 1.48 1.72 1.96 1.2704 1.0016 1.1936 1.846 1.48 1.24 1.72 0.48 第二次 1.336 1.432 1.528 1.624 1.107584 1.018496 1.003136 1.061504 1.528 1.432 1.624 0.192 第三十九页，共六十二页，2022年，8月28日 新区间的端点为： 新区间的长度为：