潮流计算的基本算法及使用方法[参照].pdfVIP

精品文档 ·可编辑版 潮流计算的基本算法及使用方法 一、 潮流计算的基本算法 1. 牛顿－拉夫逊法 1．1 概述 牛顿－拉夫逊法是目前求解非线性方程最好的一种方法。这种方法的特点就是把对非 线性方程的求解过程变成反复对相应的线性方程求解的过程，通常称为逐次线性化过程， 就是牛顿－拉夫逊法的核心。 牛顿-拉夫逊法的基本原理是在解的某一邻域内的某一初始点出发，沿着该点的一阶偏 导数——雅可比矩阵，朝减小方程的残差的方向前进一步，在新的点上再计算残差和雅可 矩阵继续前进，重复这一过程直到残差达到收敛标准，即得到了非线性方程组的解。因为 越靠近解，偏导数的方向越准，收敛速度也越快，所以牛顿法具有二阶收敛特性。而所谓 “某一邻域”是指雅可比方向均指向解的范围，否则可能走向非线性函数的其它极值点， 一般来说潮流由平电压即各母线电压(相角为0，幅值为1)启动即在此邻域内。 1．2 一般概念 对于非线性代数方程组   f x  0     即 f x ,x , ,x  0 i  1,2, n (1 －1) i 1 2 n   0 在待求量x 的某一个初始计算值x 附件，将上式展开泰勒级数并略去二阶及以上的高 阶项，得到如下的线性化的方程组 0  0 0     (1 －2) f x  f x x  0 上式称之为牛顿法的修正方程式。由此可以求得第一次迭代的修正量 0  0 1 0       x   f x f x (1 －3) 将x 0和x 0相加，得到变量的第一次改进值x 1 。接着再从x 1 出发，重复上述计 算过程。因此从一定的初值x 0 出发，应用牛顿法求解的迭代格式为  k  k  k      (1 －4) f x x  f x k 1 k  k  x  x  x (1 －5)     上两式中：f x 是函数f x 对于变量x 的一阶偏导数矩阵，即雅可比矩阵J ； 为迭代 k