定理已证明
自然三次样条的存在性与唯一性
命题陈述
给定满足 x0<x1<⋯<xn 的 n+1 个点 (x0,y0),…,(xn,yn),存在唯一函数 S,它在每个子区间 [xi,xi+1] 上是三次多项式,在整个 [x0,xn] 上二次连续可微,对每个 i 满足 S(xi)=yi,并满足自然边界条件 S′′(x0)=S′′(xn)=0。
为什么成立?
经过多个点的单一高次多项式往往会剧烈摆动,但将许多平缓的三次片段拼接起来、只要求它们在接缝处平滑衔接,就恰好提供了足够的自由度来拟合数据而不产生剧烈振荡,而自然边界条件恰好提供了使整个方程组可解所需的两个额外方程。
证明思路
第一步(未知量:各节点处的二阶导数)。设 i=0,…,n 时 Mi=S′′(xi)。自然边界条件立即固定 M0=0 与 Mn=0,剩下 n−1 个未知量 M1,…,Mn−1 待定。
第二步(由M值重构每个三次段)。在 [xi,xi+1] 上,由于 S 是三次的,故 S′′ 是线性的,必为经过 (xi,Mi) 与 (xi+1,Mi+1) 的直线。将该线性函数积分两次,并利用 S(xi)=yi 与 S(xi+1)=yi+1 固定两个积分常数,便完全确定了该段上的 S,用 Mi,Mi+1,yi,yi+1 及间距 hi=xi+1−xi 表示。因此一旦所有 Mi 已知,S 便完全确定。
第三步(斜率匹配给出线性方程组)。按构造,S 与 S′′ 在每个节点处已经连续。要求一阶导数 S′ 在每个内部节点 xi(i=1,…,n−1)处从两侧也一致,便在每个内部节点处产生一个关联三个连续未知量的线性方程:hi−1Mi−1+2(hi−1+hi)Mi+hiMi+1=6(hiyi+1−yi−hi−1yi−yi−1)。这给出关于 n−1 个未知量 M1,…,Mn−1(利用 M0=Mn=0)的 n−1 个线性方程。
第四步(方程组有唯一解)。该方程组的系数矩阵是三对角矩阵,对角元为 2(hi−1+hi),非对角元为 hi−1 与 hi;由于所有间距 hi>0,故 2(hi−1+hi)>hi−1+hi,该矩阵严格对角占优,而严格对角占优矩阵总是可逆的。因此关于 M1,…,Mn−1 的线性方程组恰有一个解,由第二步这唯一确定了一条样条 S,从而同时证明了存在性与唯一性。