计算方法 第2章 插值法均差与牛顿插值公式.pptVIP

2019/4/21 1 第二章 插值法 2.3 均差与牛顿插值公式 § 2019/4/21 2 2.3.1 均差及其性质 § ) ( x l j ? ? ? ? ? ? n j i i i j i x x x x 0 ) ( ) ( n j , , 2 , 1 , 0 ? ? 我们知道 , 拉格朗日插值多项式的插值基函数为 形式上太复杂 , 计算量很大 , 并且重复计算也很多 2019/4/21 3 拉格朗日插值公式可看作直线方程两点式的推广，若从 直线方程点斜式 1 0 1 0 0 1 0 ( ) ( ) ( ( ) ) i i i f f P x f x x f f x y x x ? ? ? ? ? ? ? 出发，将它推广到具有 n+1 个插值点的情况，可把插值 多项式表示为 ) ( ) )( ( ) )( ( ) ( ) ( 1 1 0 1 0 2 0 1 0 ? ? ? ? ? ? ? ? ? ? ? ? n n n x x x x x x a x x x x a x x a a x P ? ? 2019/4/21 4 当 0 0 0 1 0 1 1 0 1 2 0 1 2 0 2 2 0 2 1 2 ( ) ( ) ( ) ( ) ( ) ( )( ) n n n P x a f P x a a x x f P x a a x x a x x x x f ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 0 0 1 0 1 1 0 2 0 1 0 2 0 1 0 2 2 1 a f f f a x x f f f f x x x x a x x ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 依次可得到 3 4 , , , n a a a 。为写出系数的一般表达式， 现引入差商（均差）定义。 2019/4/21 5 一、差商 ( 均差 ) 定义 2. n i f x x f i i , , 1 , 0 , ) ( ? ? 处的函数值为 在互异的节点 设 称 ) 0 ( ) ( ) ( ] , [ 0 0 0 ? ? ? ? k x x x f x f x x f k k k ) ( , ) ( 0 差 商 一 阶 均 差 关 于 节 点 为 k x x x f 2019/4/21 6 2019/4/21 7 ] , , , , [ 1 1 0 k k x x x x f ? ? ? ? ? ? ? ? ? ? ? k j k j j j j j j j x x x x x x x x x f 0 1 1 0 ) ( ) )( ( ) ( ) ( ? ? 二、均差具有如下性质： 2019/4/21 8 ( ) ( ) [ , ] 0 1 0 1 0 1 f x f x f x x x x ? ? ? 0 1 0 1 1 0 f f x x x x ? ? ? ? [ , ] [ , ] [ , , ] 0 1 0 2 0 1 2 1 2 f x x f x x f x x x x x ? ? ? ( ) ( ) 0 0 1 2 1 2 0 1 1 0 1 2 0 2 2 0 f f f f 1 1 x x x x x x x x x x x x ? ? ? ? ? ? ? ? ? ? ( )( ) ( )( ) ( )( ) 0 1 2 0 1 0 2 1 0 1 2 2 0 2 1 f f f x x x x x x x x x x x x ? ? ? ? ? ? ? ? ? 例 2019/4/21 9 这个性质也表明差商与节点的排列顺序无关 ( 差商的对称性 ) 。即 0 1 1 0 2 1 2 0 [ , , , ] [ , , , , ] [ , , , , ] k k k f x x x f x x x x f x x x x ? ? ? 2019/4/21 10 性质 3 ：若 f(x) 在 [a,b] 上存在 n 阶导数，且节点 则 n 阶均差与导数关系如下： 0 , [ , ], n x x a b ? 2019/4/21 11 ) ( ) ( ) ( ) ( ) ( ) ( 4 4 3 3 2 2 1 1 0 0 x f x x f x x f x x f x x f x x f x k k 四阶差商 三阶差商 二阶差商 一阶差商 三、均差的计算方法 ( 表格法 ): ] , [ 1 0 x x f ] , [ 2 1 x x f ] , [ 3 2 x x f ]