变步长的龙格库塔法.ppt

变步长的龙格—库塔方法公式PART1产品品鉴会活动方案汇报人姓名7.4.6变步长的龙格-库塔法在微分方程的数值解中，选择适当的步长是非常重要的。单从每一步看，步长越小，截断误差就越小；但随着步长的缩小，在一定的求解区间内所要完成的步数就增加了。这样会引起计算量的增大，并且会引起舍入误差的大量积累与传播。因此微分方程数值解法也有选择步长的问题。以经典的四阶龙格-库塔法(7.20)为例。从节点xi出发，先以h为步长求出一个近似值，记为，由于局部截断误差为，故有当h值不大时，式中的系数c可近似地看作为常数。来判断所选步长是否适当添加标题这表明以作为的近似值，其误差可用步长折半前后两次计算结果的偏差添加标题由此可得添加标题这样添加标题然后将步长折半,即以为步长,从节点xi出发,跨两步到节点xi+1,再求得一个近似值,每跨一步的截断误差是,因此有添加标题当要求的数值精度为ε时：如果Δε，反复将步长折半进行计算，直至Δε为止,并取其最后一次步长的计算结果作为如果Δε，反复将步长加倍，直到Δε为止，并以上一次步长的计算结果作为。这种通过步长加倍或折半来处理步长的方法称为变步长法。表面上看，为了选择步长，每一步都要反复判断Δ，增加了计算工作量，但在方程的解y(x)变化剧烈的情况下，总的计算工作量得到减少，结果还是合算的。产品品鉴会活动方案汇报人姓名其中?i(i=1,…,m)，?i(i=2,…,m)和?ij(i=2,…,m;j=1,…,i?1)均为待定系数，确定这些系数的步骤与前面相似。§2Runge-KuttaMethod)...,(......),(),(),(]...[1122112321313312122122111--++++++=+++=++==++++=mmmmmmimiiiiiimmiihKhKhKyhxfKhKhKyhxfKhKyhxfKyxfKKKKhyybbbabbaballl?高阶Runge—KuttaMethod?Gill公式：4阶经典龙格-库塔公式的一种改进§2Runge-KuttaMethod?最常用为四级4阶经典龙格-库塔法/*ClassicalRunge-KuttaMethod*/：§2Runge-KuttaMethod注：?龙格-库塔法的主要运算在于计算Ki的值，即计算f的值。Butcher于1965年给出了计算量与可达到的最高精度阶数的关系：753可达到的最高精度642每步须算Ki的个数?由于龙格-库塔法的导出基于泰勒展开，故精度主要受解函数的光滑性影响。对于光滑性不太好的解，最好采用低阶算法而将步长h取小。