插值法

1、问题提出：生产实践中往往出现的问题是，给出一批离散样点，要求作出一条通过这些点的光滑曲线，以便满足设计要求或者进行加工，反映在数学上，即已知函数在一些点上的值，寻求它的分析表达式。因为由函数的表格形式不能直接得出表中未列点处的函数值，也不便于研究函数的性质；此外虽然有些函数有表达式，但因式复杂，不易计算和进行理论分析，也需要构造一个简单函数来近似它。
解决方法：
（1）插值法——给出函数f(x)的一些样点值，选定一个便于计算的函数形式，如多项式、分式线性函数及三角多项式等，要求它通过已知样点，由此确定函数φ(x)作为f(x)的近似。用多项式作为研究插值的工具称为代数插值。
（2）曲线（数据）拟合法——在选定近似函数的形式后，不要求近似函数过已知样点，只要求在某种意义下这些点上的总偏差最小。

2、拉格朗日（Lagrange）插值：

• 多项式插值——已知函数f(x)在区间[a,b]上n+1个不同点x_0,x_1,…,x_n处的函数值y_i=f(x_i)（i=0,1,…,n），求一个至多n次的多项式
  
  使其在给定点处与f(x)同值，即满足插值条件
  
  其中φ_n(x)称为插值多项式，x_i称为插值节点，简称节点，[a,b]称为插值区间。
  
• 误差估计——插值多项式与被插函数之间的差称为截断误差，又称插值余项。
  
  设x_0,x_1,…,x_n是区间[a,b]上的互异节点，φ_n(x)是通过这组节点的n次插值多项式，如果f(x)在[a,b]上n+1次连续可导，则对[a,b]内任意点x，插值余项为
  
• Lagrange插值多项式——
  
  优缺点——由于插值基函数只与节点有关而与函数值无关，因此当插值节点相同而函数值不同时，所有Lagrange插值基函数均不变,此时用Lagrange插值多项式方便；当新增插值节点时，用Lagrange插值多项式需要重新计算所有插值基函数，计算量大且应用不方便。

3、牛顿（Newton）插值：

• 差商——设有函数f(x)，x_0, x_1, x_2,…为一系列互不相等的点，称下式为f(x)关于点x_i, x_j的一阶差商，也称均差。
  
  f(x)关于点x_i,x_j,x_k的二阶差商，记为f[x_i, x_j, x_k]，表示如下
  
  f(x)关于点x_0,x_1,…,x_k的k阶差商，表示如下
  
  差商的性质：
  (1) 各阶差商均具有线性性；
  (2) k阶差商f[x_0, x_1, … , x_k]可表示为f(x_0)，f(x_1)，f(x_k)的线性组合；
  (3) 各阶差商均具有对称性，改变节点的位置，差商值不变；
  (4) 若f(x)是n次多项式，则一阶差商f[x,x_i]是n-1次多项式。
• Newton插值公式——
  
  其中，N_n(x)是至多n次的多项式
  优点：每增加一个节点，插值多项式只增加一项。
  
• 由插值多项式的唯一性知，n阶Newton插值多项式和n阶Lagrange插值多项式是一样的，只是表现形式不同。
  

4、Runge现象：当节点无限加密时，L_n(x)只能在很小的范围内收敛。这表明通过增加节点来提高逼近程度是不合适的。

5、分段线性插值：为提高精度，在加密节点时，可以把节点分成若干段，分段用低次多项式近似函数，这是分段插值的思想。用折线近似曲线，相当于分段用线性插值，称为分段线性插值。常用的分段低次插值有：分段线性插值、分段三次Hermite插值、三次样条插值。

• 定义——
  
• φ（x）的分段表达式——
  
  
• 误差估计——如果f(x)在[a,b]上二阶连续可微，则分段线性插值函数φ(x)的余项有如下估计
  
  其中
• 优缺点——简单易行，当节点加密时，分段线性插值的误差变小，收敛性有保证，另一方面，每个小区间上的插值函数只依赖于本段节点值，因而每个节点只影响到节点邻近的一、二个小区间，计算过程中数据误差基本上不扩大，保证了节点数增加时插值过程的稳定性；但是在节点处插值函数不可微，没有光滑性。

6、埃尔米特（Hermite）插值：对于插值函数，不仅要求它在节点处与函数同值，而且要求它与函数有相同的一阶、二阶甚至更高阶的导数值。

• 定义——设已知函数y=f(x)在n+1个互异节点x_0,x_1,…,x_n上的函数值y_i=f(x_i)，和导数值y_i’=f’(x_i)，要求一个至多2n+1次的多项式H(x)，使得
  
  满足上述条件的多项式H(x)]称为Hermite插值多项式。
• 插值余项和误差估计——
  
• 两个节点的三次Hermite插值多项式的截断误差
  

7、样条插值：许多工程技术中提出的计算问题对插值函数的光滑性有较高的要求，如飞机的机翼外形，内燃机的进、排气门的凸轮曲线，都要求曲线具有较高的光滑程度，不仅要连续，而且要有连续的曲率，即二阶导数连续。

• 定义——
  
• 三类边界条件——