logo

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

B-スプラインの性質

性質1

オーダーがm∈Nm\in \mathbb{N}のB-スプラインNmN_{m}は、次の性質を満たす。

(a) suppNm=[0,m]andNm(x)>0 for x∈(0,m)\mathrm{supp}N_{m}=[0,m] \quad \text{and} \quad N_{m}(x)>0 \text{ for } x\in(0,m)

(b) ∫−∞∞Nm(x)dx=1\displaystyle \int _{-\infty} ^{\infty} N_{m}(x)dx=1

(c) m≥2m\ge 2に対して、下の式が成立する。

∑k∈ZNm(x−k)=1,∀x∈R \begin{equation} \sum \limits_{k \in \mathbb{Z}} N_{m}(x-k)=1,\quad \forall x\in \mathbb{R} \end{equation}

(c’) m=1m=1の場合、上の式はx∈R∖Zx\in \mathbb{R}\setminus \mathbb{Z}に対して成立する。

説明

(c) とは言い換えると、{Nm(x−k)}k\left\{ N_{m}(x-k) \right\}_{k}が単位分割であるという意味だ。

証明

(b)

  • Step 1. m=1m=1の場合成立する。

    N1N_{1}の定義により自明に次の式が成立する。

    ∫−∞∞N1(x)dx=∫01dx=1 \int _{-\infty} ^{\infty} N_{1}(x)dx=\int _{0}^{1} dx = 1

  • Step 2. mmの場合成立すれば、m+1m+1の場合も成立する。

    あるm∈Nm\in \mathbb{N}に対して

    ∫−∞∞Nm(x)dx=1 \int _{-\infty} ^{\infty} N_{m}(x)dx=1

    が成立すると仮定しよう。すると定義により

    ∫Nm+1(x)dx=∫Nm∗N1(x)dx=∫x∫tNm(x−t)N1(t)dtdx=∫tN1(t)∫xNm(x−t)dxdt=∫N1(t)dt=1 \begin{align*} \int N_{m+1}(x)dx &= \int N_{m} \ast N_{1}(x)dx \\ &= \int_{x} \int_{t}N_{m}(x-t)N_{1}(t)dtdx \\ &= \int_{t}N_{1}(t) \int_{x}N_{m}(x-t)dxdt \\ &= \int N_{1}(t)dt \\ &= 1 \end{align*}

  • Step 3.

    数学的帰納法によって、任意のm∈Nm\in Nに対して次が成立する。

    ∫−∞∞Nm(x)dx=1 \int _{-\infty} ^{\infty} N_{m}(x)dx=1

■

(c)

  • Step 1. m=2m=2の場合成立する。

    N2(x)={x0≤x≤1−x+21≤x≤20otherwise N_{2}(x) = \begin{cases} x & 0\le x \le 1 \\ -x+2 & 1 \le x \le 2 \\ 0 & \text{otherwise} \end{cases}

    だから

    N2(x−k)= {x−kk≤x≤k+1−x+k+2k+1≤x≤k+20otherwiseN2(x−(k−1))= {x−k+1k−1≤x≤k−x+k+1k≤x≤k+10otherwise \begin{align*} N_{2}(x-k) =&\ \begin{cases} x-k & k\le x \le k+ 1 \\ -x+k+2 & k+1 \le x \le k+2 \\ 0 & \text{otherwise} \end{cases} \\ N_{2}(x-(k-1)) =&\ \begin{cases} x-k+1 & k-1\le x \le k \\ -x+k+1 & k \le x \le k+1 \\ 0 & \text{otherwise} \end{cases} \end{align*}

    が成立する。あるj∈Zj\in \mathbb{Z}に対してj≤x0≤j+1j \le x_{0} \le j+1としよう。すると

    N2(x0−k)=0,k∈Z∖{j,j−1} N_{2}(x_{0}-k)=0,\quad k\in \mathbb{Z}\setminus \left\{ j,j-1 \right\}

    が成立する。従って

    ∑k∈ZN2(x0−k)= N2(x0−j)+N2(x0−(j−1))= (x0−j)+(−x0+j+1)= 1 \begin{align*} \sum \limits _{k\in \mathbb{Z}} N_{2}(x_{0}-k) =&\ N_{2}(x_{0}-j)+N_{2}(x_{0}-(j-1)) \\ =&\ (x_{0}-j)+(-x_{0}+j+1) \\ =&\ 1 \end{align*}

    が成立する。これは任意のx0x_{0}に対して成立するので

    ∑k∈ZN2(x−k)=1 \sum \limits _{k\in \mathbb{Z}} N_{2}(x-k) =1

  • Step 2. mmの場合成立すれば、m+1m+1の場合も成立する。

    あるm∈Nm\in \mathbb{N}に対して(1)(1)が成立すると仮定しよう。すると定義によって次が成立する。

    ∑k∈ZNm+1(x−k)= ∑k∈Z∫01Nm(x−k−t)dt= ∫01∑k∈ZNm(x−k−t)dt= ∫01dt= 1 \begin{align*} \sum \limits _{k\in \mathbb{Z}} N_{m+1}(x-k) =&\ \sum \limits _{k\in \mathbb{Z}} \int_{0} ^{1} N_{m}(x-k-t)dt \\ =&\ \int_{0} ^{1} \sum \limits _{k\in \mathbb{Z}} N_{m}(x-k-t)dt \\ =&\ \int_{0}^{1} dt \\ =&\ 1 \end{align*}

    二番目の等式は∑k∈Z\sum _{k\in \mathbb{Z}}が**Step 1.**で示されたように有限和であるので成立する。

  • Step 3.

    数学的帰納法によって、任意のm∈Nm\in \mathbb{N}に対して次が成立する。

    ∑k∈ZNm(x−k)=1 \sum \limits _{k\in \mathbb{Z}} N_{m}(x-k) =1

■

(c')

N1N_{1}の定義によって

∑k∈ZN1(x−k)={1x∈R∖Z2x∈Z \sum \limits _{k\in \mathbb{Z}} N_{1}(x-k) =\begin{cases} 1 & x\in \mathbb{R}\setminus\mathbb{Z} \\ 2 & x \in \mathbb{Z} \end{cases}

■

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