第八章 常微分方程数值解法 - 副本.ppt

*/47 §7.2 欧拉方法 /* Convergency */ */47 例：就初值问题 考察欧拉显式格式的收敛性。 解：该问题的精确解为 欧拉公式为 对任意固定的 x = xi = i h ，有 ? ? §2 欧拉方法 */47 §7.2 欧拉方法 例：考察初值问题 在区间[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 What is wrong ??! /* Stability */ */47 §7.2 欧拉方法 一般分析时为简单起见，只考虑试验方程 /* test equation */ 常数，可以是复数 */47 §2 欧拉方法 */47 §3 龙格 —库塔方法 本节内容 一. 引言 二. 2阶龙格 —库塔格式 三. 高阶龙格 —库塔格式 返回章节目录 */47 §7.3 龙格 —库塔方法 a b 几何上的解析：寻找积分曲线 y0 一. 引言 */47 §7.3 龙格 —库塔方法 数值解法的含义：将[a,b]分成n等份，得n+1个节点 计算y(xk) (k0)的近似值yk 。 a b x1 x2 x3 积分曲线 */47 §7.3 龙格 —库塔方法 折线法法的几何构造 I a b x1 x2 x3 Euler法(折线法) 斜率近似 y’(a) y’(x1) y’(x2) y’(x3) 前进Euler法的计算公式 显式公式 */47 §7.3 龙格 —库塔方法 折线法法的几何构造II a b x1 x2 x3 斜率近似 y’(x1) y’(x2) y’(x3) y’(x4) 后退Euler法的公式 隐式公式 */47 单步法：即利用前一个节点的函数值yi，计算后一个节点的函数值yi+1。 目的：建立高精度的单步递推格式。 单步递推法的基本思想是从 ( xi , yi ) 点出发，以某一斜率沿直线达到 ( xi+1 , yi+1 ) 点。欧拉法及其各种变形所能达到的最高精度为2 阶。 §3 龙格 —库塔方法 二. 2阶龙格 —库塔格式 */47 斜率一定取K1 K2 的平均值吗？ 步长一定是一个h 吗？ §3 龙格 —库塔方法 */47 §7.3 龙格 —库塔方法 首先希望能确定系数 ?1、?2、p，使得到的算法格式有2阶精度，即在 的前提假设下，使得 Step 1: 将 K2 在 ( xi , yi ) 点作 Taylor 展开 将改进欧拉法推广为： ) , ( ) , ( ] [ 1 2 1 2 2 1 1 1 phK y ph x f K y x f K K K h y y i i i i i i + + = = + + = + l l */47 Step 2: 将 K2 代入第1式，得到 §3 龙格 —库塔方法 */47 Step 3: 将 yi+1 与 y( xi+1 ) 在 xi 点的泰勒展开作比较 要求 ，则必须有： 这里有3个未知数，2个方程。 存在无穷多个解。所有满足上式的格式统称为2阶龙格 - 库塔格式。 注意到， 就是改进的欧拉法。 §3 龙格 —库塔方法 * */47 第八章 常微分方程数值解 §1 引言 §2 欧拉方法 §3 龙格—库塔方法 §4 阿达姆斯方法 小结 作业与实验 */47 本章要求 1. 熟悉Euler显格式,梯形法及Euler预校法; 2. 熟悉局部截断误差及绝对稳定性; 3. 掌握龙格—库塔法。 */47 一. 问题提出 有一个或多个导数及其函数的方程式称为微分方程，在工程中常遇到求解微分方程的问题。 §1 引言 */47 很多微分方程的解不能用初等函数来表示，有时即使能够用解