フーリエ級数展開計算ツールは、周期関数または最大5つの区分からなる区分関数を入力し、フーリエ係数 a₀、aₙ、bₙ の計算、部分和の確認、そして元の関数とフーリエ近似を比較する再構成グラフの描画を同時に行えるオンラインツールです。
フーリエ級数の定義と計算規則
本ツールは、以下のフーリエ級数の定義に基づいて計算を行います。
f(x) ≈ a₀ / 2 + Σₙ₌₁^N [ aₙ cos(nω₀ x) + bₙ sin(nω₀ x) ]
ここで、周期は T = b - a であり、基本角周波数は ω₀ = 2π/T と定義されます。この定義に基づき、定数項および各係数の積分にはすべて 2/T という積分因子が適用されます。
- 定数項: a₀ = 2 / T∫ f(x)dx と計算され、級数における定数項は a₀/2 として代入されます。
- 高調波係数: aₙ および bₙ は、指定された高調波の数 N（1から30までの整数）の範囲で計算されます。
入力パラメータと設定方法
計算を開始するには、以下の項目を設定します。
1. 関数の定義
「1周期における関数」を直接入力するか、最大5つの「関数の区分」を定義します。「区分を追加」ボタンで区分を増やし、「区分 ‹n› を削除」ボタンで不要な区分を削除できます。各区分には「開始」と「終了」の境界値を指定します。
プリセットとして「のこぎり波」、「矩形波」、「三角波」、「整流正弦波」の例を選択することも可能です。
式 $f(x)$ の入力には、変数 x、定数 pi、e、演算子 +、−、*、/、^、括弧、および数学関数（sin, cos, tan, exp, ln, log, sqrt, abs, min, max）が使用できます。掛け算は 2*x のほか、2x のように省略して記述することも可能です。
2. 周期の設定
「始点 a」と「終点 b」を入力します。三角関数の引数はラジアンとして処理され、始点と終点は変数 x と同じ単位を使用します。
3. 級数と精度
- 高調波の数 N: 1から30までの整数を指定します。
- シンプソン分割数: 200から4000までの偶数を指定します。
- 表示桁数: 結果に表示する小数点以下の桁数を指定します。
- 指定した x で Sₙ(x) を評価（任意入力）: 周期と同じ x の単位を使用し、特定の点における部分和の値を算出します。
数値積分と安定性チェック
本ツールにおけるフーリエ係数の算出は、厳密な代数的解析ではなく、複合シンプソン公式を用いた数値積分によって行われます。
計算の信頼性を確保するため、ツールは指定された「シンプソン分割数」で積分を実行した後、自動的に2倍の解像度（分割数）で再度計算を行う安定性チェックを実行します。この2つの計算結果の間で生じた係数の最大変化量を δ として検出します。
計算実行後、以下のステータスが表示されます。
- 安定時: 「計算完了。解像度を2倍にした係数チェックは、表示された精度において安定しています。」
- 不安定時: 「計算されましたが、積分解像度を2倍にした際に係数が
‹delta›変化しました。分割数を増やすか、関数に急峻な変化がないか確認してください。」
不連続点とギブズの現象
矩形波やのこぎり波のように、関数の区分間に不連続な段差（ジャンプ）が存在する場合、再構成グラフの境界付近に特有の波打ち（振動）が発生します。これは数値積分のエラーではなく、数学的に避けられない「ギブズの現象」によるものです。
不連続点が存在する場合、ツールは「区分間に不連続点（ジャンプ）があります。境界付近の波打ちは、積分エラーではなくギブズの現象である可能性があります。」というメッセージを表示します。通常の収束条件のもとでは、不連続点におけるフーリエ級数の値は、その点における左極限値と右極限値の平均値に収束します。
エラーメッセージ一覧
入力内容や計算過程に問題がある場合、以下のエラーメッセージが表示されます。
| エラーメッセージ | 発生原因 |
|---|---|
| 「式の構文を確認してください。」 | 入力された数式に文法的な誤りがある場合。 |
| 「この周期内で関数が実数定義域から外れています。」 | 負の数の平方根や負の数の対数など、実数の範囲を超えた場合。 |
| 「この周期内で関数がゼロ除算を行っています。」 | 分母がゼロになる計算が発生した場合。 |
| 「関数または係数が、サポートされている数値範囲を超えています。」 | 計算結果がシステムの許容数値範囲をオーバーフローした場合。 |
| 「周期、高調波、および精度の範囲を確認してください。」 | N や分割数などが規定の範囲外に設定されている場合。 |
| 「関数の区分とその境界を確認してください。」 | 区分の境界値の設定に矛盾や重複がある場合。 |
プライバシーについて
入力された数式や計算結果などのデータはすべてユーザーが使用しているブラウザ内でのみ処理され、外部のサーバーに送信されることはありません。
よくある質問
係数は厳密な値ですか？
いいえ、これらは複合シンプソン公式による数値積分です。この電卓は、選択された分割数の2倍の解像度で計算を繰り返し、係数の最大変化量を報告します。変化量が小さいことは表示された桁数の妥当性を裏付けますが、正式な誤差範囲を示すものではありません。
この電卓はどのフーリエ係数の定義を使用していますか？
本ツールでは、ω₀ = 2π/T として f(x) ≈ a₀/2 + Σ [aₙ cos(nω₀ x) + bₙ sin(nω₀ x)] と定義しています。したがって、a₀、aₙ、bₙ はすべて 2/T の係数を使用します。書籍によっては平均値を直接 a₀ とする定義もあるため、係数を比較する前に表示されている数式を確認してください。
グラフが不連続点の近くで波打つのはなぜですか？
有限のフーリエ和は、不連続点の近くで振動します。高調波の数を増やすと影響を受ける領域は狭くなりますが、特徴的なギブズの現象（オーバーシュート）は解消されません。通常の収束条件のもとでは、不連続点において級数は左極限と右極限の平均値に収束します。