第五章节常微分方程数值解法11课件幻灯片.pptVIP

* * 2、 步长的自动选择 若求满足精度要求?的数值解，则 1） 若? ?，则将步长反复减半，直至??为止！此时的数值解即为所求结果； 2） 若??，则将步长反复加倍，直至??为止！此时再将步长折半一次，即为所求结果； 返回 易见，每一步计算，步长越小，截断误差就越小，但是，随着步长的减小，计算步数增加，这不但导致计算量的增加，而且可能导致舍入误差的严重积累。 因此我们要选择合适的步长，在满足精度要求的前提下，计算量尽量小。 * * 1、 收敛性简介（以显式Euler格式为例 ） 如果对任意固定的 ，数值解 当 （ ）时趋于准确解 ，则称该算法是收敛的。 为理解收敛性的含义，仅考察下列模型方程 对于Euler格式，上述问题的Euler格式具有形式 从而数值解 因此当时 , 即Euler法收敛！ 返回 这样做的根据是： 1）对试验模型不收敛的方法，不可用； 2）一般的初始问题在其解的存在区域内，可局部线性化转化为试验方程。 单步法收敛的充要条件为收敛阶p≧1。 * * 我们称算法是稳定的，如果在节点值 上有大小为 的扰动，而以后各节点值 上产生的偏差值均不超过 。 2、 稳定性简介 考察下列试验模型方程 进一步，若以后各节点产生的误差趋于零，则称算法是绝对稳定的。 1）对试验模型数值不稳定的方法，不可用； 2）一般的初始问题在其解的存在区域内，可局部线性化转化为试验方程。 * * 对于Euler格式，上述模型方程的Euler格式为 当选择h使得 时，Euler格式绝对稳定！ 其绝对稳定区域为 （1）显式欧拉法 * * 对于隐式Euler格式，上述模型方程的隐式Euler格式为 显然绝对稳定区域为 解出yk+1，有 （2）隐式欧拉法 * * （3）梯形公式 * * （4）改进欧拉法 * * （5）经典四级四阶龙格-库塔法 返回 * * 单步法小结 1、单步法计算格式简单，是自开始的算法。 2、步长h的选择原则上受到两个条件的限制: 一是要保证算法收敛，二是要保证算法稳定。 3、单步法的显格式绝对稳定区间相对比较窄，而隐格式绝对稳定区间相对比较宽，有的可是A-稳定的。 4、显式算法没有A-稳定的，高于二阶的隐式算法也没有A-稳定的。 5、利用隐格式绝对稳定区间相对比较宽的特点,可有效得到刚性问题的解。 返回 * * 线性多步法简介 返回 * * 一阶常微分方程组和刚性问题 1、一阶常微分方程组的数值方法 其中 考虑问题： 前面介绍的各种方法可以用来求解微分方程组！(以显式欧拉法为例.) 2、高阶常微分方程数值方法 可以转化为一阶常微分方程组求解！ * * 3、刚性问题及其求解 刚性方程组，是病态问题，应该使用稳定区域较大的隐式方法进行求解！ 返回 尚辅网 / * * 第五章 常微分方程数值解法 §1 问题的提出和基本概念 1、问题的提出; 2、常微分方程初值问题的数值解法; 3、一些基本概念 * * §2 单步法及其收敛性 1、Euler法，梯形公式，改进Euler法 A）设计思想： （1） 差商代替导数； （2） 数值积分法； （3） 泰勒展开法。 B）改进的Euler公式 2、龙格-库塔（Runge-Kutta）法 A) 龙格-库塔法的设计思想； B) 龙格-库塔法的设计方法； C) 常用龙格-库塔法。 C ）总结 * * §3 步长的选择 1、 精度的控制； 2、 步长的自动选择。 §4 收敛性与稳定性简介 1、 收敛性； 2、 稳定性。 单步法小结 §6 一阶常微分方程组和刚性问题 §5 线性多步法简介 * * 1、问题的提出 （1） 设初值问题（1）的解存在且唯一，求解（1）会遇到如下问题： A）（1）的解析解无法表出； B）（1）的解析解可以求出，但是计算其函数值很复杂。 如： 其解析解为 返回 * * 2、常微分方程初值问题的数值解法 给定步长h0，取节点 通过数值方法计算问题(1)的解y(x)在各个节点处的近似值 该数值方法称为初值问题(1)的数值解法，近似解yk称为数值解。 常微分方程数值解法的特点 从初值出发，顺着节点排列次序一步步向前推进，即利用已知 信息 ，来计算 。 ——步进式 返回 * * 3、一些基本概念 （1） 单步法和多步法 单步法：计算yk+1时，只用到前一节点的信息；（自开始） 多步法：计算yk+1时，用到了前面一个以上节点的信息； （2） 显式格式和隐式