Runge Kutta方法求解常微分方程

最新推荐文章于 2024-06-15 16:02:05 发布

Kevin_Opt

最新推荐文章于 2024-06-15 16:02:05 发布

阅读量1.2k

点赞数 1

分类专栏：算法文章标签：算法

本文链接：https://blog.csdn.net/qq_41523340/article/details/124887661

版权

算法专栏收录该内容

6 篇文章 0 订阅

订阅专栏

本文详细介绍了如何使用二阶和四阶Runge-Kutta方法求解常微分方程，通过实例演示了改进的Euler方法和标准四阶格式，并展示了它们在初值问题上的应用。对比分析了两种方法的精度，强调了四阶方法在提高数值解的准确性方面的优势。

摘要由CSDN通过智能技术生成

参考书籍：应用数学基础下

常微分方程

一阶常微分方程形式如下：
在这里插入图片描述

二阶Runge Kutta法

改进的Euler格式如下：
在这里插入图片描述
例题：用改进的Euler格式求解初值问题：

解：根据上面公式可以列出改进的Euler格式为：

取步长为h=0.1，初值y0=1，数值结果和解析结果如下表所示：

x_n	y_n	y(x_n)
0.1	1.09590909	1.09544511
0.2	1.18409657	1.18321596
0.3	1.26620136	1.26491106
0.4	1.34336015	1.341640798
0.5	1.41640193	1.41421356
0.6	1.48595560	1.48323970
0.7	1.55251409	1.54919334
0.8	1.61647478	1.61245155
0.9	1.67816636	1.67332005
1.0	1.73786740	1.73205081

代码如下：

clear;clc;
h=0.1;
xn=0;
yn=1;
for n = 1:10
    K1 = yn - 2*xn/yn;
    K2 = yn + h*K1 - 2*(xn + h)/(yn + h*K1);
    yn_1 = yn + h/2*(K1 + K2);
    xn = xn + h;
    yn = yn_1;
    yn_1
end

四阶Runge Kutta法

标准四阶Runge Kutta格式为：
在这里插入图片描述
例题：用改进的Euler格式求解初值问题：

解：计算格式为：

取步长为h=0.1，初值y0=1，数值结果和解析结果如下表所示：

x_n	y_n	y(x_n)
0.1	1.09544553	1.09544511
0.2	1.18321674	1.18321596
0.3	1.26491222	1.26491106
0.4	1.34164235	1.341640798
0.5	1.41421557	1.41421356
0.6	1.48324222	1.48323970
0.7	1.54919645	1.54919334
0.8	1.61245534	1.61245155
0.9	1.67332465	1.67332005
1.0	1.73205636	1.73205081

可见，四阶Runge Kutta法相比二阶Runge Kutta精度更高。
代码如下：

clear;clc;
xn=0;
yn=1;
h=0.1;
for n=1:10
    K1 = yn-2*xn/yn;
    K2 = yn + h/2*K1 - (2*xn+h)/(yn+h/2*K1);
    K3 = yn + h/2*K2 - (2*xn+h)/(yn+h/2*K2);
    K4 = yn + h*K3 - 2*(xn + h)/(yn + h*K3);
    yn_1 = yn + h/6*(K1 + 2*K2 + 2*K3 + K4);
    xn = xn + h;
    yn = yn_1;
    yn_1
    pred = sqrt(1 + 2*xn)
end