851一阶方程组与高阶方程考虑一阶常微分方程组的初值问.pptVIP

第八章常微分方程数值解法 8.5.1 一阶方程组和高阶方程 考虑一阶常微分方程组的初值问题： 那么，（8.5.1）式可以写成 （8.5.2） 若将其中的为知函数、方程的右端项都表示成向量形式： 初始条件表示成 （8.5.1） 8.5 一阶方程组的数值解法 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 可见，式（8.5.2）在形式上与一个方程的初值问题一样。关于一个方程的初值 为体的数值方法均适用于方程组。相应的理论问题也可类似地讨论。下面仅写 出两种数值方法作说明。 梯形方法： 或表达为 其中 是第 个因变量 在节点 处的近似值，相应地， 。 经典R-K方法： 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. 对于高阶方程，可把它转化为一阶方程组。例如，考察下列 m 阶微分方程。 （8.5.3） 只要引进新的变量 则可将 m 阶方程（8.5.3）化为如下的一阶方程组 （8.5.4） 因此，可用求解方程组形式的方法来求解（8.5.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 ： 这两个问题有同样的解 8.5.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方法来计算上面的两个问题，以相同的误差要求来自动选 取步长，计算从 x=0 到 x=10 。第一个问题可用相当大的步长，而第二个问题能 使用的步长小到难以接受。如果改用某种低阶隐式公式，那么这两个问题均可用 较大的步长，计算出大致符合要求的解来。上述显示出来的现象称刚性。问题2 是刚性的，问题1是非刚性的。由于这两个问题的解是相同的，因此这种现象不 是问题的解的作用，而是方程组的一种特性所引起的。基于这个事实，较为正确 的应称之为刚性方程组而不是刚性问题。 考虑方程组的通解。对于问题1，方程组的系数矩阵特征值为 和 ，其通解为 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 其中 为任意常数 其中 为任意常数。对于问题2，方程组的系数矩阵的特征值为 和 ，其通解为 数值计算中出现的现象可以用稳定性来解释。两个问题的特征值都是实的，因 此可只考虑稳定区间。经典 R-K 方法的绝对稳定区间近似为（-2.785，0）。 对于问题1，如果 或 时，可以是稳定的。对于问 题2，要求 或 ，才能保证稳定。由上 可以看出，一定精度范围内，h 完全由绝对稳定性决定。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyri