数值分析-插值法


数值分析-插值法

插值问题

插值法概述

插值法:设函数在区间上有定义,且已知在点上的值,若存在一简单函数,使 成立,就称的插值函数,点称为插值节点,包含插值节点的区间称为插值区间

多项式插值

多项式插值存在性唯一性定理:设在区间上给定n+1个点上的函数值,存在唯一次数不超过n的多项式,使

证明:由此可得关于系数的n+1元线性方程组 此方程组的系数矩阵为 由于互异,故 因此,方程组的解存在且唯一


拉格朗日插值

拉格朗日插值多项式

插值基函数:若n次多项式在n+1个节点上满足条件 就称这n+1个n次多项式为节点上的n次插值基函数

容易得到次插值基函数为

拉格朗日插值多项式:形如 的插值多项式称为拉格朗日插值多项式

插值余项与误差估计

定理:设上连续,内存在,节点是满足条件的插值多项式,则对任何,插值余项 这里且依赖于x

证明:容易得到 为了求出,构造辅助函数 可得在区间上有个零点,由罗尔定理,可知在区间上至少有个零点,以此类推,最后知在区间上存在一点,使得,即


牛顿插值

差商

差商: - 称为函数关于的一阶差商 - 称为的二阶差商 - 一般地,称 的k阶差商

  • 差商与节点的排序无关,称为差商的对称性,即
  • n阶差商与导数的关系为

牛顿插值多项式

牛顿插值多项式:形如 的插值多项式称为牛顿插值多项式

牛顿插值余项

定理:牛顿插值多项式余项满足

埃尔米特插值

插值多项式要求在插值节点上函数值相等,但有的实际问题还要求在节点上导数值相等,甚至高阶导数值也相等,将满足这种要求的插值多项式称为埃尔米特插值多项式

两点三次埃尔米特插值

两点三次埃尔米特插值:插值节点取为,插值多项式为,满足条件

:采用基函数方法,令 其中是关于节点以及的三次埃尔米特插值基函数,它们应分别满足条件 于是求得 同理可求得


利用重节点差商,还可以得到一种所谓的插值公式

插值余项与误差估计

两点三次埃尔米特插值余项:若存在,则对,两点三次Hermite插值的误差为

其中,且依赖于

分段低次插值

高次插值的病态性质

上面我们根据上给出的节点插值做插值多项式近似,一般总认为的次数n越高逼近的精度越好,但实际上并非如此。

这是因为对任意的插值节点,当时,不一定收敛于,即龙格现象

分段线性插值

分段线性插值函数:设已知节点 上的函数值 ,记,折线函数满足:

  • ;
  • ;
  • 在每个小区间上是线性函数

则称为分段线性插值函数

定理:若函数,则上一致收敛到

证明:由定义可知,分段线性插值的误差估计可以用插值余项得到 其中,由此知


分段三次埃尔米特插值函数

分段线性插值函数的导数是间断的,若在节点上除已知函数值外还给出导数值,这样就可构造一个导数连续的分段插值函数

分段三次埃尔米特插值函数:满足条件

  • ;
  • ;
  • 在每个小区间上是三次多项式

称为分段三次埃尔米特插值函数

误差估计:设在节点上的分段三次埃尔米特插值多项式,则有 其中

三次样条插值

三次样条插值函数:设函数,且在每个小区间上是三次多项式,若在节点上给定函数值,并成立 则称为三次样条插值函数

从定义知要求出,在每个小区间上要确定4个待定系数,而共有n个小区间,故应确定4n个参数。根据上二阶导数连续,在节点处应该满足连续性条件 这里共有3n-3个条件,再加上满足插值条件,共有4n-2个条件,因此还需要加上2个条件才能确定。通常,可在区间的端点上各加一个条件(称为边界条件),可根据实际问题的要求给定。常见的有以下3种:

  • 已知两端的一阶导数值,即

  • 已知两端的二阶导数,即

  • 是以为周期的周期函数时,则要求也是周期函数。这时边界条件应满足 这样确定的样条函数称为周期样条函数

定理型三次样条插值问题的解存在且唯一

证明:记,将其设为待定参数,由于上为三次多项式,故上为线性函数,于是可得的分段表达式 将上式积分两次,并利用来确定积分常数,可得

要证唯一存在,只需证明可被唯一确定

对于上述满足其自身和二阶导函数在整个区间上连续,由连接条件,还要求一阶导数在整个区间上连续,即

我们有

于是,单元右端点的左极限和左端点的右极限分别为

从而有

若记 则上式可化为 此外,在端点处应分别右连续和左连续,利用附加边界条件,可得 化简得到 综合所有条件,可以得到一个关于待定参数阶线性代数方程组,其矩阵形式为 显然方程组的系数矩阵是严格对角占优的,因此非奇异,于是可以唯一的解出,从而得到型三次样条插值问题的解存在且唯一


三次样条插值函数的基本性质

定理:若型三次样条插值函数,则对,有

证明:令,则,且满足 ,有


最佳逼近定理:若型三次样条插值函数,则对,有

证明:利用引理及Cauchy不等式,有 从而


  目录