北京科技大学计算方法10第十章高斯求积公式及常微分方程初值问题导论.ppt

北京科技大学数理学院 卫宏儒 Weihr168@ 高斯求积公式  常微分方程初值问题的数值解法 4、绝对稳定性问题 注：一阶微分方程组与高阶方程的数值解法 （1）一阶微分方程组的解法 前面介绍的单个方程的各种数值解法完全可以推广到一阶微分方程组的情形，只要将一阶微分方程组中的函数换为向量函数，得到方程组初值问题转化的类似于一阶微分方程的初值问题，按照前面的解法可给出欧拉格式或者经典的R-K方法。 （2）高阶方程的数值解法 对于高阶微分方程，可通过引入变量代换，将其化为一阶方程组，再按（1）的方法求解。 例题： 作业10： 前面已经看到,二级、四级R-K方法可分别达到最高阶数二阶、四阶，但是N级R-K方法的最高阶却不一定是N阶。N表示R-K方法的级数表示公式中计算函数值f的次数。Butcher给出了R-K方法计算函数值f的次数与阶数之间的关系表，如下： 计算f的次数 1 2 3 4 5 6 7 方法的最高阶数 1 2 3 4 4 5 6 由表可见，四级以下R-K的方法其最高阶数与计算f的次数一致，对m阶R-K公式，当m4，虽然计算f的次数增加，但是方法阶数不一定增加。因此四级四阶R-K公式是应用最为广泛的公式。 (1)Euler方法的绝对稳定性 将Euler方法应用于实验方程得到