【摘要】电力系统的暂态过程是典型的非自治动力学系统的行为，其数学模型为一组高维、非线性系统的状态方程和代数方程组，需要用各种的数值方法来求解。其中应用的最广泛的是隐式梯形积分法和龙格—库塔法。本文提出了一种基于Euler-Maclaurin积分公式的单步4阶方法，用于电力系统暂态计算，通过IEEE145节点系统测试表明，该方法有非常好的精度。
　　关键字：电力系统；暂态计算；单步法
　　0 引言
　　暂态稳定性计算是电力系统最基本的计算之一。目前，电力系统暂态稳定分析方法基本分为两种，一种是直接法，它通过检查李雅普诺夫函数的时变性来确定非线性系统的稳定性质，因此它是一种定性的方法，由于构造李雅普诺夫函数比较困难，目前电力系统暂态稳定分析的直接法仅限于比较简单的数学模型，或用暂态能量函数近似李雅普诺夫函数，因此其分析结果不能令人完全满意；另一种方法是数值积分法，其基本思想是用数值积分方法求出描述受扰运动微分方程组的解，然后用各发电机转子之间相对角度的变化判断系统的稳定性，几种常用的数值积分方法是欧拉法、改进欧拉法、Runge-Kutta法和隐式梯形积分法[1]，以及Taylor级数法[2-4]。数值积分法分析结果准确、可靠，所以得到了广泛的实际应用，并一直作为一种标准方法来考察其它分析方法的正确性和精度。
　　本文将一种基于Euler-Maclaurin数值积分方法用于电力系统暂态稳定计算中，通过仿真得到的结果表明该方法有良好的数值稳定性和精确性。
　　1 暂态稳定分析的计算方法
　　电力系统暂态稳定性计算就是分析描述系统的一组微分-代数方程式，即