【正文】 disp(39。S(39。q=B1(i,2)。 %disp(SSi(p,q))。 ZF=[39。 for i=1:nl p=B1(i,1)。 SS=S*SB。 end S(p)=E(p)*C(p)。39。)。 end %=============== 計(jì)算各輸出量 =========================== disp(39。 disp(ICT1)。 f(L)=f(L)J(k1,N)。 end if k==N0 break。 J(m,q)=X2。q=2*j11。 J(m,q)=X2。q=2*j11。 % PV 節(jié)點(diǎn)有功誤差 DV=V(i)^2V2。J(p,N)=DQ。j1~=isb %非平衡節(jié)點(diǎn) amp。J(p,N)=DQ。j1~=i %非平衡節(jié)點(diǎn) amp。%Σ (Gij*fj+Bij*ej) end P1=C(i)*e(i)+f(i)*D(i)。N=N0+1。 V(i)=B2(i,4)。 %對角元 K 側(cè) Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2。)。)。)。 %ym=input(39。e=zeros(1,n)。 0 + 1 0 0 2。請輸入由支路參數(shù)形成的矩陣： B1=39。 3 4 + 1 0。 pr=。 nl=11。 MATLAB 是此次潮流計(jì)算中最重要的工具，它是一套高性能的數(shù)學(xué)計(jì)算軟件，它集數(shù)值分析、矩陣計(jì)算、信號處理和圖形顯示于一身，構(gòu)成了一個(gè)方便的界面和友好的用戶環(huán)境，其強(qiáng)大的擴(kuò)展功能為各個(gè)領(lǐng)域的應(yīng)用提供了方便。 0 + 1 0 0 2。 5 9 + 0 1。 調(diào)節(jié)方 法 發(fā)電廠1 機(jī)端電壓 發(fā)電廠2 機(jī)端電壓 變電所1 變比 變電所2 變比 變電所3 變比 變電所4 變比 未調(diào)整 1 1 1 1 調(diào)整后 調(diào)節(jié)后的輸出結(jié)果： B1=[1 3 + 1 0。 3 4 + 1 0。 5 6 + 1 0。 電壓調(diào)節(jié)方法：通過改變發(fā)電機(jī)極端電壓、變壓器變比或增加無功功率補(bǔ)償來調(diào)整電壓。 3 4 + 1 0。 4 8 + 0 1。 B2矩陣不變。 2 6 + 1 0。 0 0 1 0 0 2。 4 5 + 1 0。 6 10 + 0 1]。 1 4 + 1 0。 0 0 3 0 0 1 0 0 2。 3 4 + 1 0。 5 6 + 1 0。 B2矩陣不變。 2 6 + 1 0。 0 0 1 0 0 2。 4 5 + 1 0。 0 + 1 0 0 2。 5 6 + 1 0。 調(diào)節(jié)方 法 發(fā)電廠1 機(jī)端電壓 發(fā)電廠2 機(jī)端電壓 變電所1 變比 變電所2 變比 變電所3 變比 變電所4 變比 未調(diào)整 1 1 1 1 調(diào)整后 調(diào)整后的輸出結(jié)果： B1=[1 3 + 1 0。 0 0 3 0 0 1 0 0 2。 0 0 3 0 0 1 0 0 2。 3 4 + 1 0。 0 + 1 0 0 2。 0 + 1 0 0 2。 5 6 + 1 0。 調(diào)節(jié)方 法 發(fā)電廠1 機(jī)端電壓 發(fā)電廠2 機(jī)端電壓 變電所1 變比 變電所2 變比 變電所3 變比 變電所4 變比 未調(diào)整 1 1 1 1 調(diào)整后 調(diào)整后的輸出結(jié)果： B1=[1 3 + 1 0。 0 0 3 0 0 1 0 0 2。 0 0 1 0 0 2。 3 7 + 0 1。利用成形的程序?qū)ο到y(tǒng)中出現(xiàn)的各種情況，例如負(fù)荷的變化以及線路上所發(fā)生的變化進(jìn)行計(jì)算，并對母線上不滿足范圍的數(shù)據(jù)進(jìn)行調(diào)控， 最終進(jìn)行合理的系統(tǒng)潮流。 三、課程設(shè)計(jì)成品基本要求： 1. 在讀懂程序的基礎(chǔ)上畫出潮流計(jì)算基本流程圖 2. 通過輸入數(shù)據(jù)，進(jìn)行潮流計(jì)算輸出結(jié)果 3. 對不同的負(fù)荷變化，分析潮流分布，寫出分析說明。 發(fā)電廠資料： 母線 1 和 2 為發(fā)電廠高壓母線，發(fā)電 廠一總裝機(jī)容量為 350MW， 母線 3 為 機(jī)壓母線，機(jī)壓母線上裝機(jī)容量為 100MW， 最大負(fù)荷和最小負(fù)荷分別為 40MW 和 25MW；發(fā) 電廠二總裝機(jī)容量為 200MW。 變電所資料： (一 ) 變電所 4 低壓母線的電壓等級分別為： 35KV 35KV 35KV 10KV (二 ) 變電所的負(fù)荷分別為： （ 1） 50MW 50MW 70MW 60MW （ 2） 60MW 70MW 50MW 45MW （ 3） 50MW 40MW 50MW 70MW 電廠一 電廠二 （ 4） 70MW 50MW 60MW 60MW （ 5） 70MW 60MW 50MW 50MW （ 6） 45MW 55MW 60MW 70MW （ 7） 70MW 45MW 50MW 60MW （ 8） 50MW 40MW 50MW 60MW （ 9） 40MW 40MW 50MW 60MW （ 10） 50MW 40MW 40MW 60MW (三 ) 每個(gè)變電所的功率因數(shù)均為 cosφ =； (四 ) 變電所 1 和變電所 2 分別配有兩臺容量為 75MVA 的變壓器，短路損耗 414KW，短路電壓（ %） =；變電所 3 和變電所 4分別配有兩臺容量為 63MVA 的變壓器，短路損耗為 245KW，短 路電壓（ %） =； 輸電線路資料： 發(fā)電廠和變電所之間的輸電線路的電壓等級及長度標(biāo)于圖中，單位長度的電阻為? ，單位長度的電抗為 ? ，單位長度的電納為 S610* 。 4. 對不同的負(fù)荷變化，進(jìn)行潮流的調(diào)節(jié)控制，并說明調(diào)節(jié)控制的方法，并列表表示調(diào)節(jié)控制的參數(shù)變化。 關(guān)鍵詞 牛頓 拉夫遜法（ NewtonRaphson） 潮流計(jì)算與 DDRT S 仿真 潮流計(jì)算結(jié)果分析與調(diào)控 利用 MATLAB 程序運(yùn)行的結(jié)果 未調(diào)整時(shí)各節(jié)點(diǎn)電壓的標(biāo)么值如下： 圖像如下： 由調(diào)壓范圍可知，變電所 4 的節(jié)點(diǎn)電壓均不滿足電壓范圍。 4 5 + 1 0。 0 0 1 0 0 2。 0 0 1 0 0 2。 1 4 + 1 0。 5 9 + 0 1。 0 + 1 0 0 2。 0 + 1 0 0 2。 3 7 + 0 1。 0 0 1 0 0 2。 0 0 1 0 0 2。 1 4 + 1 0。 5 9 + 0 1。 0 + 1 0 0 2。 4 8 + 0 1。 0 0 1 0 0 2。 3 4 + 1 0。 調(diào)整前后各節(jié)點(diǎn) 電壓 的標(biāo)幺值如下： 各支路的功率損耗如下表： 功率損耗 未調(diào)整 調(diào)整后 （ 1,4）支路 （ 2,5）支路 (2,6) 支路 (3,4) 支路 (3,7) 支路 + + (4,5) 支路 (4,8) 支路 + + (5,6) 支路 (5,9) 支路 + + (6,10) 支路 + + 各節(jié)點(diǎn)的功率 S 為 (節(jié)點(diǎn)號從小到大排列 )： + + + + 0 + 調(diào)整后的具體潮流分布 各條支路的首端功率 Si 各條支路的末端功率 Sj S(1,4)=+ S(4,1)= 節(jié)點(diǎn)電壓 1 2 3 4 5 6 7 8 9 10 調(diào)整前 0,9716 調(diào)整后 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,6)= S(6,5)= S(5,9)=+ S(9,5)=5031i S(6,10)=+ S(10,6)= 圖像如下： 斷開 4節(jié)點(diǎn)之間一回線時(shí)的潮流計(jì)算 B1矩陣變?yōu)椋? B1=[1 3 + 1 0。 5 9 + 0 1。 3 7 + 0 1。 0 0 1 0 0 2。 2 5 + 1 0。 B2矩陣不變。 4 8 + 0 1。 0 0 1 0 0 2。 3 4 + 1 0。 電壓調(diào)節(jié)方法：通過改變發(fā)電機(jī)極端電壓、變壓器變比或增加無功功率補(bǔ)償來調(diào)整電壓。 5 6 + 1 0。 3 7 + 0 1。 調(diào)節(jié)方 法 發(fā)電廠1 機(jī)端電壓 發(fā)電廠2 機(jī)端電壓 變電所1 變比 變電所2 變比 變電所3 變比 變電所4 變比 未調(diào)整 1 1 1 1 調(diào)整后 調(diào)節(jié)后的輸出結(jié)果： B1=[1 3 + 1 0。 5 9 + 0 1。 3 7 + 0 1。 1 4 + 1 0。 6 10 + 0 1]。 0 + 1 0 0 2。在潮流計(jì)算的大量數(shù)值計(jì)算過程中更顯示出其優(yōu)勢，是我們應(yīng)該掌握的一門基本技術(shù)。%input(39。%input(39。 3 7 + 0 1。)。 0 + 1 0 0 2。f=zeros(1,n)。您輸入的參數(shù)是標(biāo)么值？（若不是則輸入一個(gè)不為零的數(shù)值） 39。 YB=SB./UB./UB。 disp(B1) for i=1:n B2(i,1)=B2(i,1)./SB。 disp(B2) end % % % for i=1:nl %支路數(shù) if B1(i,6)==0 %左 節(jié)點(diǎn)處于低壓側(cè) p=B1(i,1)。 %對角元 1 側(cè) end %求導(dǎo)納矩陣 disp(39。 %PV 節(jié)點(diǎn)電壓給定模值 end for i=1:n %給定各節(jié)點(diǎn)注入功率 S(i)=B2(i,1)B2(i,2)。a=0。%節(jié)點(diǎn)功率 P 計(jì)算 eiΣ (Gij*ejBij*fj)+fiΣ (Gij*fj+Bij*ej) Q1=C(i)*f(i)e(i)*D(i)。非對角元 X1=G(i,j1)*e(i)B(i,j1)*f(i)。m=p+1。對角元 X1=C(i)G(i,i)*e(i)B(i,i)*f(i)。%擴(kuò)展列△ Q m=p+1。 % PV 節(jié)點(diǎn)電壓誤差