フーリエ級数展開シミュレーター
矩形波・三角波・のこぎり波・半波整流のフーリエ級数展開を可視化する無料ツール。項数を1〜50まで動かし、部分和が元の波形へ収束していく様子、高調波のスペクトル、ギブス現象による行き過ぎを同時に確認できます。
フーリエ級数展開とは
フーリエ級数展開は、周期をもつ関数を正弦波と余弦波の和に書き直す方法です。このツールでは矩形波・三角波・のこぎり波・半波整流の4つについて、教科書に載っている閉じた形の係数をそのまま使い、第何次までを足すかをスライダーで変えながら、部分和が元の波形へ近づいていく様子をグラフで確かめられます。数値積分は使っていないため、表示される係数は倍精度浮動小数点数の丸め誤差の範囲で厳密な値です。
このツールの見どころは、波形によって収束の速さがまったく違うところです。矩形波とのこぎり波は係数が 1/n で減衰するため、項をかなり増やしても角の近くの波打ちが消えません。一方で三角波は 1/n² で減衰するので、10項ほどでほとんど元の形に重なります。不連続点をもつ波形では、項をいくら増やしても跳びの手前で約9%行き過ぎる現象が残る点も、部分和の最大値の欄で数値として追えます。
フーリエ級数展開シミュレーターの使い方
- 波形を選ぶ 矩形波・三角波・のこぎり波・半波整流の4つから選びます。選ぶとすぐ下に、その波形の展開式と収束の特徴が表示されます。
- 項数のスライダーを動かす 項数 N を 1 から 50 まで動かせます。グラフのオレンジ色の実線が部分和、灰色の点線が元の波形です。
- 収束の指標を読む 二乗平均平方根誤差と誤差の最大値が、項を増やしたときにどれだけ小さくなるかを数値で追えます。部分和の最大値は行き過ぎの目安です。
- 係数の表とスペクトルを見る どの次数がどれだけの振幅を持つかを、棒グラフと表の両方で確認できます。矩形波なら偶数次が空になっていることが一目で分かります。
使いこなすためのヒント
- 矩形波で N を 5 → 15 → 45 と増やすと、波打ちの回数は増えるのに山の高さがほとんど変わらないことが見て取れます。これがギブス現象で、行き過ぎの量は項数によらず約9%に留まります。
- 三角波と矩形波を同じ項数で見比べると、係数の減衰が 1/n² と 1/n で違うことの効果がはっきり出ます。誤差の欄の桁が1つ以上違っているはずです。
- 半波整流だけは定数項（直流成分）を持ちます。整流回路の出力に直流分が残る理由を、展開式の 1/π という項として確認できます。
- 自分で決めた関数の係数を数値で知りたい場合は、数値積分で a₀・aₙ・bₙ を求めるフーリエ係数計算ツールをお使いください。
フーリエ級数展開シミュレーターの活用シーン
大学の講義や電気回路の教科書で習ったとき
閉じた形の式だけを見ても収束の様子は想像しづらいものです。項数を動かして波形が育っていく過程を見ることで、式の形と結果が結び付きます。
高調波歪みを扱うとき
矩形波が奇数次の高調波だけを含み、のこぎり波が偶数次も含むことは、フィルター設計やノイズ対策の議論の前提になります。スペクトルの棒グラフでその違いを一目で確認できます。
信号処理の入門でリンギングを学ぶとき
帯域を制限すると波形の角が波打つ理由は、この部分和の振る舞いそのものです。項数を制限することが帯域を制限することに対応しています。
楽器の音色を説明するとき
同じ基本周波数でも高調波の含まれ方が違えば音色は変わります。矩形波と三角波のスペクトルの差は、その最も分かりやすい例です。
フーリエ級数展開の用語
- フーリエ級数
- 周期関数を、元の周期の整数分の1の周期をもつ正弦波と余弦波の和として表したものです。理論上は無限個の項が並びます。
- 部分和
- 無限に続く級数を第 N 項で打ち切った有限個の和です。このツールのオレンジ色の実線が部分和にあたります。
- 基本波
- 級数のうち n = 1 の成分で、元の関数と同じ周期をもつ成分です。音の高さを決めるのはこの成分です。
- 高調波
- 基本波の整数倍の周波数をもつ成分です。n 次高調波は基本波の n 倍の周波数をもち、音色や波形の角の鋭さを決めます。
- 直流成分
- 周波数が 0 の成分で、1周期にわたる平均値にあたります。このツールでは定数項 a₀/2 として表示しています。
- ギブス現象
- 不連続点をもつ関数のフーリエ級数で、跳びの手前で部分和が行き過ぎる現象です。行き過ぎの量は項数を増やしても約9%残ります。
- 二乗平均平方根誤差
- 元の波形と部分和の差を1周期にわたって二乗平均し、平方根をとった値です。全体としてどれだけ合っているかの目安になります。
フーリエ級数展開についてよくある質問
余談ですが ― 「角のある関数は三角関数で書けない」と言われた時代
フーリエ級数は、1807年にジョゼフ・フーリエが熱伝導の問題を解くために持ち込んだ考え方が出発点です。当時は「角のある関数まで滑らかな三角関数の和で書けるはずがない」と受け止められ、ラグランジュをはじめとする同時代の数学者から強い反対を受けました。実際、どんな関数でも各点で収束するわけではないことが後に明らかになり、収束の条件をどう述べるかという問題が19世紀の解析学を大きく動かすことになります。ディリクレの収束条件やリーマン積分の整備は、この問いから生まれた成果です。
ギブス現象という名前は物理学者ジョサイア・ウィラード・ギブズに由来しますが、実際に先に指摘したのはイギリスの数学者ヘンリー・ウィルブラハムで、1848年のことでした。ギブズが1898年に論じたことで広く知られるようになり、名前だけが残った形です。行き過ぎの量が約8.95%というきりの悪い数値になることは、正弦積分と呼ばれる関数の最大値から導かれます。この数はどの不連続点でも、どんな波形でも同じ値になります。
この行き過ぎは今も実務で顔を出します。デジタル音声や画像で高い周波数を切り落とすと、音の立ち上がりにプリエコーと呼ばれる前置きの響きが乗り、画像の輪郭にはリンギングと呼ばれる縞が現れます。JPEG の圧縮率を上げたときに境界に見えるにじみも、根っこは同じ現象です。フィルターの設計で窓関数を工夫するのは、この行き過ぎをどこまで許すかという取引を調整する作業でもあります。項数のスライダーを動かすと、その取引の中身をそのまま目で確かめられます。