第六章常微分方程初值问题的数值解法n00907.ppt

§3 收敛性与稳定性 /* Convergency and Stability */ ? 收敛性 /* Convergency */ 定义 　　　 若某算法对于任意固定的 x = xi = x0 + i h，当 h?0 ( 同时 i ? ?) 时有 yi ? y( xi )，则称该算法是收敛的。 例：就初值问题 考察欧拉显式格式的收敛性。 解：该问题的精确解为 欧拉公式为 对任意固定的 x = xi = i h ，有 ? ? ? 稳定性 /* Stability */ 例：考察初值问题 在区间[0, 0.5]上的解。 分别用欧拉显、隐式格式和改进的欧拉格式计算数值解。 0.0 0.1 0.2 0.3 0.4 0.5 精确解 改进欧拉法 欧拉隐式 欧拉显式 节点 xi 1.0000 ?2.0000 4.0000 ?8.0000 1.6000?101 ?3.2000?101 1.0000 2.5000?10?1 6.2500?10?2 1.5625?10?2 3.9063?10?3 9.7656?10?4 1.0000 2.5000 6.2500 1.5626?101 3.9063?101 9.7656?101 1.0000 4.9787?10?2 2.4788?10?3 1.2341?10?4 6.1442?10?6 3.0590?10?7 定义 　　　若某算法在计算过程中任一步产生的误差在以后的计算中都逐步衰减，则称该算法是绝对稳定的 /*absolutely stable */。 一般分析时为简单起见，只考虑试验方程 /* test equation */ 常数，可以是复数 当步长取为 h 时，将某算法应用于上式，并假设只在初值产生误差 ，则若此误差以后逐步衰减，就称该算法相对于 绝对稳定， 的全体构成绝对稳定区域。我们称算法A 比算法B 稳定，就是指 A 的绝对稳定区域比 B 的大。 h l h = h 例：考察显式欧拉法 由此可见，要保证初始误差?0 以后逐步衰减， 必须满足： 0 - 1 - 2 Re Img 例：考察隐式欧拉法 可见绝对稳定区域为： 2 1 0 Re Img 注：一般来说，隐式欧拉法的绝对稳定性比同阶的显式法的好。 常微分方程初值问题的数值解法 第6章 引言 在实际问题中，常需要求解微分方程(如发电机转子运动方程)。只有简单的和典型的微分方程可以求出解析解，而在实际问题中的微分方程往往无法求出解析解。 常微分方程： -----------(1) -----------(2) 一阶常微分方程 -----------(3) (1),(2)式称为初值问题,(3)式称为边值问题 -----------(4) 另外,在实际应用中还经常需要求解常微分方程组: 本课程主要研究问题一阶常微分方程(1)的数值解法, 我们首先介绍初值问题(1)的解存在的条件 定理 只要 f (x, y) 连续，且关于 y 满足 Lipschitz 条件，即存在与 x, y 无关的常数 L 使 对任意定义在 [a, b] 上的 y1(x) 和 y2(x) 都成立，则初值问题（1）存在唯一解。 （通常采用等距节点） 对于问题(1) 要求它的数值解 常微分方程数值解公式的推导 求初值问题数值解的方法是步进法，即从已知的初值y0出发，通过一定的计算求y1 ，然后由y1或y0和y1求出y2 ，依次计算到yn ，即在计算出yk后计算yk+1 ，这时有 单步法：计算yk+1时，只利用yk 多步法：计算yk+1时，用到yk, yk-1, yk-2,… 常微分方程数值解公式的主要推导方法 泰勒展开 利用差商 利用数值积分法 1、泰勒展开的求解思路： 将 按泰勒级数展开 用 的近似值 代入上式右端，记所得结果为 ，则得到数值解序列的计算公式： 2、化导数为差商的求解方法思路： 若在点 处的导数用差商来近似代替，如向前差商 则微分方程初值问题化为 将近似号改为等号，精确解 改为近似解 ，得 3、数值积分的求解思路： 如果将微分方程 在各小区间 上对其两边进行积分，即 如用矩形数值积分公式可得： 以上三种方法推导出同一个数值求解公式: 这个数值公式称为欧拉(Euler)公式。 §6.1 欧拉方法 一、 欧拉格式： x0 x1 向前