直角坐标系下牛顿法潮流计算.doc

武汉理工大学《电力系统分析》课程设计说明书 直角坐标系下牛顿法潮流计算 1电力系统潮流计算 潮流计算是电力系统分析中的一种最基本的计算，它的任务是对给定的运行条件确定系统的运行状态，如母线上的电压(幅值及相角)、网络中的功率分布及功率损耗等。在电力系统规划设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性.可靠性和经济性。此外，电力系统潮流计算也是计算系统动态稳定和静态稳定的基础。 2节点导纳矩阵的形成 在图1（a）的简单电力系统中，若略去变压器的励磁功率和线路电容，负荷用阻抗表示，便可以得到一个有5个节点（包括零电位点）和7条支路的等值网络，如图1（b）所示。将接于节点1和4的电势源和阻抗的串联组合变换成等值的电流源和导纳的并联组合，变得到图1（c）的等值网络，其中 ?I?4?y40E4 ?I?1?y10E1 和 分别称为节点1和4的注入电流源。 (a) 2 4? ?4 (c) 图1 电力系统及其网络 以零电位点作为计算节点电压的参考点，根据基尔霍夫定律，可以写出4个独立节点的电流平衡方程如下： 1 武汉理工大学《电力系统分析》课程设计说明书 ??y(U??U?)?I?y10U112121?y12(U2?y23(U3?y24(U4 ? ? ?)?yU??y(U??U?)?y(U??U?)?0?U?120223232424 ? ????U2)?y34(U3?U4)?0? ??)?y(U??U?)?yU??I??U234434044? （2-1） 上述方程组经过整理可以写成 ??YU????????????????????????????I??Y11U11221 ? ??YU??YU??YU??0Y21U?1222233244 ? ??YU??YU??0?????????????Y32U?2333344??YU??YU??I???????????????Y42 （2-2） ；； Y33?y23?y34Y24?Y42??y24 式中， Y11?y10?y12 ； Y22?y20?y23?y24?y12 ；； Y44?y40?y24?y34Y34?Y43??y34 ； Y12?Y21??y12 ； Y23?Y32??y23 。 一般的，对于有n个独立节点的网络，可以列写n个节点方程 ??YU????YU??I??Y11U11221nn1 ? ??YU????YU??I??Y21U12222nn2 ? ????????????????????????? ??YU????YU??I??Yn1U1n22nnnn? （2-3） 也可以用矩阵写成 ???I???Y11???Y12???Y1n??U11 ????????Y21???Y22???Y2n?UI?2?2???? ??????????????????????????? ???Yn1????Yn2?Ynn???Un????In?? （2-4） 或缩写为 YU?I （2-5） 矩阵Y称为节点导纳矩阵。它的对角线元素 Yii 称为节点i的自导纳，其值等称为节点i、j 间的互导纳， 于接于节点i的所有支路导纳之和。非对角线元素 Yij 它等于直接接于节点i、j间的支路导纳的负值。若节点i、j间不存在直接支路，则有 Yij?0 。由此可知节点导纳矩阵是一个稀疏的对称矩阵。 2 武汉理工大学《电力系统分析》课程设计说明书 3牛顿-拉夫逊法潮流计算 3.1牛顿-拉夫逊法的基本原理 牛顿—拉夫逊法（Newton—Raphson法）是求解非线性方程代数方程组的有效迭代计算方法。在牛顿—拉夫逊法的每一次迭代过程中，对非线性方程通过线性化处理逐步近似。下面以单变量加以说明。 设有单变量非线性方程 f(x)?0 (3-1) 求解此方程时。先给出解的近似值 x? x (0) 它与真解的误差为 ?x (0) ，则 x (0) ??x (0) 将满足方程，即 f(x (0) ??x (0) )?0 (3-2) 将(3-8)式左边的函数在x f(x (0) (0) 附近展成泰勒级数，于是便得 (0) 2 ??x (0) )?f(x (0) )? f(x (n) ‘ (0) )?x (0) ? f(x ‘‘ (0) ) (?x ) 2! ?......? f (n) (x (0) ) (?x (0) ) n n! ?.... （3-3） 式中阶导数。 f ‘ (x (0) ),…… f (x (0) ) (0) 分别为函数f(x)在