logo

セントラルB-スプライン 📂フーリエ解析

セントラルB-スプライン

定義1

mNm\in \mathbb{N}に関して、中心B-スプライン BmB_{m}を次のように定義する。

Bm(x):=Tm2Nm(x)=Nm(x+12) B_{m}(x):= T_{-\frac{m}{2}}N_{m}(x)=N_{m}(x+{\textstyle \frac{1}{2}})

ここで、TTL2L^{2}空間トランスレーションである。

説明

次のように定義することもできる。

B1:=χ[1/2,1/2],Bm+1:=BmB1, mN B_{1}:= \chi_{[-1/2,1/2]},\quad B_{m+1}:=B_{m}*B_{1},\ m\in\mathbb{N}

二つの定義は実際には同じ関数を意味している。ここでのキーポイントはBmB_{m}を偶関数になるように定義したことである。B-スプラインでは、次の式が成り立つことが簡単にわかる。

Bm+1(x)=Bm(xt)B1(t)dt=1212Bm(xt)dt B_{m+1}(x)=\int _{-\infty} ^{\infty} B_{m}(x-t)B_{1}(t)dt=\int_{-{\textstyle \frac{1}{2}}}^{{\textstyle \frac{1}{2}}}B_{m}(x-t)dt

中心B-スプラインは、B-スプラインを平行移動したものに過ぎないので、B-スプラインの性質をそのまま持っている。

性質

(a) suppBm=[m2,m2]\mathrm{supp} B_{m}=[-\frac{m}{2},\frac{m}{2}]

(b) Bm(x)dx=1\displaystyle \int _{-\infty} ^{\infty} B_{m}(x)dx=1

(c) m2m\ge 2に対して、

kZBm(xk)=1,xR \sum \limits_{k\in \mathbb{Z}}B_{m}(x-k)=1,\quad \forall x\in R

(c’) m=1m=1の時、上記の式はxR{±12,±32,}x\in \mathbb{R}\setminus \left\{ \pm\frac{1}{2},\pm \frac{3}{2},\dots \right\}に対して成立する。

(d) 中心B-スプラインのフーリエ変換は次のようである。

Bm^(γ)=(eπiγeπiγ2πiγ)m=(sin(πγ)πγ)m \widehat{B_{m}}(\gamma)=\left( \frac{e^{\pi i \gamma}-e^{-\pi i \gamma}}{2\pi i \gamma} \right)^{m}=\left( \frac{\sin (\pi\gamma)}{\pi \gamma} \right)^{m}

ここで、ffのフーリエ変換 f^\widehat{f}の定義は次の通りである。

f^(γ):=f(x)e2πixγdx \widehat{f}(\gamma):=\int _{-\infty} ^{\infty} f(x)e^{-2\pi i x\gamma}dx

証明

(a)

suppNm=[m2,m2]\mathrm{supp}N_{m}=\left[ -\frac{m}{2}, \frac{m}{2} \right]であり、Bm=Tm2NmB_{m}=T_{- \frac{m}{2}}N_{m}であるため、

suppBm=[0m2,mm2]=[m2,m2] \mathrm{supp} B_{m} = \left[0-\frac{m}{2},m-\frac{m}{2} \right] = \left[ -\frac{m}{2},\frac{m}{2} \right]

(b)

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

であり、Bm=Tm2NmB_{m}=T_{- \frac{m}{2}}N_{m}であるため、

Bm(x)dx=Tm2Bm(x)dx=1 \int _{-\infty} ^{\infty} B_{m}(x)dx=\int _{-\infty} ^{\infty} T_{-\frac{m}{2}}B_{m}(x)dx=1

(c)

kZNm(xk)=1,xR \sum \limits_{k \in \mathbb{Z}} N_{m}(x-k)=1,\quad \forall x\in \mathbb{R}

であるため、

kZBm(xk)=kZTm2Nm(xk)=1,xR \sum \limits_{k\in \mathbb{Z}}B_{m}(x-k)=\sum \limits_{k \in \mathbb{Z}} T_{-\frac{m}{2}}N_{m}(x-k)=1,\quad \forall x\in \mathbb{R}

(c')

m=1m=1の時、NmN_{m}に対してxRZx\in \mathbb{R}\setminus \mathbb{Z}が成立し、Bm=Tm2NmB_{m}=T_{- \frac{m}{2}}N_{m}であるため、xR{±12,±32,}x\in \mathbb{R}\setminus \left\{ \pm\frac{1}{2},\pm \frac{3}{2},\dots \right\}に対して成立する。

(d)

B-スプラインのフーリエ変換は次の通りである。

Nm^(γ)=(1e2πiγ2πiγ)m \widehat{N_{m}}(\gamma)=\left( \frac{1-e^{-2\pi i\gamma}}{2\pi i \gamma} \right)^{m}

すると、フーリエ変換の性質によって

F[Nm(x+m2)](γ)=e2πim2γNm^(γ) \mathcal{F}\left[ N_{m}(x+\frac{m}{2}) \right] (\gamma)=e^{2\pi i \frac{ m}{2} \gamma}\widehat{N_{m}}(\gamma)

であるため、

Bm^(γ)= F[Bm(x)](γ)=F[Nm(x+m2)](γ)= e2πim2γNm^(γ)= (eπiγ)m(1e2πiγ2πiγ)m= (eπiγeπiγ2πiγ)m= (sin(πγ)πγ)m \begin{align*} \widehat{B_{m}}(\gamma) =&\ \mathcal{F}\left[ B_{m}(x) \right] (\gamma) =\mathcal{F}\left[ N_{m}(x+\textstyle \frac{m}{2}) \right] (\gamma) \\ =&\ e^{2\pi i \frac{m}{2} \gamma}\widehat{N_{m}}(\gamma) \\ =&\ \left(e^{\pi i \gamma}\right)^{m}\left( \frac{1-e^{-2\pi i\gamma}}{2\pi i \gamma} \right)^{m} \\ =&\ \left( \frac{e^{\pi i \gamma}-e^{-\pi i \gamma}}{2\pi i \gamma} \right)^{m} \\ =&\ \left( \frac{\sin (\pi\gamma)}{\pi \gamma} \right)^{m} \end{align*}


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