潮流計(jì)算的計(jì)算機(jī)算法-PQ分解法潮流計(jì)算(IEEE14)

1、高等電力系統(tǒng)分析（潮流計(jì)算的計(jì)算機(jī)算法）PQ分解法潮流計(jì)算（IEEE14）目錄一、MATLAB源程序二、對(duì)支路參數(shù)（B1）、節(jié)點(diǎn)參數(shù)（B2）的說明三、帶入數(shù)據(jù)，運(yùn)行結(jié)果MATLAB源程序clearclose alln=input(請(qǐng)輸入節(jié)點(diǎn)數(shù)：n=);n1=input(請(qǐng)輸入支路數(shù)：n1=);isb=input(請(qǐng)輸入平衡節(jié)點(diǎn)號(hào)：isb=);pr=input(請(qǐng)輸入誤差精度：pr=);B1=input(請(qǐng)輸入支路參數(shù)：B1=);B2=input(請(qǐng)輸入節(jié)點(diǎn)參數(shù)：B2=);n2=input(請(qǐng)輸入PQ節(jié)點(diǎn)個(gè)數(shù)：n2=);Y=zeros(n);for i=1:n1p=B1(i,1);q=B1(i,

2、2);Y(p,q)=Y(p,q)-1/(B1(i,3)+B1(i,4)*1j); %非對(duì)角元Y(q,p)=Y(p,q);Y(p,p)=Y(p,p)+1/(B1(i,3)+B1(i,4)*1j)+B1(i,6)*1j; %對(duì)角元Y(q,q尸丫(q,q)+1/(B1(i,3)+B1(i,4)*1j)+B1(i,6)*1j;enddisp( 導(dǎo)納矩陣Y=);disp(Y)% 下面是求P,Q,V,O矩陣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)=B2(i,3);Q(i)=B

3、2(i,4);V(i)=B2(i,5);O(i)=B2(i,6);endB3=B(1:n-1,1:n-1); %不含平衡節(jié)點(diǎn)，由節(jié)點(diǎn)導(dǎo)納虛部構(gòu)成B4=B(1:n2,1:n2);%所有 PQ點(diǎn)% 下面是求A 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)*V01)*(G(i,j1)*cos(O(i)-O(j1)+B(i,j1)*sin(O(i)-O( j1)；D(i)=D(i)+V(i)*V(j1)*(G(i,j1)*sin(O(i)-

4、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;endendendm1=m1-1;m2=m2-1;%PQDPQ=DP;DQ; %所有節(jié)點(diǎn)數(shù)節(jié)點(diǎn)數(shù)求 DP,DQV1=V(:,1:m1);V2=diag(V1);V3=inv(V2);H=V3*DP;K=-inv(B3)*H; deltO=V3*K;%-又V矩陣求逆A P/VA P/V/B3A 角=-A P/V/V/B3max1=max(abs(DP);for i=1:m1if max1prMp=0;el

5、seO(i)=O(i)+deltO(i);Mq=1;endendV4=V(:,1:m2);V5=diag(V4);V6=inv(V5);L=V6*DQ;A V=- A Q/V/BN=-inv(B4)*L;deltV=N; %max2=max(abs(DQ);for i=1:m2if max2 clearclose alln=input(請(qǐng)輸入節(jié)點(diǎn)數(shù)：n=);n1=input(請(qǐng)輸入支路數(shù)：n1=);isb=input(請(qǐng)輸入平衡節(jié)點(diǎn)號(hào)：isb=);pr=input(,請(qǐng)輸入誤差精度:pr=);B1=input(請(qǐng)輸入支路參數(shù)：B1=);B2=input(請(qǐng)輸入節(jié)點(diǎn)參數(shù)：B2=);n2=inpu

6、t(請(qǐng)輸入PQ節(jié)點(diǎn)個(gè)數(shù)：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);%非對(duì)角元Y(q,p)=Y(p,q);Y(p,p)=Y(p,p)+1/(B1(i,3)+B1(i,4)*1j)+B1(i,6)*1j;% 對(duì)角元Y(q,q尸丫(q,q)+1/(B1(i,3)+B1(i,4)*1j)+B1(i,6)*1j;enddisp(導(dǎo)納矩陣Y=);disp(Y)%下面是求P,Q,V,O矩陣V=zeros(1,n);O=zeros(1,n);P=zeros(1,n);Q=zeros(1,

7、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);%不含平衡節(jié)點(diǎn)，由節(jié)點(diǎn)導(dǎo)納虛部構(gòu)成B4=B(1:n2,1:n2);%所有 PQ 節(jié)點(diǎn)%下面是求AP,AQ矩陣DX=0;ICT=1;Mp=1;Mq=1;while ICT=0m1=1;m2=1;for i=1:n if i=isbC(i)=0；D(i)=0;for j1=1:nC(i尸C(i)+V(i)*V(j1)*(G(i,j1)*cos(O(i)-O01)+B(i,j1)*sin(O(

8、i)-O(j1);D(i尸D(i)+V(i)*V(j1)*(G(i,j1)*sin(O(i)-O01)-B(i,j1)*cos(O(i)-O01); endDP(m1)=P(i)-C(i); m1=m1+1; if B2(i,2)=1 DQ(m2)=Q(i)-D(i); m2=m2+1; end end end m1=m1-1;%所有節(jié)點(diǎn)數(shù)m2=m2-1;%PQ 節(jié)點(diǎn)數(shù)DPQ=DP;DQ;%求 DP,DQV1=V(:,1:m1); V2=diag(V1);V3=inv(V2);%對(duì)V矩陣求逆H=V3*DP;%A P/VK=-inv(B3)*H;%-A P/V/B3deltO=V3*K;% A

9、角=-A P/V/V/B3max1=max(abs(DP);for i=1:m1if max1pr Mp=0;elseO(i)=O(i)+deltO(i);Mq=1;endendV4=V(:,1:m2);V5=diag(V4);V6=inv(V5);L=V6*DQ;N=-inv(B4)*L;deltV=N;%AV=-AQ/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 Mp=0&Mq=0ICT=0;elseICT=1;endDX=DX+1;end%迭代

10、結(jié)束，開始輸出結(jié)果 disp();disp(迭代次數(shù)為：);disp(DX);for i=1:nE(i)=V(i)*cos(O(i)+1j*V(i)*sin(O(i);o(i)= 180*angle(E(i)/pi;enddisp();disp(修正后各節(jié)點(diǎn)電壓標(biāo)么值為(節(jié)點(diǎn)號(hào)從小到大排列)：);disp(V);disp();disp(修正后各節(jié)點(diǎn)電壓相角為(節(jié)點(diǎn)號(hào)從小到大排列):);disp(o);%計(jì)算各個(gè)節(jié)點(diǎn)的功率disp();disp(各節(jié)點(diǎn)的功率為:);for p=1:nc(p)=0;for q=1:nC(p尸C(p)+conj(Y(p,q)*conj(E(q);endS(p尸E(p

11、)*C(p);enddisp(S)；%計(jì)算各支路的功率for i=1:n1p=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);disp();disp(各條支路的末端功率為：)；disp(Si(q,p);DS(i尸Si(p,q)+S

12、i(q,p);disp();disp(各條支路的功率損耗為：)；disp(DS(i);end%計(jì)算平衡節(jié)點(diǎn)功率Sp=0；for i=1:nSp=Sp+V(n)*conj(Y(n,i)*conj(V(i);enddisp();disp(平衡節(jié)點(diǎn)功率為：);disp(Sp);請(qǐng)輸入節(jié)點(diǎn)數(shù)：n=14 請(qǐng)輸入支路數(shù)：n1=20 請(qǐng)輸入平衡節(jié)點(diǎn)號(hào)：isb=14 請(qǐng)輸入誤差精度：pr=0.00001請(qǐng)輸入支路參數(shù)：B1=1 20.01335 0.04211 01300.20912001400.55618001100.058110.1763200.0341110.067010.1710300.0128210

13、0.056950.1738800.0346212 00.25202 00;214 0.05403 0.22304 00.04923400.11001 00;31300.17615 00;450.031810.084500;490.127110.2703800;560.082050.1920700;6120.094980.198900;780.220920.1998800;7120.122910.2558100;890.170930.3480200;8120.066150.1302700;10110.046990.1979700.043810140.019380.0591700.0528請(qǐng)輸入節(jié)

14、點(diǎn)參數(shù):B2=1 1-0.4780.03910;21-0.076-0.01610;310010;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;10 20.18301.0450;11 2-0.94201.010;12 2-0.1120.0471.7 0;13 200.1741.90;14 0001.06 0;請(qǐng)輸入PQ節(jié)點(diǎn)個(gè)數(shù)：n2=9導(dǎo)納矩陣丫=Columns 1 through 50.0000 + 4.7819i0.0000 + 1.798

15、0i0.0000 +0.0000 + 0.0000i0.0000 + 0.0000i0.0000 +0.0000 -19.5490i0.0000 + 9.0901i0.0000 +0.0000 + 9.0901i5.3261 -24.2825i-3.90200.0000 + 0.0000i-3.9020 +10.3654i5.78290.0000 + 0.0000i0.0000 + 0.0000i-1.8809 +0.0000 + 0.0000i0.0000 + 0.0000i0.0000 +10.5130 -38.2963i0.0000i-6.8410 +21.5786i0.0000i0.0

16、000 + 4.7819i0.0000i0.0000 + 1.7980i+10.3654i0.0000 + 0.0000i-14.7683i0.0000 + 0.0000i4.4029i0.0000 + 0.0000i-6.8410 +21.5786i9.5680 -34.8916i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000iColumns 6 through 100.0000 +0.0000 +0.0000 +0.0000 +0.0000 +0.0000 +0.0

17、000 +-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 +0.0000 + 0.0000i0.0000i0.0000 + 0.0000i0.0000i-1.6860 + 5.1158i0.0000i-1.9860 + 5.0688i0.0000i0.0000 + 0.0000i0.0000i0.0000 + 0.0000i0.0000i0.0000 + 0.0000i0.0000i0.0000 + 0.0000i5.1158i0.00

18、00 + 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.0000i4.7819i-1.9550 + 4.0941i0.0000i0.0000 + 0.0000i0.0000i0.0000 + 0.0000i0.0000 +

19、 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.4890 + 2.2520i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i-1.5260 + 3.1760i0.0000 + 0.0000i

20、0.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.0000i0.0000 + 0.0000i-2.4890 + 2.2520i6.7249 -10.6697i-1.1370 + 2.3150i0.0000 + 0.0000i0.0000 + 0.0000i-3.0989 +

21、6.1028i0.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.0000 + 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.

22、0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i -4.9991+15.2631iColumns 11 through 14-12.1568 -16.1918 -14.45400-1.9860 + 5.0688i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0000i0.0000 + 0.0

23、000i0.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 + 0.0000i0.0000 + 0.0000i6.5799 -17.3407i0.00

24、00 + 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.0000 + 0.0000i-1.0259 + 4.2350i0.0000 + 0.0000

25、i0.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迭代次數(shù)為:42修正后各節(jié)點(diǎn)電壓標(biāo)么值為（節(jié)點(diǎn)號(hào)從小到大排列）：Columns 1 through 101.66541.21281.21481.56271.53641.56021.62641.67921.58091.0450Columns 11

26、 through 141.01001.70001.90001.0600修正后各節(jié)點(diǎn)電壓相角為（節(jié)點(diǎn)號(hào)從小到大排列）：Columns 1 through 10-16.3362-13.1377-11.9998-14.4540-15.0006-15.3011-15.7590-16.4954-16.0220-5.2081Columns 11 through 14各節(jié)點(diǎn)的功率為：Columns 1 through 5-0.0475 -1.4718 - 1.8800i0.9041 - 1.1423i0.0000 - 0.0000i -0.1314 + 0.1399i0.0646iColumns 6 thr

27、ough 10-1.5950 -0.0300 + 0.0035i0.0039 - 0.1228i-0.0656 - 0.0724i -0.0352 - 0.1437i1.3905iColumns 11 through 14-0.4870 - 1.2396i3.2005 + 3.0414i1.7586 + 3.1847i -2.7057 + 0.3992i各條支路的首端功率為26.1711 +87.9128i各條支路的末端功率為23.5667 +83.4812i各條支路的功率損耗為：4.9738e+01 + 1.7139e+02i各條支路的首端功率為15.2545 +65.3872i各條支路的末

28、端功率為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+02i各條支路的首端功率為18.0012 +62.4930i各條支路的末端功率為5.5097 +16.2405i各條支路的功率損耗為23.5

29、109 +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.1441i各條支路