MathLabs
定理已证明

自然三次样条的存在性与唯一性

命题陈述

给定满足 x0<x1<⋯<xnx_0<x_1<\dots<x_n 的 n+1n+1 个点 (x0,y0),…,(xn,yn)(x_0,y_0),\dots,(x_n,y_n),存在唯一函数 SS,它在每个子区间 [xi,xi+1][x_i,x_{i+1}] 上是三次多项式,在整个 [x0,xn][x_0,x_n] 上二次连续可微,对每个 ii 满足 S(xi)=yiS(x_i)=y_i,并满足自然边界条件 S′′(x0)=S′′(xn)=0S''(x_0)=S''(x_n)=0。

为什么成立?

经过多个点的单一高次多项式往往会剧烈摆动,但将许多平缓的三次片段拼接起来、只要求它们在接缝处平滑衔接,就恰好提供了足够的自由度来拟合数据而不产生剧烈振荡,而自然边界条件恰好提供了使整个方程组可解所需的两个额外方程。

证明思路

第一步(未知量:各节点处的二阶导数)。设 i=0,…,ni=0,\dots,n 时 Mi=S′′(xi)M_i = S''(x_i)。自然边界条件立即固定 M0=0M_0=0 与 Mn=0M_n=0,剩下 n−1n-1 个未知量 M1,…,Mn−1M_1,\dots,M_{n-1} 待定。

第二步(由M值重构每个三次段)。在 [xi,xi+1][x_i,x_{i+1}] 上,由于 SS 是三次的,故 S′′S'' 是线性的,必为经过 (xi,Mi)(x_i,M_i) 与 (xi+1,Mi+1)(x_{i+1},M_{i+1}) 的直线。将该线性函数积分两次,并利用 S(xi)=yiS(x_i)=y_i 与 S(xi+1)=yi+1S(x_{i+1})=y_{i+1} 固定两个积分常数,便完全确定了该段上的 SS,用 Mi,Mi+1,yi,yi+1M_i, M_{i+1}, y_i, y_{i+1} 及间距 hi=xi+1−xih_i=x_{i+1}-x_i 表示。因此一旦所有 MiM_i 已知,SS 便完全确定。

第三步(斜率匹配给出线性方程组)。按构造,SS 与 S′′S'' 在每个节点处已经连续。要求一阶导数 S′S' 在每个内部节点 xix_i(i=1,…,n−1i=1,\dots,n-1)处从两侧也一致,便在每个内部节点处产生一个关联三个连续未知量的线性方程:hi−1Mi−1+2(hi−1+hi)Mi+hiMi+1=6(yi+1−yihi−yi−yi−1hi−1)h_{i-1} M_{i-1} + 2(h_{i-1}+h_i) M_i + h_i M_{i+1} = 6\left(\dfrac{y_{i+1}-y_i}{h_i} - \dfrac{y_i-y_{i-1}}{h_{i-1}}\right)。这给出关于 n−1n-1 个未知量 M1,…,Mn−1M_1,\dots,M_{n-1}(利用 M0=Mn=0M_0=M_n=0)的 n−1n-1 个线性方程。

第四步(方程组有唯一解)。该方程组的系数矩阵是三对角矩阵,对角元为 2(hi−1+hi)2(h_{i-1}+h_i),非对角元为 hi−1h_{i-1} 与 hih_i;由于所有间距 hi>0h_i>0,故 2(hi−1+hi)>hi−1+hi2(h_{i-1}+h_i) > h_{i-1}+h_i,该矩阵严格对角占优,而严格对角占优矩阵总是可逆的。因此关于 M1,…,Mn−1M_1,\dots,M_{n-1} 的线性方程组恰有一个解,由第二步这唯一确定了一条样条 SS,从而同时证明了存在性与唯一性。

用到此定理的主题

分步证明

该定理暂无分步证明。