浅夏110 发表于 2020-6-9 14:53

常微分方程的解法 (三): 龙格—库塔(Runge—Kutta)方法 、线性多步法

§ 4 龙格—库塔(Runge—Kutta)方法
1 龙格—库塔方法的基本思想

https://img-blog.csdnimg.cn/20190430094433448.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzI5ODMxMTYz,size_16,color_FFFFFF,t_70

确定系数以提高精度

https://img-blog.csdnimg.cn/20190430094536586.png

https://img-blog.csdnimg.cn/20190430094609765.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzI5ODMxMTYz,size_16,color_FFFFFF,t_70


2 .  RK方法 :  4 阶龙格—库塔公式
要进一步提高精度,必须取更多的点,如取 4 点构造如下形式的公式:

https://img-blog.csdnimg.cn/20190430094729325.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzI5ODMxMTYz,size_16,color_FFFFFF,t_70
https://img-blog.csdnimg.cn/20190430094753706.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzI5ODMxMTYz,size_16,color_FFFFFF,t_70


这就是常用的 4 阶龙格—库塔方法(简称 RK 方法).



5  线性多步法
多步法的基本思想 、增量函数


https://img-blog.csdnimg.cn/20190430094910934.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzI5ODMxMTYz,size_16,color_FFFFFF,t_70

https://img-blog.csdnimg.cn/20190430095010460.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzI5ODMxMTYz,size_16,color_FFFFFF,t_70

https://img-blog.csdnimg.cn/20190430095044300.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzI5ODMxMTYz,size_16,color_FFFFFF,t_70
§6 一阶微分方程组与高阶微分方程的数值解法
6.1 一阶微分方程组的数值解法

https://img-blog.csdnimg.cn/20190430095312490.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzI5ODMxMTYz,size_16,color_FFFFFF,t_70

么问题(25)在 上存在唯一解 y = y(x) 。 问题(25)与(1)形式上完全相同,故对初值问题(1)所建立的各种数值解法可 全部用于求解问题(25)。

6.2 高阶微分方程的数值解法
高阶微分方程的初值问题可以通过变量代换化为一阶微分方程组初值问题。

https://img-blog.csdnimg.cn/20190430095354445.png

https://img-blog.csdnimg.cn/20190430095434426.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzI5ODMxMTYz,size_16,color_FFFFFF,t_70

刚性方程组、Stiff 方程组
最后需要指出的是,在化学工程及自动控制等领域中,所涉及的常微分方程组初值 问题常常是所谓的“刚性”问题。具体地说,对一阶线性微分方程组

https://img-blog.csdnimg.cn/2019043009553775.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzI5ODMxMTYz,size_16,color_FFFFFF,t_70
————————————————
版权声明:本文为CSDN博主「wamg潇潇」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。
原文链接:https://blog.csdn.net/qq_29831163/article/details/89703598


页: [1]
查看完整版本: 常微分方程的解法 (三): 龙格—库塔(Runge—Kutta)方法 、线性多步法