logo

B-スプラインのフーリエ変換 📂フーリエ解析

B-スプラインのフーリエ変換

1

順序がmNm \in \mathbb{N}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}

この時、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

説明

B-スプライン、フーリエ変換、畳み込みの性質を利用して、難しくなく計算できる。

証明

まず、N1N_{1}のフーリエ変換を計算すると、次のようになる。

FN1(γ)= N1(x)e2πixγdx= 01e2πixγdx= [e2πixγ2πiγ]x=01= 1e2πiγ2πiγ \begin{align*} \mathcal{F}N_{1}(\gamma) =&\ \int _{-\infty} ^{\infty} N_{1}(x)e^{-2\pi i x \gamma }dx \\ =&\ \int_{0}^{1} e^{-2\pi i x \gamma}dx \\ =&\ \left[\frac{e^{-2\pi i x\gamma} }{-2\pi i \gamma} \right]_{x=0}^{1} \\ =&\ \frac{1-e^{2\pi i \gamma}}{2\pi i \gamma} \end{align*}

B-スプラインの定義によってNm=N1N1N1mN_{m}=\overbrace{N_{1} * N_{1} * \cdots * N_{1}}^{m}であり、フーリエ変換の性質によると、

F[f1f2fn]=f1^f2^fn^ \mathcal{F}\left[ f_{1} * f_{2}*\cdots * f_{n} \right]=\hat{f_{1}} \hat{f_{2}} \cdots \hat{f_{n}}

したがって、

FNm(γ)= F[N1N1N1m](γ)= N1^(γ)N1^(γ)N1^(γ)m= (N1^(γ))m= (1e2πiγ2πiγ)m \begin{align*} \mathcal{F} N_{m}(\gamma) =&\ \mathcal{F} \left[ \overbrace{ N_{1} * N_{1} * \cdots * N_{1} }^{m} \right] (\gamma) \\ =&\ \overbrace{ \widehat{N_{1}}(\gamma)\widehat{N_{1}}(\gamma)\cdots\widehat{N_{1}}(\gamma) }^{m} \\ =&\ \left( \widehat{N_{1}}(\gamma) \right)^{m} \\ =&\ \left( \frac{1-e^{-2\pi i\gamma}}{2\pi i \gamma} \right)^{m} \end{align*}


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