编程实现戴维南诺顿等值牛顿拉夫逊潮流计算

1、高等电力网络分析编程实现戴维南诺顿等值牛顿拉夫逊潮流计算目录摘要2一、诺顿等值和戴维南等值21、诺顿等值和戴维南等值原理22、例题解析即程序结果分析2二、牛顿-拉夫逊法潮流分析21、牛顿-拉夫逊算法原理22、编程计算及结果分析2附录 程序清单2参考文献2摘要本文主要实现了两种算法，一是通过对戴维南等值和诺顿等值电路的推导过程研究，总结所需求解电力系统参数，编程实现计算机计算网络等值，使研究网络规模大大减小，提高计算速度。二是利用牛顿-拉夫逊法进行对电力系统进行潮流计算，使潮流计算的速度和收敛性得以快速提高。关键字 戴维南等值和诺顿等值 牛顿-拉夫逊法潮流计算 Matlab一、诺顿等值和戴维南等

2、值1、诺顿等值和戴维南等值原理在电网分析中，有时需要研究从网络的某一端口或多个端口看进去时该网络的表现。每个端口都是由感兴趣的一对网络节点组成的，其中一个节点还可以是公共参考节点（即地节点）。这时，可以把该电网在端口处看成一个等值的电流源或电压源，但要求等值前后端口的电气特性是相同的。这就是常规诺顿等值和戴维南等值的做法。应用诺顿等值和戴维南等值对网络进行化简，需要满足两个条件：（1）、被观察的网络是线性的；（2）、每个端口上的净流入电流为零，即要求每个端口所连接的外部电路与被观察网络没有电磁耦合，各个端口所连接的外部电路之间也没有电气耦合。基本原理：如下图所示，令原来的电力网络有N个节点，地

3、节点作为参考节点不包括在内。从中抽出m个感兴趣的端口，这m个端口分别用下标“，m”来表示，相应端口上的节点对用（p，q），（k，l）等来表示。每个端口上第一个节点的电流以流出网络为正方向，第二个节点的电流以流入网络为正方向，二者大小相等。第一二个节点和第二个节点之间的电压作为端口电压的正方向。另外不是一般性，第二个节点还可能是参考节点（即地节点）。首先引入节点-端口关联矢量和节点-端口关联矩阵的概念。以端口为例，其对应的N×1维节点-端口矢量为：P、q都不是参考点：M=0 1p -1q 0 Tq是参考点：M=0 1p 0 0 T把所有节点-端口关联矢量按列排起来，就构成了N×

4、;m维的节点-端口关联矩阵ML=M M Mm T设系统原来的网络方程是或者式中，为节点电压列矢量；为节点注入电流列矢量；Y、Z分别为节点导纳矩阵和节点阻抗矩阵。下图（c）的多端口戴维南等值电路的m×m阶等值阻抗矩阵为：戴维南等值电动势即为原网络的m个开路电压，图（b）的多端口诺顿等值电路的m×m阶等值导纳矩阵为诺顿等值电流源为（c）的网络中各端口短路时的短路电流，根据前面规定的正方向，定义端口上的电流矢量和电压矢量分别如下：从这些端口向原网络看进去，节点注入电流由两部分组成，其一是图（a）网络内部的节点注入电流，其二是与它连接的外部电路从端口注入的电流，因此，可以写出网络的

5、节点电压方程如下：由此可得两边同乘，并考虑到，及。再考虑戴维南等值阻抗矩阵，则有多端口戴维南等值电路方程为：则诺顿等值电路方程为：2、例题解析即程序结果分析高等电力网络分析P132 例5.3 如图所示的电力系统，之路阻抗和节点注入电流都标在图上，试以节点和节点为一个端口，节点和地为一个端口，建立两端口诺顿等值和戴维南等值。解：首先建立以地为参考节点的节点导纳矩阵：其逆矩阵为两个端口的关联矢量组成了关联矩阵。其中节点和 组成的端口 的关联矢量，节点与地组成的端口的关联矢量矩阵分别为： 则戴维南等值阻抗矩阵为：为求戴维南等值电动势，首先求各节点电压：求端口戴维南等值电动势：求诺顿等值导纳和诺顿等值

6、电流：程序流程见附录对教材133页例5.3验证，输出结果为：请输入节点数:n=3请输入支路数:nl=3请输入端口数m=2请输入由支路参数形成的矩阵:B1=1,1,2i;2,2,2i;3,3,2i;1,2,0.2i;1,3,0.1i;2,3,0.2i;请输入各节点参数形成的矩阵:B2=1,1;2,0.5,;3,-1.2;节点导纳矩阵 Y= 0 -15.5000i 0 + 5.0000i 0 +10.0000i 0 + 5.0000i 0 -10.5000i 0 + 5.0000i 0 +10.0000i 0 + 5.0000i 0 -15.5000i端口节点发点m1=1端口节点收点m2=2端口节

7、点发点m1=3端口节点收点m2=0关联矩阵M= 1 -1 0 0 0 1戴维南等值戴维南等值阻抗Zeq= 0 + 0.1164i 0 + 0.0127i 0 + 0.0127i 0 + 0.6970i戴维南等值电动势Veq= 0 + 0.0044i 0 + 0.1440i诺顿等值诺顿等值导纳Yeq= 0 - 8.6094i 0 + 0.1563i 0 + 0.1563i 0 - 1.4375i诺顿等值电流Ieq= 0.00230.1004分析与总结：（1）程序运行结果与理论计算相一致，验证了程序在一定范围内的正确性，因程序的设计是在通用的基础上，故在初始时刻根据所给电力系统的具体参数，对程序参

8、数进行设置。（2）矩阵B1为支路参数矩阵，前N行为N个节点的对地阻抗，第一列为支路首端序号，第二列为支路末端序号，第三列为支路阻抗。矩阵B2为节点参数矩阵，第一列为节点序号，第二类为节点注入电流。二、牛顿-拉夫逊法潮流分析1、牛顿-拉夫逊算法原理电力系统潮流计算是电力系统运行和规划中最基本和最经常的计算，其任务是在已知某些运行参数的情况下，计算出系统中全部的运行参数，一般来说，各个母线所供负荷的功率是已知的，各个节点电压是未知的（平衡节点除外），可以根据网络结构形成节点导纳矩阵，然后由节点导纳矩阵和网络拓扑结构列写功率方程，由于功率方程里功率是已知的，电压的幅值和相角是未知的，这样潮流计算的问

9、题就转化为求解非线性方程组的问题了。为了便于用迭代法解方程组，需要将上述功率方程改写成功率平衡方程，并对功率平衡方程求偏导，得出对应的雅可比矩阵，给未知节点赋电压初值，一般为额定电压，将初值带入功率平衡方程，得到功率不平衡量，这样由功率不平衡量、雅可比矩阵、节点电压不平衡量（未知的）构成了误差方程，解误差方程，得到节点电压不平衡量，节点电压加上节点电压不平衡量构成新的节点电压初值，将新的初值带入原来的功率平衡方程，并重新形成雅可比矩阵，然后计算新的电压不平衡量，这样不断迭代，不断修正，给定收敛精度最后即能收敛。牛顿拉夫逊算法的核心内容是雅可比矩阵，采用极坐标，节点电压表示为节点功率将写成式中，

10、是、两节点电压的相角差。实际上，对于每一个节点或每一个节点都可以列写一个有功功率不平衡量方程式而对于每一个节点还可以再列写一个无功功率不平衡量方程式对于上述方程式可以写出修正方程式如下式中是阶方阵，其元素为；是阶矩阵，其元素为；是阶矩阵，其元素为；是阶方阵，其元素为。可以得到雅可比矩阵元素的表达式如下当时当时2、编程计算及结果分析对如图所示系统进行潮流计算分析，结果如下：导纳矩阵 Y= 150 -50 0 -100 -50 75 -25 0 0 -25 50 -25 -100 0 -25 125雅可比矩阵 J= -150.4618 0 50.0541 0 0 0 0 149.8625 0 -5

11、0.0541 0 0 49.8625 0 -74.4929 0 24.9312 0 0 -49.8625 0 75.0946 0 -24.9312 0 0 25.0405 0 -50.2308 0 0 0 0 -25.0405 0 49.9312迭代次数 K= 3每次迭代没有达到精度要求节点个数 N= 3 3 0最后一次的不平衡量 DW= 1.0e-007 * -0.0319 0 -0.1907 0 -0.0397 0最后一次的修正量 DV= 1.0e-009 * -0.1829 0 -0.4860 0 -0.3214 014各节点电压 U: 1.0011 0.9972 1.0016 1.00

12、0014各节点的功率 S： 0.3000 -0.3000 0.1500 -0.1486分析与总结（1）利用牛顿-拉夫逊法对电力系统进行潮流分析是比较复杂的一个过程，因此未能与前程序一样对多有电力系统通用，而是在某个电力系统基础上进行编程，若有其他系统，可在此基础上做局部修改。（2）矩阵B1为支路参数矩阵，第一列为支路首端号，第二列为支路末端号，第三列为支路串联阻抗。矩阵B2为节点参数矩阵，第一列为注入发电功率，第二列为节点电压参数，第三列为节点类型标号，1为平衡节点，2为PQ节点。附录 程序清单参考文献1 张伯明,陈寿孙,严正.高等电力网络分析2 何仰赞,温增银.电力系统分析.华中科技大学出版

13、社,2002.一、戴维南诺顿等值计算clearn=input('请输入节点数:n=');nl=input('请输入支路数:nl=');m=input('请输入端口数m=');B1=input('请输入由支路参数形成的矩阵:B1=');B2=input('请输入各节点参数形成的矩阵:B2='); Y=zeros(n);Z=zeros(n);M=zeros(m,n);f=zeros(1,n);Zeq=zeros(m);V=zeros(n,1);Veq=zeros(m,1);Yeq=zeros(m);I=zeros(n,

14、1);Ieq=zeros(m,1);%-求节点导纳矩阵for i=(n+1):(n+nl) p=B1(i,1);q=B1(i,2);Y(p,q)=Y(p,q)-1./B1(i,3); Y(q,p)=Y(p,q); Y(q,q)=Y(q,q)+1./B1(i,3); Y(p,p)=Y(p,p)+1./B1(i,3); endfor p=1:n Y(p,p)=Y(p,p)+1./B1(p,3);end disp('节点导纳矩阵 Y=');disp(Y)%-求节点端口关联矩阵for j=1:m m1=input('端口节点发点m1='); m2=input('

15、端口节点收点m2='); for k=1:n if k=m1 f(1,k)=1; else if k=m2&&m2=0 f(1,k)=-1; else f(k)=0; end end end M(j,:)=f;end disp('关联矩阵M=');disp(M)%-求端口戴维南等值阻抗disp('戴维南等值')Z=inv(Y);Zeq=M*Z*M'disp('戴维南等值阻抗Zeq=');disp(Zeq)%-求各节点电压for l=1:n I(l,1)=B2(l,2);endV=Z*I;%-求端口戴维南等值电动势V

16、eq=M*V;disp('戴维南等值电动势Veq=');disp(Veq)%-求端口诺顿等值disp('诺顿等值')Yeq=inv(Zeq);Ieq=Zeq'*Veq;disp('诺顿等值导纳Yeq=');disp(Yeq);disp('诺顿等值电流Ieq=');disp(Ieq)二、牛顿-拉夫逊潮流计算clear;n=4;%节点数：n;nl=4;%支路数：nl;isb=4;%平衡母线节点号：isb;pr=0.00001;%误差精度：pr=;B1=1 2 0.02;1 4 0.01;2 3 0.04;3 4 0.04;B2

17、=0.3 1 2;-0.3 1 2;0.15 1 2;0 1 1;disp('一、潮流分析结果：');Y=zeros(n);e=zeros(1,n);f=zeros(1,n);I=zeros(1,n);%-求导纳矩阵for i=1:nlp=B1(i,1);q=B1(i,2);Y(p,q)=Y(p,q)-1./B1(i,3);Y(q,p)=Y(p,q); Y(q,q)=Y(q,q)+1./B1(i,3);Y(p,p)=Y(p,p)+1./B1(i,3);enddisp('导纳矩阵 Y=');disp(Y)%-K=0;N=1; %K迭代次数；N不满足收敛要求的节点数

18、;for i=1:ne(i)=real(B2(i,2);f(i)=imag(B2(i,2);S(i)=B2(i,1); endG=real(Y);B=imag(Y);P=real(S);Q=imag(S);%-求雅可比矩阵 while N=0 N=0;K=K+1;for i=1:nif i=isbC(i)=0;D(i)=0;for j1=1:nC(i)=C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1);%(Gij*ej-Bij*fj)D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1);%(Gij*fj+Bij*ej)endP0=C(i)*e(i)+f(i

19、)*D(i);%节点功率P=ei(Gij*ej-Bij*fj)+fi(Gij*fj+Bij*ej)Q0=C(i)*f(i)-e(i)*D(i);%节点功率Q=fi(Gij*ej-Bij*fj)-ei(Gij*fj+Bij*ej)DP=P(i)-P0;DQ=Q(i)-Q0; %-求不平衡量-v=2*i-1;w=2*i;DW(v)=DP;DW(w)=DQ; %-for j1=1:nif j1=isb&&j1=iX1=-G(i,j1)*e(i)-B(i,j1)*f(i);% X1=dP/de X2=B(i,j1)*e(i)-G(i,j1)*f(i);% X2=dP/df dP/de=

20、-dQ/dfX3=X2; % X3=dQ/de dP/df=dQ/deX4=-X1; % X4=dQ/dfp=2*i-1;q=2*j1-1;J(p,q)=X1;m=p+1; J(m,q)=X3;q=q+1; J(p,q)=X2;J(m,q)=X4; elseif j1=i&&j1=isbX1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i); X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);p=2*i-1;q=2*j1-1;J(p,q)=X1;m=p+1;J(m,q)=X3;q=q+1;J(p,q)=X2;J(m,q)=X4;endendendendDV=-J(-1)*DW'%disp('每次迭代雅可比矩阵 J=