数值分析课件之81.ppt

华长生制作 * 第八章 常微分方程数值解 8.1 Euler方法 8.2 Runge-Kutta法 8.3 线性多步法 在工程和科学技术的实际问题中，常需要求解微分方程 只有简单的和典型的微分方程可以求出解析解 而在实际问题中的微分方程往往无法求出解析解 在高等数学中我们见过以下常微分方程： -----------(1) -----------(2) -----------(3) (1),(2)式称为初值问题,(3)式称为边值问题 本课程主要研究问题(1)的数值解法 我们首先介绍初值问题(1)的解存在的条件 定理1. 对于问题(1),要求它的数值解 求微分方程数值解的方法有两类: 1. 单步法: 主要有Euler法和龙格-库塔法 计算yk时只用到yk-1 2. 多步法: 计算yk时要用到yk-1,,yk-2,…yk-m 主要有线性多步法 第一节 Euler方法 -----------(1) 对于初值问题(1) 在下列子区间上分别应用两点数值微分公式 为了讨论方便,假设以下节点为等距节点 一、欧拉公式 --------(4) (一) Euler公式 由(4)式每组的前一半可得 --------(5) --------(6) 记 其中 (5)和(6)式称为求解初值问题(1)的(前进)Euler公式和误差项 由(5)式每组的后一半可得 记 其中 --------(7) --------(8) (7)和(8)式称为求解初值问题(1)的后退Euler公式和误差项 此类公式称为隐式 从(5)或(7)式不难看出, 这种类型的方法称为单步格式或单步法 Euler方法的几何体现: 前进Euler公式 后退Euler公式 例1. 解: 由前进Euler公式 得 依此类推,有 0 1.0000 0.1000 1.1000 0.2000 1.1918 0.3000 1.2774 0.4000 1.3582 0.5000 1.4351 0.6000 1.5090 0.7000 1.5803 0.8000 1.6498 0.9000 1.7178 1.0000 1.7848 由于后退Euler公式是隐形公式,计算例1将很麻烦 事实上大多数情况下用后退Euler公式都较困难 就可得到新的Euler公式 --------(9) 此方法称为预测—校正系统 用Euler公式的预测——校正系统求解例1. 例2. 解: 由(10)式,有 依此类推,得 0 1.0000 0.1000 1.0918 0.2000 1.1763 0.3000 1.2546 0.4000 1.3278 0.5000 1.3964 0.6000 1.4609 0.7000 1.5216 0.8000 1.5786 0.9000 1.6321 1.0000 1.6819 (二) 常微分方程数值解的截断误差 评价一个微分方程求解公式的标准当然是其精度 而在求解公式 中 误差项 定义1. 因为一般情况下,求解公式的每一步都存在误差,因此有 定义2. 定义3. Euler公式的局部截断误差为 具有1阶精度 后退Euler公式的局部截断误差为 也具有1阶精度 显然一个求解公式的精度越高,计算解的精确性也就越好 从前面的分析可知,Euler法的精度并不算高 二、改进欧拉公式 -----------(1) 对于初值问题 -----------(10) 矩形求积公