数值分析龙格—库塔方法讲义.pptxVIP

会计学;7.3.1 显式龙格—库塔法的一般形式;它比欧拉法的?(xn, yn, h)=f(xn, yn), 增加了计算一个函数 f 的值，期望 p=2.若要使得到的公式阶数 p 更大, ? 就必须包含更多的 f 值. 实际上从方程(1.1)等价的积分形式(2.4) ，即;一般说来，点数 r 越多，精度越高，上式右端相当于增量函数? (xn, yn, h)，为得到便于计算的显式方法，可类似于改进欧拉法(3.1),(3.2)，将公式表示为; 注：当 r=1, ?(xn, yn, h) =f(xn, yn) 时，就是欧拉法，此时方法的阶为 p=1. 当 r=2 时，改进欧拉法(3.1)是其中一种，下面将证明其阶 p=2. 要使公式(3.4),(3.5)具有更高的阶 p，就要增加点数 r. 下面我们只就 r=2推导R-K方法. 并给出 r=3,4 时的常用公式，其推导方法与 r=2 时类似，只是计算较复杂.;补充：二元函数的泰勒公式;记号：;定理：;7.3.2 二阶显式R-K方法;其中;将以上结果代入(3.7)，则有;要使公式(3.6)具有 p=2 阶，必须使; 如取a=1/2，则c1= c2=1/2, λ2=μ21=1. 这就是改进的欧拉公式(3.1).;称为中点公式(变形的欧拉公式)，相当于数值积分的中矩形公式.也可以表示为; 对r=2的R-K公式(3.6)能否使局部误差提高到O(h4)? 为此 需把K2多展开一项，从(3.8)的 看到展开式中的项 是不能通过选择参数消掉的. 故r＝2的显式R-K方法的阶只能是 p=2，而不能得到三阶公式.;7.3.3 三阶与四阶显式R-K方法;只要K1, K2将按二元泰勒展开，使Tn+1＝O(h4)，可得待定参数满足方程;这是8个未知数6个方程的方程组，解不是唯一的. 可以得到很多公式. 满足条件(3.12)的公式(3.11)统称为三阶R-K公式. 下面只给出其中一个常见的公式. ; 继续上述过程，经过较复杂的数学演算，可以导出各种四阶R-K公式，下列经典公式是其中常用的一个：; 然而值得指出的是，龙格-库塔方法的推导基于泰勒展开方法，因而它要求所求的解具有较好的光滑性质. 反之，如果解的光滑性差，那么，使用龙格-库塔方法求得的数值解，其精度可能反而不如改进的欧拉方法. 实际计算时，我们应当针对问题的具体特点选择合适的算法.;*7.3.4 变步长的龙格-库塔方法;　　我们考察四阶R-K公式(3.13) ，从节点xn出发，先以h为步长求出一个近似值，记为 ，由于公式???局部截断误差为O(h5)，故;比较(3.14)式和(3.15)式我们看到，步长折半后，误差大约减少到1/16，即有;来判定所选的步长是否合适，具体地说，将区分以下两种情况处理：;人有了知识，就会具备各种分析能力， 明辨是非的能力。 所以我们要勤恳读书，广泛阅读， 古人说“书中自有黄金屋。 ”通过阅读科技书籍，我们能丰富知识， 培养逻辑思维能力； 通过阅读文学作品，我们能提高文学鉴赏水平， 培养文学情趣； 通过阅读报刊，我们能增长见识，扩大自己的知识面。 有许多书籍还能培养我们的道德情操， 给我们巨大的精神力量， 鼓舞我们前进。;第25页/共26页