多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值...

多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值..._第1张图片

前言

接下来的主要内容是与插值逼近相关的内容。科学计算中的许多问题需要计算函数

的函数值,但是由于计算机硬件系统只能提供加、减、乘、除四则运算和逻辑运算,所以使得绝大数函数不能直接用计算机计算。有时虽然已知
的表达式,但函数值的计算不能通过有限次的四则运算完成,又如有的函数虽然存在,却只能给出它在离散点的函数值,而无法给出其解析表达式。因此,需要为
构造一个在计算机上简单易算的近似函数,并用该近似函数代替
进行数值运算,这就是数值逼近问题。

插值方法是数值分析中一个简单而又重要的方法,利用该方法可以通过函数在有限点处的函数值求出其近似函数,进而夫算出函数在其他点处的值。插值方法在离散数据处理、函数的近似表示、数值微分、数值积分、曲线与曲面的生成等方面有重要作用。

下面我们先介绍插值的基本定义和插值格式,再介绍三种多项式插值方法,分别是Lagrange(拉格朗日)插值、Newton(牛顿)插值和Hermite插值。

一、引言

问题描述,面对给定的点和相关函数的信息,我们要估计未给出点的函数值,就需要使用到插值方法,

多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值..._第2张图片
问题描述

从几何上来理解,如下图所示

多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值..._第3张图片
插值几何角度解释

p(x)是我们计算出来的对f(x)的拟合函数,并且穿过给定的数据对

,这样如果有新的数据点
,则直接代入
来预测具体值应该是多少。

基本的思想:

  1. 简单函数类基底需要满足的条件
  2. 给出具体的基底
  3. 给出系数

多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值..._第4张图片

当我们看到

,就可以想到我们可以利用线性方程组的解法进行求解,

多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值..._第5张图片

但是我们都知道,线性方程组不一定有解,需要满足一定的条件,这个条件定义如下,

多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值..._第6张图片

我们称这个条件为解的存在的Haar条件。

接下来我们就可以定义插值基函数,

多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值..._第7张图片

上面的假设可以这样理解,插值基函数再我们此时的定义中是未知的,是这样一种假设的形式,我们希望对于每一个已知点

得到的结果与
相同,那么我们就有这样的基本认为,我们认为对于函数
来说,当k=i的时候希望
,这样我们得到的就是
,满足我们的需求,其他情况等于零即可。

二、Lagrange插值

我们之前已经得到了插值基函数的基本形式,接下来的问题就是解决如何求解

,这样我们就可以得到
的境确表达形式。

多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值..._第8张图片

多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值..._第9张图片

多项式函数显然满足Haar条件,这样我们就可以利用插值的方法来解决问题,并且解是唯一的。

多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值..._第10张图片

正如上面所说的,我们完全可以利用线性方程组来求解系数,但是非常困难,涉及很复杂的矩阵分解。这里我们就提出了Largange插值和Newton插值来求解这些系数。

那么我们先介绍Largange插值方法

多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值..._第11张图片

构造

由于

,则
,根据假设
故,上式变为

,那么我们就可以得到α的计算方式

,这样我们得到基函数
的表示。

,其中

多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值..._第12张图片

下面我们看一道例题来理解一下,

多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值..._第13张图片

Lagrange插值法计算起来非常的简单便捷,很快我们就能得到问题的精确解,但是我们也发现一个问题。

在插值问题中,为了提⾼插值精度,有时需增加插值节点个数。插值节点个数发⽣变化后,所有的Lagrange插值基函数都会发⽣变化,从⽽整个Lagrange插值多项式的结构发⽣变化,这在计算实践中是不⽅便的。

那么这种问题如何改进和解决呢?

三、Newton插值法

为了克服Lagrange插值多项式的缺点,能灵活地增加插值节点,使其具有“承袭性”,即可以充分利用已有的信息,我们引进Newton插值公式。

查找了一些资料,没有找到非常易懂的从Lagrange到Newton的推导过程,只能自己自力更生,那么这里给出我自己的推导,以供大家深入理解Newton法。

多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值..._第14张图片

商差是为了便于Newton插值的计算和表达引入的概念,便于我们迭代求解。我们给出商差的严格定义,

多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值..._第15张图片
商差

他有下面几种性质

多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值..._第16张图片

多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值..._第17张图片

第三个性质经常会在题目中出现,作为填空和选择题出现。

实际上,至此为止Newton法的基本形式和求解方式已经介绍完毕,我们可以看到实际上Newton插值和Lagrange插值完全是相同的,只是表达形式不同而已。

多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值..._第18张图片

还是通过两道例题观察一下

多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值..._第19张图片

这道例题使用到了商差性质3。

多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值..._第20张图片

考试中可以直接使用均差表来进行计算,Newton插值表达形式。

插值余项

得到Lagrange何Newton插值的求解方法,如何来衡量插值的结果的优劣程度,需要一个指标,我们引入插值余项,来衡量插值与实际真实值之间的差距。

多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值..._第21张图片
插值余项

Lagrange与Newton法使用的都是非导数信息,求解插值表达式,如果给出一些一阶导的条件如何去求出插值表达式呢?

四、Hermite插值

理论和应用中,提出的某些插值问题,要求插值函数

具有一定的光滑度,即在插值节点处满足一定的导数条件,这类插值问题称为Hermite插值问题。

那么我们定义问题,如下所示

多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值..._第22张图片

要构造满足上述插值条件的插值公式,借鉴的是Lagrange插值公式。

给出构造方法如下所示

多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值..._第23张图片
Hermite插值公式构造方式

时,Hermite退化为Lagrange插值,

我们来看两道例题帮助我们观察如何构造Hermite插值多项式。

多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值..._第24张图片

下面看一个重点要求掌握的,二点三次Hermite插值多项式的构造。

多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值..._第25张图片

多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值..._第26张图片

多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值..._第27张图片

多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值..._第28张图片

对于二点三次Hermite插值多项式四个系数的推导可以不掌握,但是需要记忆四个系数的求解公式,考试中会考察二点三次Hermite插值的构造。

应该有不少同学有和我一样的疑问,为什么二点Hermite插值结果是三次呢?这个就需要我们重新审视假设条件

多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值..._第29张图片

假设有这样一个情况

,那么每一个构造的
,这样显然就是
次数为
,所以两点三次Hermite插值就出现了。

下面简单给出二点三次Hermite的插值余项,

多项式的拟合与插值例题_矩阵与数值计算(12)——Lagrange插值、Newton插值、Hermite插值..._第30张图片

总结

综上我们介绍了多项式插值的基本内容,后续会介绍逼近相关内容。

你可能感兴趣的:(多项式的拟合与插值例题)