フーリエ変換の公式と導出

POINT

  • フーリエ変換の関係式とその導出を一覧にしました.

定義

関数$f$のフーリエ変換$\hat{f}=\mathcal{F}[f]$を
\begin{aligned}
\hat{f}(k)
&=\mathcal{F}[f](k) \\
&=\frac{1}{\sqrt{2\pi}}\int_{-\infty}^\infty f(x)e^{-ikx} \,\mathrm{d}x
\end{aligned}
で定義する.

定義の仕方は他にもあります.詳しくは後述(2πの付け方の違い).



$n$次元の場合はベクトル$k=^{t}\!\!(k_1,k_2,...,k_n)$と$x=^{t}\!\!(x_1,x_2,...,x_n)$を使って
\begin{aligned}
\hat{f}(k)
&=\mathcal{F}[f](k) \\
&=\frac{1}{(2\pi)^{n/2}}\int_{\mathbb{R}^n} f(x)e^{-ik\cdot x} \,\mathrm{d}x
\end{aligned}
となります($k\cdot x=\sum_i k_i x_i$はベクトルの内積).

逆変換

フーリエ逆変換(フーリエ反転公式)
\begin{aligned}
f(x)
&=\frac{1}{(2\pi)^{n/2}}\int_{\mathbb{R}^n} \hat{f}(k) e^{ik\cdot x} \,\mathrm{d}k
\end{aligned}
が成立します.

性質

フーリエ変換の基本的な性質として紹介される大体のものは,必要になったときに都度計算すればよいです.以下では$a$を実数とします.

$f(x+a)$

\begin{aligned}
&\frac{1}{\sqrt{2\pi}}\int_{-\infty}^\infty f(x+a)e^{-ikx} \,\mathrm{d}x \\
&=\frac{1}{\sqrt{2\pi}}\int_{-\infty}^\infty f(x)e^{-ik(x-a)} \,\mathrm{d}x \\
&=e^{ika} \hat{f}(k)
\end{aligned}

$e^{iax}f(x)$

\begin{aligned}
&\frac{1}{\sqrt{2\pi}}\int_{-\infty}^\infty [e^{iax}f(x)]e^{-ikx} \,\mathrm{d}x \\
&=\frac{1}{\sqrt{2\pi}}\int_{-\infty}^\infty f(x)e^{-i(k-a)x} \,\mathrm{d}x \\
&=\hat{f}(k-a)
\end{aligned}

$f(ax)$

$\mathrm{sgn\,}(a)=a/|a|(=|a|/a)$とするとき,
\begin{aligned}
&\frac{1}{\sqrt{2\pi}}\int_{-\infty}^\infty f(ax)e^{-ikx} \,\mathrm{d}x \\
&=\frac{1}{\sqrt{2\pi}}\int_{-\infty}^\infty f(ax)e^{-i(k/a)(ax)} \,\frac{\mathrm{d}(ax)}{a} \\
&=\frac{1}{a}\cdot\mathrm{sgn\,}(a)\hat{f}(k/a) \\
&=\frac{1}{|a|}\hat{f}(k/a).
\end{aligned}
ここで$\mathrm{sgn\,}(a)$は$a$の符号で変数変換後の積分範囲が$\displaystyle \int_{-\infty}^{\infty}$となるか$\displaystyle \int_{\infty}^{-\infty}$となるかが変わるために出てきます.わかりにくい場合は$a>0$と$a<0$に場合分けして計算すると良いでしょう.

$\bar{f}(x)$

$f$の複素共役を$\bar{f}$で表すとき,
\begin{aligned}
&\frac{1}{\sqrt{2\pi}}\int_{-\infty}^\infty \bar{f}(x)e^{-ikx} \,\mathrm{d}x \\
&=\overline{\frac{1}{\sqrt{2\pi}}\int_{-\infty}^\infty f(x)e^{ikx} \,\mathrm{d}x} \\
&=\overline{\hat{f}}(-k)
\end{aligned}

フーリエ逆変換

フーリエ逆変換(フーリエ反転公式)
\begin{aligned}
f(x)
&=\frac{1}{(2\pi)^{n/2}}\int_{\mathbb{R}^n} \hat{f}(k) e^{ik\cdot x} \,\mathrm{d}k
\end{aligned}
において$x\rightarrow -k$, $k\rightarrow x$とすれば
\begin{aligned}
f(-k)
&=\frac{1}{(2\pi)^{n/2}}\int_{\mathbb{R}^n} \hat{f}(x) e^{-ik\cdot x} \,\mathrm{d}x
\end{aligned}
となる.右辺は$ \hat{f}$のフーリエ変換であるから,
\begin{aligned}
\mathcal{F}\bigl[\mathcal{F}[f]\bigr](k)=f(-k).
\end{aligned}

偶関数

偶関数$f(-x)=f(x)$のフーリエ変換は
\begin{aligned}
&\frac{1}{\sqrt{2\pi}}\int_{-\infty}^{\infty} f(x) e^{-ikx} \,\mathrm{d}x \\
&=-\frac{1}{\sqrt{2\pi}}\int_{\infty}^{0} f(-x) e^{ikx} \,\mathrm{d}x \\
&\quad+\frac{1}{\sqrt{2\pi}}\int_{0}^{\infty} f(x) e^{-ikx} \,\mathrm{d}x \\
&=\frac{1}{\sqrt{2\pi}}\int_{0}^{\infty} f(x) (e^{ikx}+e^{-ikx}) \,\mathrm{d}x \\
&=\sqrt{\frac{2}{\pi}}\int_{0}^{\infty} f(x) \cos(kx) \,\mathrm{d}x
\end{aligned}
となる(Fourier余弦変換).

奇関数

奇関数$f(-x)=-f(x)$のフーリエ変換は
\begin{aligned}
&\frac{1}{\sqrt{2\pi}}\int_{-\infty}^{\infty} f(x) e^{-ikx} \,\mathrm{d}x \\
&=-\frac{1}{\sqrt{2\pi}}\int_{\infty}^{0} f(-x) e^{ikx} \,\mathrm{d}x \\
&\quad +\frac{1}{\sqrt{2\pi}}\int_{0}^{\infty} f(x) e^{-ikx} \,\mathrm{d}x \\
&=-\frac{1}{\sqrt{2\pi}}\int_{0}^{\infty} f(x) (e^{ikx}-e^{-ikx}) \,\mathrm{d}x \\
&=-i\sqrt{\frac{2}{\pi}} \int_{0}^{\infty} f(x) \sin(kx) \,\mathrm{d}x
\end{aligned}
となる(Fourier正弦変換).

実関数

現実のデータをフーリエ変換で解析する場合,実数値関数のフーリエ変換となる場合が多い.実関数のフーリエ変換には,下の記事のような性質がある.
実数値関数のフーリエ変換の性質 - Notes_JP

相関関数・畳み込み

相関関数と畳み込みは,フーリエ変換で計算できる.
フーリエ変換を使った相関関数と畳み込みの計算方法 - Notes_JP

フーリエ変換を使った計算方法1
\begin{aligned}
f \ast g
& = \mathcal{F}^{-1} \Bigl[\mathcal{F} [f] \cdot \mathcal{F} [g] \Bigr] \\
f \star g
& = \mathcal{F}^{-1} \Bigl[\overline{\mathcal{F} [f]} \cdot \mathcal{F} [g] \Bigr]
\end{aligned}

フーリエ変換を使った計算方法2
\begin{aligned}
(f \ast g) (\tau)
& = \frac{1}{2\pi} \mathcal{F} \Bigl[\mathcal{F} [f] \cdot \mathcal{F} [g] \Bigr] (-\tau) \\
(f \star g) (\tau)
& = \frac{1}{2\pi} \mathcal{F} \Bigl[\overline{\mathcal{F} [f]} \cdot \mathcal{F} [g] \Bigr] (-\tau)
\end{aligned}

ここで,
相互相関関数 (Cross-correlation):

\begin{aligned}
(f \star g) (\tau)
& = \int_{-\infty}^{\infty}
\bar{f} (\textcolor{red}{t}) g(\textcolor{red}{t} + \tau)
\,\mathrm{d}\textcolor{red}{t}
\end{aligned}

畳み込み積分 (Convolution):

\begin{aligned}
(f \ast g) (\tau)
& = \int_{-\infty}^{\infty}
f(\textcolor{red}{t}) g(\tau - \textcolor{red}{t})
\,\mathrm{d}\textcolor{red}{t}
\end{aligned}

フーリエ変換の具体例

基本的な関数のフーリエ変換の計算方法を整理します.より複雑な関数のフーリエ変換は,「基本的な関数のフーリエ変換」と「フーリエ変換の性質」を駆使して計算を進めることになります.

指数関数

$a>0$とする.
\begin{aligned}
\mathcal{F}[e^{-a|x|}](k)
&=\sqrt{\frac{2}{\pi}}\frac{a}{a^2+k^2}
\end{aligned}
【証明】
$f(x)=e^{-a|x|}$は偶関数なので
\begin{aligned}
&\mathcal{F}[e^{-a|x|}](k) \\
&=\sqrt{\frac{2}{\pi}}\int_{0}^\infty e^{-ax}\cos(kx)\,\mathrm{d}x
\end{aligned}
となる.ここで,$I=\displaystyle\int_{0}^\infty e^{-ax}\cos(kx)\,\mathrm{d}x$, $J=\displaystyle\int_{0}^\infty e^{-ax}\sin(kx)\,\mathrm{d}x$とすると
\begin{aligned}
I+iJ
&=\int_{0}^\infty e^{-ax}e^{ikx} \,\mathrm{d}x \\
&=-\frac{1}{a-ik}\bigl[e^{-(a-ik)x}\bigr]_0^\infty \\
&=\frac{a+ik}{a^2+k^2}.
\end{aligned}
実部を取ると$I=\displaystyle \frac{a}{a^2+k^2}$だから
\begin{aligned}
&\mathcal{F}[e^{-a|x|}](k) \\
&=\sqrt{\frac{2}{\pi}}\frac{a}{a^2+k^2}
\end{aligned}
となる.//

ガウス関数

ガウス関数(正規分布)のフーリエ変換は次の記事で紹介しています.
ガウス積分と派生公式 - Notes_JP

\begin{aligned}
\frac{1}{\sqrt{2\pi}}\int_{-\infty}^\infty e^{-x^2} e^{-ikx} \,\mathrm{d}x
&=\frac{1}{\sqrt{2}} e^{-k^2/4}
\end{aligned}

デルタ関数

デルタ関数についてはデルタ関数と公式 - Notes_JPを参照.
\begin{aligned}
\mathcal{F}[\delta](k)
&=\frac{1}{(2\pi)^{n/2}} \\
\delta(x)
&=\frac{1}{(2\pi)^{\textcolor{red}{n}}} \int_{\mathbb{R}^n} e^{ik\cdot x} \,\mathrm{d}k
\end{aligned}

【証明】
デルタ関数は$\displaystyle \int_{\mathbb{R}^n} \delta(x)f(x) \,\mathrm{d}x =f(0)$を満たすので,

\begin{aligned}
\mathcal{F}[\delta](k)
&=\frac{1}{(2\pi)^{n/2}}\int_{\mathbb{R}^n} \delta(x)e^{-ik\cdot x} \,\mathrm{d}x \\
&=\frac{1}{(2\pi)^{n/2}}
\end{aligned}
となる.したがって,逆変換により
\begin{aligned}
\delta(x)
&=\frac{1}{(2\pi)^{n/2}} \int_{\mathbb{R}^n} \mathcal{F}[\delta](k) e^{ik\cdot x} \,\mathrm{d}k \\
&=\frac{1}{(2\pi)^n} \int_{\mathbb{R}^n} e^{ik\cdot x} \,\mathrm{d}k
\end{aligned}
がわかる.//

くし型関数

くし型関数を
\begin{aligned}
\delta_{T}(t) = \sum_{n=-\infty}^{\infty} \delta(t - nT)
\end{aligned}
と表します.

くし型関数の「フーリエ変換」は

\begin{aligned}
\mathcal{F}[\delta_{T}(t)] (\omega)
&= \frac{2\pi}{T} \delta_{2\pi/T} (\omega)
\end{aligned}
と「くし型関数」となります.

導出は以下の記事を参照してください.
周期関数の「フーリエ変換」 - Notes_JP

周期関数

周期$T$の周期関数$f(t)$は
\begin{aligned}
\mathcal{F}[f](\omega)
&= \mathcal{F}[f_{0}](\omega) \cdot \mathcal{F}[\delta_{T}](\omega) \\
&= \frac{2\pi}{T} \mathcal{F}[f_{0}](\omega) \cdot \delta_{2\pi/T} (\omega) \\
f_{0}(t) &=
\begin{cases}
\, f(t) & (-T/2 \leq t < T/2) \\
\, 0 & (\text{otherwise})
\end{cases}
\end{aligned}
となり,$\mathcal{F}[f_{0}]$を周期$2\pi/T$でサンプリングした式になっている.

導出は以下の記事を参照してください.
周期関数の「フーリエ変換」 - Notes_JP

付録:2πの付け方の違い

係数の付け方には他にも流派がありますが,記法を1つ覚えておけば都度変換できます.指数部を$2\pi f t$と表したときに,$f$と$t$をフーリエ変換・逆変換で結びつける関係は$2\pi$の係数がつかないことを覚えておくのがおすすめです:
\begin{aligned}
\hat{g}(f) &= \int_{-\infty}^{\infty} g(t)e^{-i 2\pi ft} \,\mathrm{d}t \\
g(t) & =\int_{-\infty}^{\infty} \hat{g}(f) e^{i 2\pi ft} \,\mathrm{d}f
\end{aligned}


以下,$\alpha,\beta > 0$とします.

全体の定数倍

フーリエ変換を
\begin{aligned}
\hat{f}(k)
&=\mathcal{F}[f](k) \\
&=\frac{ \textcolor{red}{\alpha}^{n/2} }{(2\pi)^{n/2}}\int_{\mathbb{R}^n} f(x)e^{-ik\cdot x} \,\mathrm{d}x
\end{aligned}
で定義した場合,上の逆変換により
\begin{aligned}
\textcolor{red}{\alpha}^{n/2} f(x)
&=\frac{1}{(2\pi)^{n/2}}\int_{\mathbb{R}^n}\hat{f}(k) e^{ik\cdot x} \,\mathrm{d}k
\end{aligned}
がわかるので,
\begin{aligned}
f(x)
&=\frac{1}{(2\pi \textcolor{red}{\alpha})^{n/2}}\int_{\mathbb{R}^n}\hat{f}(k) e^{ik\cdot x} \,\mathrm{d}k
\end{aligned}
が新しい逆変換となります.

  • $\alpha=2\pi$の場合:
    • フーリエ変換:$\displaystyle\hat{f}(k) =\int_{\mathbb{R}^n} f(x)e^{-ik\cdot x} \,\mathrm{d}x$
    • フーリエ逆変換:$\displaystyle f(x)=\frac{1}{(2\pi )^{n}}\int_{\mathbb{R}^n}\hat{f}(k) e^{ik\cdot x} \,\mathrm{d}k$
  • $\alpha=1/2\pi$の場合:
    • フーリエ変換:$\displaystyle\hat{f}(k) =\frac{1}{(2\pi )^{n}}\int_{\mathbb{R}^n} f(x)e^{-ik\cdot x} \,\mathrm{d}x$
    • フーリエ逆変換:$\displaystyle f(x)=\int_{\mathbb{R}^n}\hat{f}(k) e^{ik\cdot x} \,\mathrm{d}k$

係数の定数倍

もう少し一般化しましょう.フーリエ変換を
\begin{aligned}
\hat{f}(k)
&=\mathcal{F}[f](k) \\
&=\frac{ \textcolor{red}{\alpha}^{n/2} }{(2\pi)^{n/2}}
\int_{\mathbb{R}^n} f(x)e^{-i\textcolor{green}{\beta}k\cdot x} \,\mathrm{d}x
\end{aligned}
で定義した場合,
\begin{aligned}
&=\frac{ \textcolor{red}{\alpha}^{n/2} }{(2\pi)^{n/2}}
\int_{\mathbb{R}^n} f(\textcolor{green}{\beta}x/\textcolor{green}{\beta})e^{-ik\cdot (\textcolor{green}{\beta}x)}
\frac{\,\mathrm{d}(\textcolor{green}{\beta}x)}{\textcolor{green}{\beta}^n} \\
&=\frac{\textcolor{red}{\alpha}^{n/2} }{\textcolor{green}{\beta}^n}
\frac{1}{(2\pi)^{n/2}}\int_{\mathbb{R}^n} f(x/\textcolor{green}{\beta})e^{-ik\cdot x} \,\mathrm{d}x
\end{aligned}
なので,上の逆変換により
\begin{aligned}
\frac{\textcolor{red}{\alpha}^{n/2} }{\textcolor{green}{\beta}^n}f(x/\textcolor{green}{\beta})
&=\frac{1}{(2\pi)^{n/2}}\int_{\mathbb{R}^n}\hat{f}(k) e^{ik\cdot x} \,\mathrm{d}k
\end{aligned}
となります.したがって,
\begin{aligned}
f(x)
&= \frac{\textcolor{green}{\beta}^n}{(2\pi \textcolor{red}{\alpha})^{n/2}}
\int_{\mathbb{R}^n}\hat{f}(k) e^{i\textcolor{green}{\beta}k\cdot x} \,\mathrm{d}k
\end{aligned}
が新しい逆変換となります.



$\alpha=\beta=2\pi$がよく使われます.

  • フーリエ変換:$\displaystyle\hat{f}(k) =\int_{\mathbb{R}^n} f(x)e^{-2\pi ik\cdot x} \,\mathrm{d}x$
  • フーリエ逆変換:$\displaystyle f(x)=\int_{\mathbb{R}^n}\hat{f}(k) e^{2\pi ik\cdot x} \,\mathrm{d}k$

参考文献

[1]解析入門 Ⅱ(基礎数学3)
[2]詳解物理応用数学演習

以下の書籍が詳しくておすすめです.