常微分方程及几种数值解法.docVIP

常微分方程的几种数值解法 数学与应用数学 肖振华 指导教师 张秀艳 【摘要】自然界与工程技术中的很多现象，可以归结为微分方程定解问题。其中，常微分方程求解是微分方程的重要基础内容。但是，对于许多的微分方程，往往很难得到甚至不存在精确的解析表达式，这时候，数值解提供了一个很好的解决思路。，针对于此，本文对常微分方程数值解法进行了简单研究，主要讨论了一些常用的数值解法，如欧拉法、改进的欧拉法、Runge—Kutta方法、Adams预估校正法以及勒让德谱方法等，通过具体的算例，结合MATLAB求解画图，初步给出了一般常微分方程数值解法的求解过程。同时，通过对各种方法的误差分析，让大家对各种方法的特点和适用范围有一个直观的感受。 【关键词】 常微分方程 数值解法 MATLAB 误差分析 【Abstract】 Many phenomena in nature and engineering can be attributed to the definite solution of the problem for differential equations. Among them, the ordinary differential equation solving is an important foundation for the content of the differential equations. However, many of the differential equations are often difficult to obtain accurate analytical expression .At this time, the numerical solution provides a good idea. For the Numerical Solution of Ordinary Differential Equations in this article, we focuses on some commonly used numerical solution, such as the Euler method, improved Euler method, Runge-Kutta method, Adams predictor corrector method as well as newer spectral methods. Through specific examples, combined with MATLAB solving and drawing, we initially know the solution process of general numerical solution of ordinary differential equations . At the same time, according to the error analysis of various methods , everyone has an intuitive feel of the characteristics and scope of the various methods. 【Keywords】Ordinary Differential Equations Numerical Solution MATLAB error analysis 目录 1 前言 3 1.1 常微分方程概述 3 1.2 常微分方程解与数值解法 3 2 欧拉法和改进的欧拉法 4 2.1 欧拉法 4 2.2 改进的欧拉方法 5 2.3 算例 5 3龙格-库塔法 13 3.1 龙格-库塔法与泰勒展开 13 3.2 龙格-库塔法公式与ode函数 14 3.3算例 17 4 阿达姆斯预估校正法 21 4.1 阿达姆斯（Adams）公式 21 4.2 预估校正方法 23 4.3算例 23 5 勒让德谱方法 28 5.1 谱方法介绍 28 5.2勒让德多项式与谱方法 28 5.3 算例 30 参考文献 39 1 前言 1.1 常微分方程概述 方程是一个在数学中非常熟悉的名词，在初等数学里，我们将我们要研究的问题作为一个或几个未知量，通过观察事物的规律，得出这些未知量与已知量之间的等式关系，这样就得到了一个简单的方程或方程组——当然，这只是一个很浅显粗略的定义。在数学上，物质的运动和变化规律是通过函数关系来表示的，在一些复杂的现象中，比如火箭的运动轨迹，我们要求的未知量就变成了满足特定条件的一个或一些未知函数（不再是简单的数值）。有时候，我们需要用到导数很微分的知识，即这些未知函数的导数与自变量满足某