1-3及4高斯法和牛顿法.pptVIP

任课教师：王守相; 1.3 高斯一塞德尔法潮流;缺点： 本算法的主要缺点是收敛速度很慢。 病态条件系统，计算往往会发生收敛困难 节点间相位角差很大的重负荷系统； 包含有负电抗支路(如某些三绕组变压器或线路串联电容等)的系统； 具有较长的辐射形线路的系统； 长线路与短线路接在同一节点上，而且长短线路的长度比值又很大的系统。 此外，平衡节点所在位置的不同选择，也会影响到收敛性能。 目前高斯一塞德尔法已很少使用 ;1.4 牛顿一拉夫逊法 1.4.1牛顿一拉夫逊法的一般概念; 对于非线性代数方程组 即 在待求量x的某一个初始估计值x(0)附近，将上式展开成泰勒级数并略去二阶及以上的高阶项，得到如下的经线性化的方程组 ;上式称之为牛顿法的修正方程式。由此可以求得第一次迭代的修正量 将 相加，得到变量的第一次改进值x(1)。接着就从x(1)出发，重复上述计算过程。因此从一定的初值x(0)出发，应用牛顿法求解的迭代格式为:; 上两式中：f’(x)是函数f(x)对于变量x的一阶偏导数矩阵，即雅可比矩阵J；k为迭代次数。 由上两式可见，牛顿法的核心便是反复形成并求解修正方程式。 牛顿法当初始估计值x(0)和方程的精确解足够接近时，收敛速度非常快，具有平方收敛特性。;1.4.2牛顿潮流算法的修正方程式 在将牛顿法用于求解电力系统潮流计算问题时，由于所采用f(x)的数学表达式以及复数电压变量采用的坐标形式的不同，可以形成牛顿潮流算法的不同形式。 以下讨论用得最为广泛的f(x)采用功率方程式模型，而电压变量则分别采用极坐标和直角坐标的两种形式。;(一)极坐标形式 令 则采用极坐标形式的潮流方程是： 对每个PQ节点及PV节点 对每个PQ节点;将上述方程式在某个近似解附近用泰勒级数展开，并略去二阶及以上的高阶项后，得到以矩阵形式表示的修正方程式为: 式中：n为节点总数；m为PV节点数，雅可比矩阵是(2n-m-2)阶非奇异方阵。;(二)直角坐标形式 令 在这里，潮流方程的组成与上不同，对每个节点，都有二个方程式，所以在不计入平衡节点方程式的情况下，总共有2(n-1)个方程式。 ;对每个PQ节点，根据式(1—11)和式(1—12)有：;对每个PV节点，除了有与式(1-39)相同的有功功率方程式之外，还有 采用直角坐标形式的修正方程式为; 仔细分析以上两种类型的修正方程式，可以看出两者具有以下的共同特点。 (1)修正方程式的数目分别为2(n-1)-m及2(n-1)个，在PV节点所占比例不大时，两者的方程式数目基本接近2(n-1)个。 (2)雅可比矩阵的元素都是节点电压的函数，每次迭代，雅可比矩阵都需要重新形成。; (3)分析雅可比矩阵的非对角元素的表示式可见，某个非对角元素是否为零决定于相应的节点导纳矩阵元素Yij是否为零。因此如将修正方程式按节点号的次序排列，并将雅可比矩阵分块，把每个2X2阶子阵；作为分块矩阵的元素，则按节点号顺序而构成的分块雅可比矩阵将和节点导纳矩阵具有同样的稀疏结构，是一个高度稀疏的矩阵。 (4)和节点导纳矩阵具有相同稀疏结构的分块雅可比矩阵在位置上对称，但雅可比矩阵不是对称阵。;1.4.3修正方程式的处理和求解 在本节的开头就已提到，牛顿算法的核心就是反复形成并求解修正方程式。 因此如何有效地处理修正方程式就成为提高牛顿法潮流程序计算速度并降低内存需量的关键所在。; 从算法的发展过程来看，在50年代末就已经提出了牛顿法潮流的雏形。 先是用迭代法求解修正方程式，但遇到迭代法本身不收敛的问题。 用高斯消去法等直接法求解，但如前所分析，修正方程式的数目在2(n-1)左右，如果不利用雅可比矩阵的稀疏特性，当网络节点数增加为N倍，存储雅可比矩阵的内存量将正比于N2倍,利用直接法求解修正方程的计算量将正比于N3倍地增长。 这就限制了牛顿法潮流程序的解题规模，从而使得这种方法的推广应用一度止步不前。 ; 其后正是人们注意到了雅可比矩阵高度稀疏的特点，求解修正方程式时采用了稀疏程序设计技巧，并且发展了一套在消元过程中旨在尽量保持其稀疏性、以减少内存需量并提高计算速度的有效方法(即著名的最优顺序消去法)，才使牛顿法真正得到了突破，因而在60年代中期以后被普遍采用。; 结合修正方程式的求解，目前在实用的牛顿法潮流程序中所包含的程序特点主要有以下三个方面，这些程序特点对牛顿法潮流程序性能的提高起着决定性的作用。 (1)对于稀疏矩阵，在计算机中以“压缩”方式只储存其非零元素，且只有非零元素才参加运算。; (2)修正方程