フーリエ級数 (フーリエきゅうすう、Fourier series )とは、複雑な周期関数 や周期信号を、単純な形の周期性をもつ関数の(無限の)和によって表したものである。フーリエ級数は、フランスの数学者ジョゼフ・フーリエ によって金属板の中での熱伝導に関する研究の中で導入された。
熱伝導方程式 は、偏微分方程式 として表される。フーリエの研究の前までには、一般的な形での熱伝導方程式の解法は知られておらず、熱源が単純な形である場合、例えば正弦波 などの場合の特別な解しかえられていなかった。この特別な解は現在では固有解と呼ばれる。フーリエの発想は、複雑な形をした熱源をサイン波、コサイン波の和 として考え、解を固有解の和として表すものであった。
この重ね合わせがフーリエ級数と呼ばれる。
最初の動機は熱伝導方程式を解くことであったが、数学や物理の他の問題にも同様のテクニックが使えることが分かり様々な分野に応用されている。
フーリエ級数は、電気工学 、振動 の解析、音響学 、光学 、信号処理 、量子力学 および経済学 [1] などの分野で用いられている。
フーリエ級数は、関数に対して定義されるフーリエ係数を用いて
a
0
2
+
∑
k
=
1
∞
(
a
k
cos
k
x
+
b
k
sin
k
x
)
{\displaystyle {\frac {a_{0}}{2}}+\sum _{k=1}^{\infty }(a_{k}\cos kx+b_{k}\sin kx)}
の形に表される三角級数のことである。熱方程式 を発見したフーリエ は、平衡状態 における熱方程式に注目し、適当な境界条件の下で二変数のラプラス方程式
(
∂
2
∂
x
2
+
∂
2
∂
y
2
)
ϕ
(
x
,
y
)
=
0
{\displaystyle \left({\partial ^{2} \over \partial x^{2}}+{\partial ^{2} \over \partial y^{2}}\right)\phi (x,y)=0}
に帰着させて解を求めようとした。この時、フーリエは、
∑
k
=
0
∞
(
−
1
)
k
2
k
+
1
cos
(
(
2
k
+
1
)
x
)
=
π
4
,
(
−
π
2
<
x
<
π
2
)
{\displaystyle \sum _{k=0}^{\infty }{(-1)^{k} \over 2k+1}\cos((2k+1)x)={\pi \over 4},\left(-{\pi \over 2}<x<{\pi \over 2}\right)}
という三角級数を見つけている。左辺の三角関数の一つ一つは波打っているにもかかわらず、x に依らない定数に収束しているのである。
x = 0 としたときの級数は円周率 を求めるグレゴリー級数 と同じである。x の定義域を広げるとこの三角級数は n を整数として
∑
k
=
0
∞
(
−
1
)
k
2
k
+
1
cos
(
(
2
k
+
1
)
x
)
=
(
−
1
)
n
π
4
,
(
−
π
2
+
n
π
<
x
<
π
2
+
n
π
)
{\displaystyle \sum _{k=0}^{\infty }{(-1)^{k} \over 2k+1}\cos((2k+1)x)=(-1)^{n}{\pi \over 4},\left(-{\pi \over 2}+n\pi <x<{\pi \over 2}+n\pi \right)}
という矩形波 になる。このような不連続な関数まで表せることに興味を抱いたフーリエは、さらに三角級数を詳しく調べ、1822年 に出版した著書『熱の解析的理論』の中で、全ての関数は三角級数で書けるということを主張した。
微分方程式 の解の形として、三角級数を仮定するという方法は、フーリエ以前にもダニエル・ベルヌーイ らによって行われていたが、三角級数という特別な形を仮定することによって得られる特殊な解と考えられていた。フーリエの主張は、三角級数は、そのような特別なものではなく、全ての関数が三角級数で表せると大きく出ている。
フーリエの議論は飛躍が多かったため、反論が相次ぎ、この主張は受け入れられなかった。しかし、フーリエの側にだけ非があるわけではなく、当時の数学 が、このような関数列の収束性などを扱うには未熟で、フーリエの主張の真偽を判定することは難しかったことも関係している。この後、関数がフーリエ級数で表現できるための条件などを論じるために、実数 、関数 、収束 、積分 などの概念などの見直しが行われ、フーリエ級数論は19世紀数学における解析学 の厳密化に大きな影響を与えることになった。
またフーリエ級数に始まるフーリエ解析 の研究は、フーリエ変換 などの手法を産み、画像処理 やデータ圧縮 、CT 、MRI など現代科学の基礎技術としても発展していった。
以下で述べる定義は形式的なもので、実際には f (x ) を用いた積分が存在するのかということと、f (x ) から得られたフーリエ級数が、本当に f (x ) に収束するのかといった事が問題になり、それらを解決するためには f (x ) に二乗可積分 などの制約が課される。積分 にはルベーグ積分 が用いられる。 f (x ) に収束するフーリエ級数が得られる場合、 f (x ) はフーリエ展開 できるという。
実数値関数のフーリエ級数 編集
f (x ) は、実数 x を変数 とする実数値関数で、周期 2π の周期関数 とする。
a
n
=
1
π
∫
−
π
π
f
(
t
)
cos
n
t
d
t
,
(
n
=
0
,
1
,
2
,
3
,
⋯
)
b
n
=
1
π
∫
−
π
π
f
(
t
)
sin
n
t
d
t
,
(
n
=
1
,
2
,
3
,
⋯
)
{\displaystyle {\begin{aligned}a_{n}&={1 \over \pi }\int _{-\pi }^{\pi }f\left(t\right)\cos nt\,dt,\left(n=0,1,2,3,\cdots \right)\\b_{n}&={1 \over \pi }\int _{-\pi }^{\pi }f\left(t\right)\sin nt\,dt,\left(n=1,2,3,\cdots \right)\end{aligned}}}
と置き、a n を f のフーリエ余弦係数 (Fourier cosine coefficient )、b n を f のフーリエ正弦係数 (Fourier sine coefficient ) という。これらを用いて書かれた三角級数
a
0
2
+
∑
n
=
1
∞
(
a
n
cos
n
x
+
b
n
sin
n
x
)
{\displaystyle {a_{0} \over 2}+\sum _{n=1}^{\infty }\left(a_{n}\cos nx+b_{n}\sin nx\right)}
をフーリエ級数 (Fourier series ) あるいはフーリエ級数展開 (Fourier series expansion )という。余弦項だけの
a
0
2
+
∑
n
=
1
∞
a
n
cos
n
x
{\displaystyle {a_{0} \over 2}+\sum _{n=1}^{\infty }a_{n}\cos nx}
を、フーリエ余弦級数 といい、正弦項だけの
∑
n
=
1
∞
b
n
sin
n
x
{\displaystyle \sum _{n=1}^{\infty }b_{n}\sin nx}
を、フーリエ正弦級数 という。
フーリエ係数を定める積分区間 −π < x < π に制限して f (x ) をみたときに f (x ) がフーリエ級数で表される偶関数 なら、そのフーリエ級数は余弦級数となり、f (x ) がフーリエ級数で表される奇関数 なら、そのフーリエ級数は正弦級数となる。 複素数値関数のフーリエ級数(複素フーリエ級数) 編集
オイラーの公式 を用いると、複素型のフーリエ級数を得ることができる。 f (x ) も複素数値に取ることができ
c
n
=
1
2
π
∫
−
π
π
f
(
t
)
exp
(
−
i
n
t
)
d
t
,
(
n
=
0
,
±
1
,
±
2
,
⋯
)
{\displaystyle c_{n}={\frac {1}{2\pi }}\int _{-\pi }^{\pi }f(t)\exp(-int)dt,\left(n=0,\pm 1,\pm 2,\cdots \right)}
を、 f のフーリエ係数 (Fourier coefficient )といい、これを用いて書かれた多項式
∑
n
=
−
m
m
c
n
e
i
n
x
{\displaystyle \sum _{n=-m}^{m}c_{n}e^{inx}}
を、 m 次のフーリエ多項式 (Fourier polynomial )という。この m を +∞ にした極限
∑
n
=
−
∞
∞
c
n
e
i
n
x
=
lim
m
→
+
∞
∑
n
=
−
m
m
c
n
e
i
n
x
{\displaystyle \sum _{n=-\infty }^{\infty }c_{n}e^{inx}=\lim _{m\to +\infty }\sum _{n=-m}^{m}c_{n}e^{inx}}
をフーリエ級数 という。左辺は
∑
n
=
−
∞
∞
c
n
e
i
n
x
=
lim
k
,
m
→
+
∞
∑
n
=
−
k
m
c
n
e
i
n
x
{\displaystyle \sum _{n=-\infty }^{\infty }c_{n}e^{inx}=\lim _{k,m\to +\infty }\sum _{n=-k}^{m}c_{n}e^{inx}}
の意味ではないことに注意しなければならない。
周期の変更 編集
以上に述べたフーリエ級数は、周期 2π の周期関数 f (x ) に対する定義だが、 x = (π/L ) y という変数変換により、周期 2L の周期関数 g (y ) = f ((π/L )y ) の −L ≤ y ≤ L という区間での定義に変換でき、この形で扱われることも少なくない。
a
n
=
1
L
∫
−
L
L
g
(
s
)
cos
(
n
π
s
L
)
d
s
,
(
n
=
0
,
1
,
2
,
3
,
⋯
)
b
n
=
1
L
∫
−
L
L
g
(
s
)
sin
(
n
π
s
L
)
d
s
,
(
n
=
1
,
2
,
3
,
⋯
)
g
(
y
)
=
a
0
2
+
∑
n
=
1
∞
(
a
n
cos
(
n
π
y
L
)
+
b
n
sin
(
n
π
y
L
)
)
c
n
=
1
2
L
∫
−
L
L
g
(
s
)
exp
(
−
i
n
π
s
L
)
d
s
,
(
n
=
0
,
±
1
,
±
2
,
⋯
)
g
(
y
)
=
∑
n
=
−
∞
∞
c
n
exp
(
i
n
π
y
L
)
=
lim
m
→
+
∞
∑
n
=
−
m
m
c
n
exp
(
i
n
π
y
L
)
{\displaystyle {\begin{aligned}a_{n}&={1 \over L}\int _{-L}^{L}g\left(s\right)\cos \left({n\pi s \over L}\right)ds,\left(n=0,1,2,3,\cdots \right)\\b_{n}&={1 \over L}\int _{-L}^{L}g\left(s\right)\sin \left({n\pi s \over L}\right)ds,\left(n=1,2,3,\cdots \right)\\g(y)&={a_{0} \over 2}+\sum _{n=1}^{\infty }\left(a_{n}\cos \left({n\pi y \over L}\right)+b_{n}\sin \left({n\pi y \over L}\right)\right)\\c_{n}&={1 \over 2L}\int _{-L}^{L}g\left(s\right)\exp \left(-{in\pi s \over L}\right)ds,\left(n=0,\pm 1,\pm 2,\cdots \right)\\g(y)&=\sum _{n=-\infty }^{\infty }c_{n}\exp \left({in\pi y \over L}\right)=\lim _{m\to +\infty }\sum _{n=-m}^{m}c_{n}\exp \left({in\pi y \over L}\right)\end{aligned}}}
パーセバルの等式 編集
f (x ) が二乗可積分関数ならば、パーセバルの等式 (Parseval's equality )
1
π
∫
−
π
π
|
f
(
t
)
|
2
d
t
=
a
0
2
2
+
∑
n
=
1
∞
(
a
n
2
+
b
n
2
)
=
2
∑
n
=
−
∞
∞
|
c
n
|
2
{\displaystyle {1 \over \pi }\int _{-\pi }^{\pi }|f(t)|^{2}dt={{a_{0}}^{2} \over 2}+\sum _{n=1}^{\infty }({a_{n}}^{2}+{b_{n}}^{2})=2\sum _{n=-\infty }^{\infty }|c_{n}|^{2}}
が成り立つ。
フーリエ級数の例 編集
周期関数以外の関数から周期関数を作り計算することも多い。例えば
を考え、この区間に含まれない x での値 f (x ) は、 n を整数として、周期関数の性質
によって定義し、周期関数に拡張する。今、 f (n π) の値は定義していない。この関数では、 f (π−0) ≠ f (−π+0) となっているため、この点ではどのように決めても不連続になってしまうためである。
f (x ) が区分的に連続微分可能 の時、不連続点でフーリエ級数の収束値は左右からの極限の平均を取るという性質があり、この場合であれば
f
(
π
−
0
)
+
f
(
π
+
0
)
2
=
π
−
π
2
=
0
{\displaystyle {f(\pi -0)+f(\pi +0) \over 2}={\pi -\pi \over 2}=0}
となるため、 f (x ) が、 フーリエ級数と一致するように定めるとすれば、 f (n π) = 0 となる。 −π < x < π という区間内に限ってみると、 f (x ) は、奇関数であるので、この場合は正弦級数となり
f
(
x
)
=
2
∑
n
=
1
∞
(
−
1
)
n
−
1
sin
n
x
n
{\displaystyle f(x)=2\sum _{n=1}^{\infty }(-1)^{n-1}{\sin nx \over n}}
というフーリエ級数展開になる。 x = π/2 を代入し、 2 で割れば、円周率の公式としてよく知られた級数(ライプニッツの公式 )
π
4
=
∑
k
=
0
∞
(
−
1
)
k
2
k
+
1
=
1
−
1
3
+
1
5
−
1
7
+
⋯
{\displaystyle {\pi \over 4}=\sum _{k=0}^{\infty }{(-1)^{k} \over 2k+1}=1-{1 \over 3}+{1 \over 5}-{1 \over 7}+\cdots }
が得られる。
また、
1
π
∫
−
π
π
|
t
|
2
d
t
=
2
3
π
2
{\displaystyle {1 \over \pi }\int _{-\pi }^{\pi }|t|^{2}dt={2 \over 3}\pi ^{2}}
だから、パーセバルの等式を用いることにより
2
3
π
2
=
4
∑
n
=
1
∞
1
n
2
{\displaystyle {2 \over 3}\pi ^{2}=4\sum _{n=1}^{\infty }{1 \over n^{2}}}
この両辺を 4 で割ることにより、オイラー によって得られたゼータ関数 での特殊値 ζ(2) が現れる。
の場合、π2 = (−π)2 であり、これは周期関数に拡張しても、区間の端点が不連続になることはないので、両端を含む閉区間上で定義している。 −π < x < π に制限して考えると偶関数なので、今度は余弦級数で、
f
(
x
)
=
π
2
3
+
4
∑
n
=
1
∞
(
−
1
)
n
cos
n
x
n
2
{\displaystyle f(x)={\pi ^{2} \over 3}+4\sum _{n=1}^{\infty }(-1)^{n}{\cos nx \over n^{2}}}
になる。
x 2 を微分して 2 で割ると x になるのと同じように、この右辺の級数を項別微分 して 2 で割ると、先程の f (x ) = x のフーリエ級数になる。一般には、こうはならないが、f (x ) のフーリエ級数を項別微分して得られる級数が一様収束 していれば、それは、f ′(x ) のフーリエ級数に一致する。
さらに、この級数は、 x = π を代入して整理すると
π
2
6
=
∑
n
=
1
∞
1
n
2
{\displaystyle {\pi ^{2} \over 6}=\sum _{n=1}^{\infty }{1 \over n^{2}}}
となり、ここでも ζ(2) が現れる。
三角級数の直交性 編集
フーリエ級数のようなものが考えられる背景には、関数の直交性がある。 (−π, π) 上で定義された二乗可積分関数の空間 L 2 (−π, π) を考える。 f (x ), g (x ) ∈ L 2 (−π, π) に対して、内積
⟨
f
(
x
)
,
g
(
x
)
⟩
:=
1
π
∫
−
π
π
f
(
x
)
g
(
x
)
∗
d
x
{\displaystyle \left\langle f(x),g(x)\right\rangle :={1 \over \pi }\int _{-\pi }^{\pi }f(x)g(x)^{*}dx}
g (x )* は g (x ) の複素共役 であり、実数値のときは、g (x ) と等しいを定義すると、自然数 m , n ≥ 1 に対し
⟨
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
{\displaystyle {\begin{aligned}&\langle \cos mx,\cos nx\rangle =\delta _{mn}\\&\langle \sin mx,\sin nx\rangle =\delta _{mn}\\&\langle \cos mx,\sin nx\rangle =0\\&\langle 1,1\rangle =2\\&\langle 1,\cos mx\rangle =0\\&\langle 1,\sin mx\rangle =0\end{aligned}}}
ただし、δmn はクロネッカーのデルタ で、内積の中に用いられている 1 というのは、x に依らずに 1 を値にとる定数関数 の事とする。
このような関係から
{
1
2
,
cos
x
,
sin
x
,
cos
2
x
,
sin
2
x
,
cos
3
x
,
sin
3
x
,
…
}
{\displaystyle \left\{{1 \over {\sqrt {2}}},\cos x,\sin x,\cos 2x,\sin 2x,\cos 3x,\sin 3x,\ldots \right\}}
は正規直交関数列 となり、これは L 2 (−π, π) の正規直交基底 になっている。
a
n
=
⟨
f
(
x
)
,
cos
n
x
⟩
b
n
=
⟨
f
(
x
)
,
sin
n
x
⟩
{\displaystyle {\begin{aligned}a_{n}&=\langle f(x),\cos nx\rangle \\b_{n}&=\langle f(x),\sin nx\rangle \end{aligned}}}
という計算によって、それぞれ、フーリエ級数の cos nx , sin nx の係数のみを抜き出すことができる。
また、任意の自然数 m について
⟨
f
(
x
)
,
cos
m
x
⟩
=
0
⟨
f
(
x
)
,
sin
m
x
⟩
=
0
{\displaystyle {\begin{aligned}\langle f(x),\cos mx\rangle &=0\\\langle f(x),\sin mx\rangle &=0\end{aligned}}}
が成り立てば、 f (x ) = 0 となるため、この直交関数列は完備関数列 でもあり、この内積によって、 L 2 (−π, π) は、ヒルベルト空間 になる。
複素型のフーリエ級数の場合も、整数 m , n に対して
⟨
e
i
m
x
,
e
i
n
x
⟩
=
2
π
δ
m
n
{\displaystyle \langle e^{imx},e^{inx}\rangle =2\pi \delta _{mn}}
という直交関係がなりたち、{e imx } は完備関数列になる。
ヒルベルト空間とフーリエ級数 編集
ヒルベルト空間 X と、その正規直交系 {e k } を考える。 x ∈ X に対して、その内積
⟨
x
,
e
k
⟩
{\displaystyle \langle x,e_{k}\rangle }
のことをフーリエ係数 という。この時、ベッセルの不等式
‖
x
‖
2
2
≥
∑
k
|
⟨
x
,
e
k
⟩
|
2
{\displaystyle \|x\|_{2}^{2}\geq \sum _{k}\left|\langle x,e_{k}\rangle \right|^{2}}
が成り立つ。
さらに {e k } が X の基底となっていれば、三角級数のときと同様に級数
∑
k
⟨
x
,
e
k
⟩
e
k
=
⟨
x
,
e
1
⟩
e
1
+
⟨
x
,
e
2
⟩
e
2
+
⋯
{\displaystyle \sum _{k}\langle x,e_{k}\rangle e_{k}=\left\langle x,e_{1}\right\rangle e_{1}+\left\langle x,e_{2}\right\rangle e_{2}+\cdots }
が考えられ、これも同じようにフーリエ級数 という。この級数が、元の x に等しいとき、フーリエ展開できるという。そしてこの時、プランシュレルの等式
‖
x
‖
2
2
=
∑
k
|
⟨
x
,
e
k
⟩
|
2
{\displaystyle \|x\|_{2}^{2}=\sum _{k}\left|\langle x,e_{k}\rangle \right|^{2}}
が成り立つ。
ヒルベルト空間 X について、
任意の x ∈ X がフーリエ展開できること
任意の x ∈ X に対し、プランシュレルの等式が成り立つこと
{e k } が X の正規直交基底であること の 3つは互いに同値な条件である。
^ Nerlove, Marc; Grether, David M.; Carvalho, Jose L. (1995). Analysis of Economic Time Series. Economic Theory, Econometrics, and Mathematical Economics . Elsevier. ISBN 0125157517 .