Runge_Kutta方法_数值分析.docVIP

WORD格式整理 专业资料 值得拥有 第八章　常微分方程数值解 8.3 Runge-Kutta方法8.3.1 显式 Runge-Kutta法的一般形式　　上节已给出与初值问题(8.1.1)等价的积分形式　　　　　　　　　　(8.3.1)只要对右端积分用不同的数值求积公式近似就可得到不同的求解初值问题(8.1.1)的数值方法，若用显式单步法　　　　　　　　　　(8.3.2)当，即数值求积用左矩形公式，它就是Euler法(8.2.2)，方法只有一阶，若取　　　　　　　　　(8.3.3)就是改进Euler法，这时数值求积公式是梯形公式的一种近似，计算时要用二个右端函数f的值，但方法是二阶的.若要得到更高阶的公式，则求积分时必须用更多的f值，根据数值积分公式，可将(8.3.1)右端积分表示为　　　　　注意，右端f中还不能直接得到，需要像改进Euler法(8.2.11)一样，用前面已算得的f值表示为(8.3.3)，一般情况可将(8.3.2)的表示为　　　　　　　　　　　　　　　(8.3.4)其中　　　　　　　　　　　　 这里均为待定常数，公式(8.3.2)，(8.3.4)称为r级的显式Runge-Kutta法，简称R-K方法.它每步计算r个f值(即)，而ki由前面(i-1)个已算出的表示，故公式是显式的.例