logo

Fourier Transform of B-Splines 📂Fourier Analysis

Fourier Transform of B-Splines

The Equation1

The Fourier transform of a B-spline of order mNm \in \mathbb{N} is given as follows.

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

Here, the definition of the Fourier transform of ff is as follows.

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

Explanation

Using the properties of B-splines, Fourier transforms, and convolutions, the calculation can be done without much difficulty.

Proof

First, computing the Fourier transform of N1N_{1} results in:

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*}

Since a B-spline is defined by Nm=N1N1N1mN_{m}=\overbrace{N_{1} * N_{1} * \cdots * N_{1}}^{m} and due to the properties of Fourier transforms,

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}}

It follows that,

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 ↩︎