メインコンテンツまでスキップ

フーリエ級数

フーリエ級数は、周期関数を正弦波と余弦波の線形結合で表現する数学的手法です。この概念は信号処理において極めて重要で、任意の周期信号を基本的な三角関数の組み合わせとして理解することができます。

基本概念​

数学的背景​

ベクトルの内積の視点から考えると、ベクトルを二つの基底ベクトルの線形結合で表現できるように、ある関数を正弦波と余弦波の線形結合で表現することができます。

周期関数 x(t)x(t) に対して： x(t)=x(t+T),f0=1T,ω0=2πTx(t) = x(t+T), \quad f_0 = \frac{1}{T}, \quad \omega_0 = \frac{2\pi}{T}

ここで：

  • TT：周期
  • f0f_0：基本周波数
  • ω0\omega_0：基本角周波数

フーリエ級数展開​

実関数表現​

周期関数 x(t)x(t) は以下のフーリエ級数で表現されます：

x(t)=a02+∑n=1∞(ancos⁡nω0t+bnsin⁡nω0t)x(t) = \frac{a_0}{2} + \sum_{n=1}^{\infty}\left( a_n\cos{n\omega_0t} + b_n\sin{n\omega_0t}\right)

フーリエ係数​

フーリエ係数は関数の内積として求められます：

an=2T∫0Tx(t)cos⁡nω0t dta_n = \frac{2}{T}\int_{0}^{T}x(t)\cos{n\omega_0t}\,dt

bn=2T∫0Tx(t)sin⁡nω0t dtb_n = \frac{2}{T}\int_{0}^{T}x(t)\sin{n\omega_0t}\,dt

係数の物理的意味
  • ana_n：余弦波成分の係数（関数と余弦関数の内積）
  • bnb_n：正弦波成分の係数（関数と正弦関数の内積）
  • a02\frac{a_0}{2}：直流成分（信号の平均値）

直交性​

三角関数の直交性は、フーリエ級数の数学的基礎となる重要な性質です：

2T∫0Tcos⁡mω0tcos⁡nω0t dt=δmn\frac{2}{T}\int_{0}^{T}\cos{m\omega_0t}\cos{n\omega_0t}\,dt = \delta_{mn}

2T∫0Tsin⁡mω0tsin⁡nω0t dt=δmn\frac{2}{T}\int_{0}^{T}\sin{m\omega_0t}\sin{n\omega_0t}\,dt = \delta_{mn}

2T∫0Tsin⁡mω0tcos⁡nω0t dt=0\frac{2}{T}\int_{0}^{T}\sin{m\omega_0t}\cos{n\omega_0t}\,dt = 0

クロネッカのデルタ​

δmn={1if m=n0if m≠n\delta_{mn} = \begin{cases} 1 & \text{if } m = n \\ 0 & \text{if } m \neq n \end{cases}
直交性の意味

関数の内積の特性として、自分自身との内積が1、異なる関数との内積が0となることを表しています。これにより、各周波数成分を独立して計算できます。

パーセバルの等式​

パーセバルの等式は、時間領域のエネルギーと周波数領域のエネルギーが等しいことを示します：

1T∫0Tx2(t)dt=a024+12∑n=1∞(an2+bn2)\frac{1}{T}\int_{0}^{T}x^2(t)dt = \frac{a_0^2}{4} + \frac{1}{2}\sum_{n=1}^{\infty}(a_n^2 + b_n^2)

この式は信号の平均パワーを表し、左辺は時間領域での平均パワー、右辺は各周波数成分のパワーの総和を表します。

複素数表現​

オイラーの公式​

複素数表現を用いることで、フーリエ級数をより簡潔に表現できます：

ejω0t=cos⁡ω0t+jsin⁡ω0te^{j\omega_0t} = \cos{\omega_0t} + j\sin{\omega_0t}

オイラーの公式により：

cos⁡ω0t=ejω0t+e−jω0t2\cos{\omega_0t} = \frac{e^{j\omega_0t} + e^{-j\omega_0t}}{2}

sin⁡ω0t=ejω0t−e−jω0t2j\sin{\omega_0t} = \frac{e^{j\omega_0t} - e^{-j\omega_0t}}{2j}

複素フーリエ係数​

実フーリエ係数から複素フーリエ係数への変換：

ancos⁡nω0t+bnsin⁡nω0t=an−jbn2ejnω0t+an+jbn2e−jnω0ta_n\cos{n\omega_0t} + b_n\sin{n\omega_0t} = \frac{a_n - jb_n}{2}e^{jn\omega_0t} + \frac{a_n + jb_n}{2}e^{-jn\omega_0t}

複素フーリエ係数を定義すると：

cn=an−jbn2c_n = \frac{a_n - jb_n}{2}

複素フーリエ級数​

x(t)=∑n=−∞∞cnejnω0tx(t) = \sum_{n=-\infty}^{\infty}c_n e^{jn\omega_0t}

cn=1T∫0Tx(t)e−jnω0tdtc_n = \frac{1}{T}\int_{0}^{T}x(t)e^{-jn\omega_0t}dt

複素数の直交性​

複素数版の直交性：

1T∫0Tejmω0te−jnω0tdt=δmn\frac{1}{T}\int_{0}^{T}e^{jm\omega_0t}e^{-jn\omega_0t}dt = \delta_{mn}

複素数版パーセバルの等式​

1T∫0T∣x(t)∣2dt=∑n=−∞∞∣cn∣2\frac{1}{T}\int_{0}^{T}|x(t)|^2dt = \sum_{n=-\infty}^{\infty}|c_n|^2

複素表現の利点

複素数表現を用いることで：

  • 数式が簡潔になる
  • 負の周波数成分も含めた統一的な表現
  • 数学的操作が容易になる
  • フーリエ変換への自然な拡張が可能

実用例と応用​

方形波のフーリエ級数展開​

周期的な方形波は以下のフーリエ級数で表現されます：

x(t)=4Aπ∑n=1,3,5,...∞1nsin⁡nω0tx(t) = \frac{4A}{\pi}\sum_{n=1,3,5,...}^{\infty}\frac{1}{n}\sin{n\omega_0t}

この例から分かるように：

  • 方形波は奇数次の正弦波成分のみで構成される
  • 高次成分の振幅は 1/n1/n で減少する

スペクトラム解析​

フーリエ級数係数 ∣cn∣|c_n| は、信号の振幅スペクトラムを表します：

  • 各周波数成分の振幅を示す
  • 信号の周波数特性を可視化
  • フィルタ設計の基礎となる

まとめ​

フーリエ級数は信号処理の基礎概念として重要です：

  1. 周期信号の分解：任意の周期信号を基本的な三角関数の組み合わせとして表現
  2. 直交性：各周波数成分を独立して扱える数学的基礎
  3. エネルギー保存：パーセバルの等式によるエネルギーの保存則
  4. 複素表現：より簡潔で扱いやすい数学的形式

次章では、フーリエ級数を非周期信号に拡張したフーリエ変換について学習します。

次のステップ
  • フーリエ変換：非周期信号への拡張
  • 離散フーリエ変換（DFT）：ディジタル信号処理での実装
  • 高速フーリエ変換（FFT）：効率的な計算アルゴリズム
  • 窓関数：有限長信号の処理における重要な概念