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

武汉理工大学《电力系统分析》课程设计说明书 PAGE 1 直角坐标系下牛顿法潮流计算 1电力系统潮流计算 潮流计算是电力系统分析中的一种最基本的计算，它的任务是对给定的运行条件确定系统的运行状态，如母线上的电压(幅值及相角)、网络中的功率分布及功率损耗等。在电力系统规划设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性.可靠性和经济性。此外，电力系统潮流计算也是计算系统动态稳定和静态稳定的基础。 2节点导纳矩阵的形成 在图1（a）的简单电力系统中，若略去变压器的励磁功率和线路电容，负荷用阻抗表示，便可以得到一个有5个节点（包括零电位点）和7条支路的等值网络，如图1（b）所示。将接于节点1和4的电势源和阻抗的串联组合变换成等值的电流源和导纳的并联组合，变得到图1（c）的等值网络，其中和分别称为节点1和4的注入电流源。 图1 电力系统及其网络 以零电位点作为计算节点电压的参考点，根据基尔霍夫定律，可以写出4个独立节点的电流平衡方程如下： （2-1） 上述方程组经过整理可以写成 （2-2） 式中，；；；；；；；。 一般的，对于有个独立节点的网络，可以列写个节点方程 （2-3） 也可以用矩阵写成 （2-4） 或缩写为 （2-5） 矩阵称为节点导纳矩阵。它的对角线元素称为节点的自导纳，其值等于接于节点的所有支路导纳之和。非对角线元素称为节点、 间的互导纳，它等于直接接于节点、间的支路导纳的负值。若节点、间不存在直接支路，则有。由此可知节点导纳矩阵是一个稀疏的对称矩阵。 3牛顿-拉夫逊法潮流计算 3.1牛顿-拉夫逊法的基本原理 牛顿—拉夫逊法（Newton—Raphson法）是求解非线性方程代数方程组的有效迭代计算方法。在牛顿—拉夫逊法的每一次迭代过程中，对非线性方程通过线性化处理逐步近似。下面以单变量加以说明。 设有单变量非线性方程 (3-1) 求解此方程时。先给出解的近似值它与真解的误差为，则将满足方程，即 (3-2) 将(3-8)式左边的函数在附近展成泰勒级数，于是便得 （3-3） 式中,……分别为函数在处的一阶导数，….，n阶导数。 如果差值很小，3-9式右端的二次及以上阶次的各项均可略去。于是，3-9便简化为 ＝0 (3-4) 这是对于变量的修正量的现行方程式，亦称修正方程式。解此方程可得修正量 (3-5) 用所求的去修正近似解，变得 (3-6) 由于3-10是略去高次项的简化式，因此所解出的修正量也只是近似值。修正后的近似解同真解仍然有误差。但是，这样的迭代计算可以反复进行下去，迭代计算的通式是 (3-7) 迭代过程的收敛判据为 (3-8) 或 (3-9) 式中，为预先给定的小正数。 这种解法的几何意义可以从图3－1得到说明。函数y＝f(x)为图中的曲线。f(x)＝0的解相当于曲线与x轴的交点。如果第k次迭代中得到，则过点作一切线，此切线同x轴的交点便确定了下一个近似值。由此可见，牛顿－拉夫逊法实质上就是切线法，是一种逐步线性化的方法。 应用牛顿法求解多变量非线性方程组3-1时，假定已给出各变量的初值，…. ，令，，….. 分别为各变量的修正量，使其满足方程3-2即 (3-10) 将上式中的n个多元函数在初始值附近分别展成泰勒级数,并略去含有，，……，二次及以上阶次的各项，便得 . (3-11) 方程式3-17也可以写成矩阵形式