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

第八章常微分方程数值解法 8.5.1 一阶方程组和高阶方程 考虑一阶常微分方程组的初值问题： 那么，（8.5.1）式可以写成 （8.5.2） 若将其中的为知函数、方程的右端项都表示成向量形式： 初始条件表示成 （8.5.1） 8.5 一阶方程组的数值解法 可见，式（8.5.2）在形式上与一个方程的初值问题一样。关于一个方程的初值 为体的数值方法均适用于方程组。相应的理论问题也可类似地讨论。下面仅写 出两种数值方法作说明。 梯形方法： 或表达为 其中 是第 个因变量 在节点 处的近似值，相应地， 。 经典R-K方法： 其中 或表达为 其中 对于高阶方程，可把它转化为一阶方程组。例如，考察下列 m 阶微分方程。 （8.5.3） 只要引进新的变量 则可将 m 阶方程（8.5.3）化为如下的一阶方程组 （8.5.4） 因此，可用求解方程组形式的方法来求解（8.5.4）。 先考虑两个简单的初值问题。 问题1： 问题2 ： 这两个问题有同样的解 8.5.2 刚性方程组 采用四阶经典R-K方法来计算上面的两个问题，以相同的误差要求来自动选 取步长，计算从 x=0 到 x=10 。第一个问题可用相当大的步长，而第二个问题能 使用的步长小到难以接受。如果改用某种低阶隐式公式，那么这两个问题均可用 较大的步长，计算出大致符合要求的解来。上述显示出来的现象称刚性。问题2 是刚性的，问题1是非刚性的。由于这两个问题的解是相同的，因此这种现象不 是问题的解的作用，而是方程组的一种特性所引起的。基于这个事实，较为正确 的应称之为刚性方程组而不是刚性问题。 考虑方程组的通解。对于问题1，方程组的系数矩阵特征值为 和 ，其通解为 其中 为任意常数 其中 为任意常数。对于问题2，方程组的系数矩阵的特征值为 和 ，其通解为 数值计算中出现的现象可以用稳定性来解释。两个问题的特征值都是实的，因 此可只考虑稳定区间。经典 R-K 方法的绝对稳定区间近似为（-2.785，0）。 对于问题1，如果 或 时，可以是稳定的。对于问 题2，要求 或 ，才能保证稳定。由上 可以看出，一定精度范围内，h 完全由绝对稳定性决定。 由通解（8.5.6）可见，当 时，（8.5.6）右边的第一项和第二项都趋 于零，这两项瞬态解。趋于零的快慢取决于特征值的大小。显然，第二项很快趋 于零，此项称为快瞬态解，而第一项称为慢瞬态解。（8.5.6）右边的第三项称为 稳态解。实际计算表明 ，当第二项很快趋于零以后 ，要使整个方程组的解趋于稳 态解，必须由第一项来决定计算终止与否。因此计算中，要在一个很长的区间上 处处用小步长来计算，这就是刚性现象。 由（8.5.5）和（8.5.6）可见，计算的步 数与量 有关。 一般地，考虑非齐次常系数方程组 ， （8.5.7） 其中， 为常系数矩阵。设 A 有不同的特征值， 它们对应的特征向量 ，则方程组（8.5.7）的通解为 其中 为任意常数， 为（8.5.7）的特解。 假定特征值 的实部为负，即 现设 A 的特征值按其实部的绝对值大小排列： 当 时， 趋向于零，此项称为瞬态解，而 则称为稳态解。如果 大，那么对应 项 当 x 增加快速衰减，此项称为快瞬态解。如果 小，那么对 应的