【正文】 )。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,5)^2)+B1(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)。end %求導(dǎo)納矩陣G=real(Y)。B=imag(YI)。BI=imag(Y)。fori=1:n S(i)=B2(i,1)B2(i,2)。BI(i,i)=BI(i,i)+B2(i,5)。endP=real(S)。Q=imag(S)。fori=1:n e(i)=real(B2(i,3))。f(i)=imag(B2(i,3))。V(i)=B2(i,4)。endfori=1:nif B2(i,6)==2 V(i)=sqrt(e(i)^2+f(i)^2)。 O(i)=atan(f(i)./e(i))。 endendfori=2:n ifi==n B(i,i)=1./B(i,i)。else IC1=i+1。 for j1=IC1:n B(i,j1)=B(i,j1)./B(i,i)。endB(i,i)=1./B(i,i)。 for k=i+1:n for j1=i+1:n B(k,j1)=B(k,j1)B(k,i)*B(i,j1)。endendendendp=0。q=0。fori=1:nif B2(i,6)==2 p=p+1。k=0。for j1=1:nif B2(j1,6)==2 k=k+1。 A(p,k)=BI(i,j1)。endendendendfori=1:na ifi==naA(i,i)=1./A(i,i)。 else k=i+1。for 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~=0 ICT2=0。ICT3=0。 fori=1:nifi~=isbC(i)=0。for k=1:n C(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)./V(i)。 DET=abs(DP1(i))。if DET=pr ICT2=ICT2+1。endendendNp(K)=ICT2。if ICT2~=0fori=2:nDP(i)=B(i,i)*DP(i)。ifi~=n IC1=i+1。for k=IC1:nDP(k)=DP(k)B(k,i)*DP(i)。 endelsefor LZ=3:i L=i+3LZ。 IC4=L1。for MZ=2:IC4 I=IC4+2MZ。DP(I)=DP(I)B(I,L)*DP(L)。endendendendfori=2:n O(i)=O(i)DP(i)。endkq=1。L=0。fori=1:nif B2(i,6)==2C(i)=0。L=L+1。for k=1:n C(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 =pr ICT3=ICT3+1。endendendelsekp=0。 ifkq~=0。 L=0。fori=1:n if B2(i,6)==2C(i)=0。L=L+1。for k=1:n C(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))。endendendendNq(K)=ICT3。if ICT3~=0 L=0。fori=1:naDQ(i)=A(i,i)*DQ(i)。ifi==nafor LZ=2:i L=i+2LZ。 IC4=L1。for MZ=1:IC4 I=IC4+1MZ。DQ(I)=DQ(I)A(I,L)*DQ(L)。endendelse IC1=i+1。for k=IC1:naDQ(k)=DQ(k)A(k,i)*DQ(i)。endendend L=0。fori=1:nif B2(i,6)==2 L=L+1。V(i)=V(i)DQ(L)。endendkp=1。 K=K+1。elsekq=0。ifkp~=0 K=K+1。endendfori=1:nDy(K1,i)=V(i)。endenddisp(39。迭代次數(shù)39。)。disp(K)。 disp(39。每次沒有達到精度要求的有功功率個數(shù)為39。)。disp(Np)。disp(39。每次沒有達到精度要求的無功功率個數(shù)為39。)。disp(Nq)。for k=1:nE(k)=V(k)*cos(O(k))+V(k)*sin(O(k))*j。O(k)=O(k)*180./pi。enddisp(39。各節(jié)點的電壓標(biāo)幺值E為（節(jié)點號從小到大排）：39。)。disp(E)。disp(39。各節(jié)點的電壓V大?。ü?jié)點號從小到大排）為：39。)。disp(V)。disp(39。各節(jié)點的電壓相角O（節(jié)點號從小到大排）為：39。)。disp(O)。for p=1:nC(p)=0。for q=1:n C(p)=C(p)+conj(Y(p,q))*conj(E(q))。endS(p)=E(p)*C(p)。enddisp(39。各節(jié)點的功率S（節(jié)點號從小到大排）為：39。)。disp(S)。disp(39。各條支路的首端功率Sj（順序同您輸入B1時一樣）為：39。)。 fori=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,5))conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))))。disp(Si(p,q))。enddisp(39。各條支路的末端功率Sj（順序同您輸入B1時一樣）為：39。)。 fori=1:nlif 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,5))conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))))。disp(Sj(q,p))。 enddisp(39。各條支路的功率損耗DS（順序同您輸入B1時一樣）為：39。)。 fori=1:nl if B1(i,6)==0 p=B1(i,1)。q=B1(i,2)。else p=B1(i,2)。q=B1(i,1)。 endDS(i)=Si(p,q)+Sj(q,p)。 disp(DS(i))。endfori=1:K Cs(i)=i。 for j=1:n Dy(K,j)=Dy(K1,j)。endenddisp(39。以下是每次迭代后各節(jié)點的電壓值（如圖所示）39。)。 plot(Cs,Dy),xlabel(39。迭代次數(shù)39。),ylabel(39。電壓39。),title(39。電壓迭代次數(shù)曲線39。)。47