logo

Fourier Inversion Theorem 📂Fourier Analysis

Fourier Inversion Theorem

Buildup

The process of deriving the Fourier transform also derived the definition of the inverse transform. However, this was simply explained to aid understanding, and the transformation formula was not accurately derived. The Fourier inverse transformation is as follows:

f(x)=12πf^(ξ)eiξxdξ \begin{equation} f(x) =\dfrac{1}{2\pi} \int \hat{f}(\xi) e^{i\xi x}d\xi \end{equation}

This equation implies that from ff, we can obtain f^\hat{f} and from f^\hat{f}, we can retrieve ff again. This might seem obvious or trivial, but it’s not at all. For instance, consider differentiation and integration. When we differentiate a polynomial, we lose the constant term, which cannot be retrieved through integration. However, the Fourier transform preserves information. From ff, we obtain the Fourier transform f^\hat{f}, and by applying the inverse transform to this f^\hat{f}, we can obtain ff exactly as it was.

f^(ξ):=f(x)eiξxdx \begin{equation} \hat{f} (\xi) := \int f(x) e^{-i\xi x}dx \end{equation}

When defining the Fourier transform as above, (1)(1) becomes its inverse as demonstrated below.

Fourier Inverse Transform Theorem

Let’s define the Fourier transform of ff as f^\hat{f} in (2)(2). Assume ff is integrable and piecewise continuous. At points of discontinuity, ff is defined as follows:

f(x)=12[f(x)+f(x+)] f(x)=\dfrac{1}{2} \big[ f(x-)+f(x+)\big]

Then, the following equation holds:

f(x)=limϵ012πf^(ξ)eiξxeϵ2ξ2/2dξ f(x)=\lim \limits_{\epsilon \rightarrow 0}\dfrac{1}{2\pi}\int\hat{f}(\xi)e^{i\xi x}e^{-\epsilon^2\xi^2/2}d\xi

Moreover, if we assume that f^\hat{f} \in is L1L^1, then ff is continuous and

f(x)=12πf^(ξ)eiξxdξ=F1[f^](x) f(x) =\dfrac{1}{2\pi}\int\hat{f}(\xi) e^{i\xi x}d\xi = \mathcal{F}^{-1}[\hat{f}] (x)

Corollary

If we assume f^=g^\hat{f}=\hat{g}, then f=gf=g follows.

Proof

f=F1[f^]=F1[g^] f=\mathcal{F}^{-1}[\hat{f}]=\mathcal{F}^{-1}[\hat{g}]

Therefore,

f=F1[g^]=F1F[g]=g f=\mathcal{F}^{-1}[\hat{g}]=\mathcal{F}^{-1}\mathcal{F}[g]=g

Explanation

Upon observing this corollary, one might think, ‘Isn’t this obvious?’ but it’s not at all. Consider the operation of differentiation. Suppose f(x)=x2+1f(x)=x^2+1, g(x)=x2+4g(x)=x^2+4. In this case, the fact that f(x)=2x=g(x)f^{\prime}(x)=2x=g^{\prime}(x) does not guarantee f(x)=g(x)f(x)=g(x).

Proof

Strategy:

To solve the problem easily, a cutoff function will be used. A cutoff function has the effect of cutting off outside a certain range when multiplied and taken to a limit, hence the name. It is used to maintain the original value near the origin and to converge to 00 as it moves away from the origin. It’s a concept similar to a mollifier. If this explanation is difficult to understand, it’s okay to skip it. Either way, such a cutoff function will be multiplied into the equation to derive the inverse transformation.

η(ξ)=12πex22,ηϵ(ξ)=1ϵη(ξϵ)=1ϵ2πex22ϵ2 \eta (\xi)=\dfrac{1}{\sqrt{2\pi}}e^{-\frac{x^2}{2}}, \quad \eta_\epsilon (\xi)=\dfrac{1}{\epsilon}\eta \left( \frac{\xi}{\epsilon} \right)=\dfrac{1}{\epsilon\sqrt{2\pi}}e^{-\frac{x^2}{2\epsilon^2}}

Multiplying this cutoff function yields

12πf^(ξ)eiξxeϵ2ξ2/2dξ \dfrac{1}{2\pi}\int\hat{f} (\xi) e^{i\xi x}e^{-\epsilon^2\xi^2/2}d\xi

This equation, by the definition of the Fourier transform, is as follows.

12πf(y)eiξydyeiξxeϵ2ξ2/2dξ \dfrac{1}{2\pi}\int \int f(y) e^{-i\xi y}dy e^{i\xi x}e^{-\epsilon^2\xi^2/2}d\xi

Using the definition of the Fourier transform and convolution, among others, the equation is transformed as follows:

12πf(y)eiξydyeiξxeϵ2ξ2/2dξ=12πf(y)eiξ(yx)eϵ2ξ2/2dydξ=12πF[eϵ2ξ2/2](yx)f(y)dy=12π2πϵ2e(yx)22ϵ2f(y)dy=1ϵ2πe(xy)22ϵ2f(y)dy=ηϵ(xy)f(y)dy=fηϵ(x) \begin{align*} \dfrac{1}{2\pi}\int \int f(y) e^{-i\xi y}dy e^{i\xi x}e^{-\epsilon^2\xi^2/2}d\xi &= \dfrac{1}{2\pi}\int {\color{blue}\int} f(y) {\color{blue}e^{-i\xi (y-x)}e^{-\epsilon^2\xi^2/2}}dy {\color{blue}d\xi} \\ &= \dfrac{1}{2\pi}\int {\color{blue}\mathcal{F} \left[ e^{-\epsilon^2 \xi^2 /2}\right] (y-x)} f(y) dy \\ &= \dfrac{1}{2\pi}\int {\color{blue}\sqrt{\dfrac{2\pi}{\epsilon^2}}e^{-\frac{(y-x)^2}{2\epsilon^2}}} f(y) dy \\ &= \int \dfrac{1}{\epsilon \sqrt{2\pi}} e^{-\frac{(x-y)^2}{2\epsilon^2}} f(y) dy \\ &= \int \eta_\epsilon (x-y)f(y) dy \\ &= f \ast \eta_\epsilon (x) \end{align*}

The third equality used the formula below.

Fourier Transform of a Gaussian Function

The Fourier transform of a Gaussian function f(x)=eAx2f(x)=e^{-Ax^2} is as follows:

F[f](ξ)=F[eAx2](ξ)=πAeξ24A \mathcal{F}[f] (\xi) = \mathcal{F} \left[ e^{-Ax^2} \right] (\xi)=\sqrt{\frac{\pi}{A}}e^{-\frac{\xi ^2}{4A}}

Then, assuming ff is piecewise continuous, the mollification of ff converges as follows.

limϵ012πf^(ξ)eiξxeϵ2ξ2/2dξ=limϵ0fηϵ(x)=12[f(x)+f(x+)]=f(x) \begin{align*} \lim \limits_{\epsilon \rightarrow 0} \dfrac{1}{2\pi}\int\hat{f} (\xi) e^{i\xi x}e^{-\epsilon^2\xi^2/2}d\xi &= \lim \limits_{\epsilon \rightarrow 0} f \ast \eta_\epsilon (x) \\ &= \dfrac{1}{2} \big[ f(x-)+f(x+)\big] \\ &= f(x) \end{align*}

Furthermore, assuming f^L1\hat{f} \in L^1,

eiξxeϵ2ξ2/2f^(ξ)f^(ξ) \left| e^{i\xi x}e^{-\epsilon^2|\xi|^2 /2} \hat{f}(\xi) \right| \le \left| \hat{f}(\xi) \right|

thus, by the Dominated Convergence Theorem,

f(x)=limϵ012πf^(ξ)eiξxeϵ2ξ2/2dξ=12πlimϵ0f^(ξ)eiξxeϵ2ξ2/2dξ=12πf^(ξ)eiξxdξ \begin{align*} f(x) &= \lim \limits_{\epsilon \rightarrow 0} \dfrac{1}{2\pi}\int\hat{f} (\xi) e^{i\xi x}e^{-\epsilon^2\xi^2/2}d\xi \\ &= \dfrac{1}{2\pi}\int \lim \limits_{\epsilon \rightarrow 0} \hat{f} (\xi) e^{i\xi x}e^{-\epsilon^2\xi^2/2}d\xi \\ &= \dfrac{1}{2\pi}\int \hat{f} (\xi) e^{i\xi x}d\xi \end{align*}

Therefore,

f(x)=12πf^(ξ)eiξxdξ f(x)=\dfrac{1}{2\pi}\int \hat{f} (\xi) e^{i\xi x}d\xi