【正文】 變電所4 變比 未調(diào)整 1 1 1 1 調(diào)整后 調(diào)節(jié)后的輸出結(jié)果： B1=[1 3 + 1 0。 2 5 + 1 0。 3 4 + 1 0。 4 5 + 1 0。 5 9 + 0 1。 B2=[0 0 0 1。 0 0 1 0 0 2。 0 0 1 0 0 2。 0 + 1 0 0 2。 0 + 1 0 0 2] 調(diào)整前后各節(jié)點(diǎn)的 電壓 標(biāo)幺值如下： 節(jié)點(diǎn)電壓 1 2 3 4 5 6 7 8 9 10 調(diào)整前 調(diào)整后 各支路的功率損耗如下表： 功率損耗 未調(diào)整 調(diào)整 后 （ 1,3）支路 （ 1,4）支路 （ 2,5）支路 (2,6) 支路 (3,4) 支路 (3,7) 支路 + + (4,5) 支路 (4,8) 支路 + + (5,9) 支路 + + (6,10) 支路 + + 各節(jié)點(diǎn)的功率 S 為 (節(jié)點(diǎn)號從小到大排列 )： + + + + 調(diào)整后的具體潮流分布 各條支路的首端功率 Si 各條支路的末端功率 Sj S(1,3)=+ S(3,1)= S(1,4)=+ S(4,1)= S(2,5)=+ S(5,2)= S(2,6)=+ S(6,2)= S(3,4)= S(4,3)=+ S(3,7)=+ S(7,3)=5031i S(4,5)= S(5,4)= S(4,8)=+ S(8,4)= S(5,9)=+ S(9,5)=5031i S(6,10)=+ S(10,6)= 圖像如下： 采用 DDRTS 進(jìn)行潮流分析 原始數(shù)據(jù)的潮流分布圖 4 個變電所的負(fù)荷同時以 2%的比例增大的潮流分布圖 4個變電所的負(fù)荷同時以 2%的比例減小 的潮流分布圖 1 和 4 號變電所的負(fù)荷同時以 2%的比例下降，而 2 和 3 號變電所的負(fù)荷同時以 2%的比例上升的潮流分布圖 斷開 3 節(jié)點(diǎn)之間一回線時的潮流分布圖 斷開 4 節(jié)點(diǎn)之間一回線時的潮流分布 圖 斷開 4 節(jié)點(diǎn)之間一回線時的潮流分布圖 斷開 5 節(jié)點(diǎn)之間一回線時的潮流分布圖 斷開 6 節(jié)點(diǎn)之間一回線時的潮流分布圖 斷開 6 節(jié)點(diǎn)之間一回線時的潮流分布圖 心得體會 通過這學(xué)期電力系統(tǒng)分析課程的學(xué)習(xí)和這次課程設(shè)計，我對潮流計算有了更深刻的認(rèn)識。根據(jù)系統(tǒng)給定的運(yùn)行條件、網(wǎng)絡(luò)接線及元件參數(shù)，通過潮流計算可以確定各母線的電壓（相角及幅值）、各元件中流過的功率、整個系統(tǒng)的功率損耗等。 從本次電力系統(tǒng)潮流計算課程設(shè)計中，通過計算在不同的負(fù)荷下各支路節(jié)點(diǎn)潮流分布的變化，分析負(fù)荷變動時，對電壓，功率的影響，我們掌握了常用的調(diào)整電壓的措施。 MATLAB 是此次潮流計算中最重要的工具，它是一套高性能的數(shù)學(xué)計算軟件，它集數(shù)值分析、矩陣計算、信號處理和圖形顯示于一身，構(gòu)成了一個方便的界面和友好的用戶環(huán)境，其強(qiáng)大的擴(kuò)展功能為各個領(lǐng)域的應(yīng)用提供了方便。 本次課程設(shè)計除了提高了我的專業(yè)知識水平外，也加強(qiáng)了我與同學(xué)的溝通合作能力，使我 受益匪淺。 n=10。請輸入節(jié)點(diǎn)數(shù)： n=39。 nl=11。請輸入支路數(shù)： nl=39。 isb=1。請輸入平衡母線節(jié)點(diǎn)號： isb=39。 pr=。請輸入誤差精度： pr=39。 B1=[1 3 + 1 0。 2 5 + 1 0。 3 4 + 1 0。 4 5 + 1 0。 5 6 + 1 0。 6 10 + 0 1]。請輸入由支路參數(shù)形成的矩陣： B1=39。 B2=[0 0 0 1。 0 0 1 0 0 2。 0 0 1 0 0 2。 0 + 1 0 0 2。 0 + 1 0 0 2。] %input(39。)。e=zeros(1,n)。V=zeros(1,n)。S1=zeros(nl)。 SB=100。 %ym=input(39。)。請輸 入功率基準(zhǔn)值 :SB=39。 %UB=input(39。)。 BB1=B1。 for i=1:nl B1(i,3)=B1(i,3)*YB。 end disp(39。)。 B2(i,2)=B2(i,2)./SB。 B2(i,4)=B2(i,4)./UB。 end disp(39。)。q=B1(i,2)。q=B1(i,1)。 %非對角元 Y(q,p)=Y(p,q)。 %對角元 K 側(cè) Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2。導(dǎo)納矩陣 Y=39。 disp(Y) % G=real(Y)。 %分解出導(dǎo)納陣的實部和虛部 for i=1:n %給定各節(jié)點(diǎn)初始電壓的實部和虛部 e(i)=real(B2(i,3))。 V(i)=B2(i,4)。 %i 節(jié)點(diǎn)注入功率 SGSL B(i,i)=B(i,i)+B2(i,5)。Q=imag(S)。IT2=1。N=N0+1。 while IT2~=0 IT2=0。 for i=1:n if i~=isb %非平衡節(jié)點(diǎn) C(i)=0。 for j1=1:n C(i)=C(i)+G(i,j1)*e(j1)B(i,j1)*f(j1)。%Σ (Gij*fj+Bij*ej) end P1=C(i)*e(i)+f(i)*D(i)。%節(jié)點(diǎn)功率 Q 計算 fiΣ (Gij*ejBij*fj)eiΣ (Gij*fj+Bij*ej) %求 P39。 V2=e(i)^2+f(i)^2。 %節(jié)點(diǎn)有功功率差 DQ=Q(i)Q1。j1~=i %非平衡節(jié)點(diǎn) amp。 % dP/de=dQ/df X2=B(i,j1)*e(i)G(i,j1)*f(i)。 % X2=dp/df X3=dQ/de X4=X1。q=2*j11。J(p,N)=DQ。 J(m,q)=X1。q=q+1。J(m,q)=X2。j1~=isb %非平衡節(jié)點(diǎn) amp。% dP/de X2=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)。q=2*j11。J(p,N)=DQ。 J(m,q)=X1。J(p,q)=X4。%擴(kuò)展列△ P J(m,q)=X2。 % PV 節(jié)點(diǎn)有功誤差 DV=V(i)^2V2。j1~=i %非平衡節(jié)點(diǎn) amp。 % dP/de X2=B(i,j1)*e(i)G(i,j1)*f(i)。X6=0。q=2*j11。J(p,N)=DV。 J(m,q)=X1。q=q+1。 J(m,q)=X2。j1~=isb %非平衡節(jié)點(diǎn) amp。% dP/de X2=D(i)+B(i,i)*e(i)G(i,i)*f(i)。 X6=2*f(i)。q=2*j11。J(p,N)=DV。 J(m,q)=X1。q=q+1。 J(m,q)=X2。N1=N。 % 非對角元規(guī)格化 end J(k,k)=1。 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)。 end if k==N0 break。%消去運(yùn)算 end J(k3,k)=0。%消去運(yùn)算 end J(k3,k)=0。 e(L)=e(L)J(k,N)。 f(L)=f(L)J(k1,N)。 if DET=pr %電壓偏差量是否滿足要求 IT2=IT2+1。 ICT1=ICT1+1。迭代次數(shù)： 39。 disp(ICT1)。沒有達(dá)到精度要求的個數(shù)： 39。 disp(ICT2)。 sida(k)=atan(f(k)./e(k))*180./pi。 end %=============== 計算各輸出量 =========================== disp(39。)。 EE=E*UB。 disp(39。)。各節(jié)點(diǎn)的電壓大小 V為 (節(jié)點(diǎn)號從小到大排列 )： 39。 disp(V)。 disp(VV)。39。 disp(39。)。 for p=1:n C(p)=0。 end S(p)=E(p)*C(p)。各節(jié)點(diǎn)的功率 S 為 (節(jié)點(diǎn)號從小到大排列 )： 39。 disp(S)。~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~39。 SS=S*SB。 disp(39。)。各條支路的首端功率 Si 為 (順序同您輸入 B1 時一致 )： 39。 for i=1:nl p=B1(i,1)。 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))))。 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))))。 end disp(Si(p,q))。 ZF=[39。,num2str(p),39。,num2str(q),39。,num2str(SSi(p,q))]。 %disp(SSi(p,q))。39。 end disp(39。)。q=B1(i,2)。 Sjy(i)=Sj(q,p)。 Sjy(i)=Sj(q,p)。 SSj(q,p)=Sj(q,p)*SB。S(39。,39。)=39。 disp(ZF)。 disp(39。)。各條支路的功率損耗 DS 為 (順序同您輸入 B1 時一