第四节黄金分割_二次差值法汇编.pptVIP

　构造的二次插值函数曲线通过原函数上的三个点: 解得系数 可求得 二次插值法程序框图 开 始 计算 ， Y N 给定 , 缩短区间名称置换 结 束 a1 a2 ap a3 a y 0 a* y1 y2 yp y3 a1 a2 ap a3 a 0 a* y y1 y2 yp y3 二次插值法程序编制换名方法(1) 1) a2ap y2≥yp y2yp 二次插值法程序编制换名方法(2) 2) a2ap y2 ≥ yp y2yp （1）二次插值法只要求f(x)连续，不要求其一阶可微； （2）收敛速度比黄金分割法快，但可靠性不如黄金分割 法好，程序也较长。 特点： 二次插值法（抛物线法） 例 3-3 用二次插值法求函数f(x)=3x3-4x+2的极小点， 给定 x0=0, ε=0.2。 解 1）确定初始区间:[a,b]=[0,2] 2）用二次插值法逼近极小点 相邻三点的函数值: x1=0, x2=1, x3=2; f1=2, f2=1, f3=18. 代入公式： xp＝0.555, fp=0.292 在新区间，相邻三点的函数值: x1=0, x2=0.555, x3=1; f1=2, f2=0.292, f3=1. xp=0.607, fp=0.243 由于fpf2, xpx2, 新区间[a,b]=[x2, b]=[0.555,1] |x2-xp |=|0.555-0.607|=0.0520.2, 迭代终止。 x*= 0.607, f*=0.243 由于fpf2, xp x2, 新区间[a,b]=[a,x2]=[0,1] |x2-xp |=1-0.555=0.4450.2, 应继续迭代。 本章结束 Thank You！ Y=x3-2x+5 进退法确定搜索区； 黄金分割法、二次插值法求解极小值。 * 第三章 一维搜索方法 3.3 一维搜索的试探法 3.1 搜索区间的确定 3.2 区间消去法原理 3.4 一维搜索的插值法 f(b1) f(a1) f(b1) f(a1) f(b1) a f(a1) 搜索区间确定之后，采用区间消去法逐步缩短搜索区间，找到极小点的数值近似解。 假定在搜索区间内[a,b] 任取两点a1、b1,且a1b1 f1＝f（a1）， f2＝f（b1） 一、基本思想 a1 a1 a1 b1 b a a b b b1 b1 §3.2 区间消去法原理 f1f2 f1f2 f1=f2 f(x) f(x) f(x) f(a1) f(b1) f(a1) f(b1) f(a1) f(b1) a1 a1 a1 b1 b a a b a b b1 b1 (1)f1f2, 新区间为[a,b1]; (2)f1f2, 新区间为[a1,b]; (3)如f1=f2, 新区间为[a1,b1] 对于上述缩短后的新区间，可在其内再取一个新点，然后将此点和该区间内剩下的那一点进行函数值大小的比较，以再次按照上述方法，进一步缩短区间，这样不断进行下去，直到所保留的区间缩小到给定的误差范围内，而得到近似最优解。 新区间为 a b a2 a1 a2 a a1 f(a1) f(a2) f(a2) f(a1) b 一、适用范围 适用于[a,b]区间上的任何单谷函数求极小值问题。对函数除要求“单谷”外不作其他要求，甚至可以不连续。适应面相当广。基本原理为区间消去法 §3.3 黄金分割法 黄金分割法插入两点： 黄金分割法要求在区间[a,b]中插入的两点位置是对称的，如图所示，ax1=x2b。设区间长为L，插入的两个点把区间分成较长的一段?和较短的一段(L – ?)，如图所示， ax2=bx1=? ， ax1=bx2=L-? ，这样，无论删去那一段，保留的区间长度总是? ，在每次迭代中，整个区间的长度L 与较长一段长度的比等于较长一段长度与较短长度的比。 α L L－α α L L－α 黄金分割法采用固定的缩短率0.618 λL/ (1-λ)L b x1 x2 ? ? 取点公式 原始区间 [a, b], 第一次内点为x1,x2 x1 = a x2 ? a +(1-λ)( b- a) f1=f(x1) x2 a + = λ( b- a) f2=f(x2) ? (1) 若 f1f2 f2 =f(x2) a b a= x1, b=b, x2 第k次迭代