[工学]电力系统分析2010-2.pdf

研究生学位课： 现代电力系统分析 任课教师：王继东 第六节 保留非线性潮流算法 牛顿法求解非线性潮流方程时采用了逐次线性化 的方法。70年代后期，人们开始研究这样的问 题，即如果采用更加精确的数学模型，将泰勒级 数的高阶项或非线性项也考虑进来，也许会进一 步提高算法的收敛性能及计算速度，于是便产生 了一类称之为保留非线性的潮流算法。又因为其 中大部算法主要包括了泰勒级数的前三项即取到 泰勒级数的二阶项，所以也称为二阶潮流算法， 实现这种想法的第一个尝试是在极坐标形式的牛 顿法修正方程式中增加了泰勒级数的二阶项，所 得到的算法对收敛性能略有改善，但计算速度无 显著提高。 研究表明可以进一步将上式改写为 y y x ( J )x+ Δy +x s (Δ )(0) 该式的推出促成了本算法的突破，因为可 以非常方便地计算二阶项。值得顺便指出 的是，该式是一个很重要的关系式，在研 究其它算法时将多次引用。 y y x ( J )x+ Δy +x s (Δ )(0) （二）数值计算迭代公式 上式是一个以Δx 作为变量的二次代数方程组， 从一定的初值 x （0） Δx 出发，求解满足该式的 仍 然要采用迭代的方法。上式可改写成： 1 (−0) s x J y x [ (y =Δ )y −x ( −)] + Δ 于是算法的具体迭代公式为： ( 1) 1 (0)k + − s ( )k x J y x =Δ[y ( y )−x ( −)] + Δ k表示迭代次数；J为按x=x(0)估计得到；在进行 ( =Δy ) x 0 (0) 第一次迭代时，k=0,令 ，于是和牛顿 法的第一次迭代计算完全相同。 算法的收敛判据为算法的收敛判据为 也可以采用相继二次迭代的二阶项之差，作为收 敛判据，即 (四)算法特点及性能估计 保留非线性快速潮流算法的特点可以通过 和牛顿法进行比较而得以揭示。 ( f) x ( y) x y=−0 s 设求解的方程是 ，则牛 顿法的迭代公式是 此算法的迭代公式是 下图表示了两种算法的迭代过程。 由两组公式可见，与牛顿法的在迭代过程中变化 的雅可比矩阵不同，保留非线性快速潮流算法采 用的是用初值 x （0） 计算而得的恒定雅可比矩 阵，整个计算过程只需一次形成，并三角分解构 成因子表。 下图表示了两种算法的迭代过程。 就每一次迭代所需的计算量而言，牛顿法要重新计 算y(x(k