Chapter5_4_方程组.pptVIP

§6 常微分方程组数值解法/*Linear multistep method*/ 在平衡点（1，0.0000365，0），Jacobian矩阵为 * 二、刚性(stiff)常微分方程组 一、常微分方程组数值解法 三、高阶微分方程 设一阶方程组: 如果用经典的R-K四阶方法求解,算式为: 一、常微分方程组数值解法 一般地将两个方程的情形推广到m个方程的情形 记 则方程的向量形式为： a≤t≤b 于是经典的四阶R-K方法的向量形式为： 其中： 其分量形式为： 二、刚性(stiff)常微分方程组 设常系数线性微分方程组 其中 假设A有m个互不相同的特征值λj和相应的特征 向量rj(j=1,2,…,m),则通解为: 假设A的所有特征值都有非负实部,即: 且 因此 称为瞬态解. 称为快瞬态解. 称为慢瞬态解. 称为微分方程组的刚性比. 当S1时,即认为是刚性问题. 用显式R-K法求解刚性问题,只有当步长很小时, 才能保证其稳定性. 应用隐式的R-K方法常收到 好的效果. 一个化学反应系统中提出的刚性问题的例子 A，B，C是三个化学样本，Robertson反应如下： 我们通过建模可以得到如下方程组 令 则方程组可写为： 对此Robertson系统而言，Jacobian矩阵为 特征值为 显然当t→+∞时解的各个分量yi(t)(i=1,2,3)是指数衰减的,并趋于稳态解(y1,