1、高等電力系統(tǒng)分析（潮流計算的計算機算法）PQ分解法潮流計算（IEEE14）精品文庫精品文庫歡迎下載 歡迎下載 目錄一、MATLAB源程序二、對支路參數(shù)（B1）、節(jié)點參數(shù)（B2）的說明三、帶入數(shù)據(jù)，運行結(jié)果一、MATLAB源程序clearclosealln=input(請輸入節(jié)點數(shù):n=);n1=input(請輸入支路數(shù):n1=);isb=input(請輸入平衡節(jié)點號：isb=);pr=input(請輸入誤差精度:pr=);Bl=input(請輸入支路參數(shù):B1=);B2=input(請輸入節(jié)點參數(shù)：B2=);n2=input(請輸入PQ節(jié)點個數(shù)：n2=);Y=zeros(n);fori=1:n

2、1p=Bl(i,1);q=Bl(i,2);Y(p,q)=Y(p,q)-1/(Bl(i,3)+Bl(i,4)*1j);%非對角元%對角元Y(q,p)=Y(p,q);Y(p,p)=Y(p,p)+1/(Bl(i,3)+Bl(i,4)*1j)+Bl(i,6)*1j;Y(q,q)=Y(q,q)+1/(Bl(i,3)+Bl(i,4)*1j)+Bl(i,6)*1j;enddisp(導(dǎo)納矩陣Y=);disp(Y)%下面是求P,Q,V,0矩陣V=zeros(1,n);0=zeros(1,n);P=zeros(1,n);Q=zeros(1,n);G=real(Y);B=imag(Y);fori=1:nP(i)=B

3、2(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é)點，由節(jié)點導(dǎo)納虛部構(gòu)成B4=B(1:n2,1:n2);%所有PQ節(jié)點%下面是求厶P,AQ矩陣DX=0;ICT=1;Mp=1;Mq=1;whileICT=0m1=1;m2=1;fori=1:nifi=isbC(i)=0;D(i)=0;forj1=1:nC(i)=C(i)+V(i)*V(j1)*(G(i,j1)*cos(0(i)0(j1)+B(i,j1)*sin(0(i)0(j1)；D(i)=D(i)+V(i)*V(j1)*(G(i,j1)*sin(0(

4、i)-0(j1)-B(i,j1)*cos(0(i)-0(j1)；endDP(m1)=P(i)-C(i);m1=m1+1;ifB2(i,2)=1DQ(m2)=Q(i)-D(i);m2=m2+1;endendendm1=m1-1;%所有節(jié)點數(shù)m2=m2-1;%PQ節(jié)點數(shù)DPQ=DP;DQ;%求DP,DQV1=V(:,1:m1);V2=diag(V1);V3=inv(V2);%對V矩陣求逆H=V3*DP;%P/VK=-inv(B3)*H;%-P/V/B3deltO=V3*K;%角=P/V/V/B3max1=max(abs(DP);fori=1:m1ifmax1prMp=0;elseO(i)=O(i)

5、+deltO(i);Mq=1;endendV4=V(:,1:m2);V5=diag(V4);V6=inv(V5);L=V6*DQ;N=-inv(B4)*L;deltV=N;%V=-AQ/V/Bmax2=max(abs(DQ);fori=1:m2ifmax2clearclosealln=input(請輸入節(jié)點數(shù):n=);n1=input(請輸入支路數(shù):n1=);isb=input(請輸入平衡節(jié)點號:isb=);pr=input(請輸入誤差精度:pr=);B1=input(請輸入支路參數(shù):B1=);B2=input(請輸入節(jié)點參數(shù):B2=);n2=input(請輸入PQ節(jié)點個數(shù):n2=);Y=ze

6、ros(n);fori=1:n1p=B1(i,1);q=B1(i,2);Y(p,q)=Y(p,q)-1/(B1(+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(導(dǎo)納矩陣Y=);disp(Y)%下面是求BQ,V,0矩陣V=zeros(1,n);O=zeros(1,n);P=zeros(1,n);Q=zeros(1,n);G=real(Y);B=imag(Y);fori

7、=1:nP(i)=B2(；Q(i)=B2(i,4);V(i)=B2(i,5);O(i)=B2(i,6);endB3=B(1:n-1,1:n-1);%不含平衡節(jié)點，由節(jié)點導(dǎo)納虛部構(gòu)成B4=B(1:n2,1:n2);%所有PQ節(jié)點%下面是求厶P,bQ矩陣DX=0;ICT=1;Mp=1;Mq=1;whileICT=0m1=1;m2=1;fori=1:nifi=isbC(i)=0;D(i)=0;forj1=1:nC(i)=C(i)+V(i)*V(j1)*(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

8、)*sin(O(i)-O(j1)-B(i,j1)*cos(O(i)-O(j1);endDP(m1)=P(i)-C(i);m1=m1+1;ifB2(i,2)=1DQ(m2)=Q(i)-D(i);m2=m2+1;endendendm1=m1-1;m2=m2-1;DPQ=DP;DQ;V1=V(:,1:m1);V2=diag(V1);V3=inv(V2);H=V3*DP;K=-inv(B3)*H;deltO=V3*K;%所有節(jié)點數(shù)%PQ節(jié)點數(shù)%求DP,DQ%對V矩陣求逆%P/V%-P/V/B3%角=-P/V/V/B3max1=max(abs(DP);fori=1:m1ifmax1prMp=0;else

9、O(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;%V=-AQ/V/Bmax2=max(abs(DQ);fori=1:m2ifmax2prMq=0;elseifB2(i,2)=1;V(i)=V(i)+deltV(i);Mp=1;endendendifMp=0&Mq=0ICT=0;elseICT=1;endDX=DX+1;end%迭代結(jié)束，開始輸出結(jié)果disp();disp(迭代次數(shù)為：);disp(DX);fori=1:nE(i)=V(i)*cos(0(i

10、)+1j*V(i)*sin(0(i);o(i)=180*angle(E(i)/pi;endTOC o 1-5 h zdisp();disp(修正后各節(jié)點電壓標(biāo)么值為(節(jié)點號從小到大排列):);disp(V);disp();disp(修正后各節(jié)點電壓相角為(節(jié)點號從小到大排列):);disp(o);%計算各個節(jié)點的功率disp();disp(各節(jié)點的功率為:)；forp=1:nC(p)=0;forq=1:nC(p)=C(p)+conj(Y(p,q)*conj(E(q);endS(p)=E(p)*C(p);enddisp(S);%計算各支路的功率fori=1:n1p=B1(i,1);q=B1(i,

11、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);TOC o 1-5 h zdisp();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)+Si(q,p);disp();disp(各條支路的功率損耗為:)；disp(DS(i

12、);end%計算平衡節(jié)點功率Sp=0;fori=1:nSp=Sp+V(n)*conj(Y(n,i)*conj(V(i);enddisp();disp(平衡節(jié)點功率為:)；disp(Sp);請輸入節(jié)點數(shù):n=14請輸入支路數(shù):n1=20請輸入平衡節(jié)點號:isb=14請輸入誤差精度:pr=0.00001請輸入支路參數(shù):B1=120.013350.0421100;1300.2091200;1400.5561800;1100.058110.1763200.034;1110.067010.1710300.0128;2100.056950.1738800.0346;21200.2520200;2140.0

13、54030.2230400.0492;3400.1100100;31300.1761500;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.0438;10140.019380.0591700.0528;請輸入節(jié)點參數(shù):B2=11-0.4780.03921-0.076-0.01610;310010;4

14、1-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.94201.010;122-0.1120.0471.70;13200.1741.90;140001.060;請輸入PQ節(jié)點個數(shù):n2=91導(dǎo)納矩陣Y=0；Columns1through510.5130-38.2963i-6.8410+21.5786i0.0000+4.7819i0.0000+1.7980i0.0000+O.OOOOi-6.8410+21.5

15、786i9.5680-34.8916i0.0000+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi0.0000+4.7819i0.0000+O.OOOOi0.0000-19.5490i0.0000+9.0901i0.0000+O.OOOOi0.0000+1.7980i0.0000+O.OOOOi0.0000+9.0901i5.3261-24.2825i-3.9020+10.3654i0.0000+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi-3.9020+10.3654i5.7829-14.7683i0.0000+O.OOOOi0.0000+O

16、.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi-1.8809+4.4029i0.0000+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOiO.OOOO+O.OOOOiO.OOOO+O.OOOOiO.OOOO+O.OOOOiO.OOOO+O.OOOOiO.OOOO+O.OOOOiO.OOOO+O.OOOOiO.OOOO+O.OOOOiO.OOOO+O.OOOOi-1.424O+3.O291iO.OOOO+O.OOOOi-1.686O+5.1158i-1.7O11+5.1939iO.OOOO+O.O

17、OOOiO.OOOO+O.OOOOiO.OOOO+O.OOOOi-1.986O+5.O688iO.OOOO+O.OOOOiO.OOOO+O.OOOOiO.OOOO+O.OOOOiO.OOOO+O.OOOOiO.OOOO+O.OOOOiO.OOOO+3.9679iO.OOOO+O.OOOOiO.OOOO+O.OOOOiO.OOOO+O.OOOOiO.OOOO+O.OOOOiO.OOOO+O.OOOOiO.OOOO+5.677OiO.OOOO+O.OOOOiO.OOOO+O.OOOOiO.OOOO+O.OOOOi-1.O259+4.235OiO.OOOO+O.OOOOiO.OOOO+O.OOOOi

18、O.OOOO+O.OOOOiColumns6through1OO.OOOO+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi-1.6860+5.1158iO.OOOO+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi-1.7011+5.1939iO.OOOO+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOiO.OOOO+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi-1.4240+3.0291i

19、0.0000+O.OOOOi-1.88O9+4.4O29i0.0000+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi3.8359-8.497Oi0.0000+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi4.0150-5.4279i-2.4890+2.2520i0.0000+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi-2.4890+2.2520i6.7249-10.6697i-1.1370+2.3150i0.0000+O.OOOOi0

20、.0000+O.OOOOi0.0000+O.OOOOi-1.1370+2.3150i2.5610-5.3440i0.0000+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi9.5213-30.1895i0.0000+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi-1.1350+4.7819i-1.9550+4.0941i-1.5260+3.1760i-3.0989+6.1028i0.0000+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi0

21、.0000+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi0.0000+O.OOOOi-4.9991+15.2631iColumns11through14-1.9860+5.0688i0.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.1350+4.7819i3.1210-9.7941

22、i0.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.0000+0.0000i0.0000+0.0000i0.0000+0.0000i0.0000+0.0000i0.0000+5.6770i0.0000+0.0000i0

23、.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.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.00

24、00+0.0000i0.0000+0.0000i6.0250-19.3961i迭代次數(shù)為:42修正后各節(jié)點電壓標(biāo)么值為（節(jié)點號從小到大排列）：Columns1through101.21281.21481.56271.58091.04501.53641.56021.62641.67921.6654Columns11through141.01001.70001.90001.0600修正后各節(jié)點電壓相角為（節(jié)點號從小到大排列）：Columns1through10-13.1377-11.9998-14.4540-15.0006-15.3011-15.7590-16.4954-16.3362-16.02

25、20-5.2081Columns11through14-12.1568-16.1918-14.45400各節(jié)點的功率為：Columns1through5-0.0475-1.5950-1.4718-1.8800i0.9041-1.1423i0.0000-O.OOOOi-0.1314+0.1399i0.0646iColumns6through10-0.0300+0.0035i0.0039-0.1228i-0.0656-0.0724i-0.0352-0.1437i1.3905iColumns11through14-0.4870-1.2396i3.2005+3.0414i1.7586+3.1847i-

26、2.7057+0.3992i各條支路的首端功率為:26.1711+87.9128i各條支路的末端功率為:23.5667+83.4812i各條支路的功率損耗為：4.9738e+01+1.7139e+02i各條支路的首端功率為:15.2545+65.3872i各條支路的末端功率為:0.2082+56.8008i各條支路的功率損耗為：1.5463e+01+1.2219e+02i各條支路的首端功率為:15.3538+59.6752i精品文庫精品文庫各條支路的末端功率為:11.7622+39.7706i歡迎下載 各條支路的末端功率為:11.7622+39.7706i歡迎下載 各條支路的末端功率為:12.

27、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.5109+78.7335i各條支路的首端功率為:17.0452+57.7870i各條支路的功率損耗為:28.8075+97.5576i精品文庫精品文庫各條支路的首端功率為:0.0000+64.5962

28、i歡迎下載 各條支路的首端功率為:0.0000+64.5962i歡迎下載 各條支路的首端功率為: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各條支路的功率損耗為：1.2572e+01+1.4871e+02i各條支路的末端功率為:-0.0000+37.3498i精品文庫