【導讀】(三)每個變電所的功率因數(shù)均為cosφ=；(四)變電所1和變電所2分別配有兩臺容量為75MVA的變壓器，短路損耗414KW，，單位長度的電納為S-610*。流計算，并對計算結果進行分析。以2%的比例上升；（變電所低壓母線電壓10KV要求調整范圍在之間；電壓35KV. 節(jié)控制的參數(shù)變化。潮流計算是電力系統(tǒng)最基本最常用的計算。根據(jù)系統(tǒng)給定的運行條件，網絡。接線及元件參數(shù)，通過潮流計算可以確定各母線的電壓，包括電壓的幅值和相角，將N-R法用于潮流計算是以導納矩陣為基礎的，由于利用了導。用內存，計算速度等方面的優(yōu)點都超過了阻抗法。利用成形的程序對系統(tǒng)中出現(xiàn)的各種情況，例如負荷的變化以及。行合理的系統(tǒng)潮流。由調壓范圍可知，變電所1、2、3、4的節(jié)點電壓均不滿足電壓范圍。

　　

【正文】 =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é)點 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 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。%擴展列△ Q m=p+1。 J(m,q)=X1。q=q+1。J(p,q)=X4。J(m,N)=DP。%擴展列△ P J(m,q)=X2。 end end else %=============== 下面是針對 PV節(jié)點來求取 Jacobi 矩陣的元素 =========== DP=P(i)P1。 % PV 節(jié)點有功誤差 DV=V(i)^2V2。 % PV 節(jié)點電壓誤差 for j1=1:n if j1~=isbamp。j1~=i %非平衡節(jié)點 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é)點 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 %========= 以上為求雅可比矩陣的各個元素 ===================== for k=3:N0 % N0=2*n （從第三行開始，第一、二行是平衡節(jié)點） k1=k+1。N1=N。 % N=N0+1 即 N=2*n+1 擴展列△ P、△ Q for k2=k1:N1 % 擴展列△ 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 行從第三行開始到當前行前的 k4 行消去 for k2=k1:N1 % k3 行后各行下三角元素 J(k3,k2)=J(k3,k2)J(k3,k)*J(k,k2)。%消去運算 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)。%消去運算 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)。%消去運算 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é)點電壓實部 k1=k+1。 f(L)=f(L)J(k1,N)。 %修改節(jié)點電壓虛部 end %修改節(jié)點電壓 for k=3:N0 DET=abs(J(k,N))。 if DET=pr %電壓偏差量是否滿足要求 IT2=IT2+1。 %不滿足要求的節(jié)點數(shù)加 1 end end ICT2(a)=IT2。 ICT1=ICT1+1。 end %用高斯消去法解 w=J*V disp(39。迭代次數(shù)： 39。)。 disp(ICT1)。 disp(39。沒有達到精度要求的個數(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 %=============== 計算各輸出量 =========================== disp(39。各節(jié)點的實際電壓標幺值 E 為 (節(jié)點號 從小到大排列 )： 39。)。 disp(E)。 EE=E*UB。 disp(EE)。 disp(39。39。)。 disp(39。各節(jié)點的電壓大小 V為 (節(jié)點號從小到大排列 )： 39。)。 disp(V)。 VV=V*UB。 disp(VV)。 disp(39。39。)。 disp(39。各節(jié)點的電壓相角 sida 為 (節(jié)點號從小到大排列 )： 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)。 end disp(39。各節(jié)點的功率 S 為 (節(jié)點號從小到大排列 )： 39。)。 disp(S)。 disp(39。~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~39。)。 SS=S*SB。 disp(SS)。 disp(39。39。)。 disp(39。各條支路的首端功率 Si 為 (順序同您輸入 B1 時一致 )： 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。)。 end disp(39。各條支路的末端功率 Sj 為 (順序同您輸入 B1 時一致 )： 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。)。 end disp(39。各條支路的功率損耗 DS 為 (順序同您輸入 B1 時一致