【正文】 m=p+1。 J(m,q)=X3。 J(m,N)=DQ。 q=q+1。 J(p,q)=X2。 J(m,q)=X4。 %非平衡節(jié)點(diǎn)時(shí)對角線元素 elseif k~=ph amp。 k==i X1=D(i)+B(i,i)*e(i)G(i,i)*f(i)。 X2=C(i)G(i,i)*e(i)B(i,i)*f(i)。 X3=C(i)+G(i,i)*e(i)+B(i,i)*f(i)。 X4=D(i)+B(i,i)*e(i)G(i,i)*f(i)。 p=2*i1。 q=2*k1。 J(p,q)=X1。 J(p,N)=DP。 m=p+1。 J(m,q)=X3。 J(m,N)=DQ。 q=q+1。 J(p,q)=X2。 J(m,q)=X4。 end end else DP=P(i)P1。 DV=V(i)^2V2。 for s=1:n if s~=ph amp。 s~=i X1=G(i,s)*f(i)+B(i,s)*e(i)。 X2=G(i,s)*e(i)B(i,s)*f(i)。 X5=0。 X6=0。 p=2*i1。 q=2*s1。 J(p,q)=X1。 J(p,N)=DP。 m=p+1。 J(m,q)=X5。 J(m,N)=DV。 q=q+1。 J(p,q)=X2。 J(m,q)=X6。 elseif s~=ph amp。 s==i X1=D(i)+B(i,i)*e(i)G(i,i)*f(i)。 X2=C(i)G(i,s)*e(i)B(i,s)*f(i)。 X5=2*f(i)。 X6=2*e(i)。 p=2*i1。 q=2*s1。 J(p,q)=X1。 J(p,N)=DP。 m=p+1。 J(m,q)=X5。 J(m,N)=DV。 q=q+1。 J(p,q)=X2。 J(m,q)=X6。 end end end end end disp(39。系統(tǒng)雅克比矩陣為J=39。)。 disp(J)。%對于高壓電力系統(tǒng)，最大元素出現(xiàn)在對角線位置，簡單計(jì)算如下 DJ=J(:,[1:NO2])。 DP=J([1:NO2],N)。 DU=DJ\DP。 disp(39。************************************39。)。 disp(39。 39。) disp(39。第幾次的修正值DU？39。)。 disp(Ci)。 disp(39。***********************************39。)。 disp(DU)。 for i=1:NO2 eps=abs(DP(i))。 if eps=jd。 a=a+1。 end end Ci=Ci+1。 for i=1:n1 e(i)=e(i)DU(2*i)。 f(i)=f(i)DU(2*i1)。 end E=e+j*f。 disp(39。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。39。)。 disp(39。 39。)。 disp(39。第幾次迭代后的近似計(jì)算值39。)。 disp(Ci)。 disp(39。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。amp。39。)。 for i1=1:n V(i1)=sqrt(e(i1)^2+f(i1)^2)。 ANGLE(i1)=atan(f(i1)./e(i1))*(180./pi)。 end%下列輸出結(jié)果均按節(jié)點(diǎn)號(hào)從小到大排列 disp(39。各個(gè)節(jié)點(diǎn)的電壓實(shí)際值為：E=39。)。 disp(E)。 disp(39。各個(gè)節(jié)點(diǎn)電壓幅值為：V=39。)。 disp(V)。 disp(39。各個(gè)節(jié)點(diǎn)的電壓相交為：ANGLE=39。)。 disp(ANGLE)。 disp(39。各個(gè)節(jié)點(diǎn)的功率誤差變化39。)。 disp(J([1:8],9)39。)。 MDU(Ci)=max(abs(J([1:8],9)39。))。 for p=1:n Sm(p)=0。 for q=1:n Sm(p)=Sm(p)+conj(Y(p,q))*conj(E(q))。 end S(i)=E(p)*Sm(p)。 endenddisp(39。@ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @39。)。disp(39。 39。)。disp(39。計(jì)算的最終功率近似值如下各項(xiàng)值39。)。disp(39。 39。)。disp(39。@ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @39。)。disp(39。各個(gè)節(jié)點(diǎn)功率為：S=39。)。disp(S)。disp(39。各條支路的首端功率Si為：39。)。for i=1:nl p=Br(i,1)。q=Br(i,2)。 A(i)=E(p)*[conj(E(p))conj(E(q))]*conj(1./(Br(i,3)))。 SBS(i)=[abs(E(p))]^2*conj(Br(i,4))+A(i)。 enddisp(SBS)。disp(39。各條支路的末端功率Sj為:39。)。for i=1:nl p=Br(i,1)。q=Br(i,2)。 A2=E(q)*[conj(E(q))conj(E(p))]*conj(1./(Br(i,3)))。 SBE(i)=[abs(E(q))]^2*conj(Br(i,4))+A2。enddisp(SBE)。disp(39。各條支路的功率損耗為:39。)。for i=1:nl p=Br(i,1)。q=Br(i,2)。 SBL(i)=SBS(i)+SBE(i)。enddisp(SBL)。%繪制功率誤差曲線，先用最小二乘法進(jìn)行曲線擬合for i=1:Ci CSH(i)=i。endP=polyfit(CSH,MDU,3)。plot(CSH,MDU,39。O39。)。節(jié)點(diǎn)導(dǎo)納矩陣為：Y= Columns 1 through 3 + + + + + + 0 0 + 0 0 0 0 + Columns 4 through 5 0 0 0 + 0 0 0 + 0 0 0 0 系統(tǒng)雅克比矩陣為：J= Columns 1 through 7 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Columns 8 through 9 0 0 0 0 0 0 0************************************