.-.插值法【DOC精选】.docVIP

第三章 插值法 背景问题 在现代机械工业中用计算机程序控制加工机械零件，根据设计可给出零件外形曲线的某些型值点的数据，加工时为控制每步走刀方向及步数，就要算出零件外形曲线其它点的坐标数据，才能加工出外表光滑的零件。在数学上，该问题相当于，虽然函数在某个区间上是存在的，还可以是连续的，但却不知道其解析表达式，只能给出上系列点的函数值，这只是一个函数表，如何（近似）计算在上的其它点处的函数值。有的函数虽然知道表达式，但比较复杂，计算值很不经济，通常造一个函数表，如大家熟悉的平方根表、三角函数表、对数表等，如何用简单的计算（近似）求出不在表上的函数值。因此，我们希望根据给定的函数表，寻求一个简单的函数，使它即能反映的特性：，又便于计算，用来近似。问题：选取什么函数作为近似的函数，如何求得其具体表达式，误差如何？ §1　插值法 1－1　插值问题 设函数在区间上有定义，且已知在一组互异点上的函数值，寻求一个简单的函数，使满足 　　 　 （1.1） 并用近似代替，上述问题称为插值问题。 插值条件：； 插值函数：称为的插值函数； 插值节点：； 插值区间：； 插值法：按插值条件（1.1），求函数的插值函数的方法称为插值法。 一般地，我们经常选取多项式、三角函数、有理函数等作为插值函数，尤其是代数多项式被广泛用作插值函数，此时的插值函数又称作插值多项式，相应的求插值多项式的方法称为多项式插值。下面我们重点研究多项式插值。 问题：（1）满足插值条件（1.1）的插值多项式是否存在？应该是几次多项式？（次） （2）如果满足插值条件（1.1）的多项式存在，应如何构造？ （3）用插值多项式近似代替，误差如何？ 从几何上看，插值问题即：已知平面上个不同的点，要寻找一条过这些点的多项式曲线。（不超过次） 1－2　插值多项式的存在唯一性 问题（1）即设法找出个系数，使满足插值条件：，即 系数矩阵的行列式为阶Vandermonde行列式，且互异，故行列式不等于零，所以上述关于的元方程组有唯一解。 定理1　（插值多项式存在唯一性）满足插值条件（1.1）的插值多项式存在且唯一。 1－3　插值基函数及Lagrange插值 从定理1的证明看到，要求插值多项式，通过求方程组的解即可，但这样做不但计算复杂，且难于得到的简单表达式。下面先讨论简单的情形。 时 问题1：两个点，要求一次插值多项式使得： 求得 问题2：两个点，要求一次插值多项式使得： 求得 问题3：两个点，要求一次插值多项式使得： 几何直观即直线，得点斜式； 可写成 是两个一次函数与的线性组合。称函数，为线性（一次）插值基函数。（重新验证插值条件，体会插值基函数的作用） 时 问题1：三个点，要求二次插值多项式使得： 几何直观即过三点的抛物线，利用两个零点，易得 同理可求得满足插值条件：的二次插值多项式 满足插值条件：的二次插值多项式 问题2：三个点，要求二次插值多项式使得： 易求得： 是三个二次函数、、的线性组合。称函数、、为二次插值基函数。 一般地，对于，给定个点，寻求次的多项式，使得满足条件： 根据零点易得： 定义：称次多项式为节点上的次Lagrange插值基函数。 若给定个点，容易求得过这些点的次插值多项式为： （1.7） 形如（1.7）的插值多项式称为Lagrange插值多项式。一个次Lagrange插值多项式可由个次Lagrange插值基函数组合而成，其组合系数恰好是节点上的函数值。 不超过次的多项式全体构成维线性空间，有各种各样的基函数，因而也有各种形式的插值多项式。（前面讨论的解方程组的方法求插值多项式，相当于以为基）但唯一性是肯定的。 例：给定的函数表如下： 144 169 225 12 13 15 写出二次Lagrange插值基函数，并用二次Lagrange插值多项式计算的近似值；用一次的又如何？，用一次时如何选取插值节点? 从图中可知的误差比的误差小，一般地，是否插值多项式的次数越高，其计算结果就越精确？如何估计与之间的误差？ 下面讨论与之间的误差问题。 §2　插值多项式中的误差 2－1　插值余项 已知函数在区间上的个节点处的函数值，在其余点处用Lagrange插值多项式近似代替，误差怎样估计？ 设，则，即有个零点。于是可令 构造 显然共有个零点，反复应用Rolle定理可得： 的零点 的零点 的零点 ……… ……… 的零点 的零点 与有关，使得，于是 称为次Lagrange插值多项式的余项。 定理1.1　（插值多项式余项公式）设在含节点的区间上次可微，是关