フーリエ級数

フーリエ級数（フーリエきゅうすう、テンプレート:Lang-en）とは、複雑な周期関数や周期信号を単純な形の周期性をもつ関数の無限和（級数）によって表したものである。フーリエ級数は、フランスの数学者ジョゼフ・フーリエによって金属板の中での熱伝導に関する研究の中で導入された。

熱伝導方程式は、偏微分方程式として表される。フーリエの研究の前までには、一般的な形での熱伝導方程式の解法は知られておらず、熱源が単純な形である場合、例えば正弦波などの場合の特別な解しかえられていなかった。この特別な解は現在では固有解と呼ばれる。フーリエの発想は、複雑な形をした熱源をサイン波、コサイン波の線型結合として考え、解を固有解の和として表すものであった。この重ね合わせがフーリエ級数と呼ばれる。

最初の動機は熱伝導方程式を解くことであったが、数学や物理の他の問題にも同様のテクニックが使えることが分かり様々な分野に応用されている。フーリエ級数は、電気工学、振動の解析、音響学、光学、信号処理、量子力学および経済学^[1]などの分野で用いられている。

概要

フーリエ級数は、関数に対して定義されるフーリエ係数を用いてテンプレート:Indent の形に表される三角級数のことである。熱方程式を発見したフーリエは、平衡状態における熱方程式に注目し、適当な境界条件の下で二変数のラプラス方程式

(\frac{\partial^{2}}{\partial x^{2}} + \frac{\partial^{2}}{\partial y^{2}}) ϕ (x, y) = 0

に帰着させて解を求めようとした。この時、フーリエは、

\sum_{k = 0}^{\infty} \frac{(- 1)^{k}}{2 k + 1} \cos ((2 k + 1) x) = \frac{π}{4}, (- \frac{π}{2} < x < \frac{π}{2})

という三角級数を見つけている。左辺の三角関数の一つ一つは波打っているにもかかわらず、x に依らない定数に収束しているのである。

x = 0 としたときの級数は円周率を求めるグレゴリー級数と同じである。

x の定義域を広げるとこの三角級数は n を整数として

\sum_{k = 0}^{\infty} \frac{(- 1)^{k}}{2 k + 1} \cos ((2 k + 1) x) = (- 1)^{n} \frac{π}{4}, (- \frac{π}{2} + n π < x < \frac{π}{2} + n π)

という矩形波になる。このような不連続な関数まで表せることに興味を抱いたフーリエは、さらに三角級数を詳しく調べ、1822年に出版した著書『熱の解析的理論』の中で、全ての関数は三角級数で書けるということを主張した。

微分方程式の解の形として、三角級数を仮定するという方法は、フーリエ以前にもダニエル・ベルヌーイらによって行われていたが、三角級数という特別な形を仮定することによって得られる特殊な解と考えられていた。フーリエの主張は、三角級数は、そのような特別なものではなく、全ての関数が三角級数で表せると大きく出ている。

フーリエの議論は飛躍が多かったため、反論が相次ぎ、この主張は受け入れられなかった。しかし、フーリエの側にだけ非があるわけではなく、当時の数学が、このような関数列の収束性などを扱うには未熟で、フーリエの主張の真偽を判定することは難しかったことも関係している。この後、関数がフーリエ級数で表現できるための条件などを論じるために、実数、関数、収束、積分などの概念などの見直しが行われ、フーリエ級数論は19世紀数学における解析学の厳密化に大きな影響を与えることになった。

またフーリエ級数に始まるフーリエ解析の研究は、フーリエ変換などの手法を産み、画像処理やデータ圧縮、CT、MRIなど現代科学の基礎技術としても発展していった。

定義

連続時間信号 $f$ に収束するフーリエ級数が得られるときに $f$ はフーリエ展開できるというが、 $f$ に対する形式的なフーリエ級数が収束するのか、収束するとしても本当に $f$ に収束するのかといった複雑な議論が必要で、これはフーリエ級数の収束性問題と呼ばれる。以下ではこれを考えずに形式的に述べることにする。

フーリエ係数

$n \in ℤ$ を離散変数とし、周期 $T \in ℝ^{+}$ で周期的な複素連続時間信号 $f (t)$ のテンプレート:読み仮名 $c_{n}$ は次式で定義される^[2]：

c_{n} = F [n] = \frac{1}{T} \int_{0}^{T} f (t) e^{- i n ω_{0} t} d t, ω_{0} = \frac{2 π}{T}

$ω_{0}$ は基本角周波数と解釈される^[2]。

フーリエ級数

フーリエ係数を用いて書かれた多項式

\sum_{n = - m}^{m} c_{n} e^{i n ω_{0} x}

を、テンプレート:Mvar 次のフーリエ多項式 (Fourier polynomial) という。このテンプレート:Mvar をテンプレート:Math にした極限

\sum_{n = - \infty}^{\infty} c_{n} e^{i n ω_{0} x} = \lim_{m \to + \infty} \sum_{n = - m}^{m} c_{n} e^{i n ω_{0} x}

をフーリエ級数という。左辺は

\sum_{n = - \infty}^{\infty} c_{n} e^{i n ω_{0} x} = \lim_{k, m \to + \infty} \sum_{n = - k}^{m} c_{n} e^{i n ω_{0} x}

の意味ではないことに注意しなければならない。

実フーリエ級数

実数値関数に限定したフーリエ級数は以下で定義される。

テンプレート:Mvar は、実数テンプレート:Mvar を変数とする実数値関数で、周期テンプレート:Math の周期関数であるとする。

\begin{matrix} a_{n} & = \frac{1}{π} \int_{- π}^{π} f (t) \cos n t d t, (n = 0, 1, 2, 3, \dots) \\ b_{n} & = \frac{1}{π} \int_{- π}^{π} f (t) \sin n t d t, (n = 1, 2, 3, \dots) \end{matrix}

と置き、テンプレート:Math をテンプレート:Mvar のフーリエ余弦係数 (Fourier cosine coefficient)、テンプレート:Math をテンプレート:Mvar のフーリエ正弦係数(Fourier sine coefficient) という。これらを用いて書かれた三角級数

\frac{a_{0}}{2} + \sum_{n = 1}^{\infty} (a_{n} \cos n x + b_{n} \sin n x)

をフーリエ級数(Fourier series) あるいはフーリエ級数展開（Fourier series expansion）という。余弦項だけの

\frac{a_{0}}{2} + \sum_{n = 1}^{\infty} a_{n} \cos n x

を、フーリエ余弦級数といい、正弦項だけの

\sum_{n = 1}^{\infty} b_{n} \sin n x

を、フーリエ正弦級数という。

フーリエ係数を定める積分区間テンプレート:Math に制限してテンプレート:Mvar をみたときにテンプレート:Mvar がフーリエ級数で表される偶関数なら、そのフーリエ級数は余弦級数となり、f(x) がフーリエ級数で表される奇関数なら、そのフーリエ級数は正弦級数となる。

以上に述べた実フーリエ級数は、周期テンプレート:Math の周期関数テンプレート:Mvar に対する定義だが、テンプレート:Math という変数変換により、周期テンプレート:Math の周期関数テンプレート:Math のテンプレート:Math という区間での定義に変換でき、この形で扱われることも少なくない。

\begin{matrix} a_{n} & = \frac{1}{L} \int_{- L}^{L} g (s) \cos (\frac{n π s}{L}) d s, (n = 0, 1, 2, 3, \dots) \\ b_{n} & = \frac{1}{L} \int_{- L}^{L} g (s) \sin (\frac{n π s}{L}) d s, (n = 1, 2, 3, \dots) \\ g (y) & = \frac{a_{0}}{2} + \sum_{n = 1}^{\infty} (a_{n} \cos (\frac{n π y}{L}) + b_{n} \sin (\frac{n π y}{L})) \end{matrix}

パーセバルの等式

テンプレート:Main 関数テンプレート:Math が二乗可積分ならば、以下の等式が成り立つ：

\frac{1}{π} \int_{- π}^{π} | f (t) |^{2} d t = \frac{{a_{0}}^{2}}{2} + \sum_{n = 1}^{\infty} ({a_{n}}^{2} + {b_{n}}^{2}) = 2 \sum_{n = - \infty}^{\infty} | c_{n} |^{2} .

この式はパーセヴァルの等式と呼ばれる。

フーリエ級数の例

テンプレート:Math

周期関数でない関数を周期関数へ拡張し、そのフーリエ級数を扱うことも多い。区間テンプレート:Math で定義される関数として次のような例を考える：

f (x) = x (- π < x < π) .

この関数テンプレート:Mvar を使って以下の周期関数テンプレート:Mvar を定義できる：

\begin{matrix} h (x) & = f (x) & (- π < x < π), \\ h (x + 2 π) & = h (x) & (otherwise) . \end{matrix}

この関数テンプレート:Mvar はテンプレート:Math で定義されない点に注意する。仮に定義したとして、例えば点テンプレート:Math 上では左極限テンプレート:Math と右極限テンプレート:Math が一致せず、これらの点において、値をどのように定義してもテンプレート:Mvar は不連続となる。

以降、記号を粗雑に使い、特に断りのない限り、テンプレート:Mvar の意味でテンプレート:Mvar を用いることにする。

テンプレート:Mvar が区分的に連続微分可能である場合、不連続点でフーリエ級数の収束値は左右からの極限の平均を取るという性質がある。定義した周期関数がフーリエ級数と一致することを求めるなら、テンプレート:Math での値は左右極限の平均値として定義すべきである：

f (π) : = \frac{f (π - 0) + f (π + 0)}{2} .

特に今回の場合、テンプレート:Math となる。

元の関数は奇関数なので、テンプレート:Mvar に対するフーリエ級数は正弦級数となる：

f (x) = 2 \sum_{n = 1}^{\infty} (- 1)^{n - 1} \frac{\sin n x}{n} .

上記よりテンプレート:Math について、以下の等式が得られる：

\frac{π}{2} = 2 \sum_{k = 0}^{\infty} (\frac{1}{4 k + 1} - \frac{1}{4 k + 3}) .

これはライプニッツの公式として知られる。

また、パーセバルの等式より次の関係が得られる：

\frac{1}{π} \int_{- π}^{π} | f (x) |^{2} d x = \frac{2}{3} π^{2} = 4 \sum_{n = 1}^{\infty} \frac{1}{n^{2}} .

最右辺の級数はゼータ関数の特殊値 [[バーゼル問題|テンプレート:Math]] に一致する。

テンプレート:Math

前節と同様に関数テンプレート:Mvar として以下を与える：

f (x) = x^{2} (- π \leq x \leq π) .

この場合、周期関数としてのテンプレート:Mvar は元となる関数の定義域の境界テンプレート:Math およびテンプレート:Math 上で連続となる。

元の関数は偶関数なので、テンプレート:Mvar のフーリエ級数は余弦級数のみで表される：

f (x) = \frac{π^{2}}{3} + 4 \sum_{n = 1}^{\infty} (- 1)^{n} \frac{\cos n x}{n^{2}} .

テンプレート:Indent

さらに、この級数はテンプレート:Math について以下のように整理できる：

\frac{π^{2}}{6} = \sum_{n = 1}^{\infty} \frac{1}{n^{2}} .

ここでもテンプレート:Math が現れる。

直交性

三角級数の直交性

フーリエ級数のようなものが考えられる背景には、関数の直交性がある。 (−π, π) 上で定義された二乗可積分関数の空間 L²(−π, π) を考える。 f(x), g(x) ∈ L²(−π, π) に対して、内積

⟨ f (x), g (x) ⟩ : = \frac{1}{π} \int_{- π}^{π} f (x) g (x)^{*} d x

g(x)^* は g(x) の複素共役であり、実数値のときは、g(x) と等しい

を定義すると、自然数 m, n ≥ 1 に対し

\begin{matrix} ⟨ \cos m x, \cos n x ⟩ & = δ_{m n} \\ ⟨ \sin m x, \sin n x ⟩ & = δ_{m n} \\ ⟨ \cos m x, \sin n x ⟩ & = 0 \\ ⟨ 1, 1 ⟩ & = 2 \\ ⟨ 1, \cos m x ⟩ & = 0 \\ ⟨ 1, \sin m x ⟩ & = 0 \end{matrix}

ただし、δ_mn はクロネッカーのデルタで、内積の中に用いられている 1 というのは、x に依らずに 1 を値にとる定数関数の事とする。

このような関係から

{\frac{1}{\sqrt{2}}, \cos x, \sin x, \cos 2 x, \sin 2 x, \cos 3 x, \sin 3 x, \dots}

は正規直交関数列となり、これは L²(−π, π) の正規直交基底になっている。

\begin{matrix} a_{n} & = ⟨ f (x), \cos n x ⟩ \\ b_{n} & = ⟨ f (x), \sin n x ⟩ \end{matrix}

という計算によって、それぞれ、フーリエ級数の cos nx, sin nx の係数のみを抜き出すことができる。

また、任意の自然数 m について

\begin{matrix} ⟨ f (x), \cos m x ⟩ & = 0 \\ ⟨ f (x), \sin m x ⟩ & = 0 \end{matrix}

が成り立てば、 f(x) = 0 となるため、この直交関数列は完備関数列でもあり、この内積によって、 L²(−π, π) は、ヒルベルト空間になる。

複素型のフーリエ級数の場合も、整数 m, n に対して

⟨ e^{i m x}, e^{i n x} ⟩ = 2 π δ_{m n}

という直交関係がなりたち、{e^imx} は完備関数列になる。

ヒルベルト空間とフーリエ級数

ヒルベルト空間 X と、その正規直交系 {e_k} を考える。 x ∈ X に対して、その内積 $⟨ x, e_{k} ⟩$ のことをフーリエ係数という。この時、ベッセルの不等式

‖ x ‖_{2}^{2} \geq \sum_{k} {| ⟨ x, e_{k} ⟩ |}^{2}

が成り立つ。

さらに {e_k} が X の基底となっていれば、三角級数のときと同様に級数

\sum_{k} ⟨ x, e_{k} ⟩ e_{k} = ⟨ x, e_{1} ⟩ e_{1} + ⟨ x, e_{2} ⟩ e_{2} + \dots

が考えられ、これも同じようにフーリエ級数という。この級数が、元の x に等しいとき、フーリエ展開できるという。そしてこの時、プランシュレルの等式

‖ x ‖_{2}^{2} = \sum_{k} {| ⟨ x, e_{k} ⟩ |}^{2}

が成り立つ。

ヒルベルト空間 X について、

任意の x ∈ X がフーリエ展開できること
任意の x ∈ X に対し、プランシュレルの等式が成り立つこと
{e_k} が X の正規直交基底であること

の 3つは互いに同値な条件である。

主な周期関数のフーリエ級数

主な周期関数のフーリエ級数を示す。

関数 $f (x)$ は、期間 $0 < x \leq T$ で定義されるものとする。
$a_{0}, a_{n}, b_{n}$ はそれぞれ、直流成分、フーリエ余弦係数、フーリエ正弦係数である。

時間領域 $f (x)$	周波数領域 $\begin{matrix} a_{0} \\ a_{n} for n \geq 1 \\ b_{n} for n \geq 1 \end{matrix}$	備考	出典
$f (x) = A \| \sin (\frac{2 π}{T} x) \| for 0 \leq x < T$	$\begin{matrix} a_{0} = & \frac{4 A}{π} \\ a_{n} = & {\begin{matrix} \frac{- 4 A}{π} \frac{1}{n^{2} - 1} & n even \\ 0 & n odd \end{matrix} \\ b_{n} = & 0 \end{matrix}$	正弦波の全波整流波形	^[3]テンプレート:Rp
$f (x) = {\begin{matrix} A \sin (\frac{2 π}{T} x) & for 0 \leq x < T / 2 \\ 0 & for T / 2 \leq x < T \end{matrix}$	$\begin{matrix} a_{0} = & \frac{2 A}{π} \\ a_{n} = & {\begin{matrix} \frac{- 2 A}{π} \frac{1}{1 - n^{2}} & n even \\ 0 & n odd \end{matrix} \\ b_{n} = & {\begin{matrix} \frac{A}{2} & n = 1 \\ 0 & n > 1 \end{matrix} \end{matrix}$	正弦波の半波整流波形	^[3]テンプレート:Rp
$f (x) = {\begin{matrix} A & for 0 \leq x < D \cdot T \\ 0 & for D \cdot T \leq x < T \end{matrix}$	$\begin{matrix} a_{0} = & 2 A D \\ a_{n} = & \frac{A}{n π} \sin (2 π n D) \\ b_{n} = & \frac{2 A}{n π} {(\sin (π n D))}^{2} \end{matrix}$	$0 \leq D \leq 1$ デューティ比 $D$ の矩形波
$f (x) = \frac{A x}{T} for 0 \leq x < T$	$\begin{matrix} a_{0} = & A \\ a_{n} = & 0 \\ b_{n} = & \frac{- A}{n π} \end{matrix}$	ノコギリ波（増）	^[3]テンプレート:Rp
$f (x) = A - \frac{A x}{T} for 0 \leq x < T$	$\begin{matrix} a_{0} = & A \\ a_{n} = & 0 \\ b_{n} = & \frac{A}{n π} \end{matrix}$	ノコギリ波（減）	^[3]テンプレート:Rp
$f (x) = \frac{4 A}{T^{2}} {(x - \frac{T}{2})}^{2} for 0 \leq x < T$	$\begin{matrix} a_{0} = & \frac{2 A}{3} \\ a_{n} = & \frac{4 A}{π^{2} n^{2}} \\ b_{n} = & 0 \end{matrix}$		^[3]テンプレート:Rp

出典

テンプレート:脚注ヘルプ

↑ テンプレート:Cite book
↑ ^2.0 ^2.1 テンプレート:Ilqp.1 より引用。テンプレート:Cite journal
↑ ^3.0 ^3.1 ^3.2 ^3.3 ^3.4 テンプレート:Cite book

フーリエ級数

目次

概要