数学建模-插值算法原理笔记

这里是根据清风数学建模视频课程记录的笔记，我不是清风本人。想系统学习数学建模的可以移步B站搜索相关视频

目的

比赛中常常需要根据已知的函数点进行数据、模型的处理和分析，而有时候现有的数据是极少的，不足以支撑分析的进行，这时就可以使用一些方法“模拟产生”一些新的但又比较靠谱的值，这就是插值的作用。

概念

设函数$y=f(x)$​在区间$[a,b]$​上有定义，且已知在点$a\le x_0<x_1<\dots <x_n\le b$​上的值分别为$y_0,y_1,\dots ,y_n$​，若存在一简单函数$P(x)$​使$P(x_i)=y_i(i=0,1,2,\dots ,n)$则称$P(x)$为$f(x)$的插值函数，点$x_0,x_1,\dots ,x_n$称为插值节点，包含插值节点的区间$[a,b]$称为插值区间，求插值函数$P(x)$的方法称为插值法

分类

• 插值多项式
  • $P(x)$是次数不超过$n$​的代数多项式，即$P(x)=a_0+a_{1}x+\dots +a_n{x}^n$
• 分段插值
  • $P(x)$为分段多项式
• 三角插值

一般插值多项式

定理：设有$n+1$​个互不相同的节点$(x_i,y_i),(i=0,1,2\dots ,n)$​，则存在唯一的多项式$L_n(x)=a_0+a_{1}x+a_2{x}^2+\dots +a_n{x}^n$​，使得$L_n(x_j)=y_j,(j=0,1,2\dots ,n)$​

注意：①只要$n+1$个节点互异，满足上述插值条件的多项式是唯一存在的；②如果不限制多项式次数，插值多项式并不唯一

拉格朗日插值法

• 两个点$(x_0,y_0),(x_1,y_1)$

  • $f(x)=\frac{x-x_1}{x_0-x_1}{y}_0+\frac{x-x_0}{x_1-x_0}{y}_1$​

• 三个点$(x_0,y_0),(x_1,y_1),(x_2,y_2)$

  • $f(x)=\frac{(x-x_1)(x-x_2)}{(x_0-x_1)(x_0-x_2)}{y}_0+\frac{(x-x_0)(x-x_2)}{(x_1-x_0)(x_1-x_2)}{y}_1+\frac{(x-x_0)(x-x_1)}{(x_2-x_1)(x_2-x_1)}{y}_2$

• 四个点$(x_0,y_0),(x_1,y_1),(x_2,y_2),(x_3,y_3)$

  • $f(x)=\frac{(x-x_1)(x-x_2)(x-x_3)}{(x_0-x_1)(x_0-x_2)(x_0-x_3)}{y}_0+\frac{(x-x_0)(x-x_2)(x-x_3)}{(x_1-x_0)(x_1-x_2)(x_0-x_3)}{y}_1+\frac{(x-x_0)(x-x_1)(x-x_3)}{(x_2-x_1)(x_2-x_1)(x_0-x_3)}{y}_2$

分段线性插值

分段二次插值

牛顿插值法

定义$f[x_0,x_k]=\frac{f(x_k)-f(x_0)}{x_k-x_0}$为函数$f(x)$关于点$x_0,x_k$的一阶差商，$f[x_0,x_1,x_2]=\frac{f[x_1,x_2]-f[x_0,x_1]}{x_2-x_0}$为二阶差商

牛顿插值：$f(x)=f(x_0)+f[x_0,x_1]\ast (x-x_0)+f[x_0,x_1,x_2]\ast (x-x_0)\ast (x-x_1)+\dots +f[x_0,x_1,\dots x_{n-2},x_{n-1}]\ast (x-x_0)\ast (x-x_1)\ast \dots (x-x_{n-3})(x-x_{n-2})+f[x_0,x_1,\dots x_{n-1},x_n]\ast (x-x_0)\ast (x-x_1)\ast \dots (x-x_{n-2})(x-x_{n-1})$

埃尔米特插值原理

设函数$f(x)$在区间$[a,b]$上有$n+1$个互异节点$a=x_0<x_1<x_2<\dots <x_n=b$，定义在$[a,b]$上函数$f(x)$在节点上满足$f(x_i)=y_i,f^{\prime }(x_i)=y_i^{\prime }$可唯一确定一个次数不超过$2n+1$的多项式$H_{2n+1}(x)=H(x)$满足$H(x_j)=y_j,H^{\prime }(x_j)=m_j,(j=0,1,\dots ,n)$，其余项为$R(x)=f(x)-H(x)=\frac{f^{(2n+2)}(\xi )}{(2n+2)!}{\omega }_{2n+2}(x)$

• 不但要求在节点上的函数值相等，而且还要求对应的导数值也相等

分段三次埃米尔特插值

见代码

三次样条插值

除了满足基本插值条件$S(x_i)=f_i$还应该满足

$S(x)=\{\begin{array}{l}{S}_0(x),x\in [x_0,x_1]\\ S_1(x),x\in [x_1,x_2]\\ ⋮\\ S_{n-1}(x),x\in [x_{n-1},x_n]\end{array}$

并且满足

$\{\begin{array}{l}{S}_{i-1}(x_i)=S_i(x_i),\\ S_{i-1}^{\prime }(x_i)=S_i^{\prime }(x_i),\\ S_{i-1}^{\prime \prime }(x_i)=S_i^{\prime \prime }(x_i),\end{array}$

具体见代码

注意：实际比赛中多用分段三次和三次样条