数值分析(研究生)第二章插值与拟合一.ppt

注：由上例可知，牛顿插值具有继承性，和拉格朗日插值相比节约了计算量. 另外，在估计余项时，牛顿插值对于只给出了离散点函数值的函数或是导数不存在的函数都是适用的. * 上页 下页 返回 第二章 插值、逼近与拟合（一） 第二节 拉格朗日插值 第三节 牛顿插值 第一节 实际问题的导入 插值法是广泛应用于理论研究和工程实际的重要数值方法 . 插值法的定义： 设函数 y = f (x) 在区间[a, b]上有定义，且已知在点a ≤ x0 x1 … xn ≤ b 上的值 y0, y1 , … , yn,若存在一简单函数 P(x),使 P(xi ) = yi (i = 0,1,…,n) 成立，就称P(x) 为 f (x) 的插值函数,点 x0, x1 , … , xn 称为插值节点,区间[a, b]称为插值区间，求插值函数P(x)的方法称为插值法. 若P(x)为多项式时，就称为多项式插值.同理，有有理分式插值、三角插值等. §1 实际问题的导入 如飞机、船舶的设计、制造，桥梁的设计、建造等 . 插值法的几何意义： 求曲线 y = P (x) ，使其通过给定的n+1个点(xi , yi )， i = 0,1,…,n，并用它近似已知曲线 y = f (x) . 定理1 在次数不超过n次的多项式集合Hn中，满足插值条件 P(xi)=yi (i=0,1,2,……,n) 的插值多项式存在且唯一. 一、线性插值与抛物插值 1、线性插值 §2 拉格朗日插值 2、抛物 插值 二、拉格朗日插值多项式 设节点 在[a , b]内存在, 考察截断误差 ，且 f 满足条件 , 罗尔定理： 若 充分光滑， ，则 存在 使得 。 推广：若 使得 使得 存在 使得 Rn(x) 至少有 个根 n+1 ? = - = n i i n x x x K x R 0 ) ( ) ( ) ( 任意固定 x ? xi (i = 0, …, n), 考察 ? = - ? = n i i x t x K t Rn t 0 ) ( ) ( ) ( ) ( j ?(x)有 n+2 个不同的根 x0 … xn x ! ) 1 ( ) ( ) ( ) 1 ( + - + n x K R x n n x 注意这里是对 t 求导 = + - - + + ! ) 1 )( ( ) ( ) ( ) 1 ( ) 1 ( n x K L f x n n x n x x ! ) 1 ( ) ( ) ( ) 1 ( + = + n f x K x n x 三、插值余项与误差估计 (也称为插值多项式的余项). 当 f(x) 为任一个次数? n 的多项式时， , 可知 ，即插值多项式对于次数? n 的多项式是精确的. 例1 求经过A(0, 1)，B(1, 2)，C(2, 3)三个插值点的插值多项式. 解：三个插值节点及对应的函数值为 注：本题还有更适用的方法. 例3 已知 分别利用 sin x 的1次、2次拉格朗日插值计算 sin 50?， 并估计误差. 解： n = 1 分别利用x0, x1 以及 x1, x2 计算 ?利用 这里 而 ?sin 50? = 0.7660444… ) 18 5 ( 50 sin 1 0 ? ? p L 0.77614 外推 的实际误差 ? ?0.01001 ?利用 sin 50? ? 0.76008, 内插 的实际误差 ? 0.00596 内插通常优于外推. 选择要计算的 x 所在的区间的端点，插值效果较好. n = 2 ) 18 5 ( 50 sin 2 0 ? ? p L 0.76543 ?sin 50? = 0.7660444… 2次插值的实际误差 ? 0.00061 高次插值通常优于低次插值 但绝对不是次数越高就越好，嘿嘿…… 四、拉格朗日插值公式的优缺点 一、均差及其性质 §2 均差与牛顿插值公式 1 2 … … … … n?1 1 + (x ? x0) ? 2 + … … + (x ? x0)…(x ? xn?1) ? n?1 Nn(x) Rn(x)