Cht2插值法.ppt

湖南科技大学 §1 引 言 一、线性插值与抛物线插值 四个点逐次插值，列表 例如三次Aitken插值 §6 埃尔米特插值 在实际问题中，对所构造的插值多项式，不仅要求通过点(xi, yi) (i = 0, 1, 2,…, n)，同时还要求在xi上,若干阶导数值也与被插函数的导数值相等。把此类插值称为Hermite插值，插值多项式称为Hermite插值多项式，记为H(x)。 其求法如同不带导数的插值多项式的求法， 一是借助Lagrange插值的基函数方法，可得出Lagrange型Hermite插值多项式， 二是通过构造重节点差商表，得出Newton型的Hermite插值多项式。 差分表 求解方法见差分表 (4.9)的系数 (4.12)的系数 2.5 逐次线性插值法 Lagrange插值法和Newton插值法分别是选用不同的两组基函数，得到的插值多项式。这两种方法有一个缺点： 我们先从简单情形开始，已知三点P0(x0,y0),P1(x1,y1),P2(x2,y2). 求过此三点的插值函数。采取分步的策略。 第一步：过点P0,P1作线性插值，有 为克服这一缺点，采用逐次线性插值方法。 不能确定到底用几次多项式逼近， ◎如果次数太低，多项式不能很好地逼近函数f(x)； ◎如果次数太高，多项式的值会产生振动。 第二步：过点P1,P2作线性插值，有 第三步：将(x0,P01),(x2,P12) 也理解为“点”，过此两点作“直线” P012有何特性？一是它是二次函数，二是P012(xi)=yi(i=0,1,2) P012就是所求的插值多项式，它是按步通过三个线性插值获得的，所以称为逐次线性插值。 一般的，利用两个n-1次插值多项式 P01...n-1和P12...n，再作线性插值，结果得到n次插值多项式。即 Neville逐次线性插值公式。 Aitken逐次线性插值公式。 也可利用两个n-1次插值多项式 P01...n-1和P012...(n-2),n，再作线性插值，结果得到n次插值多项式。即 从表上看每增加一个节点就计算一行，斜线上是1次到3次插值多项式的值，如精度不满足要求，再增加一个节点，前面计算完全有效，这个算法适用于计算机上计算，且具有自动选节点并逐步比较精度的特点，程序也比较简单。算例见教材(略）。 问题的一般提法是： (2.6.1) 带导数的两点插值（重要特例) 结论： 函数表及导数表 问题:已知 求3次多项式 使满足插值条件: 存在且唯一，表达式为 其中 证明思路与不带导数插值类似。 一般地，已知 x0 , …, xn 处有 y0 , …, yn 和 y0, …, yn ，求 H2n+1(x) 满足 H2n+1(xi) = yi ， H2n+1(xi) = yi。 解：设 ? + = n i ) ( ) ( ) ( = 0 i i x H x h yi x H2n+1 ? n = 0 i yi 其中 hi(xj) = ?ij , hi(xj) = 0, (xj) = 0, (xj) = ?ij Hi Hi hi(x) 有根 x0 , …, xi , …, xn且都是2重根 ? ? ) ( ) ( ) ( 2 x l B x A x h i i i i + = 由余下条件 hi(xi) = 1 和 hi(xi) = 0 可解Ai 和 Bi ? (x) Hi 有根 x0 , …, xn, 除了xi 外都是2重根 ? Hi ) ( ) ( i i li2(x) x x C x - = 又: Hi (xi) = 1 ? Ci = 1 Hi ) ( x ) ( i li2(x) x x - = 设 则 这样的Hermite 插值唯一 方法：基函数法，待定参数法，重节点差商表 例：在[?5, 5]上考察 的Ln(x)。取 - 5 - 4 - 3 - 2 - 1 0 1 2 3 4 5 - 0.5 0 0.5 1 1.5 2 2.5 n 越大， 端点附近抖动 越大，称为 Runge 现象 Ln(x) ? f (x) ? Runge现象 问题：高次插值过程的收敛性如何？ ？ 结论：高次插值产生龙格(Rnuge)现象。 注:实际应用时取 取等距节点 举例：对 采用10次插值，则有 Rnuge现象： 等距节点高次插值所产生的小区间内逼近很差的 现象。 §7 分段低次插值 一、高次插值的病态性质 1. 分段线性插值 Ih(x)是平面上以点(xi,yi