数值分析学习笔记——三次样条插值

2.6 三次样条插值

1.三次样条函数

（1）定义：若函数$S(x)\in C^2[a,b]$，且每个小区间$[x_j,x_{j+1}]$上为三次多项式，其中$a=x_0<x_1<\cdots <x_n=b$为给定节点，则称$S(x)$为节点$x_0,x_1,\cdots ,x_n$的三次样条函数

​ 若在节点$x_j$上给定函数值$y_j=f(x_j)(j=0,1,\cdots ,n)$并满足
$$S(x_j)=y_j,j=0,1,\cdots ,n\text{\hspace{1em}(1.46)}$$
称$S(x)$为三次样条插值函数

（2）定值条件：要求$S(x)$，在每个小区间$[x_j,x_{j+1}]$上要确定4个待定系数，共有$n$个区间，故应确定$4n$个参数。根据$S(x)$在$[a,b]$上二阶导数连续，在节点$x_j(j=1,2,\cdots ,n-1)$满足连续性条件
$$S(x_j-0)=S(x_j+0),S^{\prime }(x_j-0)=S^{\prime }(x_j+0),S^{\prime \prime }(x_j-0)=S^{\prime \prime }(x_j+0)\text{\hspace{1em}(1.47)}$$
这里共有$3n-3$个条件，加上$S(x)$满足（1.46）这$n+1$个条件，再加上边界条件区间$[a,b]$左右端点$a=x_0,b=x_n$，刚好满足$4n$个条件，三次样条函数问题可解，边界条件的三种常见情况：

​ 1）已知区间两端的一阶导数值，即
$$S^{\prime }(x_0)=f_0^{\prime },S^{\prime }(x_n)=f_n^{\prime }\text{\hspace{1em}(1.48)}$$
​ 2）已知区间两端的二阶导数值，即
$$S^{\prime \prime }(x_0)=f_0^{\prime \prime },S^{\prime \prime }(x_n)=f_n^{\prime \prime }\text{\hspace{1em}(1.49)}$$
​ 特殊情况（自然边界条件）
$$S^{\prime \prime }(x_0)=S^{\prime \prime }(x_n)=0\text{\hspace{1em}(1.49-C)}$$
​ 3）当$f(x)$以$x_n-x_0$为周期的周期函数时，则要求$S(x)$也是周期函数，满足
$$\{\begin{array}{l}S(x_0+0)=S(x_n-0),S^{\prime }(x_0+0)=S^{\prime }(x_n-0)\\ S^{\prime \prime }(x_0+0)=S^{\prime \prime }(x_n-0)\end{array}\text{\hspace{1em}(1.50)}$$
​ 此时满足（1.46）式中的$y_0=y_n$，称这类样条函数$S(x)$为周期样条函数

2.三次样条插值函数

推导建立插值函数的过程：

​ 由二阶导数值$S^{\prime \prime }(x_j)=M_j,j=0,1,\cdots ,n$表达$S(x)$.由$S(x)$在区间$[x_j,x_{j+1}]$为三次多项式得$S^{\prime \prime }(x)$为线性函数，写为
$$S^{\prime \prime }(x)=M_j\frac{x_{j+1}-x}{h_j}+M_{j+1}\frac{x-x_j}{h_j}\text{\hspace{1em}(1.51)}$$
对$S(x)$进行积分并利用$S(x_j)=y_j,S(x_{j+1})=y_{j+1}$条件得到三次样条表达式
$$\begin{gathered} S(x)=M_j\frac{(x_{j+1}-x{)}^3}{6h_j}+M_{j+1}\frac{(x-x_j{)}^3}{6h_j}+(y_j-\frac{M_j{h}_j^2}{6})\frac{x_{j+1}-x}{h_j} \\ +(y_{j+1}-\frac{M_{j+1}{h}_j^2}{6})\frac{x-x_j}{h_j},j=0,1,\cdots ,n-1\text{\hspace{1em}(1.52)} \end{gathered}$$
求导得
$$S^{\prime }(x)=-M_j\frac{(x_{j+1}-x{)}^2}{2h_j}+M_{j+1}\frac{(x-x_j{)}^2}{2h_j}+\frac{y_{j+1}-y_j}{2h_j}-\frac{M_{j+1}-M_j}{6}{h}_j\text{\hspace{1em}(1.53)}$$
由此得到
$$\begin{array}{rl}{S}^{\prime }(x_j+0)& =-\frac{h_j}{3}{M}_j-\frac{h_j}{6}{M}_{j+1}+\frac{y_{j+1}-y_j}{h_j}\\ S^{\prime }(x_j-0)& =\frac{h_{j-1}}{6}{M}_{j-1}+\frac{h_{j-1}}{3}{M}_j+\frac{y_j-y_{j-1}}{h_{j-1}}\end{array}$$
由$S^{\prime }(x_j-0)=S^{\prime }(x_j+0)$得
$${\mu }_j{M}_{j-1}+2M_j+{\lambda }_j{M}_{j+1}=d_j,j=1,2,\cdots ,n-1\text{\hspace{1em}(1.54)}$$
其中
$$\begin{array}{rl}{\mu }_j& =\frac{h_{j-1}}{h_{j-1}+h_j},{\lambda }_j=\frac{h_j}{h_{j-1}+h_j}\\ d_j& =\frac{6}{h_{j-1}+h_j}(f[x_j,x_{j+1}]-f[x_j,x_{j-1}])=f[x_{j-1},x_j,x_{j+1}],j=1,2,\cdots ,n-1\end{array}\text{\hspace{1em}(1.55)}$$
最后再代入边界条件即可得到矩阵形式的方程组，解出$M_k(k=j-1,j,j+1)$即可得到三次样条插值函数$S(x)$​的表达式.

比较常见的是第一类边界条件的运用，将（1.48）式代入（1.53）式可以推出两个方程
$$\begin{array}{l}2M_0+M_1=\frac{6}{h_0}(f[x_0,x_1]-f_0^{\prime })\\ M_{n-1}+2M_n=\frac{6}{h_{n-1}}(f_n^{\prime }-f[x_{n-1},x_n])\end{array}\}\text{\hspace{1em}(1.56)}$$
若令
$$\begin{array}{rl}{\lambda }_0& =1;\\ d_0& =\frac{6}{h_0}(f[x_0,x_1]-f_0^{\prime })=\frac{6}{h_0}(f[x_0,x_1]-P_0^{\prime });\\ {\mu }_n& =1;\\ d_n& =\frac{6}{h_{n-1}}(f_n^{\prime }-f[x_{n-1},x_n])=\frac{6}{h_{n-1}}(P_n^{\prime }-f[x_{n-1},x_n]).\end{array}$$
写成矩阵形式
$$\left[\begin{array}{ccccc}2& {\lambda }_0& & & \\ {\mu }_1& 2& {\lambda }_1& & \\ & \ddots & \ddots & \ddots & \\ & & {\mu }_{n-1}& 2& {\lambda }_{n-1}\\ & & & {\mu }_n& 2\end{array}\right]\left[\begin{array}{c}{M}_0\\ M_1\\ ⋮\\ M_{n-1}\\ M_n\end{array}\right]=\left[\begin{array}{c}{d}_0\\ d_1\\ ⋮\\ d_{n-1}\\ d_n\end{array}\right]\text{\hspace{1em}(1.57)}$$
利用追赶法可解方程组.