【正文】 elseif(i=m+1 amp。amp。 i=n)%如果i號(hào)節(jié)點(diǎn)是PV節(jié)點(diǎn) dPQV2(2*(i1)1)=P_PV(im)real(A)。%PV節(jié)點(diǎn)有功功率不平衡量 dPQV2(2*(i1))=U_PV(im)^2abs(U(i))^2。%PV節(jié)點(diǎn)電壓平方不平衡量 end end end i=1。 %節(jié)點(diǎn)的不平衡量的dPQV2計(jì)算結(jié)束 for i=1:n%計(jì)算節(jié)點(diǎn)注入電流I； if(i=s1) for a=1:n I(i)=I(i)+Y(i,a)*U(a)。 end a=1。 elseif(i=s+1) for a=1:n I(i)=I(i)+Y(i,a)*U(a)。 end a=1。 end end i=1。 for i=1:n %計(jì)算雅克比矩陣 if(i=s1)%如果i號(hào)節(jié)點(diǎn)是平衡節(jié)點(diǎn)之前的PQ節(jié)點(diǎn) for a=1:n if(a~=i) if(a=s1)%如果a號(hào)節(jié)點(diǎn)是平衡節(jié)點(diǎn)之前的PQ節(jié)點(diǎn) Ykb(2*i1,2*a1)=imag(Y(i,a))*real(U(i))real(Y(i,a))*imag(U(i))。 %PQ節(jié)點(diǎn)Hia參數(shù) Ykb(2*i1,2*a)=real(Y(i,a))*real(U(i))imag(Y(i,a))*imag(U(i))。 %PQ節(jié)點(diǎn)Nia參數(shù)Ykb(2*i,2*a1)=Ykb(2*i1,2*a)。%PQ節(jié)點(diǎn)的Jia參數(shù)=NiaYkb(2*i,2*a)=Ykb(2*i1,2*a1)。%PQ節(jié)點(diǎn)的Lia參數(shù)=Hia else if(a=s+1 amp。amp。 a=m)%如果a號(hào)節(jié)點(diǎn)是平衡節(jié)點(diǎn)之后的PQ節(jié)點(diǎn),Ykb中的aa1。i不變 Ykb(2*i1,2*(a1)1)=imag(Y(i,a))*real(U(i))real(Y(i,a))*imag(U(i))。%PQ節(jié)點(diǎn)Hia參數(shù) Ykb(2*i1,2*(a1))=real(Y(i,a))*real(U(i))imag(Y(i,a))*imag(U(i))。%PQ節(jié)點(diǎn)Nia參數(shù) Ykb(2*i,2*(a1)1)=Ykb(2*i1,2*(a1))。%PQ節(jié)點(diǎn)的Jia參數(shù)=Nia Ykb(2*i,2*(a1))=Ykb(2*i1,2*(a1)1)。%PQ節(jié)點(diǎn)的Lia參數(shù)=Hia elseif(a=m+1 amp。amp。 a=n)%Ykb中的aa1。i不變 Ykb(2*i1,2*(a1)1)=imag(Y(i,a))*real(U(i))real(Y(i,a))*imag(U(i))。%PQ節(jié)點(diǎn)Hia參數(shù) Ykb(2*i1,2*(a1))=real(Y(i,a))*real(U(i))imag(Y(i,a))*imag(U(i))。%PQ節(jié)點(diǎn)Nia參數(shù)Ykb(2*i,2*(a1)1)=Ykb(2*i1,2*(a1))。%PV節(jié)點(diǎn)Jia參數(shù)=NiaYkb(2*i1,2*(a1))=Ykb(2*i1,2*(a1)1)。%PV節(jié)點(diǎn)Lia參數(shù)=Hia end end elseif(a==i) if(a=s1)%如果a號(hào)節(jié)點(diǎn)是平衡節(jié)點(diǎn)之前的PQ節(jié)點(diǎn) Ykb(2*i1,2*a1)=imag(Y(i,a))*real(U(i))real(Y(i,a))*imag(U(i))imag(I(i))。%PQ節(jié)點(diǎn)Hia參數(shù) Ykb(2*i1,2*a)=real(Y(i,a))*real(U(i))imag(Y(i,a))*imag(U(i))real(I(i))。%PQ節(jié)點(diǎn)Nia參數(shù)Ykb(2*i,2*a1)=Ykb(2*i1,2*a)2*real(I(i))。%PQ節(jié)點(diǎn)的Jia參數(shù)Ykb(2*i,2*a)=Ykb(2*i1,2*a1)+2*imag(I(i))。%PQ節(jié)點(diǎn)的Lia參數(shù) else if(a=s+1 amp。amp。 a=m)%如果a號(hào)節(jié)點(diǎn)是平衡節(jié)點(diǎn)之后的PQ節(jié)點(diǎn) Ykb(2*i1,2*(a1)1)=imag(Y(i,a))*real(U(i))real(Y(i,a))*imag(U(i))imag(I(i))。%PQ節(jié)點(diǎn)Hia參數(shù) Ykb(2*i1,2*(a1))=real(Y(i,a))*real(U(i))imag(Y(i,a))*imag(U(i))real(I(i))。%PQ節(jié)點(diǎn)Nia參數(shù) Ykb(2*i,2*(a1)1)=Ykb(2*i1,2*(a1))2*real(I(i))。%PQ節(jié)點(diǎn)的Jia參 數(shù)=Nia2real（I(i)）； Ykb(2*i,2*(a1))=Ykb(2*i1,2*(a1)1)+2*imag(I(i))。%PQ節(jié)點(diǎn)的Lia參數(shù)=Hia+2imag(I(i))。 elseif(a=m+1 amp。amp。 a=n) Ykb(2*i1,2*(a1)1)=imag(Y(i,a))*real(U(i))real(Y(i,a))*imag(U(i))imag(I(i))。%PQ節(jié)點(diǎn)Hia參數(shù) Ykb(2*i1,2*(a1))=real(Y(i,a))*real(U(i))imag(Y(i,a))*imag(U(i))real(I(i))。%PQ節(jié)點(diǎn)Nia參數(shù) Ykb(2*i,2*(a1)1)=Ykb(2*i1,2*(a1))2*real(I(i))。%PQ節(jié)點(diǎn)的Jia參數(shù)=Nia2real（I(i)）； Ykb(2*i1,2*(a1))=Ykb(2*i1,2*(a1)1)+2*imag(I(i))。%PQ節(jié)點(diǎn)的Lia參數(shù)=Hia+2imag(I(i))。 end end end end a=1。 else if(i=s+1 amp。amp。 i=m)%如果i號(hào)節(jié)點(diǎn)是平衡節(jié)點(diǎn)之后的PQ節(jié)點(diǎn)。Ybk中的ii1。 for a=1:n if(a~=i) if(a=s1)%如果a號(hào)節(jié)點(diǎn)是平衡節(jié)點(diǎn)之前的PQ節(jié)點(diǎn)。Ybk中的ii1。a不變 Ykb(2*(i1)1,2*a1)=imag(Y(i,a))*real(U(i))real(Y(i,a))*imag(U(i))。%PQ節(jié)點(diǎn)Hia參數(shù) Ykb(2*(i1)1,2*a)=real(Y(i,a))*real(U(i))imag(Y(i,a))*imag(U(i))。%PQ節(jié)點(diǎn)Nia參數(shù) Ykb(2*(i1),2*a1)=Ykb(2*(i1)1,2*a)。%PQ節(jié)點(diǎn)的Jia參數(shù)=Nia Ykb(2*(i1),2*a)=Ykb(2*(i1)1,2*a1)。%PQ節(jié)點(diǎn)的Lia參數(shù)Hia else if(a=s+1 amp。amp。 a=m)%如果a號(hào)節(jié)點(diǎn)是平衡節(jié)點(diǎn)之后的PQ節(jié)點(diǎn)。Ybk中的ii1。aa1。 Ykb(2*(i1)1,2*(a1)1)=imag(Y(i,a))*real(U(i))real(Y(i,a))*imag(U(i))。%PQ節(jié)點(diǎn)Hia參數(shù) Ykb(2*(i1)1,2*(a1))=real(Y(i,a))*real(U(i))imag(Y(i,a))*imag(U(i))。%PQ節(jié)點(diǎn)Nia參數(shù) Ykb(2*(i1),2*(a1)1)=Ykb(2*(i1)1,2*(a1))。%PQ節(jié)點(diǎn)的Jia參數(shù)=Nia Ykb(2*(i1),2*(a1))=Ykb(2*(i1)1,2*(a1)1)。%PQ節(jié)點(diǎn)的Lia參數(shù)=Hia elseif(a=m+1 amp。amp。 a=n) Ykb(2*(i1)1,2*(a1)1)=imag(Y(i,a))*real(U(i))real(Y(i,a))*imag(U(i))。%PQ節(jié)點(diǎn)Hia參數(shù) Ykb(2*(i1)1,2*(a1))=real(Y(i,a))*real(U(i))imag(Y(i,a))*imag(U(i))。%PQ節(jié)點(diǎn)Nia參數(shù) Ykb(2*(i1),2*(a1)1)=Ykb(2*(i1)1,2*(a1))。%PQ節(jié)點(diǎn)的Jia參數(shù)=Nia Ykb(2*(i1)1,2*(a1))=Ykb(2*(i1)1,2*(a1)1)。%PQ節(jié)點(diǎn)的Lia參數(shù)=Hia end end elseif(a==i) if(a=s1)%如果a號(hào)節(jié)點(diǎn)是平衡節(jié)點(diǎn)之前的PQ節(jié)點(diǎn)。Ybk中的ii1。a不變 Ykb(2*(i1)1,2*a1)=imag(Y(i,a))*real(U(i))real(Y(i,a))*imag(U(i))imag(I(i))。%PQ節(jié)點(diǎn)Hia參數(shù) Ykb(2*(i1)1,2*a)=real(Y(i,a))*real(U(i))imag(Y(i,a))*imag(U(i))real(I(i))。%PQ節(jié)點(diǎn)Nia參數(shù) Ykb(2*(i1),2*a1)=Ykb(2*(i1)1,2*a)。%PQ節(jié)點(diǎn)的Jia參數(shù) Ykb(2*(i1),2*a)=Ykb(2*(i1)1,2*a1)+2*imag(I(i))。%PQ節(jié)點(diǎn)的Lia參數(shù) else if(a=s+1 amp。amp。 a=m)%如果a號(hào)節(jié)點(diǎn)是平衡節(jié)點(diǎn)之后的PQ節(jié)點(diǎn) Ykb(2*(i1)1,2*(a1)1)=imag(Y(i,a))*real(U(i))real(Y(i,a))*imag(U(i))imag(I(i))。%PQ節(jié)點(diǎn)Hia參數(shù) Ykb(2*(i1)1,2*(a1))=real(Y(i,a))*real(U(i))imag(Y(i,a))*imag(U(i))real(I(i))。%PQ節(jié)點(diǎn)Nia參數(shù) Ykb(2*(i1),2*(a1)1)=Ykb(2*(i1)1,2*(a1))2*real(I(i))。%PQ節(jié)點(diǎn)的Jia參數(shù)=Nia2real(I(i)) Ykb(2*(i1),2*(a1))=Ykb(2*(i1)1,2*(a1)1)+2*imag(I(i))。%PQ節(jié)點(diǎn)的Lia參數(shù)=Hia+2*imag(I(i)) elseif(a=m+1 amp。amp。 a=