第四章 数值插值方法解析.ppt

第五章 插值方法 第五章 插值方法 插值的基本概念 Lagrange插值 分段低次插值 均差与Newton插值 Hermite插值 三次样条插值 5.1 代数插值问题 例. 某地区某年夏季时节间隔 30 天的日出日落时间为 日照时间的变化设为 y(x)=a0+ a1x + a2x2, 设 Pn (x)=a0 + a1x + a2x2 + ……+ anxn 是y=f(x)在[a,b]上的n+1个互异节点x0,x1,…,xn的插值多项式，则求Pn (x)问题归结为求系数a0,a1,…,an。 故Pn (x)存在且唯一。 5.2 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次多项式本身。 这说明截断误差很小。 设节点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 二、三次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= 0,1,2,···,n) 则称S(x)为三次样条插值函数。（即全部通