ホイン: O(h²) 2次
RK4: O(h⁴) 4次
刻み幅を½にすると
RK4の誤差は1/16倍
オイラー・ホイン・RK4法を同一問題で比較。刻み幅による精度・安定性の違いを可視化。ロジスティック成長・ファン・デル・ポール振動子・ローレンツ方程式のプリセット付き。
機械・自動車工学:エンジンのピストン運動やサスペンションの振動解析に利用されます。リアルタイム制御が必要なECU(エンジン制御ユニット)内の計算では計算コストの低いオイラー法が、設計段階での詳細な性能評価にはRK4やRK45がよく用いられます。
電気・電子回路シミュレーション:トランジスタやオペアンプを含む非線形回路の過渡応答を計算するために必須です。スイッチングの瞬間など電圧・電流が急激に変化する部分では、RK45のような刻み幅制御機能を持つ手法が威力を発揮します。
化学反応・生物増殖モデリング:複数の化学物質が関与する反応速度論や、人口増加のロジスティック方程式の予測に使われます。反応速度が速い物質と遅い物質が混在する「硬い方程式」では、解法の選択が計算の成否を分けます。
気象・流体シミュレーション(CFD):ナビエ-ストークス方程式を時間方向に離散化する際の時間積分スキームとして応用されます。乱流の計算などでは膨大な計算コストがかかるため、精度と効率を両立させる高度な解法が研究されています。
まず、「高次=常に正義」ではないという点を押さえよう。確かにRK4はオイラー法より圧倒的に精度が高い。しかし、例えば制御システムのリアルタイムシミュレーションのように、計算速度が最優先で、多少の誤差がフィードバックループで吸収される場合は、オイラー法が採用されることもある。このツールで「ロジスティック成長」を選び、刻み幅h=1.0で各手法を比較してみてほしい。RK4とオイラー法の結果は大きく違うが、h=0.1まで細かくすると、オイラー法でも実用に足る解に近づく。つまり、「刻み幅」と「次数」はトレードオフの関係にあるんだ。
次に、「安定性」を見落とさないでほしい。例えば「ファン・デル・ポール振動子」で、刻み幅をh=2.0など非常に大きく設定してオイラー法を実行すると、解が発散してグラフが吹き飛ぶことがある。これは数値的不安定性と呼ばれる現象で、手法によって許容できる刻み幅の上限が違う。実務では、計算が破綻しない安定した刻み幅を事前に見積もることが必須だ。
最後に、「刻み幅さえ小さくすれば安心」ではないこと。本ツールは固定刻み幅のEuler・Heun・RK4を扱う。hを小さくすれば誤差は減るが計算ステップ数は反比例で増える。実務ではRK45(Dormand-Prince)のような適応刻み幅法で、許容誤差に応じてhを自動調整するのが定石だが、それでも「どの程度の精度が必要か」という要求仕様から目標精度を決める必要がある。本ツールでは刻み幅スライダーを動かし、各手法が真の解に収束する様子を比較して感覚を養おう。
ローレンツ方程式(既定 σ=10, ρ=28, β=2.67, 初期値[1,0,0], T=20)で固定刻み幅 h=0.01 を選ぶと、Euler・Heun・RK4 はいずれも約2000ステップを実行します。Euler法は刻み幅が粗いとカオス軌道がすぐ真の解から離れますが、RK4は同じhでも軌道をはるかに長く追従します。T=50・h=0.005 では各手法とも約10000ステップとなり、ステップ数が4倍の精度(O(h⁴))に効くRK4の優位が解の比較グラフで確認できます。