第5章_常微分方程数值解法精读.ppt

§3 Convergency and Stability §3 Convergency and Stability §3 Convergency and Stability 例：隐式龙格-库塔法 而显式 1~ 4 阶方法的绝对稳定区域为 其中2阶方法 的绝对稳定区域为 0 Re Img k = 1 k = 2 k = 3 k = 4 - 1 - 2 - 3 - - - 1 2 3 Re Img 无条件稳定 注：一般来说，隐式欧拉法的绝对稳定性比同阶的显式法的好。 小结 尤拉两步公式与尤拉单步公式相比，使用两个节点上的已知信息将精度提高一阶。可以设想，计算y(xn+1)时，充分利用前面已经求出的节点上的 y 及 y’ 值的线性组合来近似y(xn+1)，精度会大大提高。----线性多步法的基本思想。 构造线性多步法有多种途径，这里介绍两种： 基于数值积分的构造方法； 基于Taylor展开的构造方法。 §4 线性多步法 /* Multistep Method */ 例 HW: p.201 #1-5 证明中点法和梯形公式的精度为2阶 §2 龙格 - 库塔法 /* Runge-Kutta Method */ 建立高精度的单步递推格式:在改进尤拉法和尤拉两步法预测-校正系统中,预测公式都是单步法,如果预测误差很小,则通过校正后得到的近似值误差会更小,因此需要研究高精度的单步法. ? 1. Taylor级数法 IVP: 设其解为y=y(x) 由Taylor展开，有 ? ? ? (1) §2 Runge-Kutta Method (2) §2 Runge-Kutta Method 要使公式具有p阶精度，则在(1)式中截取前p+1项，用(2)式计算各阶导数，即得下面Taylor公式： 局部截断误差 (3) §2 Runge-Kutta Method 2.Taylor公式(3)表面上看形式简单，但具体构造时往往很困难，因为按(2)式求导，这一过程可能很复杂。因此通常不直接用Taylor公式，而借鉴其思想提出其它公式。 1. 由此看出，一种方法具有p阶精度?公式对不超过p次的多项式准确成立（局部截断误差为0）。 这一等价条件也可以用来判断一种方法的精度。 §2 Runge-Kutta Method 单步递推法的基本思想是从 ( xn , yn ) 点出发，以某一斜率沿直线达到 ( xn+1 , yn+1 ) 点。欧拉法及其各种变形所能达到的最高精度为2阶。 ? 2. Runge—Kutta Method 由微分中值定理，有 k*称为区间[xn, xn+1]上的平均斜率，只要知道平均斜率，就可计算y(xn+1).因此只要对平均斜率提供一种近似算法，则由(4)式可导出一种相应的求解公式。 (4) §2 Runge-Kutta Method 例 §2 Runge-Kutta Method 由此看出，改进的尤拉公式用xn与xn+1两个节点的斜率的算术平均作为平均斜率， xn+1点的斜率通过已知信息yn来预测。 ? 考察改进的欧拉法，可以将其改写为： 斜率 一定取K1， K2 的平均值吗？ 步长一定是h 吗？即第二个节点一定是xn+1吗？ §2 Runge-Kutta Method §2 Runge-Kutta Method 首先希望能确定系数 ?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 2: 将 K2 代入第1式，得到 ? 2阶Runge—Kutta Method §2 Runge-Kutta Method Step 3: 将 yi+1 与 y( xi+1 ) 在 xi 点的泰勒展开作比较 要求 ，则必须有： 这里有 个未知数， 个方程。 3 2 存在无穷多个解。所有满足上式的格式统称为2阶龙格 - 库塔格式。 注意到， 就是改进的欧拉法； p=1/2, ?1=0, ?2=1, 变形尤拉公式。 Q: 为获得更高的精度，应该如何进一步推广？ 改进的Euler 公式推广