袁海燕，曲绍平，贺 丹

（黑龙江工程学院 数学系，黑龙江 哈尔滨150050）

近年来，有关延迟微分方程（DDEs）数值方法稳定性的研究越来越多［1－3］，这些方程广泛出现于科学工程领域，如电路分析、计算机辅助设计、光学控制等。同期出现了有关模糊延迟微分方程组的结论，例如，多芯片的互连问题［4］。Zhu和Petzold用θ方法、Runge－Kutta方法、BDF方法和线性多步法［5］验证过模糊延迟微分方程的渐近稳定性，Yu用多步法研究了一般中立型延迟微分代数方程（GNDDAEs）［6］。 最 近，人 们 对 延 迟 积 分 微 分 方 程（DIDEs）的研究越来越多。

本文主要探讨延迟积分微分方程（DIDEs）的二步Runge－Kutta方法和它的稳定区域，给出并证明DIDEs的二步Runge－Kutta方法的渐近稳定性结果。

1 二步龙格－库塔方法和它的稳定域

考虑二步龙格－库塔方法（TSRK）具有如下形式：

其中：tj＝tn＋cjh，~tj＝tn－1＋cjh，0≤θ≤1；ui是u（ti）的一个近似；h是固定步长；θ，＾bj，~bj和~ai，j，cj是系数。这些方法形成了一类一般线性方法［7］，也可以看作二步分裂方法。这里，将式（1）～（3）表示为

将式（1）～（3）应用到实验方程

为了研究式（1）～（3）的稳定性，必须分析式（5）解的渐近性，由下面的特征多项式的根决定

φ（z）＝z2－R（a，θ）z－S（a，θ）.

二步龙格－库塔方法的稳定区域是所有使得φ（z）的根在单位圆内或圆上的α的集合，其中φ（z）在单位圆上的根是单根。如果对于任意的α，Reα＜0，φ（z）是一个schur多项式，二步龙格－库塔方法的稳定区域包含负半平面，则称此方法是A－稳定的。

2 NMDIDAEs的TSRK 方法的渐近稳定性

这里考虑方程

假定τq＝qτ，τ＝Mh，M 是一个正整数，q＝1，2，…，m。

定义［8］如果一个数值方法用到渐近稳定系统（6）上的数值解满足＝0，则该方法就是渐近稳定的。

将式（1）～（3）代入到式（6）中，得

其中 Kn，i＝，i＝1，2，…，s是乘h 的级导数。……