分段三次样条

一、插值和样条？

给定若干已知点位，连接、逼近这些点位的方式有很多种，我们既可以给定一个足够阶数的多项式一次性穿过、逼近这些点位，这些方法可以是：牛顿插值法、拉格朗日插值法；也可以按段给出每个区间的表达式，比如说三次样条。

上面展示的就是按段给出的表达式，按段一条条给出的，所以叫做样条。当然你也可以用一个多项式穿过这些点：

红色的线就是用一个多项式连接的方法。用一个多项式拟合所有的点位不都是好的，尤其是均匀分布的插值点两端因为它会产生龙格现象：

从上图可以看出，左右两端的拟合对于原函数完全不匹配。缓解方法有很多种，可以使用切比雪夫或者分段三次样条进行处理。也就是文章的重点，分段三次多项式样条。

二、分段三次样条（自然边界条件）

样条的表达式是通过多个分段函数给出的，约束条件有：

• 满足节点自变量递增关系；
• 相邻多项式之间一二阶导数连续；
• 多项式的次数是3次；

如果每相邻两个点之间的连线是三次多项式，那么对应的整段组合而成的曲线就叫做三次样条函数。假设点数是 $n+1$，那么段数为：$n$

• 三次多项式有四个未知数，一共有 $n$ 个，所以待求的方程组个数为 $4n$；
• 一共有 $n+1$ 个点，我已有 $(n+1)$ 个方程组，这样一来我就剩下 $4n-(n+1)=3n-1$；
• 一共有 $n-1$ 个衔接点，衔接点处函数值连续、一阶导连续，二阶导连续，提供 $3(n-1)$个方程，剩余$3n-1-3(n-1)=2$

根据不同的约束可以将三次样条分为三类：

• 端点处二阶导函数值约束，如果都是0，则为自然样条；
• 端点处一阶导数约束；
• 周期约束，略；

给定以下插值点：

$x$	$t_0$	$t_1$	$\cdots$	$t_n$
$y$	$y_0$	$y_1$	$\cdots$	$y_n$

根据这个表格结合三次样条多项式的定义：
$$S(x)=\{\begin{array}{rl}& S_0(x)x\in [t_0,t_1]\\ & S_1(x)x\in [t_1,t_2]\\ & ⋮\\ & S_{n-1}(x)x\in [t_{n-1},t_n]\end{array}$$

在此之前我们先来推导区间 $[t_i,t_{i+1}]$ 上 $S_i(x)$ 的表达式。首先我们定义一组数 $z_i=S^{\prime \prime }(t_i)$，因为 $S^{\prime \prime }$ 在每个内结点上连续，所以显然，对于 $0\le i\le n$ 存在 $z_i$ 满足：
$$\underset{x\to t_i^{+}}{\lim }{S}^{\prime \prime }(x)=\underset{x\to t_i^{-}}{\lim }{S}^{\prime \prime }(x)$$
因为 $S_i$ 是在 $[t_i,t_{i+1}]$ 上的三次多项式，因此是满足 $S_i^{\prime \prime }(t_i)=z_i$、 $S_i^{\prime \prime }(t_{i+1})=z_{i+1}$ 的线性函数，根据根据拉格朗日插值法有：
$$S_i^{\prime \prime }(x)=\frac{z_i}{h_i}(t_{i+1}-x)+\frac{z_{i+1}}{h_i}(x-t_i)$$

其中 $h_i=t_{i+1}-$