学业论文龙格库塔法求解微分方程组.doc

龙格—库塔求解微分方程组 数学111班：何荣航 指导老师：于鹏 （陕西科技大学理学院 陕西 西安 710021） 摘 要：在工程实践当中，经常会遇到求解微分方程组的问题，一般情形下，求微分方程组的精确解是困难的，为了解决此问题需要求微分方程组的数值解。龙格—库塔法就是求解微分方程组的一种有效方法。一个龙格—库塔法可由一个合适的泰勒方法推导而来，省去高阶导数计算，大大简化了微分方程组的求解。本文就是以N=4的龙格—库塔方法为基础，求解了一个具体的微分方程组。 关键字：龙格—库塔法,微分方程组,欧拉方法 Runge - Kutta for Solving Differential Equations Abstract：In engineering practice, often encounter the problem of solving differentialequations, under normal circumstances, the exact solution of the differentialequations is difficult, in order to solve this problem requires the numericalsolution of differential equations. An effective method of Runge Kutta method is the solution of differential equations. Can be made of a suitable method ofTaylor derived a Runge Kutta method, calculation of high order derivative is omitted, greatly simplifies the solution of differential equations. This paper is based on Runge Kutta method based N=4, solving a specific differentialequations. Keywords：Runge Kutta method ; Differential equations Formula; Euler method 龙格—库塔法解微分方程组的原理： 在求解微分方程组时我们可以考虑微分 ， （1） 可得出微分方程组 , 在区间上的数值解。 易得求解该方程组的欧拉方法，将差分式，和代入式，得 （2） 将区间分成M个子区间，宽度为,网格点为。把它们代入式中可得欧拉方法的递归公式：. （3） 要得到合理的精确度，应该使用更高阶的方法来缩小误差。而我们通常所说的龙格-库塔法是指四阶而言的也是在工程中应用广泛的经典龙格－库塔算法 解：第一个点为，计算和所需的中间计算为 用它们得到最终的计算结果： 数值解在每一步都有一定的误差。上面的例子中误差递增，并在处达到最大值： 参考文献： 1 马修斯.数值方法（MATLAB版）第四版.北京：电子工业出版社，2010 2 李庆扬.数值计算原理.北京：清华大学出版社，2000 3 李庆扬.数值分析第五版.北京：清华大学出版社，2008 1