常微分方程的Euler解法.doc

. PAGE . . . 毕 业 论 文 题 目： 常微分方程的Euler解法 及其程序设计 学 院： 数学与信息科学学院 专 业： 数学与应用数学 毕业年限： 2011年6月 学生姓名： 学 号： 指导教师： 常微分方程的Euler解法及其程序设计 摘要 本文总结了常微分方程的Euler解法，对各种格式给出了误差估计，设计了这些格式的计算程序． 关键词 常微分方程；Euler解法；误差分析；程序设计 Euler Method of Ordinary Differential Equation and Its Programming Abstract Euler method of ordinary differential equation is summarized，the error of each format is analyzed and its programming is designed in this paper. Keywords Ordinary differential equation; Euler method; Error analysis; Programming 科学技术中常常需要求解常微分方程的定解问题，这类问题最简单的形式，即为微分方程 （1） 的初值问题 （2） 定理 （存在与唯一性定理）如果方程（1）的右端函数在闭矩形域 上满足如下条件： （1）在上连续； （2）在上关于变量满足利普希茨（Lipschitz）条件，即存在常数，使 对于上任何一对点和有不等式： , 则初值问题（2）在区间上存在唯一解 , 其中． 根据存在与唯一性定理，只要关于满足Lipschitz条件 ， 即可保证其解存在并唯一． 然而解析方法只能用来求解少数较简单和典型的常微分方程，例如线性常系数微分方程等，对于变系数常微分方程的解析求解就比较困难，而一般的非线性常微分方程就更不用说，因此，在大多数情况下，实际问题中归结出来的微分方程主要靠数值解法求解． 所谓数值解法，就是寻找在一系列离散节点上的近似值．相邻两个节点的间距称为步长，假定为定数，节点为. Euler解法 （1）Euler格式 Euler格式的计算公式为 . （3） 下面用4种方法推导公式（3） a 泰勒展开法 在处展开有 , （4） 略去余项，得 . 用近似值代替，把上式右端所得结果记为，得 . b 数值差商 由导数定义知 ， 所以 . 用代替，把上式右端所得结果记为，即得公式（3）. c 数值积分 在区间上对微分方程（1）进行积分，得 , （5） 利用左矩形公式得 . 用代替，把上式右端所得结果记为，即得公式（3）. d 几何方法 在平面中，微分方程 { 的解称为它的积分曲线．积分曲线上的一点的切线斜率等于函数的值．如果按函数在平面上建立一个方向场，那么，积分曲线上的每一点的切线方向均与方向场在该点的方向一致． 基于上述对微分方程解的几何解释，从初始点出发，依方向场在该 点的方向上推进到上一点，然后再从依方向场的方向推进到上 的一点．循环前进，可作出一条折线（如图1） 图1 一般地，设已作出该折线的极点，过依方向场的方向再推进到，显然，两个极点，的坐标有下列关系： , . 若初值已知，则依上式可逐步算出 ， ， . 即得Euler格式的计算公式（3）. （2）梯形格式 如果用梯形公式计算（5）式右端的积分，得 . 对于上式右端，用代替，把上式右端所得结果记为，即得梯形公式 . （6） 梯形公式（6）是关于的隐式方程，因此称梯形法为隐式方法，而公式（3）是的显式形式，因此称欧拉法为显式方法． （3）改进的Euler格式 先用Euler格式求得一个初步的近似值，称之为预测值，预测值的精度可能很差，再用梯形公式将它校正一次，得，这个结果称之为校正值．而这样建立的预测—校正