Advances in Applied Mathematics
Vol.
08
No.
12
(
2019
), Article ID:
33528
,
5
pages
10.12677/AAM.2019.812235
Numerical Algorithms for Ordinary Differential Equations and Application
Sheng Huang, Zihui Jiang, Jiahui Li, Xiang Chi
College of Science, Shenyang Aerospace University, Shenyang Liaoning
Received: Nov. 23rd, 2019; accepted: Dec. 12th, 2019; published: Dec. 19th, 2019
ABSTRACT
In this paper, we mainly use MATLAB software to program the numerical solution of ordinary differential equation, such as Euler method and classical fourth-order Runge Kutta algorithm. Through the results of corresponding numerical example, we analyze the advantages and disadvantages of these two algorithms. Finally, we apply those algorithms to a concrete example; we solve it and compare the results.
Keywords:Ordinary Differential Equations, Numerical Algorithm
常微分方程常用的数值算法及其应用
黄 盛，姜紫惠，李佳辉，迟祥
沈阳航空航天大学理学院，辽宁 沈阳
收稿日期：2019年11月23日；录用日期：2019年12月12日；发布日期：2019年12月19日
摘 要
本文主要利用了MATLAB软件对常用的常微分方程数值解法编程实现，如欧拉法和经典的四阶龙格库塔算法，通过相应的数值算例的结果，分析了两种算法的优缺点。最后我们把算法应用到了一个具体的实例，分别用欧拉法和经典的四阶龙格库塔法进行求解并对计算结果做了比较。
关键词 :常微分方程，数值算法
Copyright © 2019 by author(s) and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
1. 常用的数值算法
常微分方程的理论和方法不仅在自然科学领域中应用非常广泛，而且在社会科学的各个领域有着越来越多的研究。但是很多时候，我们都没办法求出常微分方程的解析解，因而需要用数值算法去求出相应的数值解。常微分方程一般所用的求解初值问题的数值解法简单来说分为两类：单步法和多步法 [1] - [6]。本文主要讨论了单步法中的两类方法。
考虑如下的初值问题：
数值算法一般是：取步长
，将
分成n等份，得到离散点
，进而求解
在这些离散点处的近似值
。
1.1. 算法介绍
欧拉法的迭代格式为：
。在已知初值
的情况下，我们便可以依次推算出：
经典的四阶龙格库塔算法的迭代格式为：
1.2. 数值算例
我们分别用欧拉法和四阶的龙格库塔算法，利用MATLAB软件编程来计算如下的初值问题：
得到结果如表1所示：
Table 1. Comparison of calculation results
表1. 计算结果对比表
从表1可以看出，四阶龙格库塔法在的精确度是要高于欧拉算法的，但是欧拉算法较四阶的龙格库塔算法来说，要更简单。
2. 算法的实际应用
在文献 [7] 中，研究了三维的微生物连续发酵生产1,3丙二醇的微分方程动力系统如下：
其中：
，
，
。
参数：
，
，
，
，
，
，
，
，
，
，
，
，
，
(具体的参数意义可以参见相关文献)。
我们用欧拉法来求解该系统，建立迭代格式如下：
参数为：
，
，
，
，
，步长h取0.01，计算结果如表2所示：
Table 2. Calculation results of Euler method
表2. 欧拉法计算结果
我们用四阶的龙格库塔算法来求解该系统，建立迭代格式如下：
核心代码如下：
计算结果如表3所示：
Table 3. Calculation results of fourth-order classical Runge Kutta algorithm
表3. 四阶经典的龙格库塔算法的计算结果
从两组计算结果可以看出，欧拉法和经典的四阶龙格库塔算法都能够对该系统进行求解，两组数据相差不大。因为系统本身和实验数据都是带有误差的，而且该系统没有解析解，所以没有对两组数据的精确度作对比。
3. 结论
本文主要研究了利用MATLAB软件对欧拉法和经典的四阶龙格库塔算法进行了编程并上机实现，通过数值算例的结果分析了两种数值算法的特点。最后，分别用欧拉算法和经典的龙格库塔算法求解了一个三阶的微分方程动力系统，从计算结果来看，两种方法都适合对该系统的求解。常微分方程的数值算法还有很多，以后会再讨论。
基金项目
沈阳航空航天大学大学生创新创业训练项目，项目编号：201810143271。
文章引用
黄 盛,姜紫惠,李佳辉,迟 祥. 常微分方程常用的数值算法及其应用
Numerical Algorithms for Ordinary Differential Equations and Application[J]. 应用数学进展, 2019, 08(12): 2045-2049. https://doi.org/10.12677/AAM.2019.812235
参考文献
- 1. 钟巍. 一种新的常微分方程初值问题数值解法[J]. 重庆文理学院学报(自然科学版), 2011(4): 15-18.
- 2. 林爽, 张杰. 常微分方程初值问题的基本数值解法分析[J]. 海南师范大学学报(自然科学版), 2012, 25(2): 119-121.
- 3. 霍晓程. 常微分方程数值解法的研究[J]. 佳木斯教育学院报, 2011(5): 167-167.
- 4. 孙美玲. 常微分方程数值解法的Matlab计算与可视比较[J]. 高教学刊, 2016(19): 60-61.
- 5. 余丽. 常微分方程数值解法的若干求解格式[J]. 鄂州大学学报, 2013(S1): 123-125.
- 6. 孙燕新, 李娟. 关于几种常用一阶常微分方程数值解法的导出[J]. 常州工程学院报, 2010(6): 71-74.
- 7. 李艳杰. 非线性时滞脉冲动力系统的辨识模型, 优化算法及应用[D]: [硕士学位论文]. 大连: 大连理工大学, 2005, 12.