数值分析第五章常微分方程数值解法讲解.pptVIP

) ... ( ... 1 1 0 1 1 1 1 0 1 r n r n n n r n r n n n f f f f h y y y y ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 线性多步法的通式可写为： ) , ( j j j y x f f ? 当 ? ? 1 ? 0 时，为 隐式公 式 ; ? ? 1 =0 则为 显式公式 。 ? 基于数值积分的构造法 将 在 上积分，得到 ) , ( y x f y ? ? 1 [ , ] n n x x ? 1 1 ( ) ( ) ( , ( )) n n x n n x y x y x f x y x dx ? ? ? ? ? 只要 近似地算出右边的积分 ，则可通 过 近似 y ( x n +1 ) 。而 选用不同近似式 I ，可得到不 同的计算公式 。例如利用左矩形积分公式得到尤拉公式； 梯形积分公式得到梯形公式。 1 ( , ( )) n n x x I f x y x dx ? ? ? 1 n n y y I ? ? ? 一般地，利用插值原理所建立的一系列数值积分 方法也可以导出解微分方程的一系列计算公式。 运用插值方法的关键在于选取合适的插值节点。 假设已构造出 f ( x , y ( x )) 的插值多项式 P r ( x ), 则 1 1 1 1 ( , ( )) ( ) ( ) n n n n n n x x r x x x n n r x f x y x dx P x dx y y P x dx ? ? ? ? ? ? ? ? ? ? ? 2.Taylor 公式 (3) 表面上看形式简单，但具体构造时 往往很困难，因为按 (2) 式求导，这一过程可能很复 杂。因此通常不直接用 Taylor 公式，而借鉴其思想 提出其它公式。 1. 由此看出，一种方法具有 p 阶精度 ? 公式对不 超过 p 次的多项式准确成立（局部截断误差为 0 ）。 这一等价条件也可以用来判断一种方法的精度。 § 2 Runge-Kutta Method 单步递推法的 基本思想 是从 ( x n , y n ) 点出发，以 某一 斜率 沿直线达到 ( x n +1 , y n +1 ) 点。欧拉法及其各种变形 所能达到的最高精度为 2 阶。 ? 2. Runge — Kutta Method 由微分中值定理，有 ? ? ? ? * 1 * 1 ( ) ( ) ( ) , , 0 1 ( ) ( ) n n n n n n n y x y x y x h f x h y x h k h y x y x hk ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? k * 称为区间 [ x n , x n +1 ] 上的 平均斜率 ，只要知道平均斜 率，就可计算 y ( x n +1 ). 因此只要对平均斜率提供一种 近似算法，则由 (4) 式可导出一种相应的求解公式。 (4) § 2 Runge-Kutta Method ? ? ? ? ? ? ? ? * 1 * 1 1 2 * 1 2 1. 0, , Eulers formula 2. 1, , Eulers implicit formula 3. Trapezoid formula 2 n n n n k f x y x k k f x y x k k k k ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 例 ? ? 1 2 1 * 1 2 ( , ) 4. ( , ) 2 modified Eulers method n n n n k f x y k f x h y hk k k k ? ? ? ? ? ? ? ? ? ? ? ? ? ? § 2 Runge-Kutta Method 由此看出，改进的尤拉公式用 x n 与 x n +1 两个节点 的斜率的算术平均作为平均斜率， x n +1 点的斜率通