插值法

一、拉格朗日插值法

①原理

对于平面上已知的n个点，能够找到一个n-1次多项式y=a₀+a₁x+a₂x²+…+a_n-1x^n-1，使得该多项式曲线通过这n个点。对于缺失的值，将对应的x值代入多项式便可获得近似值L(x)。

②公式推算

已知：

将n个点的坐标(x_i,y_i)代入多项式，得：

解得拉格朗日插值多项式为：

其中l_i(x)称为插值基函数：

 

拉格朗日插值公式是在节点上给出节点基函数，而后作基函数的线性组合，组合系数为节点函数值。

 

存在性证实：

对于给定的n个点,拉格朗日插值法的思路是找到一个在点(x_i,y_i)取值为1，而在其余点取值都是0的多项式l_i(x)，这样，多项式y_il_i(x)在点x_i处取值为y_i，而在其余点取值为0。

惟一性证实：

对任意两个次数不超过n-1的多项式，的差值在全部n个点的取值都为0：

因此：

当φ(x)≠0时，其次数必定不小于n，可是次数不超过n-1，φ(x)次数不该该超过n-1。因此有：

几何性质：

基函数l_i(x)可当作由次数不超过n-1的多项式组成的线性空间的一组基底。

若是存在一组系数λ_i，使得：

 

则有：

所以：

证实l_i(x)线性无关。

缺点：

一、在计算中，当插值点增长或减小一个时，所对应的基本多项式就须要所有从新计算，因而整个公式都会变化，很是繁琐。计复杂度为O(n²)。

二、当插值点比较多的时候，拉格朗日插值多项式的次数可能会很高，所以具备数值不稳定的特色，也就是说尽管在已知的几个点取到给定的数值，但在附近却会和“实际上”的值之间有很大的误差（如右下图）[6]。这类现象也被称为龙格现象，解决的办法是分段用较低次数的插值多项式。

③改进——重心拉格朗日插值

不妨设：

其中w_i称为重心权，则改进后：

优势：

一、当数据点的个数增长一个x_n+1时，可在O(n)复杂度内更新每一个w_i(w_i←w_i/(x_i-x_n+1))，以及计算新的重心权w_n+1，并在O(n)时间内计算获得新的L(x)。

二、牛顿插值法

①原理

牛顿（Newton）插值公式是代数插值方法的一种形式。牛顿差值引入了差商的概念，使其在差值节点增长时便于计算。

②公式推演

求已知的n个点对(x_i,y_i)的全部阶的差商公式：

变换获得：

而后其中1式+(x-x1)·2式+···+(x-x₁)···(x-x_n-1)·n式，获得：

其中：

P(x)称为牛顿插值逼近函数，R(x)称为偏差函数。

一、对于P(x)，每当新增一个插值数据点，前面已有的各项不变，计算只增长一项。每项系数均为各阶差商。

二、仍没有改变拉格朗日的插值曲线在节点处有尖点，不光滑，插值多项式在节点处不可导等缺点。

三、Hermite插值(埃尔米特插值)

①原理

埃尔米特插值是另外一类插值问题，这类插值在给定的节点处，不但要求插值多项式的函数值与原函数值相同。同时还要求在节点处，插值多项式的一阶直至指定阶的导数值，也与被插函数的相应阶导数值相等。

若在某节点x_i，要求插值函数多项式的函数值、一阶导数值、直至m_i-1阶导数值均与被插函数的函数值相同及相应的导数值相等。咱们称x_i为m_i重插值点节。

设原函数为f(x)，x_i的函数值为y_i，x_i的插值多项式为H(x)。

②相关定理

一、惟一性定理： f(x)关于结点x_i的二重Hermite插值多项式存在且惟一。

二、偏差定理

③二重Hermite插值多项式推导

设在[A,B]范围内，一共有x₀、x₁、···、x_n共n+1个结点，不妨假设x₀<x₁<···<x_n。令设：

根据Hermite插值的要求：

这里一共给出了2(n+1)个条件，可惟一肯定一个次数不超过2n+1的多项式。列出多项式为：

采用拉格朗日插值多项式的基函数法进行推导这2n+1个系数。

先求解插值基函数α_i(x)和β_i(x)，共有2n+1个，且每一个基函数都是2n+1次多项式，且知足：

H_2n+1(x)可用插值基函数表示：

接下来就是肯定基函数。

设：

令·：

根据条件可知：

易知l_i(x_i)=1，整理得：

 

解出：

从l_i(x)定义式，对两边求导，获得：

代入获得α_i(x)表达式：

同理，因为β_i(x_i)=0以及 ，可设：

因为l_i(x_i)=1，可知c=1。所以：

 

④证实Hermite插值多项式惟一

使用反证法，假设另存在一个插值多项式为,则以及知足Hermite插值条件，有：

在每一个结点x_i上均有二重根，则φ(x)有2n+2重根。

但因为φ(x)是不高于2n+1的多项式，故φ(x)=0。惟一性得证。

⑤计算余项----偏差定理

仿照拉格朗日插值余项的证实方法，若f(x)在[A,B]范围内的2n+2阶导数存在，则：

ξ的值与x有关：

⑥三次Hermite插值多项式

当只有2个结点，也就是n=1时，插值多项式记为H₃(x),不妨设这2个节点为x_k、x_k+1。

H₃(x)知足：

相应的插值基函数为：

所以：

其他项为:

⑦通常Hermite插值多项式

四、分段插值

分段差值就是对每个分段区间[x_i,x_i+1]分别进行插值，则最后所得插值函数为一分段函数。

①分段线性插值

对每个分段区间[x_i,x_i+1]用1阶多项式逼近f(x)，相似拉格朗日插值：

②分段三次Hermite插值：

五、样条插值(待补充)

样条插值法是一种以可变样条来做出一条通过一系列点的光滑曲线的数学方法。插值样条是由一些多项式组成的，每个多项式都是由相邻的两个数据点决定的，这样，任意的两个相邻的多项式以及它们的导数(不包括仇阶导数)在链接点处都是连续的。该方法不须要已知导数。

①样条函数定义

设在区间[A,B]上有n+1个结点x₀<x₁<···<x_n，若是存在正整数k使得在[A,B]上的分段函数S(x)知足：

• 在[A,B]上直到k-1阶都有连续导数

• 在每一个小区间[x_i,x_i+1]上为次数不大于k的多项式

则称S(x)为k次样条函数。

②三次样条插值

一、三次样条函数定义

在区间[A,B]上有n+1个结点x₀<x₁<···<x_n，已知y_i=f(x_i)，若是分段函数：

知足：

• S(x)在[A,B]为三次多项式
• 在区间[A,B]上连续

 

• S(x_i)=y_i

则称S(x)为三次样条插值函数。

二、公式推导

易知一共有n个三次多项式，须要待定4n个系数。以[x_i,x_i+1]为例：

插值条件：

连续性条件：

上述条件一共有4n-2个方程，但未定系数有4n个，剩余条件经过边界给出：

第一边界条件（固定边界条件）：

第二边界条件：

当二阶导数恒为0时，称之天然边界条件。

周期边界条件：

三、三转角法——待定一阶系数

四、三弯矩法——待定二阶系数

五、B样条基函数法