コンテンツにスキップ

第8回 フーリエ変換の導入 講義ノート

学習目標

この授業では、以下の概念と技能を習得することを目指します: - フーリエ級数からフーリエ変換への移行を理解する - フーリエ変換の定義と基本的な性質を理解する - 具体的な関数のフーリエ変換を計算できる - フーリエ変換の応用について理解する

1. フーリエ級数からフーリエ変換へ

1.1 フーリエ級数の復習

周期 \(2\pi\) の関数 \(f(x)\) のフーリエ級数展開: $$ f(x) \sim \frac{a_0}{2} + \sum_{n=1}^{\infty} (a_n \cos nx + b_n \sin nx) $$

周期 \(2L\) の関数への拡張: $$ f(x) \sim \frac{a_0}{2} + \sum_{n=1}^{\infty} \left(a_n \cos \frac{n\pi x}{L} + b_n \sin \frac{n\pi x}{L}\right) $$

1.2 フーリエ変換への移行

オイラーの公式からの導出: 1. オイラーの公式を用いて三角関数を複素指数関数で表現:

\[ \cos \theta = \frac{e^{i\theta} + e^{-i\theta}}{2}, \quad \sin \theta = \frac{e^{i\theta} - e^{-i\theta}}{2i} \]
  1. 周期 \(2L\) のフーリエ級数を複素指数関数で表現:
\[ \begin{aligned} f(x) &\sim \frac{a_0}{2} + \sum_{n=1}^{\infty} \left(a_n \cos \frac{n\pi x}{L} + b_n \sin \frac{n\pi x}{L}\right) \\ &= \frac{a_0}{2} + \sum_{n=1}^{\infty} \left(a_n \frac{e^{i\frac{n\pi x}{L}} + e^{-i\frac{n\pi x}{L}}}{2} + b_n \frac{e^{i\frac{n\pi x}{L}} - e^{-i\frac{n\pi x}{L}}}{2i}\right) \\ &= \frac{a_0}{2} + \sum_{n=1}^{\infty} \left(\frac{a_n - ib_n}{2} e^{i\frac{n\pi x}{L}} + \frac{a_n + ib_n}{2} e^{-i\frac{n\pi x}{L}}\right) \end{aligned} \]
  1. 係数を整理して複素フーリエ級数に:
  2. \(c_0 = \frac{a_0}{2}\)
  3. \(c_n = \frac{a_n - ib_n}{2}\)\(n > 0\)
  4. \(c_{-n} = \frac{a_n + ib_n}{2}\)\(n > 0\)

とおくと:

\[ f(x) \sim \sum_{n=-\infty}^{\infty} c_n e^{i\frac{n\pi x}{L}} \]
  1. 係数 \(c_n\) の計算: 複素指数関数の直交性:
\[ \int_{-L}^{L} e^{i\frac{m\pi x}{L}} e^{-i\frac{n\pi x}{L}} dx = \begin{cases} 2L & (m = n) \\ 0 & (m \neq n) \end{cases} \]

を用いて、\(f(x)\) の両辺に \(e^{-i\frac{n\pi x}{L}}\) を掛けて積分すると:

\[ \int_{-L}^{L} f(x) e^{-i\frac{n\pi x}{L}} dx = c_n \cdot 2L \]

したがって:

\[ c_n = \frac{1}{2L} \int_{-L}^{L} f(x) e^{-i\frac{n\pi x}{L}} dx \]
  1. 周期 \(L\) を無限大に近づける:
  2. \(\Delta \xi = \frac{\pi}{L}\) とおくと、\(n\Delta \xi = \frac{n\pi}{L}\)
  3. \(L \to \infty\) のとき \(\Delta \xi \to 0\) となり、離散的な和が積分に
  4. \(n\Delta \xi \to \xi\) となり、離散的な周波数が連続的な周波数に

  5. 極限への移行:

\[ \begin{aligned} f(x) &\sim \sum_{n=-\infty}^{\infty} c_n e^{i\frac{n\pi x}{L}} \\ &= \sum_{n=-\infty}^{\infty} \frac{1}{2L} \left(\int_{-L}^{L} f(x) e^{-i\frac{n\pi x}{L}} dx\right) e^{i\frac{n\pi x}{L}} \\ &= \sum_{n=-\infty}^{\infty} \frac{\Delta \xi}{2\pi} \left(\int_{-L}^{L} f(x) e^{-i\xi x} dx\right) e^{i\xi x} \\ &\to \frac{1}{2\pi} \int_{-\infty}^{\infty} \left(\int_{-\infty}^{\infty} f(x) e^{-i\xi x} dx\right) e^{i\xi x} d\xi \end{aligned} \]
  1. 最終的に:
\[ f(x) = \frac{1}{2\pi} \int_{-\infty}^{\infty} \hat{f}(\xi) e^{i\xi x} d\xi \]

ここで \(\hat{f}(\xi)\) はフーリエ変換:

\[ \hat{f}(\xi) = \int_{-\infty}^{\infty} f(x) e^{-i\xi x} dx \]

2. フーリエ変換の定義と性質

2.1 フーリエ変換の定義

\[ \mathcal{F}[f](\xi)= \hat{f}(\xi) = \int_{-\infty}^{\infty} f(x) e^{-i\xi x} dx \]

2.2 逆フーリエ変換

\[ \mathcal{F}^{-1}[\hat{f}](x) = \frac{1}{2\pi} \int_{-\infty}^{\infty} \hat{f}(\xi) e^{i\xi x} d\xi \]

意味: - 周波数成分 \(\hat{f}(\xi)\) から元の関数 \(f(x)\) を再構築 - フーリエ変換と逆フーリエ変換は互いに逆変換の関係

2.3 基本的な性質

  1. 線形性
\[ \mathcal{F}[af(x) + bg(x)] = a\mathcal{F}[f(x)] + b\mathcal{F}[g(x)] \]
  1. 平行移動
\[ \mathcal{F}[f(x-a)] = e^{-i\xi a} \hat{f}(\xi) \]
  1. スケーリング
\[ \mathcal{F}[f(ax)] = \frac{1}{|a|} \hat{f}\left(\frac{\xi}{a}\right) \]

導出: - \(a > 0\) の場合: $$ \mathcal{F}[f(ax)] = \int_{-\infty}^{\infty} f(ax) e^{-i\xi x} dx $$ 変数変換 \(y = ax\) を行うと: $$ = \frac{1}{a} \int_{-\infty}^{\infty} f(y) e^{-i\frac{\xi}{a}y} dy = \frac{1}{a} \hat{f}\left(\frac{\xi}{a}\right) $$

  • \(a < 0\) の場合: 変数変換 \(y = ax\) を行うと:

    • \(x \to -\infty\) のとき \(y \to +\infty\)\(a < 0\) なので)
    • \(x \to +\infty\) のとき \(y \to -\infty\)
    • \(dx = \frac{1}{a} dy\)\(a < 0\) なので負の値)

    したがって:

$$ \begin{aligned} \mathcal{F}[f(ax)] &= \int_{-\infty}^{\infty} f(ax) e^{-i\xi x} dx \ &= \int_{+\infty}^{-\infty} f(y) e^{-i\frac{\xi}{a}y} \frac{1}{a} dy \ &= -\frac{1}{a} \int_{-\infty}^{\infty} f(y) e^{-i\frac{\xi}{a}y} dy \ &= \frac{1}{|a|} \hat{f}\left(\frac{\xi}{a}\right) \end{aligned} $$ ここで、積分の上下限が逆転し、\(a\) が負の値であるため、最終的に \(\frac{1}{|a|}\) となる。

3. 具体的な関数のフーリエ変換

3.1 指数関数

\[ f(x) = e^{-a|x|} \quad (a > 0) \]

のフーリエ変換:

計算過程

\[ \begin{aligned} \hat{f}(\xi) &= \int_{-\infty}^{\infty} e^{-a|x|} e^{-i\xi x} dx \\ &= \int_{-\infty}^{0} e^{ax} e^{-i\xi x} dx + \int_{0}^{\infty} e^{-ax} e^{-i\xi x} dx \\ &= \int_{-\infty}^{0} e^{(a-i\xi)x} dx + \int_{0}^{\infty} e^{-(a+i\xi)x} dx \\ &= \left[\frac{1}{a-i\xi} e^{(a-i\xi)x}\right]_{-\infty}^{0} + \left[\frac{-1}{a+i\xi} e^{-(a+i\xi)x}\right]_{0}^{\infty} \\ &= \frac{1}{a-i\xi} + \frac{1}{a+i\xi} \\ &= \frac{2a}{a^2 + \xi^2} \end{aligned} \]

別の計算過程(オイラーの公式を使用)

オイラーの公式 \(e^{-i\xi x} = \cos(\xi x) - i\sin(\xi x)\) を用いて計算することもできます:

\[ \begin{aligned} \hat{f}(\xi) &= \int_{-\infty}^{\infty} e^{-a|x|} e^{-i\xi x} dx \\ &= \int_{-\infty}^{\infty} e^{-a|x|} \cos(\xi x) dx - i \int_{-\infty}^{\infty} e^{-a|x|} \sin(\xi x) dx \end{aligned} \]

ここで、\(e^{-a|x|}\) は偶関数、\(\cos(\xi x)\) は偶関数、\(\sin(\xi x)\) は奇関数であることを利用すると:

\[ \begin{aligned} \hat{f}(\xi) &= 2\int_{0}^{\infty} e^{-ax} \cos(\xi x) dx - i \cdot 0 \\ &= 2\int_{0}^{\infty} e^{-ax} \cos(\xi x) dx \end{aligned} \]

この積分は部分積分を2回行うことで計算できます:

\[ \begin{aligned} \int e^{-ax} \cos(\xi x) dx &= \frac{1}{a} e^{-ax} \cos(\xi x) + \frac{\xi}{a} \int e^{-ax} \sin(\xi x) dx \\ &= \frac{1}{a} e^{-ax} \cos(\xi x) + \frac{\xi}{a} \left(-\frac{1}{a} e^{-ax} \sin(\xi x) + \frac{\xi}{a} \int e^{-ax} \cos(\xi x) dx\right) \end{aligned} \]

これを整理すると:

\[ \left(1 + \frac{\xi^2}{a^2}\right) \int e^{-ax} \cos(\xi x) dx = \frac{1}{a} e^{-ax} \cos(\xi x) - \frac{\xi}{a^2} e^{-ax} \sin(\xi x) \]

したがって:

\[ \int_{0}^{\infty} e^{-ax} \cos(\xi x) dx = \frac{a}{a^2 + \xi^2} \]

よって:

\[ \hat{f}(\xi) = 2 \cdot \frac{a}{a^2 + \xi^2} = \frac{2a}{a^2 + \xi^2} \]

このように、オイラーの公式を用いてsinとcosに分けて計算しても、同じ結果が得られます。

3.2 矩形関数

\[ f(x) = \begin{cases} 1 & (|x| \leq a) \\ 0 & (|x| > a) \end{cases} \]

のフーリエ変換:

計算過程

\[ \begin{aligned} \hat{f}(\xi) &= \int_{-\infty}^{\infty} f(x) e^{-i\xi x} dx \\ &= \int_{-a}^{a} e^{-i\xi x} dx \\ &= \left[\frac{-1}{i\xi} e^{-i\xi x}\right]_{-a}^{a} \\ &= \frac{-1}{i\xi} (e^{-i\xi a} - e^{i\xi a}) \\ &= \frac{2\sin(a\xi)}{\xi} \end{aligned} \]

別の計算過程(オイラーの公式を使用)

オイラーの公式 \(e^{-i\xi x} = \cos(\xi x) - i\sin(\xi x)\) を用いて計算することもできます:

\[ \begin{aligned} \hat{f}(\xi) &= \int_{-a}^{a} e^{-i\xi x} dx \\ &= \int_{-a}^{a} \cos(\xi x) dx - i \int_{-a}^{a} \sin(\xi x) dx \end{aligned} \]

ここで、\(\cos(\xi x)\) は偶関数、\(\sin(\xi x)\) は奇関数であることを利用すると:

\[ \begin{aligned} \hat{f}(\xi) &= 2\int_{0}^{a} \cos(\xi x) dx - i \cdot 0 \\ &= 2\left[\frac{1}{\xi} \sin(\xi x)\right]_{0}^{a} \\ &= \frac{2}{\xi} \sin(a\xi) \end{aligned} \]

このように、オイラーの公式を用いてsinとcosに分けて計算しても、同じ結果が得られます。

3.3 ガウス関数

\[ f(x) = e^{-ax^2} \quad (a > 0) \]

のフーリエ変換:

計算過程

\[ \begin{aligned} \hat{f}(\xi) &= \int_{-\infty}^{\infty} e^{-ax^2} e^{-i\xi x} dx \\ &= \int_{-\infty}^{\infty} e^{-a(x^2 + \frac{i\xi}{a}x)} dx \\ &= e^{-\frac{\xi^2}{4a}} \int_{-\infty}^{\infty} e^{-a(x + \frac{i\xi}{2a})^2} dx \\ &= e^{-\frac{\xi^2}{4a}} \sqrt{\frac{\pi}{a}} \end{aligned} \]

補足説明

最後の積分 \(\int_{-\infty}^{\infty} e^{-a(x + \frac{i\xi}{2a})^2} dx\)\(\sqrt{\frac{\pi}{a}}\) となる理由について説明します。

  1. 複素平面での積分経路を考える:
  2. 実軸上の積分 \(\int_{-\infty}^{\infty} e^{-ax^2} dx = \sqrt{\frac{\pi}{a}}\) は既知
  3. 複素平面で \(x\)\(x + \frac{i\xi}{2a}\) に平行移動した積分を考える

  4. コーシーの積分定理の適用:

  5. 被積分関数 \(e^{-az^2}\) は複素平面全体で正則
  6. 実軸から虚軸方向に \(\frac{i\xi}{2a}\) だけ平行移動した経路での積分値は同じ

  7. 積分経路の評価:

  8. 実軸から虚軸方向への平行移動による積分値の変化は0
  9. 無限遠での寄与も0になる
  10. したがって、積分値は \(\sqrt{\frac{\pi}{a}}\) のまま

このように、複素数が絡んでも積分値が変わらないのは、複素解析の重要な性質によるものです。

別の計算過程(オイラーの公式と微分方程式を使用)

オイラーの公式 \(e^{-i\xi x} = \cos(\xi x) - i\sin(\xi x)\) を用いて計算することもできます:

\[ \begin{aligned} \hat{f}(\xi) &= \int_{-\infty}^{\infty} e^{-ax^2} e^{-i\xi x} dx \\ &= \int_{-\infty}^{\infty} e^{-ax^2} \cos(\xi x) dx - i \int_{-\infty}^{\infty} e^{-ax^2} \sin(\xi x) dx \end{aligned} \]

ここで、\(e^{-ax^2}\) は偶関数、\(\cos(\xi x)\) は偶関数、\(\sin(\xi x)\) は奇関数であることを利用すると:

\[ \begin{aligned} \hat{f}(\xi) &= 2\int_{0}^{\infty} e^{-ax^2} \cos(\xi x) dx - i \cdot 0 \\ &= 2\int_{0}^{\infty} e^{-ax^2} \cos(\xi x) dx \end{aligned} \]

この積分を \(I = \int_{0}^{\infty} e^{-ax^2} \cos(\xi x) dx\) とおいて、\(I\) が満たすべき微分方程式を導きます。\(I\)\(\xi\) で微分すると:

\[ \frac{dI}{d\xi} = -\int_{0}^{\infty} x e^{-ax^2} \sin(\xi x) dx \]

部分積分を行うと:

\[ \frac{dI}{d\xi} = -\frac{\xi}{2a} I \]

この微分方程式を解くと:

\[ I = C e^{-\frac{\xi^2}{4a}} \]

ここで、\(C\) は定数で、\(\xi = 0\) のときの積分値から:

\[ C = \int_{0}^{\infty} e^{-ax^2} dx = \frac{1}{2} \sqrt{\frac{\pi}{a}} \]

したがって:

\[ I = \frac{1}{2} e^{-\frac{\xi^2}{4a}} \sqrt{\frac{\pi}{a}} \]

よって:

\[ \hat{f}(\xi) = 2I = e^{-\frac{\xi^2}{4a}} \sqrt{\frac{\pi}{a}} \]

このように、オイラーの公式から始めて微分方程式を解く方法でも、同じ結果が得られます。

4. フーリエ変換の収束性と応用

4.1 収束性の条件

  • 絶対可積分性\(\int_{-\infty}^{\infty} |f(x)| dx < \infty\)
  • \(L^2\) 空間\(\int_{-\infty}^{\infty} |f(x)|^2 dx < \infty\)

4.2 リーマン・ルベーグの補題

  • 絶対可積分関数のフーリエ変換は \(\xi \to \infty\) で 0 に収束
  • 高周波成分は急速に減衰

4.3 応用例

  1. 微分方程式の解法
  2. 微分方程式をフーリエ変換すると代数方程式に
  3. 熱方程式、波動方程式の解法に応用

  4. 畳み込み積分

\[ \mathcal{F}[f * g] = \hat{f}(\xi) \cdot \hat{g}(\xi) \]
  • 信号処理や画像処理で重要

  • パーセバルの等式 $$ \int_{-\infty}^{\infty} |f(x)|^2 dx = \frac{1}{2\pi} \int_{-\infty}^{\infty} |\hat{f}(\xi)|^2 d\xi $$

5. 重要ポイントのまとめ

  • フーリエ変換は、関数を周波数成分に分解する強力なツール
  • 基本的な性質(線形性、平行移動、スケーリング)を理解することが重要
  • 具体的な関数のフーリエ変換を計算できるようになる
  • 収束性の条件と応用について理解を深める