数值教案解析.ppt

除了上述算子外，常用的还有不变算子和移位算子定义如下： ?32? 结束 （1）各阶差分可以表示成函数值的线性组合，如： （2）各阶差分之间可以互化，如： （3）可用差分表示差商，如： 差分性质： ?33? 结束 （4）可用各阶差分表示函数值。例如可用 向前差分表示 一般地有： ?34? 结束 例5 给出函数y=?(x)的函数表 解 差商表如下 3 2 1 0 i 21 17 3 5 ?(xi) 2 1 -1 -2 xi 写出函数y=?(x)的差商表. -2 -1 1 2 xi -1 3 -1 -2 7 4 5 3 17 21 0 1 2 3 三阶差商 二阶差商 一阶差商 ?(xi) i ?35? 结束 §3.2 Newton插值多项式及其余项 ?(x)=?(x0)+(x-x0)?[x0,x] 由差商的定义可得 所以有 ?[x0,x]=?[x0,x1]+(x-x1)?[x0,x1,x] ?[x0,x1,…,xn-1,x]=?[x0,x1,…,xn] +(x-xn)?[x0,x1,…,xn,x] ? ?[x0,x1,x]=?[x0,x1,x2]+(x-x2)?[x0,x1,x2,x] ?36? 结束 ?(x)=?(x0)+(x-x0)?[x0,x1] +(x-x0)(x-x1)?[x0,x1,x2]+… +(x-x0)(x-x1)…(x-xn-1)?[x0,x1,…,xn] +(x-x0)(x-x1)…(x-xn)?[x0,x1,…,xn,x] 记 Nn(x)=?(x0)+(x-x0)?[x0,x1] (2.8) +(x-x0)(x-x1)?[x0,x1,x2]+… +(x-x0)(x-x1)…(x-xn-1)?[x0,x1,…,xn] Rn(x)=(x-x0)(x-x1)…(x-xn)?[x0,x1,…,xn,x] 则有 ?(x)=Nn(x)+Rn(x) ?37? 结束 而且Nn(x)是n次多项式,且满足 称Nn(x)为n次Newton插值多项式,称Rn(x)为n次Newton插值余项.由插值多项式的唯一性有 Nn(xi)=?(xi)(i=0,1,…,n) Nk+1(x)=Nk(x)+?k+1(x)?[x0,x1,…,xk+1] k=1,2,…,n-1 由(2.8)式易见 ?38? 结束 解 由例4的差商表知?[x0,x1]=-2,?[x0,x1,x2]=3, ?[x0,x1,x2,x3]=-1,于是有 N1(x)=5-2(x+2)=1-2x N2(x)=1-2x+3(x+2)(x+1)=3x2+7x+7 N3(x)=3x2+7x+7-(x+2)(x+1)(x-1)=-x3+x2+8x+9 例6 对例5中的 ?(x),求节点为 x0,x1 的一次插值 x0,x1,x2的二次插值和x0,x1,x2,x3的三次插值 多项式. ?39? 结束 §4 Hermite 插值多项式 要求函数值重合，而且要求若干阶导数也重合。 即：要求插值函数 P (x) 满足 p (xi) = f (xi), P’ (xi) = f ’ (xi), …, P(m) (xi) = f (m) (xi). 在实际问题中，对所构造的插值多项式，不仅 把此类插值多项式称为埃米尔特（Hermite） 插值多项式或称带导数的插值多项式，记为H (x)。 注： ? N 个条件可以确定 阶多项式。 ?要求在1个节点 x0 处直到m0 阶导数都重合的插值多项式即为Taylor多项式 其余项为 N ? 1 例：设 x0 ? x1 ? x2, 已知 f(x0)、 f(x1)、 f(x2) 和 f ’(x1), 求多项式 P(x) 模仿 Lagrange 多项式的思想，设 解：首先，P 的阶数 = 3 ?  + = 2 1 3 ) ( ) ( ) ( ) ( ) ( = 0 i i i x h x1 f ’ x h x f x P ? h0(x) 有根 x1, x2，且 h0’(x1) = 0 ? x1 是重根。 ) ( ) ( ) ( 2 2 1 0 0 x x x x C x h - - = 又: h0(x0) = 1 ? C0 h2(x) h1(x) 有根 x0, x2 ? ) )( )( ( ) ( 2 0 1 x x x x B Ax x h - - + = 由余下条件 h1(x1) = 1 和 h1’(x1) = 0 可解。 与h0(x) 完全