【正文】 ejBij*fj)eiΣ(Gij*fj+Bij*ej)求P39。,Q39。 V2=e(i)^2+f(i)^2。 電壓模平方 ========= 以下針對(duì)非PV節(jié)點(diǎn)來(lái)求取功率差及Jacobi矩陣元素 ========= if B2(i,6)~=3 非PV節(jié)點(diǎn) DP=P(i)P1。 節(jié)點(diǎn)有功功率差 DQ=Q(i)Q1。 節(jié)點(diǎn)無(wú)功功率差 =============== 以上為除平衡節(jié)點(diǎn)外其它節(jié)點(diǎn)的功率計(jì)算 ================================== 求取Jacobi矩陣 =================== for j1=1:n if j1~=isbamp。j1~=i 非平衡節(jié)點(diǎn)amp。非對(duì)角元 X1=G(i,j1)*e(i)B(i,j1)*f(i)。 % dP/de=dQ/df X2=B(i,j1)*e(i)G(i,j1)*f(i)。 % dP/df=dQ/de X3=X2。 % X2=dp/df X3=dQ/de X4=X1。 % X1=dP/de X4=dQ/df p=2*i1。q=2*j11。 J(p,q)=X3。J(p,N)=DQ。m=p+1。 J(m,q)=X1。J(m,N)=DP。q=q+1。 J(p,q)=X4。J(m,q)=X2。 elseif j1==iamp。j1~=isb 非平衡節(jié)點(diǎn)amp。對(duì)角元 X1=C(i)G(i,i)*e(i)B(i,i)*f(i)。 dP/de X2=D(i)+B(i,i)*e(i)G(i,i)*f(i)。 dP/df X3=D(i)+B(i,i)*e(i)G(i,i)*f(i)。 dQ/de X4=C(i)+G(i,i)*e(i)+B(i,i)*f(i)。 dQ/df p=2*i1。q=2*j11。J(p,q)=X3。J(p,N)=DQ。 擴(kuò)展列△Q m=p+1。 J(m,q)=X1。q=q+1。J(p,q)=X4。J(m,N)=DP。 擴(kuò)展列△P J(m,q)=X2。 end end else =============== 下面是針對(duì)PV節(jié)點(diǎn)來(lái)求取Jacobi矩陣的元素 =========== DP=P(i)P1。 PV節(jié)點(diǎn)有功誤差 DV=V(i)^2V2。 PV節(jié)點(diǎn)電壓誤差 for j1=1:n if j1~=isbamp。j1~=i 非平衡節(jié)點(diǎn)amp。非對(duì)角元 X1=G(i,j1)*e(i)B(i,j1)*f(i)。 dP/de X2=B(i,j1)*e(i)G(i,j1)*f(i)。 dP/df X5=0。X6=0。 p=2*i1。q=2*j11。J(p,q)=X5。J(p,N)=DV。 m=p+1。 J(m,q)=X1。J(m,N)=DP。q=q+1。J(p,q)=X6。 J(m,q)=X2。 elseif j1==iamp。j1~=isb 非平衡節(jié)點(diǎn)amp。對(duì)角元 X1=C(i)G(i,i)*e(i)B(i,i)*f(i)。 dP/de X2=D(i)+B(i,i)*e(i)G(i,i)*f(i)。 dP/df X5=2*e(i)。 X6=2*f(i)。 p=2*i1。q=2*j11。J(p,q)=X5。J(p,N)=DV。 m=p+1。 J(m,q)=X1。J(m,N)=DP。q=q+1。J(p,q)=X6。 J(m,q)=X2。 end end end end end========= 以上為求雅可比矩陣的各個(gè)元素 ===================== for k=3:N0 N0=2*n （從第三行開始，第一、二行是平衡節(jié)點(diǎn)） k1=k+1。N1=N。 N=N0+1 即 N=2*n+1擴(kuò)展列△P、△Q for k2=k1:N1 擴(kuò)展列△P、△Q J(k,k2)=J(k,k2)./J(k,k)。 非對(duì)角元規(guī)格化 end J(k,k)=1。 對(duì)角元規(guī)格化 if k~=3 不是第三行 ============================================================ k4=k1。 for k3=3:k4 用k3行從第三行開始到當(dāng)前行前的k4行消去 for k2=k1:N1 k3行后各行下三角元素 J(k3,k2)=J(k3,k2)J(k3,k)*J(k,k2)。%消去運(yùn)算 end J(k3,k)=0。 end if k==N0 break。 end========================================== for k3=k1:N0 for k2=k1:N1 J(k3,k2)=J(k3,k2)J(k3,k)*J(k,k2)。%消去運(yùn)算 end J(k3,k)=0。 end else for k3=k1:N0 for k2=k1:N1 J(k3,k2)=J(k3,k2)J(k3,k)*J(k,k2)。%消去運(yùn)算 end J(k3,k)=0。 end end end ====上面是用線性變換方式將Jacobi矩陣化成單位矩陣===== for k=3:2:N01 L=(k+1)./2。 e(L)=e(L)J(k,N)。 修改節(jié)點(diǎn)電壓實(shí)部 k1=k+1。 f(L)=f(L)J(k1,N)。 修改節(jié)點(diǎn)電壓虛部 end 修改節(jié)點(diǎn)電壓 for k=3:N0 DET=abs(J(k,N))。 if DET=pr 電壓偏差量是否滿足要求 IT2=IT2+1。 不滿足要求的節(jié)點(diǎn)數(shù)加1 end end ICT2(a)=IT2。 ICT1=ICT1+1。end用高斯消去法解w=J*Vdisp(39。迭代次數(shù)：39。)。disp(ICT1)。disp(39。沒有達(dá)到精度要求的個(gè)數(shù)：39。)。disp(ICT2)。for k=1:n V(k)=sqrt(e(k)^2+f(k)^2)。 sida(k)=atan(f(k)./e(k))*180./pi。 E(k)=e(k)+f(k)*j。end%=============== 計(jì)算各輸出量 ===========================disp(39。各節(jié)點(diǎn)的實(shí)際電壓標(biāo)幺值E為(節(jié)點(diǎn)號(hào)從小到大排列)：39。)。disp(E)。EE=E*UB。disp(EE)。disp(39。39。)。disp(39。各節(jié)點(diǎn)的電壓大小V為(節(jié)點(diǎn)號(hào)從小到大排列)：39。)。disp(V)。VV=V*UB。disp(VV)。disp(39。39。)。disp(39。各節(jié)點(diǎn)的電壓相角sida為(節(jié)點(diǎn)號(hào)從小到大排列)：39。)。disp(sida)。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)。enddisp(39。各節(jié)點(diǎn)的功率S為(節(jié)點(diǎn)號(hào)從小到大排列)：39。)。disp(S)。disp(39。~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~39。)。SS=S*SB。disp(SS)。disp(39。39。)。disp(39。各條支路的首端功率Si為(順序同您輸入B1時(shí)一致)：39。)。for i=1:nl p=B1(i,1)。q=B1(i,2)。 if B1(i,6)==0 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))))。 Siz(i)=Si(p,q)。 else 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))))。 Siz(i)=Si(p,q)。 end disp(Si(p,q))。 SSi(p,q)=Si(p,q)*SB。 ZF=[39。S(39。,num2str(p),39。,39。,num2str(q),39。)=39。,num2str(SSi(p,q))]。 disp(ZF)。 %disp(SSi(p,q))。 disp(39。