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=1