【文章內(nèi)容簡介】 元 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。 % PV 節(jié)點 電壓誤差 J(p,N)=DV。 m=p+1。 J(m,q)=X1。q=q+1。J(p,q)=X6。 % PV 節(jié)點有功誤差 J(m,N)=DP。 J(m,q)=X2。 end end end %(if B2(i,5)~=3 else) end %(if i~=isb) end %(for i=1:n)n 個節(jié)點 2n行 (每節(jié)點兩個方程 P和 Q或 U) JZ0=[39。形成的第 (39。,num2str(a),39。)次 Jacobi 矩陣： 39。]。 disp(JZ0)。disp(J)。 %=============================== 以 上 為 形 成 完 整 的 Jacobi 矩陣 ================================ %====下面用高斯消去法對由 Jacobi 矩陣形成的修正方程 進行求解 (按列消去、回代 ) ========== for k=3:N0 % N0=2*n （從第三行開始，第一、二行是平衡節(jié)點） for k1=k+1:N1 % 從 k+1 列的 Jacobi 元素到擴展列的△ P、△ Q 或 △ U J(k,k1)=J(k,k1)./J(k,k)。% 用 K行 K列對角元素去除 K行 K列后的非對角元素進行規(guī)格化 end J(k,k)=1。 % 對角元規(guī)格化 K行 K列對角元素賦 1 %================== 按列消去運算 ================================== for k2=k+1:N0 % 從 k+1 行到 2*n 最后一行 for k3=k+1:N1 %從 k2+1列到擴展列消去 k+1行后各行下三角元素 J(k2,k3)=J(k2,k3)J(k2,k)*J(k,k3)。%消去運算 end %用當前行 K3列元素減去當前行 K列元素乘以第 k行 K3列元素 J(k2,k)=0。 %當前行第 k 列元素已消為 0 end end JZ=[39。Jacobi 矩陣第 (39。,num2str(a),39。) 次消去運算 39。]。JZ1=[39。Jacobi 矩陣第(39。,num2str(a),39。)次回代運算 39。]。 disp(JZ)。disp(J)。 %==================== 按列回代運算 ======================================= for k=N0:1:3 for k1=k1:1:3 J(k1,N1)=J(k1,N1)J(k1,k)*J(k,N1)。 J(k1,k)=0。 end end for m=1:N0 JJN1(m)=J(m,N1)。 end disp(JZ1)。disp(JJN1)。%disp(J)。 %修改節(jié)點電壓 for k=3:2:N01 L=(k+1)./2。 e(L)=e(L)J(k,N1)。 %修改節(jié)點電壓實部 k1=k+1。 f(L)=f(L)J(k1,N1)。 %修改節(jié)點電壓虛 部 U(L)=sqrt(e(L)^2+f(L)^2)。 end disp(39。各個節(jié)點電壓模 39。)。disp(U)。 %============================== 結(jié)束一次迭代 ============================== end %********************** 下面為迭代計算結(jié)束后的有關(guān)輸出過程 ***************** disp(39。迭代次數(shù)： 39。)。 disp(ICT11)。 XV disp(39。沒有達到 精度要求的個數(shù)： 39。)。 disp(ICT2)。 for k=1:n V(k)=sqrt(e(k)^2+f(k)^2)。 %計算各節(jié)點電壓的模值 sida(k)=atan(f(k)./e(k))*180./pi。 %計算各節(jié)點電壓的角度 E(k)=e(k)+f(k)*j。 %將各節(jié)點電壓用復(fù)數(shù)表示 end %=============== 計算各輸出量 =========================== disp(39。各節(jié)點的電壓復(fù)數(shù)值 E為 (節(jié)點號從小到大 排列 )： 39。)。 disp(E)。 %顯示各節(jié)點的實際電壓值 E用復(fù)數(shù)表示 disp(39。39。)。 disp(39。各節(jié)點的電壓模值大小 V為 (節(jié)點號從小到大排列 )： 39。)。 disp(V)。 %顯示各節(jié)點的電壓大小 V的模值 disp(39。39。)。 disp(39。各節(jié)點的電壓相角 sida為 (節(jié)點號從小到大排列 )： 39。)。 disp(sida)。 %顯示各節(jié)點的電壓相角 for p=1:n C(p)=0。 for q=1:n C(p)=C(p)+conj(Y(p,q))*conj(E(q))。 %計算各節(jié)點注入電流的共軛值 end S(p)=E(p)*C(p)。 %計算各節(jié)點的功率 S = 電壓 X 注入電流的共軛值 end disp(39。各節(jié)點的功率 S為 (節(jié)點號從小到大排列 )： 39。)。 disp(S)。 %顯示各節(jié)點的注入功率 disp(39。39。)。 disp(39。各條支路的首端功率 Si為 (順序同您輸入 B1時一致 )： 39。)。 for i=1:nl p=B1(i,1)。q=B1(i,2)。 if B1(i,7)==0 Si(p,q)=E(p)*conj(E(p)*B1(i,4)./2+(E(p)E(q))./B1(i,3))。 Siz(i)=Si(p,q)。 else if B1(i,6)==0 Si(p,q)=E(p)*(conj(E(p)*B1(i,4)... +(E(p)*B1(i,5)E(q))*(1./(B1(i,3)*B1(i,5)))))。 Siz(i)=Si(p,q)。 else Si(p,q)=E(p)*conj((E(p)E(q)*B1(i,5))*(1./(B1(i,3)*B1(i,5)^2)))。 Siz(i)=Si(p,q)。 end end ZF=[39。S(39。,num2str(p),39。,39。,num2str(q),39。)=39。,num2str(Si(p,q))]。 disp(ZF)。 disp(39。39。)。 end disp(39。各條支路的末端功率 Sj為 (順序同您輸入 B1時一致 )： 39。)。 for i=1:nl p=B1(i,1)。q=B1(i,2)。 if B1(i,7)==0 Sj(q,p)=E(q)*conj(E(q)*B1(i,4)./2+(E(q)E(p))./B1(i,3))。 Sjy(i)=Sj(q,p)。 else if B1(i,6)==0 Sj(q,p)=E(q)*conj((E(q)E(p)*B1(i,5))*(1./(B1(i,3)*B1(i,5)^2)))。 Sjy(i)=Sj(q,p)。 else Sj(q,p)=E(q)*(conj(E(q)*B1(i,4)... +(E(q)*B1(i,5)E(p))*(1./(B1(i,3)*B1(i,5)))))。 Sjy(i)=Sj(q,p)。 end end ZF=[39。S(39。,num2str(q),39。,39。,num2str(p),39。)=39。,num2str(Sj(q,p))]。 disp(ZF)。 disp(39。39。)。 end disp(39。各條支路的功率損耗 DS為 (順序同您輸入 B1時一致 )： 39。)。 for i=1:nl XVII p=B1(i,1)。q=B1(i,2)。 DS(i)=Si(p,q)+Sj(q,p)。 ZF=[39。DS(39。,num2str(p),39。,39。,num2str(q),39。)=39。,num2str(DS(i))]。 disp(ZF)。 disp(39。39。)。 end zws=0。JDDY=[]。JDP=[]。JDQ=[]。JDDYJD=[]。 for i=1:n %總網(wǎng)損為所有節(jié)點注入功率的代數(shù)和 zws=zws+S(i)。 JDDYJD=strcat(JDDYJD,num2str(i),39。(39。,num2str(sida(i)),39。),39。)。 JDDY=strcat(JDDY,num2str(i),39。(39。,num2str(V(i)),39。),39。)。 JDP=strcat(JDP,num2str(i),39。(39。,num2str(real(S(i))),39。),39。)。 JDQ=strcat(JDQ,num2str(i),39。(39。,num2str(imag(S(i))),39。),39。)。 end JDDYJD=strcat(39。節(jié)點電壓角度： 39。,JDDYJD)。 JDDY=strcat(39。節(jié)點電壓模值： 39。,JDDY)。 JDP=strcat(39。節(jié)點有功： 39。,JDP)。 JDQ=strcat(39。節(jié)點無功： 39。,JDQ)。 ZF subplot(4,1,2)。 plot(sida)。 title(JDDYJD)。 ylabel(39。節(jié)點電壓角 度 39。)。 grid on。 subplot(4,1,3)。 P=real(S)。Q=imag(S)。 bar(P)。 title(JDP)。 ylabel(39。節(jié)點注入有功 39。)。 grid on。 subplot(4,1,4)。 bar(Q)。 title(JDQ)。 xlabel(ZF)。ylabel(39。節(jié)點注入無功 39。)。 grid on。