第八章_常微分方程的数值解法.pptVIP

四阶龙格 - 库塔法(经典R-K法) /* Classical Runge-Kutta Method */ e(h)=O(h5) 4阶精度 ?优缺点： 优点是：(1)都是一步法，因此只要给定一个初始值就 可以一直计算下去； (2)精度相对较高，如经典R-K法为四阶精度 缺点是：(特别是三阶和四阶法)计算量较大。 ?其他问题讨论： R-K法推导基于Taylor展式，因而要求y(x)有较好的光滑性（即有高阶导数）。 最常用的是四阶公式，它适用于一般的问题，准确、稳定、易于编程。 步长h减小，局部误差O(h5)减小 但步数增加，舍入误差积累增加 h要适当， 总误差才最小 四阶龙格 - 库塔法(经典R-K法) /* Classical Runge-Kutta Method */ 例：取h=0.2，用四阶龙格－库塔法求解初值问题 举例 解： 这里 ，经典的四阶龙格-库塔公式为 举例 表中列出了计算结果，同时列出了相应的精确解.比较本章第一个例子的计算结果，显然龙格-库塔方法的精度高. (Matlab) 习题八 3，7 第八章 作 业 龙格－库塔法 （Runge－Kutta Method ） * 提问学生：数值解法的步进式特点，让学生回答这个特点是怎么回事，然后让学生按照递归的特点，得出