フーリエ級数可視化ツール 戻る
信号処理・数学

フーリエ級数可視化ツール

矩形波・鋸波・三角波を回転する円(フェーザー)の重ね合わせで合成。高調波の追加につれて信号が近似される様子をリアルタイムアニメーション。ギブス現象も体験。

波形選択
パラメータ
高調波数 N
アニメーション速度
×
表示切替
計算結果
7
使用高調波数
f₁
基本周波数
最大振幅
近似誤差
フェーザーアニメーション + 合成波形
周波数スペクトル(高調波振幅)
理論・主要公式

矩形波:

$$x(t) = \frac{4}{\pi}\sum_{n=1,3,5,...}^{N}\frac{\sin(n\omega t)}{n}$$

奇数次高調波のみ。
N→∞ で完全な矩形波に収束。

フーリエ級数可視化ツールとは

🙋
このシミュレーターで、複雑な波が単純なサイン波の足し算で作れるって本当ですか?
🎓
本当だよ。大まかに言うと、どんな周期的な形の波も、基本となるサイン波と、その整数倍の速さで振動する「高調波」というサイン波を足し合わせることで作り出せるんだ。このツールの「高調波数 N」のスライダーを動かしてみて。N=1から順に増やしていくと、三角波や矩形波がサイン波の積み木のように組み上がっていく様子がわかるよ。
🙋
画面でくるくる回っている矢印(フェーザー)は何を表しているんですか?
🎓
あれは、それぞれのサイン波を複素平面上の回転ベクトルとして表したものだ。各矢印の長さが振幅、回る速さが周波数に対応している。この矢印たちを、先端から次へとつなげていった、最後の矢印の先端の高さ(y座標)が、今見ている合成波形の瞬間の値になるんだ。「アニメーション速度」を変えると、この回転の様子がゆっくり観察できるぞ。
🙋
矩形波でNを大きくしても、角のところがピョコンと飛び出してしまうのはなぜですか?
🎓
あれが有名な「ギブス現象」だ。不連続な点(矩形波の角)を滑らかなサイン波だけで無理やり近似しようとすると、どうしても約9%のオーバーシュート(飛び出し)が起きてしまうんだ。Nをどこまで大きくしてもこの飛び出しは消えない、数学的に面白い現象さ。シミュレーターでNを最大まで上げて、その様子を確かめてみよう。

よくある質問

それはギブス現象です。フーリエ級数で不連続点を近似する際、高調波をいくら増やしても約9%のオーバーシュートが残る現象です。このツールで高調波数を変えながら、その振動が収束せずに尖ったまま残る様子を確認できます。
矩形波は奇数次正弦波のみ(1,3,5…次)、鋸波は全次数の正弦波(1,2,3…次)、三角波は奇数次正弦波で係数が1/n²と急速に減衰します。この違いにより、三角波は少ない高調波でも滑らかに近似されるのが特徴です。
各円の回転が各高調波のcos/sin成分に対応します。円の半径がフーリエ係数、回転速度が周波数nωです。すべての円の中心の軌跡を時間方向に投影すると、それらの和として元の周期波形が再現される仕組みです。
高調波の打ち切り次数と波形の近似精度の関係を視覚的に理解できます。例えば、矩形波は多くの高調波が必要で帯域制限時に歪みやすい一方、三角波は少ない次数でも実用可能です。音色設計やフィルタ設計の直感的な参考になります。

実世界での応用

音声合成・電子音楽:三角波、矩形波、鋸波はシンセサイザーの基本波形です。フーリエ級数の原理を用いて、これらの倍音(高調波)成分を制御することで、豊かな音色を作り出しています。

振動・騒音解析:エンジンやモーターの回転から生じる周期的な振動や騒音を測定し、フーリエ級数(フーリエ変換)で周波数成分に分解します。どの周波数の振動が問題かを特定し、対策を講じるために使われます。

交流回路・電力工学:実際の交流電源の波形は完全な正弦波ではなく歪みを含みます。この歪み波形をフーリエ級数で分解し、高調波成分が機器に与える影響を評価・対策します(高調波抑制)。

画像処理(2次元拡張):フーリエ級数の考え方は2次元に拡張され、JPEGなどの画像圧縮技術の基礎となっています。画像の濃淡パターンを異なる周波数と方向の波の重ね合わせとして表現します。

よくある誤解と注意点

このツールを使い始めるとき、いくつか勘違いしやすいポイントがあるから気をつけてね。まず、「高調波数Nを無限にすれば完全な波形になる」と思いがちだけど、矩形波のように急峻な不連続点を持つ波形では、ギブス現象で約9%のオーバーシュートが残るんだ。これは数学的に消えないもので、N=100にしても1000にしても変わらない。実務では、この「完全には再現できない」という特性を頭に入れて設計する必要があるよ。

次に、「基本波の周波数」と「波形の周期」の関係を混同しないで。例えば、基本周波数 $f_1$ が 50Hz なら、その周期 $T_1$ は $1/50 = 0.02$秒だよね。でも、合成される矩形波の周期も同じ0.02秒になる。基本波1周期分の時間で、合成波形も1周期を完了するんだ。ツールで「基本周波数」スライダーを動かすと、フェーザーの回転速度と波形の繰り返し速度が連動して変わることを確認してみよう。

最後に、実践的な落とし穴として「打ち切り誤差」がある。無限級数を有限のNで打ち切るから、どうしても誤差が生じる。例えば、N=9の矩形波は、理論式 $$x_{\text{square}}(t) = \frac{4}{\pi}\sum_{n=1,3,5,...}^{\infty}\frac{\sin(n\omega t)}{n}$$ のn=17以降を無視している状態だ。この誤差は、滑らかな三角波では小さく早く収束するけど、矩形波では振動(リンギング)として残りやすい。実務の信号処理では、必要な精度と計算コストのトレードオフを考えてNを決めるんだ。

使い方ガイド

  1. 波形タイプ(矩形波・鋸波・三角波)を選択し、基本周波数を1~50Hzの範囲で設定します
  2. スライダーで使用高調波数を1~50次まで調整し、フェーザーの重ね合わせアニメーションを確認します
  3. 再生速度(0.5~2.0倍)でアニメーション表示を制御し、各高調波成分の位相関係を観察します
  4. 左側のグラフで時間領域波形、右側でスペクトラムを同時表示して収束過程を検証します

具体的な計算例

矩形波(振幅1V、基本周波数10Hz)を近似する場合、1次のみでは実効値0.707V、9次まで加算すると0.998Vに収束します。鋸波の場合、N=20次で近似誤差は2.1%まで低下し、各高調波の振幅は1/Nに比例して減少します。三角波では奇数次高調波のみが存在(5次:0.04V、7次:0.020V)し、12次までで99.5%の波形再現精度が得られます。

実務での注意点

  1. 電源の高調波歪み測定(矩形波近似)では、電流計の帯域幅がサンプリング周波数の3倍以上必要です。50Hz基本波から250Hz(5次)まで対応する計器を選定してください
  2. ギブス現象により不連続点で±9%のオーバーシュートが発生するため、PWM制御などの急峻な遷移を伴う装置では充分な次数(N≧20)を確保してください
  3. アナログフィルタ設計時は、ナイキスト周波数(最大周波数の2倍)を超える高調波成分を除去し、時間領域でのリンギング現象を防止します