電力系統分析課程設計--PQ分解法潮流計算

1、電力系統分析課程設計之P-Q 分解法潮流計算程序編程及相應數據、圖形和結果編程:n=5;nl=5;isb=1;pr=0.00001;B1=1 2 0.03i 0 1.05 0;2 3 0.08+0.3i 0.5i 1 0;2 4 0.1+0.35i 0 1 0;3 4 0.04+0.25i 0.5i 1 0;3 5 0.015i 0 1.05 1;B2=0 0 1.05 1.05 0 1;0 3.7+1.3i 1.05 0 0 2;0 2+1i 1.05 0 0 2;0 1.6+0.8i 1.05 0 0 2;5 0 1.05 1.05 0 3;X=1 0;2 0;3 0;4 0;5 0;na

2、=3;Y=zeros(n;YI=zeros(n;e=zeros(1,n;f=zeros(1,n;V=zeros(1,n;O=zeros(1,n;for i=1:nif X(i,2=0;p=X(i,1;Y(p,p=1./X(i,2;endendfor i=1:nlif B1(i,6=0p=B1(i,1;q=B1(i,2;else p=B1(i,2;q=B1(i,1;endY(p,q=Y(p,q-1./(B1(i,3*B1(i,5;YI(p,q=YI(p,q-1./B1(i,3;Y(q,p=Y(p,q;YI(q,p=YI(p,q;Y(q,q=Y(q,q+1./(B1(i,3*B1(i,52+B1(

3、i,4./2;YI(q,q=YI(q,q+1./B1(i,3;Y(p,p=Y(p,p+1./B1(i,3+B1(i,4./2;YI(p,p=YI(p,p+1./B1(i,3;endG=real(Y;B=imag(YI;BI=imag(Y;for i=1:nS(i=B2(i,1-B2(i,2;BI(i,i=BI(i,i+B2(i,5;endP=real(S;Q=imag(S;for i=1:ne(i=real(B2(i,3;f(i=imag(B2(i,3;V(i=B2(i,4;endfor i=1:nif B2(i,6=2V(i=sqrt(e(i2+f(i2;O(i=atan(f(i./e(i;

4、endendfor i=2:nif i=nB(i,i=1./B(i,i;else IC1=i+1;for j1=IC1:nB(i,j1=B(i,j1./B(i,i;endB(i,i=1./B(i,i;for k=i+1:nfor j1=i+1:nB(k,j1=B(k,j1-B(k,i*B(i,j1; endendendendp=0;q=0;for i=1:nif B2(i,6=2p=p+1;k=0;for j1=1:nif B2(j1,6=2k=k+1;A(p,k=BI(i,j1;endendendendfor i=1:naif i=naA(i,i=1./A(i,i;else k=i+1;fo

5、r j1=k:naA(i,j1=A(i,j1./A(i,i;endA(i,i=1./A(i,i;for k=i+1:nafor j1=i+1:naA(k,j1=A(k,j1-A(k,i*A(i,j1;endendendendICT2=1;ICT1=0;kp=1;kq=1;K=1;DET=0;ICT3=1;while ICT2=0|ICT3=0ICT2=0;ICT3=0;for i=1:nif i=isbC(i=0;for k=1:nC(i=C(i+V(k*(G(i,k*cos(O(i-O(k+BI(i,k*sin(O(i-O(k; endDP1(i=P(i-V(i*C(i;DP(i=DP1(i

6、./V(i;DET=abs(DP1(i;if DET>=prICT2=ICT2+1;endendendNp(K=ICT2;if ICT2=0for i=2:nDP(i=B(i,i*DP(i;if i=nIC1=i+1;for k=IC1:nDP(k=DP(k-B(k,i*DP(i;endelsefor LZ=3:iL=i+3-LZ;IC4=L-1;for MZ=2:IC4I=IC4+2-MZ;DP(I=DP(I-B(I,L*DP(L;endendendendfor i=2:nO(i=O(i-DP(i;endkq=1;L=0;for i=1:nif B2(i,6=2C(i=0;L=L+1;

7、for k=1:nC(i=C(i+V(k*(G(i,k*sin(O(i-O(k-BI(i,k*cos(O(i-O(k; endDQ1(i=Q(i-V(i*C(i;DQ(L=DQ1(i./V(i;DET=abs(DQ1(i;if DET>=prICT3=ICT3+1;endendendelse kp=0;if kq=0;L=0;for i=1:nif B2(i,6=2C(i=0;L=L+1;for k=1:nC(i=C(i+V(k*(G(i,k*sin(O(i-O(k-BI(i,k*cos(O(i-O(k; endDQ1(i=Q(i-V(i*C(i;DQ(L=DQ1(i./V(i;DET=

8、abs(DQ1(i;endendendendNq(K=ICT3;if ICT3=0L=0;for i=1:naDQ(i=A(i,i*DQ(i;if i=nafor LZ=2:iL=i+2-LZ;IC4=L-1;for MZ=1:IC4I=IC4+1-MZ;DQ(I=DQ(I-A(I,L*DQ(L; endendelseIC1=i+1;for k=IC1:naDQ(k=DQ(k-A(k,i*DQ(i; endendendL=0;for i=1:nif B2(i,6=2L=L+1;V(i=V(i-DQ(L;endendkp=1;K=K+1;elsekq=0;if kp=0K=K+1;endendf

9、or i=1:nDy(K-1,i=V(i;endenddisp('迭代次數 'disp(K;disp('每次沒有達到精度要求的有功功率個數為 ' disp(Np;disp('每次沒有達到精度要求的無功功率個數為 ' disp(Nq;for k=1:nE(k=V(k*cos(O(k+V(k*sin(O(k*j; O(k=O(k*180./pi; end disp('各節(jié)點的電壓標么值 E 為' disp(E; disp('各節(jié)點的電壓 V 大小' disp(V; disp('各節(jié)點的電壓相角 O' d

10、isp(O; for p=1:n C(p=0; for q=1:n C(p=C(p+conj(Y(p,q*conj(E(q; end S(p=E(p*C(p; end disp('各節(jié)點的功率為' disp(S; disp('各條支路的首端功率為' for i=1:nl if B1(i,6=0 p=B1(i,1;q=B1(i,2; else p=B1(i,2;q=B1(i,1; end Si(p,q=E(p*(conj(E(p*conj(B1(i,4./2+(conj(E(p*B1(i,5conj(E(q*conj(1./(B1(i,3*B1(i,5; disp

11、(Si(p,q; end disp('各條支路的末端功率為' for i=1:nl if B1(i,6=0 p=B1(i,1;q=B1(i,2; else p=B1(i,2;q=B1(i,1; end Sj(q,p=E(q*(conj(E(q*conj(B1(i,4./2+(conj(E(q./B1(i,5conj(E(p*conj(1./(B1(i,3*B1(i,5; disp(Sj(q,p; end disp('各條支路的功率損耗為' for i=1:nl if B1(i,6=0 p=B1(i,1;q=B1(i,2; else p=B1(i,2;q=B1(i

12、,1; end DS(i=Si(p,q+Sj(q,p; disp(DS(i; end for i=1:K Cs(i=i; for j=1:n Dy(K,j=Dy(K-1,j; end end disp('每次迭代后各節(jié)點的電壓值如圖所示' plot(Cs,Dy xlabel('迭代次數' ylabel('電壓' title('電壓迭代次數曲線' 運行結果： 運行結果： 迭代次數 10 每次沒有達到精度要求的有功功率個數為 4 4 4 4 4 4 4 4 3 0 每次沒有達到精度要求的無功功率個數為 3 3 3 3 3 3 3 3

13、1 0 各節(jié)點的電壓標么值 E 為 1.0500 0.6466 2.6784 各節(jié)點的電壓 V 大小 1.0500 1.0364 各節(jié)點的電壓相角 O 0 -4.2819 各節(jié)點的功率為 0 +15.1965i +73.6194i 各條支路的首端功率為 0 +15.1965i -1.0902 - 4.1929i 0.5001 2.9280 1.0779 0.8622 1.0500 17.8535 -4.7785 21.8433 -1.0187 -12.8552i 8.1568 -28.0336i -0.7350 - 4.5044i 0 0.0715 + 0.2503i 3.6408 +20.9613i 0 +73.6194i 各條支