フーリエ解析 · フーリエ級数

フーリエ級数

フーリエ級数は周期関数を正弦と余弦の和、a₀/2 + Σ (aₖ cos kx + bₖ sin kx)で表します。各係数は独立に測れます。関数に1つの波を掛け、その面積を求めます。

この可視化で見るフーリエ級数

上のグラフには、周期2π2\piの関数ffを黒で、その部分和SNS_Nを黄色で描きます。部分和には定数項と、周波数NNまでのすべての余弦・正弦を含めます。中央の帯ではffに1つのプローブ波を掛け、積と軸の間の符号付き面積を色で示します。下の帯には、こうして求めた係数を、余弦は上、正弦は下に並べます。

係数は面積から求まる
ak=1π∫−ππf(x)cos⁡kx dx,bk=1π∫−ππf(x)sin⁡kx dxa_k = \frac{1}{\pi}\int_{-\pi}^{\pi} f(x)\cos kx\,dx, \qquad b_k = \frac{1}{\pi}\int_{-\pi}^{\pi} f(x)\sin kx\,dx

ffにはすべての波が含まれているのに、なぜ1つの面積で1つの係数を取り出せるのでしょうか。1周期全体では、1,cos⁡x,sin⁡x,cos⁡2x,sin⁡2x,…1, \cos x, \sin x, \cos 2x, \sin 2x, \dotsの異なる2つの波の積は符号付き面積が0になります。一方、正弦や余弦をそれ自身と掛けたときの面積はπ\piです。したがって、これらの波の和にsin⁡kx\sin kxを掛けて積分すると、bksin⁡kxb_k \sin kx以外の項がすべて消え、πbk\pi b_kが残ります。2つの正弦の例では、帯の中で積を2つに分け、異なる周波数の成分が打ち消し合う様子を示します。波は互いに垂直な軸のように振る舞い、各係数はその軸への射影に当たります。そのため、項を加えても、すでに求めた係数は変わりません。

十分な項を加えると、部分和は跳びのそばで元の関数を越えます。NNが増えるにつれてピークは細くなり、跳びに近づきますが、オーバーシュートは消えません。矩形波では、上端を越える量が跳びの高さの約8.949%に近づきます。各SNS_Nは連続なので、短い区間を使って跳びを横切る必要があり、跳びを含むどの区間でも最大誤差は0に近づけません。ただし、連続性だけで超過量が決まるわけではありません。8.949%という値は部分和を計算すると得られ、区分的に滑らかな関数のどの跳びでも、最終的に同じ割合になります。

滑らかな関数ほど、少ない項で近づけます。ここでffは無限に繰り返す関数なので、±π\pm\piの両端が一致しなければ、そこも跳びになります。部分積分をすると、波に掛かっていた微分がffへ移り、因子1/k1/kが現れます。跳びがあると、最初の段階で因子1/k1/kを含む境界項が残ります。矩形波の係数はkkが奇数ならbk=4/(kπ)b_k = 4/(k\pi)、kkが偶数なら0です。ffが連続で傾きだけが跳ぶ場合、たとえば三角波の角では、最初の段階はその境界項を残さずに進み、2段階目で項が残るため、係数は1/k21/k^2のように減衰します。対数尺度で、棒を破線1/k1/k、1/k21/k^2と比べてください。

エネルギー表示では、同じ直交性を平方に使います。f2f^2の積分は交差項を含まず、係数の平方に分かれます。1π∫−ππf(x)2 dx=a022+∑k≥1(ak2+bk2).\frac{1}{\pi}\int_{-\pi}^{\pi} f(x)^2\,dx = \frac{a_0^2}{2} + \sum_{k \ge 1} \left(a_k^2 + b_k^2\right). 表示されるのは、この合計のうちSNS_Nがすでに持つ割合です。残りは同じ積分で測った二乗誤差に当たります。矩形波ではS1S_1が約81%、S9S_9が約96%を持ちますが、オーバーシュートはまだ残っています。細いピークが加える面積は小さいからです。部分和は平均的にはffに近くても、跳びのそばではずれることがあります。

注目する点

  • 定数項をa0/2a_0/2と書くのは、余弦係数の式をk=0k = 0にも適用できるようにするためです。定数プローブはffの平均値a0/2a_0/2を直接測ります。
  • 矩形波のような奇関数には正弦項だけが、∣sin⁡x∣|\sin x|のような偶関数には余弦項だけが現れるため、スペクトルの片方の段は空になります。

操作方法

  • スペクトルの棒をクリックすると、そこへプローブが移動します。
  • 「描く」を選び、上のグラフの1周期分をドラッグしてください。描いた点は線分でつながれ、左右に繰り返されます。

関連項目