计算物理lecture1102.pptVIP

关于Langrange插值的几点说明 仅与已知数据 有关，与 的原来形式无关，但余式与 密切相关。 若 本身是一个不超过n次多项式，则 在实际应用中 ，常是等距节点情况，即 这里h0为常数，称为步长，这时Newton插值公式就可以简化，为此我们引入差分概念。 例题 Lagrange插值公式所求得L(x)保证了节点处的函数值相等，也就是保证了函数的连续性，但不少实际问题还需要插值得光滑，也就是还要求在节点处的导数值也相等，导数的阶数越高则光滑度越高。现代的仿生学就是一个典型的例子。在设计交通工具的外形，就是参照海豚的标本上已知点及已知点的导数，做插值在计算机上模拟海豚的外形制成飞机、汽车等外形。 构造 Hermite 插值多项式H(x) 分段线性插值 分段线性插值 Hermite 插值 三点三次Hermite 插值 插值举例 插值举例 插商与差分的关系 在等距节点条件下有： 0 2 2 0 1 0 2 1 0 2 1 2 1 0 ! 2 1 2 1 1 ] , [ ] , [ ] , , [ f h h f h f h x x x x f x x f x x x f D = D - D = - - = 0 1 0 ! 1 ] ,..., , [ f h n x x x f n n n D = 0 2 2 2 1 0 ! 2 1 ] , , [ f h x x x f D = 一般有 ] ,... , , [ ) )( )...( )( ( ] ,... , [ ) )...( )( ( ... ] , , [ ) )( ( ] , [ ) ( ) ( ) ( 1 0 1 1 0 1 1 0 2 1 0 2 1 0 1 0 1 0 0 0 n n n n n x x x x f x x x x x x x x x x x f x x x x x x x x x f x x x x x x f x x x f x f - - - - + - - - + + - - + - + = - - - ) , ( ), ( ) )...( 1 ( )! 1 ( ) ( ! ) 1 )...( 1 ( ... ! 2 ) 1 ( ! 1 ) ( ) ( 0 0 ) 1 ( 1 0 0 2 0 0 0 h x x f n t t t n h x R f n n t t t f t t f t f th x N x N n n n n n n + ? - - + = D + - - + + D - + D + = + = + + x x , 0 Newton th x x + = 式 插值公式和余式具有形 则 若令 ] ,... , , [ ) )( )...( )( ( ] ,... , [ ) )...( )( ( ... ] , , [ ) )( ( ] , [ ) ( ) ( ) ( 1 0 1 1 0 1 1 0 2 1 0 2 1 0 1 0 1 0 0 0 n n n n n x x x x f x x x x x x x x x x x f x x x x x x x x x f x x x x x x f x x x f x f - - - - + - - - + + - - + - + = - - - 上面介绍的差分叫向前差分。鉴于计算实践的需要，我们再介绍向后差分概念。 设函数 f(x) 在等距节点 上值为 fi =f(xi )，则 为f(x)在点 处的一阶向后差分 1 - - = ? i i i f f f 为f(x)在点 xi 处的n阶向后差分。 i i i f f f - = D + 1 插商与向后差分的关系 用后插表示N(x) 用后插表示N(x) 对于给定的等距结点的数据，我们可以灵活运用插值余项极小化原则，给出适应具体要求的插值公式。一般说来，在左端点附近进行插值，宜用Newton向前插值公式；在右端点附近插值，宜用Newton向后插值公式。 Hermite 插值 所谓 Hermite 插值是插值函数与被插值函数 不仅在节点处函数值相等, 而且导数相等. 设函数 f ( x ) 在区间 [ a , b ] 上具有一阶连续导数, 如果有不超过 2n+1 次多项式 H2n+1(x) 满足 ) ,... 2 , 1 , 0 (i n = ) ( 1 2 x f H i n = + ), ( ) ( 1 2 x f x H i i n = + 则称 为 f ( x