MathLabs
定理証明済み

自然3次スプラインの存在と一意性

内容

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) が与えられたとき、各小区間 [xi,xi+1][x_i,x_{i+1}] 上で3次多項式であり、[x0,xn][x_0,x_n] 全体で2回連続微分可能であり、すべての ii について S(xi)=yiS(x_i)=y_i を満たし、自然境界条件 S′′(x0)=S′′(xn)=0S''(x_0)=S''(x_n)=0 を満たす唯一の関数 SS が存在する。

なぜ正しいのか?

多くの点を通る単一の高次多項式は波打つ傾向があるが、多くの緩やかな3次片をつなぎ合わせ、継ぎ目で滑らかに一致することだけを要求すると、激しい振動なしにデータに適合するのにちょうど十分な自由度が得られ、自然境界条件は系全体を解けるようにするために必要なちょうど2つの追加の方程式を与える。

証明の概略

ステップ1(未知数:各節点での2階導関数)。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} を決定すればよい。

ステップ2(各3次片をM値から再構成する)。[xi,xi+1][x_i,x_{i+1}] 上では、SS が3次なので S′′S'' は線形であり、(xi,Mi)(x_i,M_i) と (xi+1,Mi+1)(x_{i+1},M_{i+1}) を通る直線でなければならない。この線形関数を2回積分し、S(xi)=yiS(x_i)=y_i と S(xi+1)=yi+1S(x_{i+1})=y_{i+1} を用いて2つの積分定数を固定すると、その区間上の 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 は完全に決まる。

ステップ3(傾きを一致させると線形方程式系が得られる)。構成により SS と S′′S'' はすでに各節点で連続である。各内部節点 xix_i(i=1,…,n−1i=1,\dots,n-1)で1階導関数 S′S' も両側から一致することを要求すると、内部節点ごとに連続する3つの未知数を関連づける線形方程式が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 個の線形方程式が得られる。

ステップ4(系は唯一の解を持つ)。この系の係数行列は三重対角であり、対角成分は 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} についての線形系はちょうど1つの解を持ち、ステップ2によりこれはちょうど1つのスプライン SS を定めるので、存在と一意性の両方が証明される。

この定理を使うトピック

ステップごとの証明

この定理のステップごとの証明はまだありません。