第五章节_数值插值方法.pptVIP

第五章 插值方法 第五章 插值方法 插值的基本概念 Lagrange插值 分段低次插值 均差与Newton插值 Hermite插值 三次样条插值 Lagrange插值基函数构造法 给定插值节点 x0,x1， y0=f(x0)，y1=f(x1). 求线性插值多项式L1 (x)=a0+ a1x，使满足: L1(x0)=y0 , L1(x1)=y1. 由两点式可以看出， L1 (x)是由两个线性函数 l0(x0)=1 , l0(x1)=0 , l0(x2)=0. l1(x0)=0 , l1(x1)=1 , l1(x2)=0. l2(x0)=0 , l2(x1)=0 , l2(x2)=1. 满足上式的插值基函数很容易求出。如求l0(x), 因x1, x2 为其零点，故可表为 显然 L(x)=l0(x)y0+l1(x)y1+l2(x)y2 满足条件 L2(xj)=yj (j=0,1,2) 取x0=4,y0=2，x1=9, y1=3 ，x2=16, y2=4. 设有n+1个互异节点x0 x1…xn，且 yi=f(xi) (i=0,1,2…,n) 构造Ln (x)，使 Ln (xj)= yj (j = 0,1,2,…,n) 由n=1,2时的讨论可得 定理 设 f(x)在[a,b]上具有n阶连续导数, 且 f (n+1)(x) 存在,节点a ≤ x0 x1…xn≤b， Ln (x)是满足条件Ln (xj)= yj (j = 0,1,2,…,n)的插值多项式，则对任何x?[a,b]，插值余项 当 f(x) 是n次的多项式时, Ln(x)= f(x)。即n次多项式的n次插值函数即为该n次多项式本身。 取xk=-5+k 计算: f(xk) (k=0,1,…,10) 构造L10(x). 取:tk=-5+0.05k (k=0,1,…,200),计算: L10(tk) 这说明截断误差很小。 设节点x0 x1…xn，分段插值函数Hn (x)在两个相邻节点构成的小区间[xj , xj+1 ] (j=0,1,…n-1)上满足条件： 其中 (1)S(x)在每一个子区间[xj-1 , xj ] ( j= 0,1,2,···,n)上是一个三次多项式； (2) S(x)在每一个内接点xj ( j= 1,2,···,n-1)上具有直到二阶的连续导数； 例 已知f(x):f(-1)=1,f(0)=0,f(1)=1，求f(x)在[-1,1]上的三次自然样条插值函数。 P137 习题五：1 ，2，3，4，7，8 H3(x)应用四个插值基函数表示。 设H3(x)的插值基函数为α0(x)， α1(x) ，β0(x)， β1(x)， 则 其中 可知x1是α0(x)的二重零点，即可假设 由 可得 Lagrange 插值基函数 同理可得 将以上结果代入 得两个节点的三次Hermite插值公式： 二、三次Hermite插值的余项 定理： 设f(x)在区间[a,b]上有定义，f(x)在(a,b)内有4阶导数，H3(x)是满足插值条件 (j=0,1) 的三次Hermite插值函数，则对任意的x∈[a,b]，H(x)的插值余项为 证明： 由 (i=0,1) 可知，x0，x1均为R3(x)的二重零点，因此可设 其中K(x)待定 构造辅助函数 i=0,1 因此φ(t)至少有5个零点。 连续4次使用Rolle定理可得，至少存在一点ξ∈[x0,x1]，使得 即 所以，两点三次Hermite插值的余项为 例1. 解: 作为多项式插值，三次已是较高的次数，次数再高就有 可能发生Runge现象，因此，对有n+1个节点的插值问题， 我们可以使用分段两点三次Hermite插值。 三、分段三次Hermite插值 用三次Hermite插值，当x? [xj , xj+1 ]时，有 5.6 三次样条插值 样条：是 指飞机或轮船等的制造过程中为描绘出光滑的外形曲线(放样)所用的工具。 样条本质上是一段一段的三次多项式拼合而成的曲线，在拼接处，不仅函数是连续的，且一阶和二阶导数也是连续的。 1946年，Schoenberg将样条引入数学，即所谓的样条函数。 因分段线性插值导数不连续，埃尔米特插值导数连续但需要已知，故引入样条插值概念。 一、三次样条插值函数的定义 定义： 给定区间[a,b]上的一个划分：a = x0 x1…xn=b，已知函数f(x)在点xj上的函数值为 f (xj) = yj, ( j= 0,1,2,···,n)如果存在分段函数 满足下述条件： 则称S(x)为节点x0,x1,…, xn 上的三次样条函数。 若S(x)在节点x0,x1,…, xn 上还满足插值条件： (3) S (xj)= yj ( j