【正文】 =X2。 end end else %=============== 下面是針對PV節(jié)點(diǎn)來求取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。非對角元 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。對角元 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)。 % 非對角元規(guī)格化 end J(k,k)=1。 % 對角元規(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)號從小到大排列)：39。)。disp(E)。EE=E*UB。disp(EE)。disp(39。39。)。disp(39。各節(jié)點(diǎn)的電壓大小V為(節(jié)點(diǎn)號從小到大排列)：39。)。disp(V)。VV=V*UB。disp(VV)。disp(39。39。)。disp(39。各節(jié)點(diǎn)的電壓相角sida為(節(jié)點(diǎn)號從小到大排列)：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)號從小到大排列)：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。39。)。enddisp(39。各條支路的末端功率Sj為(順序同您輸入B1時(shí)一致)：39。)。for i=1:nl p=B1(i,1)。q=B1(i,2)。 if B1(i,6)==0 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))))。 Sjy(i)=Sj(q,p)。 else 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))))。 Sjy(i)=Sj(q,p)。 end disp(Sj(q,p))。 SSj(q,p)=Sj(q,p)*SB。 ZF=[39。S(39。,num2str(q),39。,39。,num2str(p),39。)=39。,num2str(SSj(q,p))]。 disp(ZF)。 %disp(SSj(q,p))。 disp(39。39。)。enddisp(39。各條支路的功率損耗DS為(順序同您輸入B1時(shí)一致)：39。)。for i=1:nl p=B1(i,1)。q=B1(i,2)。 DS(i)=Si(p,q)+Sj(q,p)。 disp(DS(i))。 DDS(i)=DS(i)*SB。 ZF=[39。DS(39。,num2str(p),39。,39。,num2str(q),39。)=39。,num2str(DDS(i))]。 disp(ZF)。 %disp(DDS(i))。 disp(39。39。)。endfigure(1)。subplot(2,2,1)。plot(V)。xlabel(39。節(jié)點(diǎn)號39。)。ylabel(39。電壓標(biāo)幺值39。)。grid on。subplot(2,2,2)。plot(sida)。xlabel(39。節(jié)點(diǎn)號39。)。ylabel(39。電壓角度39。)。grid on。subplot(2,2,3)。bar(real(S))。ylabel(39。節(jié)點(diǎn)注入有功39。)。grid on。subplot(2,2,4)。bar(Siz)。ylabel(39。支路首端無功39。)。grid on。