2020变步长的龙格库塔法.pptVIP

( ) 1 ( ) 5 1 1 1 5 ( 2) 1 5 ( 2) 1 1 ( 2) 1 1 ( ) 1 1 , ( ) , 2 , , 2 ( ) 2 , 2 ( ) 1 . ( ) 16 h n h n n n n h n h n n h n n h n n h y y x y ch h x x h y c h y x y c y x y y x y ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? n ／ ／ ／ 以经典四阶龙格 — 库塔公式为例。从节点 x 出发，以 为 步长求一近似值 将步长折半，即取 为步长从 跨两步到 ，求一近似 值 每跨一步的截断误差是 因此有 由上两式 ( 2) ( 2) ( ) 1 1 1 1 1 ( ) [ ]. 15 h h h n n n n y x y y y ? ? ? ? ? ? ? ? ／ ／ 变步长的龙格 — 库塔方法 1 1 1 1 1 RK ( , ) ( , ) ( 2, 3 , ) ( , ) ( ) p n n i i i n n i i n i n ij j j n n n y y h c K K f x y K f x a h y h b K i p x y Taylor y x x Taylor ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? L 一般地， 方法设近似公式为 确定原则是使近似公式在 处的 展开式与 在 处的 展开式的前面项尽可能多地重合。 1 1 1 1 1 2 1 2 1 ( , ) 1 1 ( ) 2 2 ( , ) ( , ) n n n n n n n n n n y y hK Euler K f x y y y h K K Euler K f x y K f x h y hK ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 公式 改进 公式 公 式 7.4.6 变步长的龙格 - 库塔法 在微分方程的数值解中，选择适当的步长是非常 重要的。单从每一步看，步长越小，截断误差就越 小；但随着步长的缩小，在一定的求解区间内所要 完成的步数就增加了。这样会引起计算量的增大， 并且会引起舍入误差的大量积累与传播。因此微分 方程数值解法也有选择步长的问题。 以经典的四阶龙格 - 库塔法 ( 7.20) 为例。从节点 x i 出发，先以 h 为步长求出一个近似值，记为 ， 由于局部截断误差为 ，故有 ) ( 1 h i y ? ) ( 5 h O 5 ) ( 1 1 ) ( ch y x y h i i ? ? ? ? 当 h 值不大时，式中的系数 c 可近似地看作为常数。 然后将步长折半 , 即以为 步长 , 从节点 x i 出发 , 跨 两步到节点 x i+1 , 再求得一个近似值 , 每跨一步的 截断误差是 , 因此有 2 h ) 2 ( 1 h i y ? 5 2 ? ? ? ? ? ? h c 5 ) 2 ( 1 1 2 2 ) ( ) ( ? ? ? ? ? ? ? ? ? ? h c x y x y h i i 这样 16 1 ) ( ) ( ) ( 1 1 ) 2 ( 1 1 ? ? ? ? ? ? ? h i i h i i y x y y x y ) ( 15 1 ) ( ) ( 1 ) 2 ( 1 ) 2 ( 1 1 h i h i h i i y y y x y ? ? ? ? ? ? ? 由此可得 这表明以 作为 的近似值，其误差可用步 长折半前后两次计算结果的偏差 ) 2 ( 1 h i y ? ) ( 1 ? i x y ) ( 1 ) 2 ( 1 h i h i y y ? ? ? ? ? 来判断所选步长是否适当 当要求的数值精度为ε时： （ 1 ）如果Δε，反复将步长折半进行计算，直 至Δε为止 , 并取其最后一次步长的计算结果作为 （ 2 ）如果Δε，反复将步长加倍，直到Δε 为止，并以上一次步长的计算结果作为 。 这种通过步长加倍或折半来处理步长的方法称为 变步长法。表面上看，为了选择步长，每一步都要 反复判断Δ，增加了计算工作量，但在方程的解 y(x) 变化剧烈的情况下，总的计算工作量得到减少，结 果还是合算的。 1 ? i y 1 ? i y 其中 ? i ( i = 1, …, m ) ， ? i ( i = 2, …, m ) 和 ? ij ( i = 2, …, m ; j = 1, …, i ? 1 ) 均为待定系数，确 定这些系数的步骤与前面相似。 § 2 Runge-Kutta Method ) ...