潮流计算的计算机算法

1、高等电力系统分析（潮流计算的计算机算法）PQ分解法潮流计算（IEEE14）目录一、MATLAB源程序二、对支路参数（B1）、节点参数（B2）的说明三、带入数据，运行结果、MATLAB源程序clearclose alln=in put(请输入节点数:n=);n仁in put(请输入支路数：n 1=)Jisb=in put(请输入平衡节点号:isb=)pr=i nput(请输入误差精度:pr=);B1=i nput(请输入支路参数：B1=);B2=i nput(请输入节点参数:B2=);n2=in put(请输入PQ节点个数：:n 2=);Y=zeros (n);for i=1: n1P=B1(i

2、,1)；q=B1(i,2);Y(p,q)=Y(p,q)-1/(B1(i,3)+B1(i,4)*1j); %Y(q,p)=Y(p,q);Y(p,p)=Y(p,p)+1/(B1(i,3)+B1(i,4)*1j)+B1(i,6)*1j; %Y(q,q)=Y(q,q)+1/(B1(i,3)+B1(i,4)*1j)+B1(i,6)*1j;end非对角元对角元disp(导纳矩阵Y=);disp(Y)% 下面是求P,Q,V,0矩阵V=zeros(1, n);O=zeros(1, n);P=zeros(1, n);Q=zeros(1, n); G=real(Y);B=imag(Y);for i=1: nP(i

3、)=B2(i,3);Q(i)=B2(i,4);V(i)=B2(i,5);O(i)=B2(i,6);endB3=B(1: n-1,1: n-1); %B4=B(1: n2,1: n2);%不含平衡节点，由节点导纳虚部构成 所有PQ节点% 下面是求 P, Q矩阵DX=0;ICT=1;Mp=1;Mq=1;while ICT=0m1=1;m2=1;for i=1: n所有节点数节点数求 DP,DQ对V矩阵求逆 P/V P/V/B3角=- P/V/V/B3 V=- Q/V/Bif i=isbC(i)=O;D(i)=O;for j1=1:nC(i)=C(i)+V(i)*V(j1)*(G(i,j1)*cos

4、(O(i)-O(j1)+B(i,j1)*si n( O(i)-O(j1)；D(i)=D(i)+V(i)*V(j1)*(G(i,j1)*si n(O(i)-O(j1)-B(i,j1)*cos(O(i)-O( j1)；endDP(m1)=P(i)-C(i); m1=m1+1;if B2(i,2)=1 DQ(m2)=Q(i)-D(i); m2=m2+1;endendendm1=m1-1;%m2=m2-1; %PQ DPQ=DP;DQ; %V仁 V(:,1:m1);V2=diag(V1);V3=i nv(V2);%H=V3*DP;%K=-i nv(B3)*H;%-deltO=V3*K; % max1=

5、max(abs(DP);for i=1:m1if max1prMp=0;elseO(i)=O(i)+deltO(i); Mq=1;endendV4=V(:,1:m2);V5=diag(V4);V6=i nv(V5);L=V6*DQ;N=-i nv(B4)*L; deltV=N; % max2=max(abs(DQ);for i=1:m2if max2 clearclose alln=input（请输入节点数：n=）;n1=input（请输入支路数：n仁）;isb=i nput（请输入平衡节点号：isb=）;pr=input（请输入误差精度:pr=）;B1= input（请输入支路参数:B仁）；

6、B2=input（请输入节点参数:B2=）;n2=input（请输入PQ节点个数:n2=）;Y=zeros （n）;for i=1: n1p=B1(i,1);q=B1(i,2);Y(p,q)=Y(p,q)-1/(B1(i,3)+B1(i,4)*1j);%非对角元Y(q,p)=Y(p,q);Y(p,p)=Y(p,p)+1/(B1(i,3)+B1(i,4)*1j)+B1(i,6)*1j;% 对角元Y(q,q)=Y(q,q)+1/(B1(i,3)+B1(i,4)*1j)+B1(i,6)*1j;enddisp(导纳矩阵Y=);disp(Y)%下面是求PQ,V,0矩阵V=zeros(1, n);O=ze

7、ros(1, n);P=zeros(1, n);Q=zeros(1, n); G=real(Y);B=imag(Y);for i=1: nP(i)=B2(i,3);Q(i)=B2(i,4);V(i)=B2(i,5);O(i)=B2(i,6);endB3=B（1: n-1,1: n-1）;%不含平衡节点，由节点导纳虚部构成B4=B（1:n2,1:n2）;%所有 PQ 节点%下面是求 P,A Q矩阵DX=0;ICT=1;Mp=1;Mq=1;while ICT=0m1=1;m2=1;for i=1: nif i=isbC(i)=0;D(i)=0;for j1=1: nC(i)=C(i)+V(i)*V

8、(j1)*(G(i,j1)*cos(O(i)-O(j1)+B(i,j1)*si n(O(i)-O(j1); D(i)=D(i)+V(i)*V(j1)*(G(i,j1)*si n(O(i)-O(j1)-B(i,j1)*cos(O(i)-O(j1);endDP(m1)=P(i)-C(i); m1=m1+1;if B2(i,2)=1DQ(m2)=Q(i)-D(i); m2=m2+1;endendend m1=m1-1; m2=m2-1;DPQ=DP;DQ;V1=V(:,1:m1);V2=diag(V1);V3=i nv(V2); H=V3*DP;K=-i nv(B3)*H; deltO=V3*K;%

9、 所有节点数%PQ节点数%求 DPQQ%对V矩阵求逆%A P/V%-A P/V/B3%A 角=-Amax1=max(abs(DP);for i=1:m1if max1prMp=0;elseO(i)=O(i)+deltO(i);Mq=1;endendV4=V(:,1:m2);V5=diag(V4); V6=i nv(V5);L=V6*DQ;N=-i nv(B4)*L;deltV=N;% V=-A Q/V/Bmax2=max(abs(DQ);for i=1:m2if max2prMq=0;elseif B2(i,2)=1;V(i)=V(i)+deltV(i)： Mp=1;endendendif M

10、p=0&M q=0ICT=0;elseICT=1;endDX=DX+1;end%迭代结束，开始输出结果disp();disp(迭代次数为：);disp(DX);for i=1: nE(i)=V(i)*cos(O(i)+1j*V(i)*si n(O (i);o(i)= 180*a ngle(E(i)/pi;enddisp();disp(修正后各节点电压标么值为(节点号从小到大排列):);disp(V);disp();disp(修正后各节点电压相角为(节点号从小到大排列):);disp(o);%计算各个节点的功率 disp();disp(各节点的功率为:);for p=1: nC(p)=0;for

11、 q=1: nC(p)=C(p)+conj(Y(p,q)*conj(E(q);endS(p)=E(p)*C(p); end disp(S);%计算各支路的功率for i=1: n1 p=B1(i,1);q=B1(i,2);Si(p,q)=E(p)*(conj(E(p)*conj(Y(p,p)-Y(p,q)+(conj(E(p)-conj(E(q)*conj(Y(p,q);disp();disp(各条支路的首端功率为:);disp(Si(p,q);Si(q,p)=E(q)*(conj(E(q)*conj(Y(q,q)-Y(p,q)+(conj(E(q)-conj(E(p)*conj(Y(p,q)

12、;disp();disp(各条支路的末端功率为:);disp(Si(q,p);DS(i)=Si(p,q)+Si(q,p);disp();disp(各条支路的功率损耗为:);disp(DS(i);end%计算平衡节点功率Sp=0;for i=1: nSp=Sp+V( n)*conj(Y( n,i)*conj(V(i);enddisp();disp(平衡节点功率为:);disp(Sp);请输入节点数：n=14请输入支路数：n 1=20请输入平衡节点号:isb=14请输入误差精度：pr=0.00001请输入支路参数:B1=1 20.01335 0.04211 00;1300.20912 005140

13、0.55618 0051100.058110.1763200.0341110.067010.1710300.01282100.056950.1738800.034621200.25202 0052140.054030.2230400.0492340 0.11001 00 ;31300.17615 005450.031810.084500 ;490.127110.2703800 ;560.082050.1920700 ;6120.094980.198900 ;780.220920.1998800 ;7120.122910.2558100 ;890.170930.3480200 ;8120.066

14、150.1302700 ;10110.046990.1979700.0438;10140.019380.0591700.0528;：请输入节点参数：B2=1 1-0.4780.03921-0.076-0.01610;310 01 0;41-0.295-0.16610;51-0.09-0.05810;61-0.035-0.01810;71-0.061-0.01610;81-0.135-0.05810;91-0.149-0.0510;1020.18301.0450;112-0.9420 1.010;122-0.1120.0471.70;13200.1741.90;1400 01.06 0;请输入P

15、Q节点个数:n2=9导纳矩阵Y=Columns 1 through 50；10.5130 -38.2963i-6.8410 +21.5786i0.0000 + 4.7819i0.0000 + 1.7980i0.0000 +0.0000i-6.8410 +21.5786i9.5680 -34.8916i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 +0.0000i0.0000 + 4.7819i0.0000 + 0.0000i0.0000 -19.5490i0.0000 + 9.0901i0.0000 +0.0000i0.0000 + 1.7980i0.0000

16、+ 0.0000i0.0000 + 9.0901i5.3261 -24.2825i-3.9020+10.3654i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i-3.9020 +10.3654i5.7829-14.7683i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i-1.8809 +4.4029i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 +O.OOOOiColu

17、mns 6 through 100.0000 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000 + -1.6860 + -1.7011 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000 +9.5213 -1.1350 + 0.0000 + 0.0000 +-4.99910.0000 + O.OOOOiO.OOOOi0.0000 + 0.0000i0.0000i-1.6860 + 5.1158i0.0000i-1.9860 + 5.0688i0.0000i0.0000 + 0.0000i0

18、.0000i0.0000 + 0.0000i0.0000i0.0000 + 0.0000i0.0000i0.0000 + 0.0000i5.115 8i0.0000 + 0.0000i5.1939i0.0000 + 0.0000i0.0000i0.0000 + 0.0000i0.0000i-1.8809 + 4.4029i0.0000i3.8359 - 8.4970i0.0000i0.0000 + 0.0000i0.0000i0.0000 + 0.0000i0.0000i0.0000 + 0.0000i0.0000i0.0000 + 0.0000i-30.1895i0.0000 + 0.000

19、0i4.7819i-1.9550 + 4.0941i0.0000i0.0000 + 0.0000i0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i-1.7011 + 5.1939i0.0000 + 0.0000i0.0000 + 3.9679i0.0000 + 0.0000i-1.0259 + 4.2350i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i4.0150 - 5.4279i-2

20、.4890 + 2.2520i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i-1.5260 + 3.1760i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 5.6770i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0

21、000i0.0000 + 0.0000i-2.4890 + 2.2520i6.7249 -10.6697i-1.1370 + 2.3150i0.0000 + 0.0000i0.0000 + 0.0000i-3.0989 + 6.1028i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i-1.4240 + 3.0291i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.00

22、00 + 0.0000i-1.4240 + 3.0291i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i-1.1370 + 2.3150i2.5610 - 5.3440i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i+15.2631iColu mns 11 through 14-1.9860 + 5.0688i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.000

23、0i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i-1.1350 + 4.7819i3.1210 - 9.7941i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 3.9679i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i-1.9550 + 4.0941i-1.5260 + 3.1760i-3.0989 + 6.1028i0.0000 + 0.0000i0.0000

24、 + 0.0000i0.0000 + 0.0000i6.5799 -17.3407i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 5.6770i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 - 5.6770i0.0000 + 0.0000i0.

25、0000 + 0.0000i-1.0259 + 4.2350i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i-4.9991 +15.2631i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i6.0250 -19.3961i迭代次数为:42修正后各节点电压标么值为Colu mns 1 through 10（节点号从小到大排列）1.21281.21481.58091.0450

26、1.56271.53641.56021.62641.67921.6654Colu mns 11 through 141.01001.70001.90001.0600修正后各节点电压相角为（节点号从小到大排列）Colu mns 1 through 10-13.1377-11.9998-14.4540-15.0006-15.3011-15.7590-16.4954-16.3362-16.0220-5.2081Colu mns 11 through 14-12.1568-16.1918-14.4540各节点的功率为：Columns 1 through 5-0.0475 -1.4718 - 1.880

27、0i0.9041 - 1.1423i0.0000 - 0.0000i-0.1314 + 0.1399i0.0646iColumns 6 through 10-1.5950 -0.0300 + 0.0035i0.0039 - 0.1228i-0.0656 - 0.0724i -0.0352 - 0.1437i1.3905iColu mns 11 through 14-0.4870 - 1.2396i3.2005 + 3.0414i1.7586 + 3.1847i-2.7057 + 0.3992i各条支路的首端功率为26.1711 +87.9128i各条支路的末端功率为23.5667 +83.48

28、12i各条支路的功率损耗为：4.9738e+01 + 1.7139e+02i各条支路的首端功率为15.2545 +65.3872i各条支路的末端功率为0.2082 +56.8008i各条支路的功率损耗为：1.5463e+01 + 1.2219e+02i各条支路的首端功率为15.3538 +59.6752i各条支路的末端功率为12.6813 +60.6685i各条支路的功率损耗为：2.8035e+01 + 1.2034e+02i各条支路的首端功率为18.4735 +62.4536i各条支路的末端功率为11.6194 +39.6840i各条支路的功率损耗为：3.0093e+01 + 1.0214e

29、+02i各条支路的首端功率为18.0012 +62.4930i各条支路的末端功率为5.5097 +16.2405i各条支路的功率损耗为23.5109 +78.7335i各条支路的首端功率为17.0452 +57.7870i各条支路的末端功率为11.7622 +39.7706i各条支路的功率损耗为28.8075 +97.5576i各条支路的首端功率为13.5220 +59.6677i各条支路的末端功率为19.6150 +58.2875i各条支路的功率损耗为：3.3137e+01 + 1.1796e+02i各条支路的首端功率为16.5471 +56.5546i各条支路的末端功率为6.9281 +27.4025i各条支路的功率损耗为23.4752 +83.9571i各条支路的首端功率为-0.2082 +69.5644i各条支路的末端功率为12.7806 +79.1