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

第五章 常微分方程的差分方法 一、教学目标及基本要求 通过对本节课的学习，使学生掌握常微分方程、常微分方程方程组的数值解法。 二、教学内容及学时分配 本节课主要介绍常微分方程的数值解法。具体内容如下： 讲授内容：欧拉公式、改进的欧拉公式。 三、教学重点难点 1．教学重点：改进的欧拉公式、龙格库塔方法、收敛性与稳定性。 2. 教学难点：收敛性与稳定性。 四、教学中应注意的问题 多媒体课堂教学为主。适当提问，加深学生对概念的理解。 五、正文 基于数值积分的求解公式：欧拉公式、改进的欧拉公式 引 言 1．主要考虑如下的一阶常微分方程初值问题的求解 : y ( x ) f ( x, y) y (x ) y 0 0 微分方程的解就是求一个函数 y=y(x) ，该函数满足微分方程并且符合初值条件。 2. 例如微分方程 : xy-2y=4x ；初始条件 : y(1)=-3 。 于是可得一阶常微分方程的初始问题 2 y y 4 x y (1) 3 。 2 显然函数 y(x)=x -4x 满足以上条件 ,因而是该初始问题的微分方程的解。 3. 但是 ,只有一些特殊类型的微分方程问题能够得到用解析表达式表示的函数 解，而大量的微分方程问题很难得到其解析解， 有的甚至无法用解析表达式来表 示。因此，只能依赖于数值方法去获得微分方程的数值解。 4 ． 微分方程的数值解： 设微分方程问题的解 y(x) 的存在区间是 [a,b] ，初始点 x0=a，将[a,b]进行划分得一 系列节点 x , x ,...,x ，其中 a= x x … x =b 。y(x) 的解析表达式不容易得到 0 1 n 0 1 n 或根本无法得到，我们用数值方法求得 y(x) 在每个节点 xk 的近似值 y(xk)，即 y ≈y(x)，这样 y , y ,...,y 称为 微分方程的数值解 。 k 0 1 n 如果计算 yn 时，只利用 yn-1，称这种方法为 单步法 ；如果在计算 yn 时不仅利用 yn-1 ，而且还要利用 yn-2, yn-3,…, yn-r,则称这种方法为 r 步方法，也称 多步法 。 §5.1 欧拉方法 §5.1.1 欧拉格式 y (x ) y (x ) n 1 n 方程 y (xn ) f (xn , yn ) 中， y (x ) 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 称为解一阶常微分方程初值问题的欧拉公式，也称显示欧拉公式。 欧拉公式的几何意义非常明显， 因为微分方程的解在 xoy 平面上表示一族积分曲 线。用欧拉公式求数值解的几何意义如图： p2 p1 y1 y0 x0 x1 x2 容易验证，该折线各个顶点的纵坐标 yn (n 1,2...) 就是欧拉公式算得的近似值解， 所以，欧拉方法又称为折线法。 算例： P98