物理を解説 ♪
ツイッター用のシェアボタン フェイスブック用のシェアボタン はてなブックマーク用のシェアボタン ライン用のシェアボタン
ツイッター用のシェアボタン フェイスブック用のシェアボタン はてなブックマーク用のシェアボタン ライン用のシェアボタン

フーリエ級数の基本

ウェブアプリで遊んで行ってちょうだい。
作成:2012/12/19

驚くべき公式!

まずは0≦x≦2πの範囲で定義された連続な関数f(x)を考える.この関数がどんな形をしていようとも三角関数の足し合わせで表現できそうだという驚くべき内容をフランスの学者フーリエが論文中で使い,それが本当なのかどうかを巡って議論が沸き起こったのであった.19 世紀初めのことである. 数式 波長が2πsin波とcos,その1/2の波長のsin波とcos,1/3の波長のsin波とcos,・・・というように,どんどん細かく上下するようになる波を次々と色んな振幅で重ね合わせていくのである.

そんなことで本当に「どんな形でも」表せるのだろうか

いや,そうはいくまい.sin波もcos波も上下に同じだけ振動していて平均すれば 0 なので,そのようなものをどれだけ重ね合わせたとしても平均は 0 だろう.だから平均が 0 になるような形の関数しか表せないことになる.しかしそのような弱点を補うために (1) 式には平均値であるcを入れておいた.それならどうだろう

それが本当であることを実感してもらえるようにウェブアプリを用意してみた.

フーリエ級数のウェブアプリのサムネイル

自由に波形を描いてみて欲しい.どんな形でも最終的にはかなり正確に再現してくれるはずだ.


係数を求める方法

なぜこのようなことが可能なのかという証明は放っておくことにしよう.それよりも (1) 式に出てくる係数a_nb_nをどのように決めたら (1) 式が成り立つように出来るのかを説明したい.

その前に,cは関数f(x)の平均値なので次のように計算すれば良いことは分かるはずだ. 数式 係数a_nb_nもこれに少し似ていて,次のようにして求めるのである. 数式 なぜπで割る2πで割るのではないのなぜsincosを掛けて積分する色んな疑問が出るかも知れないが,徐々に解決してゆこう.

この計算を見ていると,例えばa_nを求めるときにはcos(nx)f(x)を掛けたものを積分している.つまり,0≦x≦2πの範囲内でf(x)cos(nx)と似た動きをしていれば結果は大きめに出て,合わない動き方をしていれば,結果は打ち消されて小さめに出てきそうだと想像できる.要するにこれは,f(x)の中からcos(nx)に似た成分がどれだけあるかを抜き出してくる操作なのであろう.そのために0≦x≦2πの範囲に渡って積分したので,それを平均するために2πで割るというのなら何となく意味は繋がる気がするのだが,なぜかπだけで割っている.結果を 2 倍せねばならぬ事情がありそうだ.b_nについても同じようなものである.

(2) 式と (3) 式は形式が似ている.(3) 式のa_nの式でn=0とすれば,cos0=1であるので積分のところは同じ形になる.するとa_{0}cとは係数が違うだけであり,c=a_{0}/2だと言えそうだ.だから (1) 式を次のように表しておけば (2) 式は不要になるだろう. 数式 意味は分かりにくくなるが,式の数を一つ減らせて,公式を書くためのスペースと手間を節約できるという利点がある.

関数を (1) 式や (1') 式のように無限に続く三角関数の和の形で表したものを「フーリエ級数」と呼ぶ.関数の形によっては有限項で終わる場合もあり,その場合でもフーリエ級数と呼んで構わない.


公式が成り立つ理由

しかし (3) 式で係数が求められるというのはなぜだろうか.この辺りのことを理解するために,次のような公式を知っていると助けになる. 数式 ただしmnは正の整数であるとする.右辺のδ_{mn}は「クロネッカーのデルタ」というもので,mnが等しければ 1 で,それ以外は 0 であることを意味している.この公式は三角関数の積和の公式を使えば簡単に導けるので説明を省略したいところだが,m=nとなる場合とm≠nとなる場合とで状況が異なることに気付かないと混乱する可能性があるので一つだけ例を示しておこう. 数式 この計算はm≠nの場合には問題ないが,m=nでは分母が 0 になってしまうところがあって正しくない.実はm=nの場合には積分する前にcos(m-n)x=cos0=1となっている.そのことに気付けばこの問題は回避できて,違った結果が得られることになるだろう.

(4) 式はとても重要なことに気付かせてくれる.\{cosx,cos2x,cos3x,…,sinx,sin2x,sin3x,…\}という関数は,互いに掛け合わせて積分した時,どの組み合わせを取ってみても 0 にしかならないただ自分自身と掛け合わせた時に限ってπになるのである

この知識を利用しよう.(3) 式の右辺にはf(x)が入っている.f(x)は (1) 式のように表されるというのを仮定だと考えてやって,これを (3) 式の右辺に代入してやると,その計算結果はどうなるだろうか(4) 式を利用してやれば,ほとんどの項は消え去ることが分かるだろう.残る項は一つだけであって,その係数部分しか残らない.このようにして (3) 式が正しいことが示されることになる.先ほどの「全体を2πで割るべきところがπで割られているのはなぜか」という疑問はあまり意味がなくて,ただ (4) 式がそういう形になっているから,というだけの事だったようだ.

何か騙されたような気がするかもしれないし,循環論法的に感じるかも知れない.これではどうも説明になっていない感じがする.実際その通りだ.今のところ,関数f(x)が (1) 式のように表せると仮定すれば,そこで使われている係数は (3) 式のようであるべきだということを説明しただけであって,どんな関数の場合にでも (1) 式のように等式が成り立つという点についてはまだ解決していない.この点については昔の学者たちもすぐには認めることができなかったのである.

本当にこんなものであらゆる関数を表すことができるのだろうか


周期を変えて公式を作ってみる

ここまでは0≦x≦2πの範囲だけで考えていたが,sin関数もcos関数も周期関数なのでこの範囲外であっても全く同じ振る舞いを何度も繰り返すだけである.(1) 式のように表された関数f(x)についても周期2πで同じ動きを繰り返すのである.

ここまでに出てきた公式では全て0≦x≦2πの範囲で積分していたのだが,一つの周期に渡って積分すれば結果は同じなのだから,例えば-π≦x≦πのような範囲で積分しても同じことである.

しかし周期が2πに限られているのはどうにも不自由さを感じる.周期を好きに設定できるように公式を改造できないだろうか.例えば (1) 式を次のように変更すれば,周期がLで繰り返すようにできそうだ.x=Lの時にどうなるかを考えてみれば納得が行くだろう. 数式 (2) 式の代わりには次のようなものを計算すればいいだろう. 数式 ではa_nb_nはどうなるだろうかそれを探るために,(4) 式に代わるものを計算してみよう.やることは大して変わらないので結果だけ書くことにする. 数式 なるほど,先ほどの話と比べてほとんど変更はない.そんなに難しいことを考える必要は無さそうだ.係数a_nb_nを次のように決めておけば話が合うだろう. 数式 cはやはりc=a_{0}/2とすることで (6) 式に吸収できそうである.(5) 式を書き直しておこう. 数式 この (5') 式と (6) 式が,周期がLになるように拡張したフーリエ級数の公式である.

教科書によっては[-L,L]の範囲で積分してあるものがあるが,その場合,周期は2Lになるので上の公式のL2Lに置き換えれば同じ形になり,話は合うだろう.積分範囲については周期と同じ幅になっていればどう選んだって構わないのである.

フーリエ級数は,積分した範囲のf(x)の形と同じ形を周期Lで何度も何度も繰り返すような関数を再現してくれることになる.


偶関数と奇関数

sin関数は奇関数であり,cos関数は偶関数である.もしどんな関数でもフーリエ級数のように表せるとしたならば,どんな関数でも,偶関数と奇関数に分けて表せるということになる.

まぁ,それについてはフーリエ級数に頼らなくてもいつでも言えることではある. 数式 としておけば,g(-x)=-g(x)となるのでg(x)は奇関数だし,h(-x)=h(x)となるのでh(x)は偶関数だし, 数式 なので,f(x)は偶関数と奇関数に分けて表せたことになるからである.

本当に言いたいのはそのことではないのだった.f(x)が偶関数ならcos関数だけの項で表せるし,f(x)が奇関数ならsin関数だけの和で表せるだろうということを記憶に留めておいてもらいたいのである.f(x)が偶関数なら全てのb_nは 0 になるし,f(x)が奇関数なら全てのa_nは 0 になる.なぜちゃんとそんなことになるのかを考えるのは読者に任せよう.偶関数と奇関数の積は奇関数になるとか,奇関数と奇関数の積は偶関数になるだとかはちゃんと知ってるだろうかその辺りを使えばいい.

a_nが全て 0 でsin関数ばかりの項で出来たフーリエ級数のことを「フーリエ正弦級数」と呼び,b_nが全て 0 で,定数a_{0}cos関数ばかりの項で出来たフーリエ級数のことを「フーリエ余弦級数」と呼ぶ.

基礎知識として知っておけばいいことはだいたいこれくらいだろうと思う.



趣味の物理学書店の案内バナー

EMAN物理note出張所の案内バナー