算法-----龙格－库塔法（转）

最新推荐文章于 2024-10-13 12:14:17 发布

知耻而后勇的蜗牛

最新推荐文章于 2024-10-13 12:14:17 发布

阅读量7.9k

点赞数

分类专栏：算法文章标签：算法框架工具 bi hp c

本文链接：https://blog.csdn.net/yf210yf/article/details/7638334

版权

算法专栏收录该内容

3 篇文章 2 订阅

订阅专栏

数值分析中，龙格－库塔法（Runge-Kutta）是用于模拟常微分方程的解的重要的一类隐式或显式迭代法。这些技术由数学家卡尔·龙格和马丁·威尔海姆·库塔于1900年左右发明。

龙格库塔法的家族中的一个成员如此常用，以至于经常被称为“RK4”或者就是“龙格库塔法”。该方法主要是在已知方程导数和初值信息，利用计算机仿真时应用，省去求解微分方程的复杂过程。经典四阶龙格库塔法

令初值问题表述如下。

$y' = f(t, y), \quad y(t_0) = y_0$

则，对于该问题的RK4由如下方程给出：

$y_{n+1} = y_n + {h \over 6} (k_1 + 2k_2 + 2k_3 + k_4)$

其中

$k_1 = f \left( t_n, y_n \right)$

$k_2 = f \left( t_n + {h \over 2}, y_n + {h \over 2} k_1 \right)$

$k_3 = f \left( t_n + {h \over 2}, y_n + {h \over 2} k_2 \right)$

$k_4 = f \left( t_n + h, y_n + hk_3 \right)$

这样，下一个值(yn+1)由现在的值(yn)加上时间间隔(h)和一个估算的斜率的乘积决定。该斜率是以下斜率的加权平均：
k1是时间段开始时的斜率；
k2是时间段中点的斜率，通过欧拉法采用斜率k1来决定y在点tn + h/2的值；
k3也是中点的斜率，但是这次采用斜率k2决定y值；
k4是时间段终点的斜率，其y值用k3决定。
当四个斜率取平均时，中点的斜率有更大的权值：

$\mbox{slope} = \frac{k_1 + 2k_2 + 2k_3 + k_4}{6}.$