ルンゲ・クッタ法は、コンピュータを使って常微分方程式を解くための基本的かつ非常に重要な数値計算手法の一つです。特に、精度の高い4次のルンゲ・クッタ法は、工学部の基礎教育はもちろん、企業の研究・開発現場でも活用されています。
本記事では、常微分方程式を数値的に解く手法としてのルンゲ・クッタ法の数値計算を中心に、基本的な理論、オイラー法との違い、Excelでの実装例、さらに応用問題としてロトカ・ヴォルテラ方程式の解析まで、順を追って丁寧に解説していきます。
A. ルンゲ・クッタ法 関連資料ダウンロード一覧
以下のリンクから、ルンゲ・クッタ法に関するExcelファイルや練習用資料などをダウンロード可能です。
| 番号 | 内容 | ダウンロードリンク |
|---|---|---|
| 1. | ルンゲ・クッタ法による核崩壊モデルのExcelファイル | ダウンロードはこちら |
| 2. | ロトカ・ヴォルテラ方程式（連立微分方程式）実装例 | ダウンロードはこちら |
| 3. | オイラー法との比較用Excelテンプレート | ダウンロードはこちら |
| 4. | 微分方程式解法マクロ付きVBAソースコード | ダウンロードはこちら |
| 5. | 解説付きグラフ描画Excelテンプレート | ダウンロードはこちら |
B. 数値計算の基本とルンゲ・クッタ法の概要
微分方程式を数値で解くとは？
微分方程式をコンピュータで解く場合、「連続関数」を離散化して四則演算に落とし込む必要があります。これが数値計算の基本です。関数 y = f(x) の微分を用いた近似値を段階的に求め、点の集合として関数を再構築します。
オイラー法の限界と問題点
最も基本的な数値解法であるオイラー法は、下記の形で次の値 yi+1 を求めます。
\[ y_{i+1} = y_i + h f(x_i, y_i) \]
- 簡単だが精度が低く、刻み幅 hh を小さくしても誤差の蓄積が避けられません。
- 急激な変化を持つ関数や曲率の大きいグラフでは、著しい誤差が生じます。
4次のルンゲ・クッタ法のアルゴリズム
オイラー法の進化形がルンゲ・クッタ法です。特に4次のルンゲ・クッタ法では、4つの異なる傾きを計算し、その加重平均により次の点を求めます。
手順は以下の通りです。
\[ \begin{aligned} k_1 &= h f(x_i, y_i) \\ k_2 &= h f\left(x_i + \frac{h}{2}, y_i + \frac{k_1}{2}\right) \\ k_3 &= h f\left(x_i + \frac{h}{2}, y_i + \frac{k_2}{2}\right) \\ k_4 &= h f(x_i + h, y_i + k_3) \end{aligned} \]
そして、
\[ y_{i+1} = y_i + \frac{1}{6}(k_1 + 2k_2 + 2k_3 + k_4) \]
この計算により、高精度な予測が可能になります。
Excelでの実装例（例題1：核崩壊）
以下の常微分方程式をルンゲ・クッタ法で解いてみましょう。
\[ \frac{dN}{dt} = -\lambda N \quad (\lambda = 0.08664, \, N(0) = 100) \]
目的：Nが50（半減）になるまでの時間を求める。
Excelの構成
- h = 0.1
- 各ステップで k1,k2,k3,k4 を算出
- yi+1 を反復計算
結果
MATCH関数により50を下回る時間を算出すると、およそ8日であることが分かります（理論値と一致）。
オイラー法との比較表
| 項目 | オイラー法 | ルンゲ・クッタ法（4次） |
|---|---|---|
| 精度 | 低い | 高い |
| 実装の容易さ | 非常に簡単 | やや複雑 |
| 計算負荷 | 軽い | 中程度 |
| 曲線への追従性 | 弱い | 強い |
| 実務での利用 | 限定的 | 多用される |
応用例（例題2：ロトカ・ヴォルテラ方程式）
- xx：ウサギの個体数（初期値 1000）
- yy：キツネの個体数（初期値 100）
- a=0.01,b=0.05,c=0.0001
- 時間：1000日分をシミュレーション
結果
Excelによるルンゲ・クッタ法での解析により、ウサギとキツネの個体数が時間と共に周期的に変動する様子が視覚的に確認できます。
まとめ：ルンゲ・クッタ法の数値計算は、微分方程式の解法として極めて有効
ルンゲ・クッタ法の数値計算は、精度が高く、現実の複雑な現象を数値的に再現できる信頼性の高い手法です。本記事では、基本的な原理とオイラー法との比較、Excelによる実装例、さらに連立方程式への応用までを丁寧に解説しました。
特に、単純な微分方程式のみならず、ロトカ・ヴォルテラ方程式のような実用的なモデルにおいても、ルンゲ・クッタ法は高い効果を発揮します。数値計算の基礎をしっかりと理解し、状況に応じて最適な刻み幅やアルゴリズムを選択することが、正確なシミュレーションには欠かせません。
数値解析のスキルは、工学や物理、生態学などさまざまな分野で活用可能です。今後もぜひ、ルンゲ・クッタ法の数値計算を活かして、実践的な課題解決に挑戦してみてください。
