【正文】 OrgS=zeros(2*n2,1)。 C=zeros(n,1)。 D=zeros(n,1)。 h=0。 %P的處理 for s=1:n if s~=isbamp。B2(s,3)==2 h=h+1。 for j=1:n C(s,1)=C(s,1)+G(s,j)*e(j)B(s,j)*f(j)。 D(s,1)=D(s,1)+G(s,j)*f(j)+B(s,j)*e(j)。 end OrgS(h,1)=C(s,1)*e(s)+f(s)*D(s,1)。 end end for s=1:n if s~=isbamp。B2(s,3)==3 h=h+1。 for j=1:n C(s,1)=C(s,1)+G(s,j)*e(j)B(s,j)*f(j)。 D(s,1)=D(s,1)+G(s,j)*f(j)+B(s,j)*e(j)。 end OrgS(h,1)=C(s,1)*e(s)+f(s)*D(s,1)。 end end %Q的處理 for s=1:n if s~=isbamp。B2(s,3)==2 h=h+1。 OrgS(h,1)=f(s)*C(s,1)D(s,1)*e(s)。 end end %V的處理 for s=1:n if s~=isbamp。B2(s,3)==3 h=h+1。 OrgS(h,1)=B2(s,4)。 end end %創(chuàng)建DetaS,用于存儲不平衡量 DetaS=zeros(2*n2,1)。 h=0。 r=0。 %對P的處理 for s=1:n if s~=isbamp。B2(s,3)==2 h=h+1。 r=r+1。 DetaS(h,1)=real(B2(s,1))OrgS(h,1)。 end end for s=1:n if s~=isbamp。B2(s,3)==3 h=h+1。 r=r+1。 DetaS(h,1)=real(B2(s,1))OrgS(h,1)。 end end %對Q的處理 for s=1:n if s~=isbamp。B2(s,3)==2 h=h+1。 DetaS(h,1)=imag(B2(s,1))OrgS(h,1)。 end end %對V的處理 for s=1:n if s~=isbamp。B2(s,3)==3 h=h+1。 DetaS(h,1)=real(OrgS(h))^2+imag(OrgS(h))^2real(B2(s,2))^2imag(B2(s,2))^2。 end end %創(chuàng)建雅克比矩陣Jacobi Jacobi=zeros(2*n2)。 h=0。 k=0。 %對P的處理 for i=1:n if B2(i,3)==2 h=h+1。 for j=1:n if j~=isb k=k+1。 if i==j Jacobi(h,k)=G(i,j)*e(i)+B(i,j)*f(i)C(h,1)。 Jacobi(h,r+k)=G(i,j)*f(i)+B(i,j)*e(i)D(h,1)。 Jacobi(r+h,k)=G(i,j)*f(i)+B(i,j)*e(i)+D(h,1)。 Jacobi(r+h,r+k)=G(i,j)*e(i)+B(i,j)*f(i)C(h,1)。 else Jacobi(h,k)=G(i,j)*e(i)B(i,j)*f(i)。 Jacobi(h,r+k)=G(i,j)*f(i)+B(i,j)*e(i)。 Jacobi(r+h,k)=G(i,j)*f(i)+B(i,j)*e(i)。 Jacobi(r+h,r+k)=G(i,j)*e(i)+B(i,j)*f(i)。 end if k==(n1) k=0。 end end end end end k=0。 for i=1:n if B2(i,3)==3 h=h+1。 for j=1:n if j~=isb k=k+1。 if i==j Jacobi(h,k)=G(i,j)*e(i)+B(i,j)*f(i)C(h,1)。 %H(i,i) Jacobi(h,r+k)=G(i,j)*f(i)+B(i,j)*e(i)D(h,1)。 Jacobi(r+h,k)=2*e(h)。 Jacobi(r+h,r+k)=2*f(h)。 else Jacobi(h,k)=G(i,j)*e(i)B(i,j)*f(i)。 Jacobi(h,r+k)=G(i,j)*f(i)+B(i,j)*e(i)。 Jacobi(r+h,k)=0。 Jacobi(r+h,r+k)=0。 end if k==(n1) k=0。 end end end end end Jacobi DetaU=zeros(2*n2,1)。 %電壓偏移量 U=zeros(1,n)。 DetaU=inv(Jacobi)*DetaS。 DetaU。 %修正節(jié)點電壓 h=0。 for j=1:n if j~=isb h=h+1。 B2(j,2)=B2(j,2)(DetaU(h,1)+DetaU(r+h,1)*sqrt(1))。 end end %節(jié)點電壓 for j=1:n U(1,j)=B2(j,2)。 end disp(39。節(jié)點電壓U39。)。 U end%計算支路節(jié)點功率y=zeros(n)。for s=1:m if B1(s,6)==0 %不含變壓器的支路 p=B1(s,1)。 q=B1(s,2)。 y(p,q)=B1(s,4)。 y(q,p)=y(p,q)。 else %含有變壓器的支路 p=B1(s,1)。 q=B1(s,2)。 y(p,q)=(B1(s,5)^2B1(s,5))/B1(s,3)。 y(q,p)=(1B1(s,5))/B1(s,3)。 end endy。sum=0。for j=1:n sum=sum+conj(Y(m,j))*conj(U(1,j))。endS=zeros(n)。for s=1:n if B1(s,6)==0 %不含變壓器的支路 p=B1(s,1)。 q=B1(s,2)。 S(p,q)=U(1,p)*(conj(U(1,p))*conj(y(p,q))+(conj(U(1,p))conj(U(1,q)))*conj(inv(B1(s,3))))。 S(q,p)=U(1,q)*(conj(U(1,q))*conj(y(q,p))+(conj(U(1,q))conj(U(1,p)))*conj(inv(B1(s,3))))。 else p=B1(s,1)。 q=B1(s,2)。 S(p,q)=U(1,p)*(conj(U(1,p))*conj(y(p,q))+(conj(U(1,p))conj(U(1,q)))*conj(inv(B1(s,3)/B1(s,5))))。 S(q,p)=U(1,q)*(conj(U(1,q))*conj(y(q,p))+(conj(U(1,q))conj(U(1,p)))*conj(inv(B1(s,3)/B1(s,5))))。 endenddisp(39。各條支路的功率為：39。)。S袁節(jié)膅薂羄肅蒃薁蚃芀荿薀螆肅芅蕿袈羋膁蚈羀肁蒀蚇蝕襖莆蚇螂肀莂蚆羅袂羋蚅蚄膈膄蚄螇羈蒂蚃衿膆莈螞羈罿芄螁蟻膄膀螁螃羇葿螀裊膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃螞肂莈蒂螄羋芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羈莀蒈羃膇芆蕆蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃蠆羆艿薃袁節(jié)膅薂羄肅蒃薁蚃芀荿薀螆肅芅蕿袈羋膁蚈羀肁蒀蚇蝕襖莆蚇螂肀莂蚆羅袂羋蚅蚄膈膄蚄螇羈蒂蚃衿膆莈螞羈罿芄螁蟻膄膀螁螃羇葿螀裊膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃螞肂莈蒂螄羋芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羈莀蒈羃膇芆蕆蚃羀