第六章—常微分方程的数值解法.pptVIP

第六章常微分方程的数值解法 本章内容 §6.1 引言 §6.2 欧拉方法 §6.3 龙格—库塔方法 §6.4 边值问题的数值方法 §6.1 引言 一. 问题提出 有一个或多个导数及其函数的方程式称为微分方程，在工程中常遇到求解微分方程的问题。 §6.1 引言 二. 两类定解问题 ◆常微分方程的定解问题有两种基本类型类:初值问题和边值问题 ◆定解指已知因变量和/或其导数在某些点上是已知的(约束条件)。 ●1. 边值问题 约束条件为已知，在自变量的任一非初值上，已知函数值和/或其导数值，如 例如,受连续分布横向荷载的变截面简支梁弯曲问题 §6.1 引言 §6.1 引言 §6.1 引言 §6.2 欧拉方法及其改进 Euler’s Method 内容 一. 欧拉格式 二. Euler预估—校正法 三. 误差估计、收敛性和稳定性 §6.2 欧拉方法及其改进 §6.2 欧拉方法及其改进 §7.2 欧拉方法 §6.2 欧拉方法及其改进 §7.2 欧拉方法 §7.2 欧拉方法 §6.2 欧拉方法及其改进 §6.2 欧拉方法及其改进 6.2.4 梯形格式 §6.2 欧拉方法及其改进 §7.2 欧拉方法及其改进 §7.2 欧拉方法 §7.2 欧拉方法 §6.3 龙格 —库塔方法 内容 一. 2阶龙格 —库塔格式 三. 高阶龙格 —库塔格式 §6.3 龙格 —库塔方法 二. 2阶龙格 —库塔格式 §6.3 龙格 —库塔方法 §7.3 龙格 —库塔方法 §6.3 龙格 —库塔方法 三. 高阶龙格 —库塔格式 §6.3 龙格 —库塔方法 §6.3 龙格 —库塔方法 §6.3 龙格 —库塔方法 §7.3 龙格 —库塔方法 §7.3 龙格 —库塔方法 §7.3 龙格 —库塔方法 §7.3 龙格 —库塔方法 §7.2 欧拉方法 §6.5 收敛性与稳定性 §7.2 欧拉方法 §7.2 欧拉方法 §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 Step 3: 将 y( xi+1 ) 在 xi 点的泰勒展开并与yi+1作比较 要求 ，则必须有： 这里有3个未知数，2个方程。 存在无穷多个解。所有满足上式的格式统称为2阶龙格 - 库塔格式。 注意到， 就是改进的欧拉法。 Step 2: 将 K2 代入第1式，得到 其中?i ( i = 1, …, m )，?i ( i = 2, …, m ) 和 ?ij ( i = 2, …, m; j = 1, …, i?1 ) 均为待定系数，确定这些系数的步骤与前面相似。 ) ... , ( ... ... ) , ( ) , ( ) , ( ] ... [ 1 1 2 2 1 1 2 32 1 31 3 3 1 21 2 2 1 2 2 1 1 1 - - + + + + + + = + + + = + + = = + + + + = m m m m m m i m i i i i i i m m i i hK hK hK y h x f K hK hK y h x f K hK y h x f K y x f K K K K h y y b b b a b b a b a l l l 问题: 为获得更高的精度，应该如何进一步推广？ 3阶龙格-库塔法 最常用为四级4阶经典龙格-库塔法 注： ? 龙格-库塔法的主要运算在于计算 Ki 的值，即计算 f 的值。Butcher 于1965年给出了计算量与可达到的最高精度阶数的关系： 7 5 3 可达到的最高精度 6 4 2 每步须算Ki 的个数 ? 由于龙格-库塔法的导出基于泰勒展开，故精度主要受解函数的光滑性影响。对于光滑性不太好的解，最好采用低阶算法而将步长h 取小。 例：使用高阶R-K方法计算初值问题 解: （1）使用三阶R-K方法(公式) §6.3 龙格 —库塔方法 其余结果如下: i xi K1 K2