CH5+常微分方程数值解法—5.1+引言（基本求解公式） .ppt

* 第五章 常微分方程 数值解法 张红梅 自动化学院 2010年4月 5.1 引言 5.1 引言(基本求解公式) 基于数值微分的求解（Euler公式） 基于数值积分的求解（梯形公式 Simpson 公式） 5.2 Runge-Kutta法 本章主要研究基于微积分数值解法的常微分方程数值解，主要研究内容： ? 本章要点 5.1 引言 (基本求解公式) 工程和科学技术的实际问题中,常需要求解微分方程： ----------- (1) 常微分方程 初值问题: ----------- (2) 积分方程: 问:初值问题(1)的解是否存在？如何判断解的存在性？ 本章重点研究形如(1)和(2)问题的数值解法. ? 对于上述问题，可以用解析法求解。但实际问题中的很多常微分方程，解析解很难求得或不存在。 定理1. 如果 f(x,y) 满足 Lipschitz 条件,即 正数 L ,使得 ,均有 则初值问题(1)的解存在且唯一. 初值问题(1)的数值解,即未知函数 y(x) 在区间[a, b]上一系列离散点（节点）： 上函数值 的近似值 ,而 就是初值问题(1)的数值解. 初值问题解的存在唯一性定理 导数 y’(x) 的数值计算或积分 的数值计算 数值积分问题 知，常微分方程初值问题的数值解问题涉及： 由 数值微分问题 常微分方程数值解问题的实质：计算函数 y(x) 在离散点 xk 处函数值 y(xk) 的近似。 “步进法”求解常微分方程 …… 步进法： 单步法： 多步法： 1-1 基于数值微分的求解公式 两点数值微分公式： --------(3) 以下涉及的节点均为等距, 即: 对于初值问题(1),有： --------(4) -------- (5) -------- (6) 得近似解及误差： 其中: (前进)Euler 公式，显示公式，直接计算 yj+1的公式 (4)式中区间 上第一等式的一般格式： 即 -------- (7) -------- (8) (4)式中区间 上第二等式的一般格式： 其中: 得近似解及误差 即 后退 Euler 公式，隐示公式，右端函数含有未知量 yj+1 (5)式和(7)式特点： 称这类方法为单步法／单步格式 Euler方法的几何体现—— 前进 Euler 公式 在计算 yj+1 时只用到之前一个值 yj 后退 Euler 公式 例1. 解: 用前进Euler 公式求解初值问题： 显然 由前进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 --------(9) 预测—校正系统 先由前进 Euler 公式得到 的预测值 ,然后将 代入后退Euler 公式计算 ,就可得到新的Euler 公式: 预测 校正 求解过程： 后退 Euler 公式 中, 为了把 从函数公式中求解出来很难。 实际上，隐式公式公式能获得更高精度的解。 为了避免求解函数方程，采用“显式+隐式”的方式—— 用 Euler 公式的预测——校正系统求解例1. 例2. 解: 由(9)式, 有 依此类推, 得 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中前进Euler公式的结果以及精确值比较发现,用预测——校正系统的结果与精确解更接近. 常微分方程求解公式中,虽然显式比隐式方便,但由隐式可获得更高精度. 评价一个微分方程求解 公式好坏的标准是什么？ 只能表示求解公式第 j+1 步的误差。 在求解公式