常微分方程数值法.pptVIP

常微分方程数值法 第一页，共三十四页，2022年，8月28日 华长生制作 * 9.1 引言 9.3 龙格-库塔方法 9.2 欧拉方法 本章要点： 本章主要研究基于微积分数值解法的常 微分方程数值解，主要方法有 线性单步法中的Euler方法、Runge-Kutta方法 第二页，共三十四页，2022年，8月28日 华长生制作 * 9.1 引言 在工程和科学技术的实际问题中，常需要求解微分方程 只有简单的和典型的微分方程可以求出解析解 而在实际问题中的微分方程往往无法求出解析解 在高等数学中我们见过以下常微分方程： -----------(1) -----------(2) 第三页，共三十四页，2022年，8月28日 华长生制作 * -----------(3) (1),(2)式称为初值问题,(3)式称为边值问题 -----------(4) 另外,在实际应用中还经常需要求解常微分方程组: 本课程主要研究问题(1)的数值解法,对(2)~(4)只作简单介绍 我们首先介绍初值问题(1)的解存在的条件 第四页，共三十四页，2022年，8月28日 华长生制作 * 定理1. 对于问题(1),要求它的数值解 第五页，共三十四页，2022年，8月28日 华长生制作 * -----------(1) 从(1)的表达式 可以看出,求它的数值解的关键在于 而数值微分或数值积分问题我们都已经学习过 第六页，共三十四页，2022年，8月28日 华长生制作 * 9.2 Euler 方法 　　考虑一阶常微分方程初值的问题： 设f（x,y）是连续函数，对y满足Lipschitz条件，这样初值问题的解是存在唯一的，而且连续依赖于初始条件。 为了求得离散点上的函数值，将微分方程的连续问题进行离散化。一般是引入点列{ }，这里　　　　　　　　　　　　　　　为步长，经常考虑定长的情形，即 　　　。 记 为初始问题（1）的问题准确解 在 处的值，用均差近似代替（1）的导数得 (1) 第七页，共三十四页，2022年，8月28日 华长生制作 * 令 为 的近似值，将上面两个近似写成等式，整理后得 （2） （3） 从 处的初值 开始，按（2）可逐步计算以后各点上的值。称 （2）式为显式Euler。由于（3）式的右端隐含有待求函数值 ， 不能逐步显式计算，称（3 ）式为隐式Euler公式或后退Euler公式。如果 将（2）和（3）两式作算术平均，就得梯形公式。 第八页，共三十四页，2022年，8月28日 华长生制作 * 梯形公式也是隐式公式。以上公式都是由 去计算 ，故称它们为单步法。 例1 取h=0.1，用Euler方法、隐式Euler方法和梯形方法解 解 本题有 如果用Euler方法，由（2） 并代入h=0.1得 同理，用隐式Euler方法有 （4） 第九页，共三十四页，2022年，8月28日 华长生制作 * 同样可得梯形公式的显式为 三种方法及准确解 的数值结果如表1所示。从表中看 到，在 处，Euler方法和隐式Euler方法的误差 分 别是 和 ，而梯形方法的误差却是 。