第三章插值与拟合.pptVIP

第三章 插值与拟合 在科学研究与工程技术中，常常需要从一组测量数据出发，寻找变量之间函数关系的近似表达式，这就是函数逼近问题，插值和拟合都是函数逼近的重要方法。 应用举例:胡克定律---弹簧在力F的作用下伸长x，在一定范围内服从胡克定律，F与x成正比，即F=kx，k为倔强系数，现在实验中测得以下一组x，F的数据，并在(x，F)坐标下作图。可以看出，当F达到一定数值后，就不服从胡克定律。 试由数据确定 k，并给出不服从胡克定律时的近似公式。 第一节 插值问题 在生产实践的许多领域里，例如机械工业、造船、汽车制造，常常有这样的问题：给了一批离散样点，要求作出一条光滑曲线（乃至于曲面），使其通过或尽可能的靠近这些样点，以满足设计要求或者据此进行机械加工－－这就是插值问题。 插值法历史悠久。据考证在公元六世纪时，我国的刘焯（zhuo）已把等距二次插值应用于天文计算。十七世纪时，Newton和Gregory（格雷格里）建立的等距节点上的一般插值公式，十八世纪时，Lagrange给出了更一般的非等距节点插值公式。 插值问题的基本概念 设函数y=f (x)在区间[a,b]上有定义， 且它在该区间上的n+1个互异点 的函数值已知，即 如果选取简单函数P(x)作为y=f (x)的 近似表达式，并满足以下条件： （3-1） 称为插值条件；满足插值条件的近似函数P(x)就称为f (x)的插值函数，而y=f (x)称为被插值函数；互异点x0,x1,...,xn被称为插值节点(简称节点)，x称为插值点；区间[a,b]称为插值区间。 构造插值函数需要关心下列问题： 插值多项式存在的唯一性 定理2-1 在n+1个互异节点(xi，f (xi))上满足 插值条件 的次数不高于 n次的插值多项式 由于节点互异，即 ，所以 　　　 ，由克莱姆法则可知方程组（3-3）有唯一的一组解a0 , a1 , ...an ，也就是插值多项式（3-2）存在且唯一。 第二节 拉格朗日插值多项式 拉格朗日插值多项式是一种在形式上不同于（3-2）的插值多项式，这种多项式的明显优点是无需求解方程组，只要给出了n+1个互异节点及对应的函数值，便能直接写出这种形式的插值多项式。 拉格朗日插值多项式的构造方法 已知函数y =f (x)在 n+1个节点 处的函数值为 ，现 要作一个n次插值多项式 ，并使 在节点xi处满足 将构成n+2个方程的联立方程组，而未知量 只有n+1个，根据非齐次线性方程组有解的 充要条件必有 利用范德蒙行列式的计算公式，经整理后可得到的表达形式如下： 其中 叫做拉格朗日插值多项式的插值基函数，显然它满足 两个节点 带公式得 线性插值也可以直接列方程组求解，即 二次插值多项式 插值余项的估计： (3-10)式称为多项式插值的余项公式，由于无法求出公式中 的准确值，所以实际计算中用它来估计误差仍有困难。因此，只在理论分析中应用。实际应用中通常采用(3-11)式来作为多项式插值的余项 估计式。 余项 线性插值：选11,12两点带公式计算 例2:已知e-x在x=1,2,3的值如下，试用二次 Lagrange插值求e-2.1的近似值，并进行误差估计。 因 ，故误差的余项 第三节　牛顿插值多项式 　　 当使用拉格朗日插值多项式为了提高精度增加插值节点时，插值多项式就得重新构造，整个公式改变之后，之前的计算结果不能继续使用，计算工作必须全部从头做起。为了克服这一缺点，下面介绍另一种插