第7单元插值法.pptVIP

第7章 插值法 7.1 引言 在科学研究和许多实际工程问题中，常常需要求函数 的函数值、导数值、零点、极值或积分值。但有许多函数没有明显的表达式，而是根据实验或观测得到的一系列离散的数据来表示函数关系的；即使有些函数有表达式，但表达式比较复杂，直接使用它感到不方便。针对这种情况，人们希望构造一个简单的连续函数 来近似地代替所要考察的函数 ，这样使得针对该函数的问题求解得到简化。插值法就是建立这种近似表达式的一种基本方法。 根据实际需要，可以用各种不同的函数来近似原来的函数，比如选代数多项式、三角多项式、有理函数等。由于代数多项式最简单，计算其值只需用到加、减、乘运算，且积分和微分都很方便，所以常用它近似表达复杂的函数或表格函数，这样的插值方法叫做代数插值法，简称插值法。代数插值的数学定义：对于给定的列表函数，求一个 次（或次数不高于 的）多项式 (7.1) 使它在已知点 取给定函数值 ，即满足： (7.2) 则称 为列表函数的插值多项式。 点 称为插值节点，包含所有插值节点的区间 称为插值区间。 插值法的几何意义就是通过 个点 ， 作一条代数曲线 ，使其近似于曲线 。如图7.1所示。 一般地，显然过 的 个给定点的曲线存在，且不惟一。那么，插值多项式是否存在并且惟一呢？ 次多项式 中有 个待求系数 ，由插值条件(7.2)式应用下列方程组： (7.3) 成立。 此方程的系数 的行列式 是一范德蒙行列式。 因为 互不相同，故 ，因此方程组存在惟一的解 ，这说明 存在并惟一。 上述过程说明，只要通过求解方程组(7.3)得出各 的值，便可以确定 了。但由于 较大时，不仅解线性方程组的计算量大，而且方程组往往是病态的，因此这种方法一般不宜采用。通常我们采用的是构造型的方法，即构造一个满足插值多项式(7.2)的 次多项式，而这个多项式是容易求得的。本章介绍直接构造 的拉格朗日插值法和牛顿插值法法。 令 (7.5) 显然 (7.6) 于是由(7.2)式 多项式 是 个 次多项式 的线性组合，因而 的次数不高于 。形如(7.4)式 的称为拉格朗日插值多项式。为了以后便于区别，常用 代替 ，以突出表示这是由拉格朗日插值所得到的插值多项式，即 (7.7) 拉格朗日插值多项式(7.7)可以写成较简单的形式： 其中 特别地，当 时，即得到 的线性插值多项式 这就是通过两点 、 的直线。见图7.2所示。 当 时，即得到 的二次插值多项式 如图7.3所示，这是一条过点 、 、 的抛物线。这说明用二次插值多项式近似函数，其几何意义是用过曲线上三点的抛物线近似曲线，故三点二次插值又称为抛物线插值。 7.2.2 实现拉格朗日插值法的基本步骤 （1）输入插值次数 、插值点 以及计算插值多项式的 值； （2）对 ，计算： （3）输出 ； （4）结束。 例7.1 已知函数 的函数表如下：