第8章常微分方程数值解法本章主要内容：1.欧拉法、改进欧拉法2.pdf

第 8 章 常微分方程数值解法 本章 主要内容 ： 1 ．欧拉法、改进欧拉法 . 2 ．龙格 -库塔法。 3 ．单步法的收敛性与稳定性。 重点、难点 一、微分方程的数值解法 在工程技术或自然科学中，我们会遇到的许多微分方程的问题，而我们只能对其中具 有较简单形式的微分方程才能够求出它们的精确解。对于大量的微分方程问题我们需要考 虑求它们的满足一定精度要求的近似解的方法，称为微分方程的数值解法。本章我们主要 dy f ( x, y ) 讨论常微分方程初值问题 dx 的数值解法。 y (x0 ) y0 数值解法的基本思想是：在常微分方程初值问题解的存在区间 [a,b] 内，取 n+1 个节点 a=x0 ＜x 1＜ … ＜xN=b （其中差 hn= x n –xn-1 称为步长，一般取 h 为常数，即等步长），在这 些节点上把常微分方程的初值问题离散化为差分方程的相应问题，再求出这些点的上的差 分方程值作为相应的微分方程的近似值（满足精度要求）。 二、欧拉法与改进欧拉法 欧拉法与改进欧拉法是用数值积分方法对微分方程进行离散化的一种方法。 xn 1 将常微分方程 y f (x, y) 变为 y( x ) y( x ) f (t , y(t ))dt n 1 n x n 1 ．欧拉法（欧拉折线法） 欧拉法是求解常微分方程初值问题的一种最简单的数值解法。 欧拉法的基本思想：用左矩阵公式计算（＊）式右端积分，则得欧拉法的计算公式 b a 为： y y hf (x , y ) (n 0,1,..., N 1) h n 1 n n n N 欧拉法局部截断误差 2 h Rn 1 y ( n 1 ) x n n 1 xn 1 或简记为 O （h2 ）。 2 1 我们在计算时应注意欧拉法是一阶方法，计算误差较大。 欧拉法的几何意义：过点 A 0 0 0 1 1 1 n n n （x ，y ）， A （x ，y ）， …，A （x ，y ），斜率分 别为 f （x 0，y0 ）， f （x 1，y 1）， … ，f （x n ，y n ）所连接的一条折线，所以欧拉法亦称为 欧拉折线法。 例 1 用欧拉法解初值问题 dy 2xy (0 x 1) dx y(0) 1 在 x ＝0 (0.2) 1 处的近似解。（计算过程保留 4 位小数）。 【思路】 用欧拉法求解常微分方程的初值问题时，首先熟练掌握欧拉公式的一般形式， 根据具体题目写出找出欧拉公式的迭代式，并根据初始条件和所给步长进行迭代求解。 解 ∵ f(x ，y) ＝－ 2xy ，h ＝0.2， 欧拉公式为： y y hf ( x , y ) y 0