【正文】 N)=DP。 end end end end end%========= 以上為求雅可比矩陣的各個元素 ===================== for k=3:N0 % N0=2*n （從第三行開始，第一、二行是平衡節(jié)點） k1=k+1。 % 對角元規(guī)格化 if k~=3 % 不是第三行 %============================================================ k4=k1。 end%========================================== for k3=k1:N0 for k2=k1:N1 J(k3,k2)=J(k3,k2)J(k3,k)*J(k,k2)。 end end end %====上面是用線性變換方式將Jacobi矩陣化成單位矩陣===== for k=3:2:N01 L=(k+1)./2。 %修改節(jié)點電壓虛部 end %修改節(jié)點電壓 for k=3:N0 DET=abs(J(k,N))。end%用高斯消去法解w=J*Vdisp(39。disp(39。for k=1:n V(k)=sqrt(e(k)^2+f(k)^2)。各節(jié)點的實際電壓標(biāo)幺值E為(節(jié)點號從小到大排列)：39。disp(EE)。disp(39。VV=V*UB。)。disp(sida)。enddisp(39。disp(39。disp(SS)。disp(39。q=B1(i,2)。 Siz(i)=Si(p,q)。S(39。)=39。 disp(39。各條支路的末端功率Sj為(順序同您輸入B1時一致)：39。 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))))。 end disp(Sj(q,p))。,num2str(q),39。,num2str(SSj(q,p))]。39。)。 disp(DS(i))。,num2str(p),39。,num2str(DDS(i))]。39。plot(V)。ylabel(39。subplot(2,2,2)。)。grid on。節(jié)點注入有功39。bar(Siz)。grid on。支路首端無功39。grid on。bar(real(S))。電壓角度39。xlabel(39。)。節(jié)點號39。endfigure(1)。 %disp(DDS(i))。,num2str(q),39。 ZF=[39。q=B1(i,2)。enddisp(39。 %disp(SSj(q,p))。,num2str(p),39。 ZF=[39。 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))))。for i=1:nl p=B1(i,1)。)。 disp(ZF)。,39。 SSi(p,q)=Si(p,q)*SB。 Siz(i)=Si(p,q)。)。39。)。)。 for q=1:n C(p)=C(p)+conj(Y(p,q))*conj(E(q))。各節(jié)點的電壓相角sida為(節(jié)點號從小到大排列)：39。disp(39。)。39。disp(E)。 E(k)=e(k)+f(k)*j。)。)。 %不滿足要求的節(jié)點數(shù)加1 end end ICT2(a)=IT2。 %修改節(jié)點電壓實部 k1=k+1。 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。 % N=N0+1 即 N=2*n+1擴展列△P、△Q for k2=k1:N1 % 擴展列△P、△Q J(k,k2)=J(k,k2)./J(k,k)。J(p,q)=X6。 m=p+1。 p=2*i1。對角元 X1=C(i)G(i,i)*e(i)B(i,i)*f(i)。J(p,q)=X6。 m=p+1。 p=2*i1。非對角元 X1=G(i,j1)*e(i)B(i,j1)*f(i)。 end end else %=============== 下面是針對PV節(jié)點來求取Jacobi矩陣的元素 =========== DP=P(i)P1。q=q+1。J(p,q)=X3。% dP/df X3=D(i)+B(i,i)*e(i)G(i,i)*f(i)。 elseif j1==iamp。J(m,N)=DP。 J(p,q)=X3。 % dP/df=dQ/de X3=X2。 %節(jié)點無功功率差 %=============== 以上為除平衡節(jié)點外其它節(jié)點的功率計算 =================%================= 求取Jacobi矩陣 =================== for j1=1:n if j1~=isbamp。,Q39。%Σ(Gij*ejBij*fj) D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1)。a=a+1。N0=2*n。 %i節(jié)點無功補償量 end%===================================================================P=real(S)。 f(i)=imag(B2(i,3))。)。 Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2。 else p=B1(i,2)。B2矩陣B2=39。 B2(i,3)=B2(i,3)./UB。B1矩陣B1=39。 BB2=B2。請輸入電壓基準(zhǔn)值:UB=39。if ym~=0 %SB=input(39。UB=220。sida=zeros(1,n)。Y=zeros(n)。 0 + 1 0 0 2。 0 0 1 0 0 2。 0 0 3。%input(39。 4 8 + 0 1。 2 6 + 1 0。)。)。)。)。非常感謝老師給我們提供這次寶貴的實踐機會！ 附錄%本程序的功能是用牛頓——拉夫遜法進行潮流計算% B1矩陣：支路首端號；末端號；支路阻抗；支路對地電納% 支路的變比；支路首端處于K側(cè)為1，1側(cè)為0% B2矩陣：該節(jié)點發(fā)電機功率；該節(jié)點負(fù)荷功率；節(jié)點電壓初始值% PV節(jié)點電壓V的給定值；節(jié)點所接的無功補償設(shè)備的容量% 節(jié)點分類標(biāo)號clear。在本次課程設(shè)計過程當(dāng)中，鍛煉了自己實際操作分析能力，理論聯(lián)系實際，采用三種調(diào)整電壓的措施：⑴改變機端電壓調(diào)壓，一般情況下首先考慮的就是改變機端電壓調(diào)壓，通常采用逆調(diào)壓方式，尤其在整體電壓偏高或偏低時；⑵改變變壓器的變比調(diào)壓，因為改變變壓器的變比不需要二次投資，所以該種調(diào)壓方式也很重要；⑶加無功補償設(shè)備調(diào)壓，當(dāng)調(diào)節(jié)變壓器的變比后仍然不能滿足電壓質(zhì)量的要求，或網(wǎng)絡(luò)中無功不足時就需要就地補償無功，減小支路損耗，改善電壓質(zhì)量，通常我們將調(diào)節(jié)變壓器變比和無功補償一起考慮，盡量做到電容器投入容量最小。潮流計算時電力系統(tǒng)的最基本、最常用的計算。 0 + 1 0 0 2。 0 0 3 0 0 1 0 0 2。 4 8 + 0 1。 2 6 + 1 0。電壓調(diào)節(jié)方法：通過改變發(fā)電機極端電壓、變壓器變比或增加無功功率補償來調(diào)整電壓。 4 8 + 0 1。 2 6 + 1 0。B2矩陣不變。 4 8 + 0 1。 2 5 + 1 0。B2矩陣不變。 4 8 + 0 1。 2 5 + 1 0。6 10 + 0 1]。 4 5 + 1 0。 1 4 + 1 0。6 10 + 0 1]。 4 5 + 1 0。 1 4 + 1 0。 0 + 1 0 0 2。 0 0 1 0 0 2。 5 9 + 0 1。 3 7 + 0 1。調(diào)節(jié)方法發(fā)電廠1機端電壓發(fā)電廠2機端電壓變電所1變比變電所2變比變電所3變比變電所4變比未調(diào)整 1 1 1 1調(diào)整后 調(diào)整后的輸出結(jié)果：B1=[1 3 + 1 0。 5 9 + 0 1。 3 7 + 0 1。0 + 1 0 0 2]調(diào)整前后各節(jié)點的電壓標(biāo)幺值如下：節(jié)點電壓12345678910調(diào)整前調(diào)整后各支路的功率損耗如下表：功率損耗未調(diào)整調(diào)整后（1,3）支路 （2,5）支路 (2,6) 支路 (3,4) 支路 (3,7) 支路 + +(4,5) 支路 (4,8) 支路 + + (5,6) 支路 (5,9) 支路 + + (6,10) 支路 + +各節(jié)點的功率S為(節(jié)點號從小到大排列)： + + + 調(diào)整后的具體潮流分布各條支路的首端功率Si各條支路的末端功率Sj S(1,3)=+ S(3,1)= S(2,5)=+ S(5,2)= S(2,6)=+