计算机数值方法第5章节常微分方程数值解法lz.pptVIP

三阶龙格 - 库塔法 三、四阶R-K法思想： 利用各种显化方法，对以上公式显化 Simpson公式求积 xi xi +h xi +h/2 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 三阶龙格 - 库塔法 K1 K2 K3 显Euler P177 要计算3个函数 可以证明(Taylor展开K1，K2，K3） e(h)=O(h4) 具有3阶精度 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. P180 例：取h=0.1，用三阶、 四阶R-K法求解初值问题 举例 解：用三阶R-K法 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 取n=1 计算K1， K2， K3， = y1 取n=2 计算K1， K2， K3， = y2 取n=3 计算K1， K2， K3， = y3 n K1 K2 K3 y1 1 … … … … 2 … … … … 列表计算 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 四阶龙格 - 库塔法(经典R-K法) /* Classical Runge-Kutta Method */ e(h)=O(h5) 4阶精度 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. ?优缺点： 优点是：(1)都是一步法，因此只要给定一个初始值就 可以一直计算下去； (2)精度相对较高，如经典R-K法为四阶精度 缺点是：(特别是三阶和四阶法)计算量较大。 ?其他问题讨论： R-K法推导基于Taylor展式，因而要求y(x)有较好的光滑性（即有高阶导数）。 最常用的是四阶公式，它适用于一般的问题，准确、稳定、易于编程。 步长h减小，局部误差O(h5)减小 但步数增加，舍入误差积累增加 h要适当， 总误差才最小 四阶龙格 - 库塔法(经典R-K法) /* Classical Runge-Kutta Method */ Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. P180 例：取h=0.1，用三阶、 四阶R-K法求解初值问题 解：用三阶R-K法 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 取n=1 计算K1， K2， K3， K4 = y1 取n=2 计算K1， K2， K3， K4 = y2 取n=3 计算K1， K2， K3， K4 = y3 n K1 K2 K3 K4 y1 1 … … … … … 列表计算 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 例：取h=0.2，用四阶龙格－库塔法求解初值问题 举例 解： 这里 ，经典的四阶龙格-库塔公式为 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 举例 表中列出了计算结果，