更新时间：2024-09-16 14:05
插值法又称“内插法”，是利用函数f (x)在某区间中已知的若干点的函数值，作出适当的特定函数，在区间的其他点上用这特定函数的值作为函数f (x)的近似值，这种方法称为插值法。如果这特定函数是多项式，就称它为多项式插值。常用的几种多项式插值法有：直接法、拉格朗日插值法和牛顿插值法。
给定n+1个点 （称为插值点），所谓多项式插值就是找到一个多项式（称为插值多项式）
使得它满足条件
其中，i=0,1,...,n。也就是说，多项式y=P(x)的图像要经过给定的n+1个点。
在实际应用中，这些插值点可能来自某次实验测量所得的数据，也可能来自某个复杂函数 的值。通过计算插值多项式，我们可以找到这些实验数据间的规律，或者使用简单的多项式函数 来近似复杂的函数 。
定理一：
给定n+1个点 ,若 两两不同，则存在唯一一个次数不超过n的多项式 ，使得 成立。
证明：利用范德蒙德矩阵和代数学基本定理即得。
当 的值来自某个函数 ，且f(x)具有n+1阶连续导数时，下面的定理可以用来计算多项式插值的（截断）误差。
定理二：
给定n+1个点 ,其中 ，进一步假设函数f(x)具有n+1阶连续导数，则插值多项式P(x)的误差R(x)为
其中，
给定n+1个点, 计算插值多项式的主要方法有：直接法、拉格朗日多项式插值和牛顿多项式插值。下面我们分别介绍这三种方法。
（注意，根据定理一，这三种方法得到的插值多项式在理论上说应该是一致的，而且误差也相同。）
直接法
根据定理一，假设插值多项式为
由插值条件 ，我们得到关于系数 , ,…, , 的线性方程组
通过求解这个线性方程组，即得到插值多项式。
优点：直接，性质一目了然。
缺点：待求解的线性方程组的系数矩阵为范德蒙德(Vandermonde)矩阵，它是一个病态矩阵，这使得在实际求解方程组时将产生很大的误差。
拉格朗日多项式插值
拉格朗日(Lagrange)多项式插值的原理是：先构造一组拉格朗日基函数 ，这些基函数为次数不超过n的多项式，且具有性质
然后将这些基函数做线性组合，得到拉格朗日插值多项式
容易验证，多项式L(x)满足插值条件
拉格朗日基函数 的构造如下：
由基函数的性质，当 时， ，即 为 的零点，可以假设
其中，K为待定系数。再由 ，得到
从而得到
因此，基函数
令 ，则 还可以表示为
下面的定理说明 称为基函数的原因：
定理三：令 为全体次数不超过n的多项式构成的集合，则 是线性空间 的一组基。
Matlab 实现
均差与牛顿多项式插值
牛顿多项式插值是基于均差的计算。首先定义均差如下：
函数f(x)关于点的一阶均差（或差商）为
一阶均差反映了函数在区间的平均变化率
用递归的方式，我们定义二阶均差为
同理，k阶均差为
特别地，0阶均差定义为 。
根据均差的定义，构造均差表如下：
如果将x也看作一个点，由均差的定义可以得到
其中，
称为牛顿插值多项式。
为插值余项。
由定理一和定理二得到均差和导数的关系如下：
其中，
Matlab实现
拉格朗日多项式插值的计算量大于牛顿多项式插值的计算量。
特别地，当新增一个插值点时，拉格朗日插值需要重新计算全部的基函数，而牛顿插值只需计算均差表中新的一行的值即可。