入门级理解插值方法

插值方法

1 Hermite插值

1.1 Hermite插值的定义

设在节点$a\le x_0<x_1<...<x_n\le b(i=1,2,...,n)$上给出插值条件$$y_k=f(x_k),y_k^{‘}=f^{‘}(x_k),k=0,1,...,n$$此时，在小区间$[x_{i-1},x_i]$上有四个插值条件，所以能构造一个三次多项式$H_3^{(i)}(x)$，只要令$$H_3^{(i)}(x)={\varphi }_{i-1}(x)y_{i-1}+{\varphi }_i(x)y_i+{\psi }_{i-1}(x)y_{i-1}^{‘}+{\psi }_i(x)y_i^{‘}$$其中${\varphi }_{i-1}(x),{\varphi }_i(x),{\psi }_{i-1},{\psi }_i(x)$都是三次多项式，而且满足
$${\varphi }_{i-1}(x_{i-1})=1,{\varphi }_{i-1}(x_i)=0,{\varphi }_{i-1}^{‘}(x_{i-1})=0,{\varphi }_{i-1}^{‘}(x_i)=0$$$${\varphi }_i(x_{i-1})=0,{\varphi }_i(x_i)=1,{\varphi }_i^{‘}(x_{i-1})=0,{\varphi }_i^{‘}(x_i)=0$$$${\psi }_{i-1}(x_{i-1})=0,{\psi }_{i-1}(x_i)=0,{\psi }_{i-1}^{‘}(x_{i-1})=1,{\psi }_{i-1}^{‘}(x_i)=0$$$${\psi }_i(x_{i-1})=0,{\psi }_i(x_i)=0,{\psi }_i^{‘}(x_{i-1})=0,{\psi }_i^{‘}(x_i)=1$$把${\varphi }_{i-1}(x),{\varphi }_i(x),{\psi }_{i-1}(x),{\psi }_i(x)$称为三次Hermite插值基函数。

1.2 Hermite插值的求法

因为${\varphi }_{i-1}(x_i)=0,{\varphi }_{i-1}^{‘}(x_i)=0$，所以可以将${\varphi }_{i-1}(x)$写成$${\varphi }_{i-1}(x)=(ax+b)(x-x_i{)}^2$$将${\varphi }_{i-1}(x_{i-1})=1,{\varphi }_{i-1}^{‘}(x_{i-1})=0$带入可得：$${\varphi }_{i-1}(x)=\frac{1}{h_i^3}(2x+x_i-3x_{i-1})(x-x_i{)}^2$$$${\varphi }_i(x)=\frac{1}{h_i^3}(3x_i-2x-x_{i-1})(x-x_{i-1}{)}^2$$

因为${\psi }_{i-1}(x_{i-1})=0,{\psi }_{i-1}(x_i)=0,{\psi }_{i-1}^{‘}(x_i)=0$所以有$${\psi }_{i-1}(x)=C(x-x_{i-1})(x-x_i{)}^2$$将${\psi }_{i-1}^{‘}(x_{i-1})=1$带入得$C{h}_i^2=1$，即$C=h_i^{-2}$。所以$${\psi }_{i-1}(x)=\frac{1}{h_i^2}(x-x_{i-1})(x-x_i{)}^2$$$${\psi }_i(x)=\frac{1}{h_i^2}(x-x_{i-1}{)}^2(x-x_i)$$

2 三次样条插值

2.1 三次样条插值的定义

定义：给定节点$a=x_0<x_1<...<x_n=b$，及其上的函数值$y_k=f(x_k),k=0,1,...,n$.如果函数满足：
（1）S(x)是一个分段的三次多项式且$S(x_k)=y_k$；
（2）$S(x)\in C^2[a,b]$.
则称S(x)是区间[a,b]上的三次样条插值函数。
S(x)在区间$[x_{i-1},x_i]$上是三次多项式，$S(x)=a_i{x}^3+b_i{x}^2+c_{i}x+d_i$有四个待定系数，所以在整个区间上总的来说，$S(x)$共有4n个待定系数。
$$S(x_j-0)=S(x_j+0)$$$$S^{‘}(x_j-0)=S^{‘}(x_j+0)$$$$S^{‘’}(x_j-0)=S^{‘’}(x_j+0)$$共有$3n-3$个条件，再加上$S(x)$满足$S(x_j)=f_j,j=0,1,2,...n$这又有$n+1$个条件。目前，共有$4n-2$个条件。
因此，还需要两个条件才能确定$S(x)$。
通常可在$[a,b]$端点$a=x_0,b=x_n$上各加上一个条件，称为边界条件。

2.2 边界条件

常见的边界条件有：
1）给定两端一阶导数值，$S‘(x_0)=y{’}_0$，$S‘(x_n)=y{’}_n$；
2）给定两端二阶导数值，$S’’(x_0)=y’{’}_0$，$S’’(x_n)=y’{’}_n$。特殊的有$S’’(x_0)=S’’(x_n)=0$，称为自然边界条件；
3）假设$f(x)$是以$b-a$为周期的周期函数，这时要求$$S(x_0+0)=S(x_n-0)$$$$S’(x_0+0)=S’(x_n-0)$$$$S’’(x_0+0)=S’’(x_n-0)$$此时，$y_n=y_0$，这样确定的$S(x)$为周期样条函数。

2.3 三次样条函数的求法

2.3.1 三转角方法

当给定左右端点的一阶导数值时，假定$S’(x_j)=m_j(j=0,1,2,...,n)$，再利用三次Hermite插值的一半表达式可知：$$S(x)={\Sigma }_{j=0}^n[y_j{\alpha }_j(x)+m_j{\beta }_j(x)]$$，其中${\alpha }_j(x),{\beta }_j(x)$为分段三次Hermite插值基函数，再利用三次样条函数的定义及其相应的边界条件$$S’’(x_0)=f{’}_0,S’(x_n)=f{’}_n$$当$x\in [x_{i-1},x_i]$时，有$$\begin{gathered} S(x)=\frac{1}{h_i^3}[(x_i+2x-3x_{i-1})(x-x_i{)}^2{y}_{i-1}+(3x_i-2x-x_i)(x-x_{i-1}{)}^2{y}_i] \\ +\frac{1}{h_i^2}[(x-x_{i-1})(x-x_i{)}^2{m}_{i-1}+(x-x_{i-1}{)}^2(x-x_i)m_i] \end{gathered}$$其中$h_i=x_i-x_{i-1}$，为了确定$S(x)$，只需要确定$m_i,i=0,1,...,n$可利用三次样条函数二阶导数连续的性质，即$S’’(x_i-0)=S’’(x_i+0)$。
当$x\in [x_{i-1},x_i]$时，由于 S ’ ( x ) = 2 h i 3 { [ ( x − x i ) 2 + ( x i + 2 x − 3 x i − 1 ) ( x − x i ) ] y i − 1 + [ ( 3 x i − 2 x − x i ) ( x − x i − 1 ) − 2 ( x − x i − 1 ) 2 ] y i } + 1 h i 2 { [ ( x − x i ) 2 + 2 ( x − x i − 1 ) ( x − x i ) ] m i − 1 + [ ( x − x i − 1 ) 2 + 2 ( x − x i − 1 ) ( x − x i ) ] m i } S’(x)=\frac{2}{h_i^3}\{[(x-x_i)^2+(x_i+2x-3x_{i-1})(x-x_i)]y_{i-1} \\ +[(3x_i-2x-x_i)(x-x_{i-1})-2(x-x_{i-1})^2]y_i\} \\ +\frac{1}{h_i^2}\{[(x-x_i)^2+2(x-x_{i-1})(x-x_i)]m_{i-1} \\ +[(x-x_{i-1})^2+2(x-x_{i-1})(x-x_i)]m_i\} S’(x)=hi3​2​{[(x−xi​)2+(xi​+2x−3xi−1​)(x−xi​)]yi−1​+[(3xi​−2x−xi​)(x−xi−1​)−2(x−xi−1​)2]yi​}+hi2​1​{[(x−xi​)2+2(x−xi−1​)(x−xi​)]mi−1​+[(x−xi−1​)2+2(x−xi−1​)(x−xi​)]mi​}
$$\begin{gathered} S’’(x)=\frac{6}{h_i^3}(2x-x_{i-1}-x_i)(y_{i-1}-y_i) \\ +\frac{2}{h_i^2}[(3x-x_{i-1}-2x_i)m_{i-1}+(3x-2x_{i-1}-x_i)m_i] \end{gathered}$$可以看出$S’’(x)$是一个线性函数的形式。
在$x_i$点处考虑它的左极限和右极限：
$$S’’(x_i-0)=\frac{6}{h_i^2}(y_{i-1}-y_i)+\frac{2}{h_i}(m_{i-1}+2m_i)$$$$S’’(x_i+0)=\frac{6}{h_{i+1}^2}(y_{i+1}-y_i)-\frac{2}{h_{i+1}}(2m_i+m_{i+1})$$由连续性条件$S’’(x_i-0)=S’’(x_i+0)$可得：$$\frac{1}{h_i}{m}_{i-1}+2\left(\frac{1}{h_i}+\frac{1}{h_{i+1}}\right)m_i+\frac{1}{h_{i+1}}{m}_{i+1}=3\left(\frac{y_{i+1}-y_i}{h_{i+1}^2}+\frac{y_i-y_{i-1}}{h_i^2}\right)$$两侧同时除以$\left(\frac{1}{h_i}+\frac{1}{h_{i+1}}\right)$，并记$\frac{h_{i+1}}{h_i+h_{i+1}}={\lambda }_i$，$\frac{h_i}{h_i+h_{i+1}}=1-{\lambda }_i={\mu }_i$，$3({\lambda }_{i}f[x_{i-1},x_i]+{\mu }_{i}f[x_i,x_{i+1}])=g_i$，则有$${\lambda }_i{m}_{i-1}+2m_i+{\mu }_i{m}_{i+1}=g_i,i=1,2,...,n-1.$$这个式子中未知的是$m_i,m_i,m_{i+1}$
  若边界条件为：$m_0=y{’}_0,m_n=y{’}_n$，带入上面的式子可以得到
可以看出这个形式就变成了$n-1$个未知数和$n-1$个方程。左边是一个稀疏的三对角矩阵，求解方式有很多。
  若边界条件为：$S’’(x_0)=y’{’}_0,S’’(x_n)=y’{’}_n$则有：$$2m_0+m_1=3f[x_0,x_1]-\frac{1}{2}{h}_{1}y’{’}_0=g_0$$$$m_{n-1}+2m_n=3f[x_{n-1},x_n]+\frac{1}{2}{h}_{n}y’{’}_n=g_n$$将这两个方程和上式放在一起就得到了$n+1$个未知数和$n+1$个方程。

  若边界条件为周期性边界条件，由$S’(x_0+0)=S’(x_n-0)$，$S’’(x_0+0)=S’’(x_n-0)$，有$m_0=m_n$和${\lambda }_n{m}_{n-1}+2m_n+{\mu }_n{m}_1=g_n$，${\lambda }_n=\frac{h_1}{h_1+h_n}$，${\mu }_n=1-{\lambda }_n=\frac{h_n}{h_1+h_n}$，$g_n=3({\lambda }_{n}f[x_0,x_1]+{\mu }_{n}f[x_{n-1},x_n])$

所以是一个n个未知数n个方程的线性方程组。
在给定第一类边界条件的情况下，三对角方法很好。8

2.3.2 三弯矩算法

  由于三次样条函数在每个小区间上都是一个三次多项式，所以它的二阶导数在整个区间上是一个分段的线性函数，在每个小区间上是一个线性函数。设$S’’(x_i)=M_i,i=0,1,...,n$，则对$x\in [x_{i-1},x_i]$有：$$S’’(x)=\frac{x-x_i}{x_{i-1}-x_i}{M}_{i-1}+\frac{x-x_{i-1}}{x_i-x_{i-1}}{M}_i$$连续积分两次，并利用$S(x_{i-1})=y_{i-1}$，$S(x_i)=y_i$确定积分常数，可得$$\begin{gathered} S(x)=\frac{1}{6h_i}[(x_i-x{)}^3{M}_{i-1}+(x-x_{i-1}{)}^3{M}_i] \\ +\left(\frac{y_{i-1}}{h_i}-\frac{h_i{M}_{i-1}}{6}\right)(x_i-x)+\left(\frac{y_i}{h_i}-\frac{h_i{M}_i}{6}\right)(x-x_{i-1}) \end{gathered}$$其中$h_i=x_i-x_{i-1}$，为了确定$S(x)$，只需要确定$M(i),i=0,1,...,n$，可利用$S’(x_i-0)=S’(x_i+0)$来求出$M_i$。则有：$$S’(x)=\frac{1}{2h_i}[(x-x_{i-1}{)}^2{M}_i-(x-x_i{)}^2{M}_{i-1}]+f[x_{i-1},x_i]+\frac{h_i}{6}(M_{i-1}-M_i)$$于是有$$S’(x_i-0)=\frac{h_i}{6}(M_{i-1}+2M_i)+f[x_{i-1},x_i]$$$$S’(x_i+0)=-\frac{h_{i+1}}{6}(2M_i+M_{i+1})+f[x_i,x_{i+1}]$$$$\frac{h_i}{6}{M}_{i-1}+\frac{h_i+h_{i+1}}{3}{M}_i+\frac{h_{i+1}}{6}{M}_{i+1}=f[x_i,x_{i+1}]-f[x_{i-1},x_i]$$引入记号：
${\lambda }_i=\frac{h_{i+1}}{h_i+h_{i+1}},{\mu }_i=\frac{h_i}{h_i-h_{i+1}},d_i=6f[x_{i-1},x_i,x_{i+1}]$$${\mu }_i{M}_{i-1}+2M_i+{\lambda }_i{M}_{i+1}=d_i,i=1,2,...,n-1$$
假如我们有边界条件为：$M_0=y’{’}_0,M_n=y’{’}_n$，可得

假如我们有边界条件为：$S’(x_0)=y{’}_0,S’(x_n)=y{’}_n$，则有$$2M_0+M_1=\frac{6}{h_1}(f[x_0,x_1]-y{’}_0)=d_0$$$$M_{n-1}+2M_n=\frac{6}{h_n}(y{’}_n-f[x_{n-1},x_n])=d_n$$

假如边界条件为周期性边界条件，由$S’(x_0+0)=S’(x_n-0)$和$S’’(x_0+0)=S’’(x_n-0)$，有$$M_0=M_n$$和$${\lambda }_n{M}_1+{\mu }_n{M}_{n-1}+2M_n=d_n$$其中${\lambda }_n=\frac{h_1}{h_1+h_n},{\mu }_n=1-{\lambda }_n=\frac{h_n}{h_1+h_n},d_n=6f[x_0,x_1.x_{n-1}]$。
所以针对给定第二类边界条件的情况下，三弯矩方法很好。