常微分方程的差分方法.docx

第五章 常微分方程的差分方法 一、教学目标及基本要求 通过对本节课的学习，使学生掌握常微分方程、常微分方程方程组的数值解法。 二、教学内容及学时分配 本节课主要介绍常微分方程的数值解法。具体内容如下： 讲授内容：欧拉公式、改进的欧拉公式。 三、教学重点难点 教学重点：改进的欧拉公式、龙格库塔方法、收敛性与稳定性。 教学难点：收敛性与稳定性。 四、教学中应注意的问题 多媒体课堂教学为主。适当提问，加深学生对概念的理解。 五、正文 基于数值积分的求解公式：欧拉公式、改进的欧拉公式 引 言 主要考虑如下的一阶常微分方程初值问题的求解: ? y?(x) ? f (x, y) ? ?y(x ? 0 ) ? y 0 微分方程的解就是求一个函数 y=y(x)，该函数满足微分方程并且符合初值条件。 例如微分方程: xy-2y=4x ；初始条件: y(1)=-3。 于是可得一阶常微分方程的初始问题 ?? y? ? ? ? 2 y ? 4 x ?? y(1)? ?3 。 显然函数 y(x)=x2-4x 满足以上条件,因而是该初始问题的微分方程的解。 但是,只有一些特殊类型的微分方程问题能够得到用解析表达式表示的函数 解，而大量的微分方程问题很难得到其解析解，有的甚至无法用解析表达式来表示。因此，只能依赖于数值方法去获得微分方程的数值解。 微分方程的数值解： 设微分方程问题的解 y(x)的存在区间是[a,b]，初始点 x0=a，将[a,b]进行划分得一系列节点 x0 , x1 ,...,xn，其中 a= x0 x1… xn =b。y(x)的解析表达式不容易得到或根本无法得到，我们用数值方法求得 y(x)在每个节点 xk 的近似值 y(xk)，即y≈y(xk)，这样 y0 , y1 ,...,yn 称为微分方程的数值解。 如果计算 yn 时，只利用 yn-1，称这种方法为单步法；如果在计算 yn 时不仅利用 yn-1，而且还要利用 yn-2, yn-3,…, yn-r,则称这种方法为 r 步方法，也称多步法。 §5.1 欧拉方法 §5.1.1 欧拉格式 ? y(x ) ? y(x ) 方程 y?(x n ) ? f (x , y n n ) 中， y (x ) ? n n?1 n h y(x ) ? y(x ) ? hf (x , y(x )) ? y ? y ? hf (x , y ) n?1 n n n n?1 n n n 称为解一阶常微分方程初值问题的欧拉公式，也称显示欧拉公式。 p2y1 p 2 y 1 p 1 y 0 x 0 x 1 x 2 线。用欧拉公式求数值解的几何意义如图： 容易验证，该折线各个顶点的纵坐标 y n 所以，欧拉方法又称为折线法。 算例：P98 (n ? 1,2...) 就是欧拉公式算得的近似值解， 可以看出误差随着计算在积累。 Euler 法的特点和误差 ?特点：(1)单步方法；(2) 显式格式；(3)局部截断误差为 ?  ? ?h 2? 局部截断误差：当yn ? y ?x n ? 时，由 y ?x n  。按照欧拉方法计算来的 y 。  n ?1  的误差称为 局部截断误差。即， y(x ) ? y 是局部截断误差。 n ?1 n ?1 按泰勒展开 y ?x n?1 ?? y ?x n ?? hy??x n ?? 1 h2 y ??? ? 2 y ? y ? hf ?x , y ? 欧拉法得： n?1 n n n ? h 2 ? h 2 因此，局部截断误差是 。 如果局部截断误差是O(hp?1 ) ，称这种数值方法是p 阶的。欧拉方法显然是一阶的。一阶的含义： 常微分方程的解为 y(x) ? ax ? b 时，若第 n 步精确，即 y ? y(x ) 前提下，欧拉公式能准 n n 确求解 y(x )。 n?1 y ? y n ?1 n hf (x , y n n ) ? y(x n ) ? hf (x n , y(x )) n ? y(x n ) ? hy (x n ) ? ax n ?b ? ah ? a(x n h) ? b ? ax  n ?1 b ? y(x ) n ?1 §5.1.2 隐式欧拉格式 在对微分方程初值问题进行离散化时，如果用向后差商  y(x n?1  ) ? y(x ) n h  代替方程 y?(x ) ? f (x , y(x )) 中的 y?(x ) ,并用近似值 y 表示 y(x ) , y 表示 y(x ) ,得： n?1 n?1 n?1 n ?1 n?1 n ?1 n n y ? y ? hf (x , y ) n?1 n n?1 n?1 称为隐式欧拉公式，或后退的欧拉公式，它是关于 y 的一个函数方程，其计算远比显式