计算方法讲义课件 五 插值.docVIP

第五章 插值 插值在科学计算和工程技术中有广泛应用。例如由实验得到一系列点x0, x1,…, xn对应的值y0, yi,…, yn，要构造函数y = f (x)，使yi=f(xi)，这就是简单的插值问题。插值核心问题是：存在性、唯一性、表示方法以及误差分析。插值和逼近有广泛应用，例如构造曲线曲面等。 5.1 代数插值 用代数多项式作为工具来研究插值的方法叫做代数插值。 插值 插值问题就是根据已知数据来构造函数y = f (x)的近似表达式。常用方法就是利用多项式Pn(x)，使 ，作为f(x)的近似。多项式求值方便，且有导数。称Pn(x)为f(x)的一个插值函数，称x0, x1,…, xn为插值节点。用代数多项式作为工具来研究插值的方法叫做代数插值。设x0 x1 … xn，记a = x0, b = xn，则[a, b]为插值区间。 设所要构造的插值多项式为：，由插值条件 。得到如下线性代数方程组：。该线性方程组的系数行列式为，为范得蒙行列式。当，时，D ? 0，所以Pn(x)由a0, a1,…, an唯一确定。 5.2 Lagrange插值 已知y = f (x)在给定点x0, x1上的值为y0，y1。线性插值就是构造一个一次多项式 P1(x) = ax + b，使它满足条件P1 (x0) = y0，P1 (x1) = y1。几何解释就是一条直线。由解析几何， 或。 例 用线性插值求（x* = 10.723805）。 解：设，取x0 = 100，x1 = 121， 则y0 = 10，y1 = 11，从而。 用简单的曲线近似地代替复杂的曲线，最简单的曲线是二次曲线。设函数y=f (x) 在给定互异的自变量值x0, x1, x2上对应的函数值为y0, y1, y2，二次插值就是构造一个二次多项式，使之满足，称抛物插值。令，须有：，可以求解：，，。 例 用抛物插值求（x* = 10.7238）。解：设，取x0=100,x1=121,x2=144,。 例 用抛物插值求25的立方根（x* = 2.924）。 设连续函数y = f（x）在[a, b]上对给定n + 1个不同结点：x0, x1,…,xn分别取函数值 y0, y1,…,yn其中y-i = f (xi),i = 0, 1, 2,…, n，构造次数不超过n的插值多项式，使之满足条件，i = 0, 1, 2,…, n。为此令，则。令，则Pn(x)满足，i = 0,1,…, n，Pn(x)称为拉格朗日插多项式。 把差f (x)－Pn(x)称为用插值多项式Pn(x)代替f (x)的余项，记为： 。设x是 [a, b]中任意固定的数，若x是插值节点xi，则，如x不是节点，考虑如下函数：，其中 。由插值条件知，i = 0, 1,…, n，且。所以，?(t)在[a, b]上有n+2个互异零点。应用罗尔定理，?’(t)在?(t)的每两个零点间有一个零点，即?’(t) 在[a, b]上至少有n+1个零点，进一步对?’(t)应用罗尔定理可知，?? (t)在[a, b]上至少有n个零点，继续上述讨论就可推得?(n+1)(t)在[a, b]内至少有一个零点，记之为?，即。因Pn (t)为不高于n次的多项式。所以于是，整理可得。得以下余项定理： 定理：设f (n)(x)在区间[a, b]上连续，f (n + 1)(x)在[a, b]上存在，x0, x1,…, xn是[a, b]上互异的数，记插值余项为，那么当x ? [a, b]时，有如下估计，其中，。 5.3 Hermite插值多项式 在实际问题中，有时要求通过点（xi, yi）i = 0, 1, 2,…, n，还要求在xi上，阶导数值也与被插函数的导数值相等。把此插值多项式称为埃米尔特（Hermite）插值多项式，记为H (x)。 要求在1个节点 x0处直到m0 阶导数都重合的插值多项式即为Taylor展开式：。 依照Lagrange插值的做法，先确定多项式插值空间的维数。 注意到插值条件共有2(n+1)个条件，所以次数为2n+1 。设（Hermite插值多项式），问题变为求函数，满足：。 【例】构造H(x)使在x1,x2取值为f1,f2,导数为f1’,f2’。令H(x)=h1(x)*f1+h2(x)*f2+g1(x)*f1’+g2(x)*f2’。令x1=0,x2 =1则可以简化问题。用待定系数法求解，可以得到以下结果：h1(x)=(1+2x)(x-1)2 , h2(x)=(3-2x)x2 , g1(x)=x(x-1)2 , g2(x)= (x-1)x2 。 图 h1,h2,g1,g2的几何意义 根据h1,h2,g1,g2的几何意义，可以构造出无穷多个h1,h2,g1,g2函数，满足表中所列出的性质