2013计算方法第八章.pptVIP

1.5 绝对稳定性 一个常微分方程数值解法称为绝对稳定，是指在某一步（如第一步）产生的误差（如计算机的存储误差），在计算中会逐步减小。 称方程 为试验方程，设计算步长为 h ，设在计算开始时产生误差（存储误差），此误差在以后会逐步减弱，我们称该算法相对于 是绝对稳定的， 的全体称为该算法的绝对稳定域。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 欧拉法的绝对稳定区域 后退欧拉法的绝对稳定区域 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 1.6 局部截断误差的实用估计 （1）用两种阶数相同的算法求解，计算出n+1步的近似值，从而得到局部截断误差估计。 （2）用同样的公式，用不同步长计算出n+1步的近似值，从而得到局部截断误差估计。 1.7 隐式法 隐式法具有较好的绝对稳定性！ 只不过在使用隐式法的时候，需要进行迭代，或者使用预测－校正计算格式。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 第二节 泰勒级数法与龙格－库塔法 对于方程： 取计算步长为 h ，则 ，将函数进行泰勒展开 如函数 y( x )有p+1阶导数，容易得到p阶泰勒级数展开法： Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 公式中的导数用下面公式计算： 例子： Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 步长0.1 步长0.01 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 龙格－库塔法： 对于常微分方程的数值解法，一个关键在于选择精度高的算法计算下面公式中的积分。 要高精度的计算积分，常用的方法是适当增加计算节点，不妨 设用 m 个节点计算上面积分，节点为 则积分为 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 将积分改写为： 则得公式： 这样的公式称为显式龙格－库塔公式。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 确定二阶 R－K 法: 系数满足： 此为四个未知数的三个方程，任意取 ，得 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 特别，取 ，得到 通常也称为变形欧拉法，也常写为 它具有二阶精度，也称为二阶二级 R－K 方法。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright