logo

B-スプラインの正則性 📂フーリエ解析

B-スプラインの正則性

定理1

m=2,3,m=2,3,\dotsに対して、B-スプライン NmN_{m}は以下の性質を持つ。

(a) NmCm2(R)N_{m}\in C^{m-2}(\mathbb{R})

(b) kZk\in \mathbb{Z}において、各区間 [k,k+1][k,k+1]で、NmN_{m}は最大でもm1m-1の次数の多項式である。

B-スプラインの明示的な式

Nm(x)=1(m1)!j=0m(1)j(mj)(xj)+m1,xR N_{m}(x) = \frac{1}{(m-1)!}\sum \limits_{j=0}^{m} \left( -1 \right)^{j}\binom{m}{j}\left( x-j \right)_{+}^{m-1},\quad x\in \mathbb{R}

この時

f(x)+:=max(0,f(x))&f(x)+n:=(f(x)+)n f(x)_{+}:=\max \left( 0,f(x) \right) \quad \& \quad f(x)_{+}^{n}:=\left( f(x)_{+} \right)^{n}

補助定理

m=2,3,m=2,3,\cdotsに対して、x+m1x_{+}^{m-1}m2m-2回まで微分可能で、m2m-2回目の導関数は連続である。

証明

m=2m=2に対して、

x+1=max(0,x)={0if x0xif x0 x_{+}^{1}=\max(0,x)=\begin{cases} 0 & \text{if}\ x\le0 \\ x & \text{if}\ x\ge0 \end{cases}

であるため、全ての点で連続であり、x=0x=0を除いて微分可能である。m=3m=3に対して、

x+2=(max(0,x))2={0if x0x2if x0 x_{+}^{2}=\left( \max(0,x) \right)^{2}=\begin{cases} 0 & \text{if}\ x\le0 \\ x^{2} & \text{if}\ x\ge0 \end{cases}

であるため、全ての点で一度微分可能である。導関数は2x+12x_{+}^{1}であるため、全ての点で連続であり、x=0x=0を除く所で微分可能である。

証明

(a)

B-スプラインの明示的な式により、NmN_{m}x+m1x_{+}^{m-1}の平行移動の線形結合であることがわかる。従って、下の補助定理により、NmN_{m}m2m-2回微分可能で、各導関数は連続である。

(b)

j=0,1,,mj=0,1,\dots,mに対して、以下の式が成立する。

(xj)+m1=(max(0,xj))m1={0if xj(xj)m1if x>j (x-j)_{+}^{m-1}=\left( \max \left( 0,x-j \right) \right)^{m-1}=\begin{cases} 0 &\text{if}\ x\le j \\ (x-j)^{m-1}& \text{if}\ x>j \end{cases}

NmN_{m}はこのような関数の線形結合であるため、自明に成立する。


  1. Ole Christensen, Functions, Spaces, and Expansions: Mathematical Tools in Physics and Engineering (2010), p208 ↩︎