两类两步Runge-Kutta方法代数稳定性定量误差分析.pdf

摘 要 刚性问题广泛出现于控制系统、生物学、电子网络、物理及化学等重要的科 学和工程领域。二十年多来，刚性问题算法理论的研究迅速发展，现已形成一套 相对成熟的理论—B-理论。但该理论只适用于单边 Lipschitz 常数具有适度大小的 刚性微分方程。 奇异摄动问题是一类典型而又很特殊的刚性问题，由于它的经典 Lipschitz 常 数和单边 Lipschitz 常数都具有 −1 的量级，故不能被已有的经典收敛 O(ε )(0 ε 1) 理论和 B-理论覆盖。因此有必要对其数值方法的稳定性和收敛性作专门的研究。 两步 Runge −Kutta 方法和 Rosenbrock 方法是两类常用的数值方法。用 Rosenbrock 方法（含两步W -方法）求解刚性初值问题可以大大简化计算过程而且 也很容易实现，因此它的计算工作量比一般同级隐式 Runge −Kutta 方法小得多而 仅与 BDF 方法相当。两步 Runge −Kutta 与级数相同的一般单步Runge −Kutta 方法 相比，它对于求解刚性问题具有比较好的稳定性且可达到更高的阶。 本文分成两部分。在第一部分，我们研究了并行两步 W -方法和并行分裂两步 W -方法关于两参数奇异摄动问题的定量收敛性，获得了变步长情况下的局部和 整体误差估计。在第二部分，我们研究了一类两步 Runge −Kutta 方法的代数稳定 性，并提出了几类代数稳定的高阶两步 Runge −Kutta 方法。这将文献中关于两步 Runge −Kutta 方法的线性稳定性讨论拓展到非线性稳定性情形。此外，每部分均 有数值例子证实相应的理论结果。 关键词：刚性问题；两参数奇异摄动问题；并行两步W -方法； 两步Runge −Kutta 方法；代数稳定 I Abstract Stiff problems arise widely in many scientific and engineering fields such as control systems、biology 、electronic network、physics and chemistry, etc. In recent more than twenty years, the study on the theory of computational methods for stiff problems developed rapidly，and we have a relative mature theory—B-theory. But this theory only fits stiff differential equations with moderate-size one-sided Lipschitz constants. Singular perturbation problems are typical classes of stiff problems. But they can’t be covered by the known classical convergence theory and B-theory, due to their classical Lipschitz constants and one-sided Lipschitz constants are in general of size −1 . Therefore, it’s necessary to do special study about stability properties O(ε )(0 ε 1) and convergence properties of their numerical methods. Rosenbrock methods