【正文】 J(m,N)=DP。%擴(kuò)展列△ Q m=p+1。% dQ/df p=2*i1。對(duì)角元 X1=C(i)G(i,i)*e(i)B(i,i)*f(i)。 J(p,q)=X4。m=p+1。 % X1=dP/de X4=dQ/df p=2*i1。非對(duì)角元 X1=G(i,j1)*e(i)B(i,j1)*f(i)。 %電壓模平方 %========= 以下針對(duì)非 PV 節(jié)點(diǎn)來求取功率差及 Jacobi 矩陣元素 ========= if B2(i,6)~=3 %非 PV 節(jié)點(diǎn) DP=P(i)P1。% 節(jié)點(diǎn)功率 P 計(jì)算 ei Σ(Gij*ejBij*fj)+fiΣ (Gij*fj+Bij*ej) Q1=C(i)*f(i)e(i)*D(i)。D(i)=0。a=0。 ICT1=0。 %PV 節(jié)點(diǎn)電壓給定模值 end for i=1:n %給定各節(jié)點(diǎn)注入功率 S(i)=B2(i,1)B2(i,2)。B=imag(Y)。 %對(duì)角元 1 側(cè) end %求導(dǎo)納矩陣 disp(39。 end Y(p,q)=Y(p,q)1./(B1(i,3)*B1(i,5))。 disp(B2) end % % % for i=1:nl %支路數(shù) if B1(i,6)==0 %左節(jié)點(diǎn)處于低壓側(cè) p=B1(i,1)。 B2(i,5)=B2(i,5)./SB。 disp(B1) for i=1:n B2(i,1)=B2(i,1)./SB。 B1(i,4)=B1(i,4)./YB。 YB=SB./UB./UB。)。您輸入的參數(shù)是標(biāo)么值？（若不是則輸入一個(gè)不為零的數(shù)值）39。 %－－－－－－－修改部分－－－－－－－－－－－－ ym=0。f=zeros(1,n)。請(qǐng)輸入各節(jié)點(diǎn)參數(shù)形成的矩陣： B2=39。 0 0 220 0 2。 0 0 220 0 2。 7 9 + 1 0 0。 2 9 + 1 0 0。)。)。)。)。 附錄 1：潮流分布圖 附錄 2： MATLAB 程序 %本程序的功能是用牛頓 —— 拉夫遜法進(jìn)行潮流計(jì)算 % B1 矩陣： 支路首端號(hào)； 末端號(hào)； 支路阻抗； 支路對(duì)地電納 % 支路的變比； 支路首端處于 K側(cè)為 1， 1側(cè)為 0 % B2 矩陣： 該節(jié)點(diǎn)發(fā)電機(jī)功率； 該節(jié)點(diǎn)負(fù)荷功率； 節(jié)點(diǎn)電壓初始值 % PV 節(jié)點(diǎn)電壓 V 的給定值； 節(jié)點(diǎn)所接的無功補(bǔ)償設(shè)備的容量 % 節(jié)點(diǎn)分類標(biāo)號(hào) clear。我們應(yīng)該努力學(xué)好專業(yè)課，打好專業(yè)基礎(chǔ)，在以后的工作崗位上有所成就。這個(gè)方案充分利用了前 兩種調(diào)壓措施，使電壓都滿足了要求，保證了調(diào)壓的經(jīng)濟(jì)性。 七、心得體會(huì) 從本次電力系統(tǒng)潮流計(jì)算課程設(shè)計(jì)中，通過計(jì)算在不同的負(fù)荷下各支路節(jié)點(diǎn)潮流分布的變化，分析負(fù)荷變動(dòng)時(shí)，對(duì)電壓，功率的影響，我們掌握了常用的調(diào)整電壓的措施 ， 鍛煉了自己實(shí)際操作分析能力 。有上述的數(shù)據(jù)分析， 當(dāng)調(diào)節(jié)機(jī)端電壓不能滿足各個(gè)變電所的電壓要求時(shí)， 若 個(gè)別 變電所電壓偏高， 則升高變壓器的變比；若 個(gè)別 變電所電壓偏低，則降低變壓器的變比， 通常可以達(dá)到調(diào)壓要求。 對(duì)于電力系統(tǒng)，調(diào)節(jié)電壓首先要考慮改變機(jī)端電壓，其是最經(jīng)濟(jì)最直接的調(diào)壓措施。電壓過高或過低，都將對(duì)人身及其用電設(shè)備產(chǎn)生重大的影響。 當(dāng)斷開雙回線中的一條時(shí)，由于對(duì)應(yīng)支路阻抗的增大， 與該支路相連線路輸送的功率會(huì)有所改變， 有 功損耗會(huì)隨之增加 ，其余支路的潮流幾乎不發(fā)生變化 。 在該網(wǎng)絡(luò)中，由于變電所 1和 2離發(fā)電廠一比較近，變電所 4離發(fā)電廠二比較近，在選擇 PV 節(jié)點(diǎn)的功率時(shí)，要考慮功率輸送的經(jīng)濟(jì)性，所以發(fā)電廠二的功率選擇了 70MW， 可將網(wǎng)損降到很小。 電壓 1 2 3 4 5 6 7 8 9 10 DDRTS MATLAB 由上述比較可觀察到 : 應(yīng)用 DDRTS 仿真和 MATLAB 計(jì)算的到的結(jié)果大致相同，可以認(rèn)為此潮流計(jì)算是正確的 。具體的潮流分布如下： 各條支路的首端功率 Si 各條支路的末端功率 Sj S(1,3)=+ S(3,1)= S(1,5)=+ S(5,1)= S(1,7)=+ S(7,1)= S(2,9)=70+ S(9,2)= S(3,4)=+ S(4,3)= S(5,6)=+ S(6,5)= S(7,8)=+ S(8,7)= S(7,9)= S(9,7)= S(9,10)=+ S(10,9)= 各節(jié)點(diǎn)的功率 S為 (節(jié)點(diǎn)號(hào)從小到大排列 )： + + + 斷開 9 支路的一條線 當(dāng)斷開 9節(jié)點(diǎn)間的雙回線中的 一條 時(shí) ， B1矩陣 的改變， B1 矩陣變?yōu)? 1 3 + 1 0 0 1 5 + 1 0 0 1 7 + 1 0 0 2 9 17+ 1 0 0 3 4 + 0 220/ 1 1 5 6 + 0 220/11 1 1 7 8 + 0 220/ 1 1 7 9 + 1 0 0 9 10 + 0 220/11 1 1 用修改后的 B2 矩陣帶入程序，用 MATLAB 進(jìn)行潮流計(jì)算并進(jìn)行調(diào)整， 對(duì)應(yīng)的調(diào)節(jié)方法和結(jié)果見下表： 調(diào)節(jié)方法： 調(diào)節(jié)方法 發(fā)電廠一 發(fā)電廠二 分接頭一 分接頭二 分接頭三 分接頭四 未調(diào)整 235 235 220 220 220 220 調(diào)整后 225 225 231 231 231 調(diào)整結(jié)果： 電壓 1 2 3 4 5 6 7 8 9 10 未調(diào)整 調(diào)整后 由上表觀察到， 調(diào)節(jié)后 的電壓均在題目允許的范圍內(nèi)，對(duì)線路損耗進(jìn)行分析 ， 統(tǒng)計(jì)每次調(diào)整后各個(gè)支路的有功損耗 ， 記錄于下表 ： 功率損耗 未調(diào)整 調(diào)整后 (1,3)支路 (1,5)支路 (1,7)支路 (2,9)支路 (3,4)支路 + + (5,6)支路 + + (7,8)支路 + + (7,9)支路 (9,10)支路 + + 總損耗 由電壓和功率損耗可以觀察到，最后一次的電壓和有功損耗均符合題目要求，可認(rèn)為是合理的。 將 DDRTS 仿真得到的潮流圖中個(gè)點(diǎn)電壓的有名值與 MATLAB 計(jì)算得到的各節(jié)點(diǎn)電壓進(jìn)行比較： 電壓 1 2 3 4 5 6 7 8 9 10 DDRTS MATLAB 由上述比較可觀察到 : 應(yīng)用 DDRTS仿真和 MATLAB 計(jì)算的到的 電壓 結(jié)果大致相同， 各個(gè)支路的功率分布也大概相同， 可以認(rèn)為此潮流計(jì)算是正確的。 各節(jié)點(diǎn)的功率 S為 (節(jié)點(diǎn)號(hào)從小到大排列 )： + + + 4 號(hào)變電所負(fù)荷成 2%下降， 3 號(hào)變電負(fù)荷成 2%上升 此時(shí)的各負(fù)荷變?yōu)椋? 變電所 1 LS =+ 變電所 2 LS = 51+ 變電所 3 LS =+ 變電所 4 LS =+ 相應(yīng)的， B2矩陣變?yōu)椋? 0 0 235 0 1 70 0 235 0 3 0 0 220 0 2 0 + 35 0 2 0 0 220 0 2 0 51+ 10 0 2 0 0 220 0 2 0 + 35 0 2 0 0 220 0 2 0 + 10 0 2 將程序中的 B2矩陣替換后，運(yùn)行并進(jìn)行調(diào)節(jié)，此系統(tǒng)經(jīng)過調(diào)節(jié)后達(dá)到要求，對(duì)應(yīng)的調(diào)節(jié)方法和結(jié)果見下表： 調(diào)節(jié)方法： 調(diào)節(jié)方法 發(fā)電廠一 發(fā)電廠二 分接頭一 分接頭二 分接頭三 分接頭 四 未調(diào)整 235 235 220 220 220 220 調(diào)整后 225 225 231 231 231 調(diào)整結(jié)果： 電壓 1 2 3 4 5 6 7 8 9 10 未調(diào)整 調(diào)整后 由上表觀察到， 調(diào)節(jié)后 的電壓均在題目允許的范圍內(nèi)， 負(fù)荷增加的變電所電壓有B2= 所下降，負(fù)荷減小電壓升高。 具體的潮流分布如下： 各條支路的首端功率 Si 各條支路的末端功率 Sj S(1,3)=+ S(3,1)= S(1,5)=+ S(5,1)= S(1,7)=+ S(7,1)= S(2,9)= S(9,2)= S(3,4)=+ S(4,3)= S(5,6)=+ S(6,5)= S(7,8)=+ S(8,7)= S(7,9)=+ S(9,7)= S(9,10)=+ S(10,9)= 由表可以看出，各支路輸送的功率較正常時(shí)略有 增加。調(diào)整后的各支路的功率如下表所示： 各條支路的首端功率 Si 各條支路的末端功率 Sj S(1,3)=+ S(3,1)= S(1,5)=+ S(5,1)= S(1,7)=+ S(7,1)= S(2,9)= S(9,2)= S(3,4)=+ S(4,3)= S(5,6)=+ S(6,5)= S(7,8)=+ S(8,7)= S(7,9)=+ S(9,7)= S(9,10)=+ S(10,9)= 支路的末端功率加上支路損耗等于該支路的首段功率，符合實(shí)際情況。 圖 等值電路圖 60+ 發(fā)電廠一 1:k21:k3 1 3 5 7 4 6 8 δ =0 U=235 P=70MW +2j 5+15j 5+18.09j 5455+1.65024j 5455+1.65024j 0373+0.10084j 1:k1 50+ 40+ j + j 8.5+2j 發(fā)電廠二 U=235 60+ 利用牛頓拉夫遜法進(jìn)行求解，用 MATLAB 軟件編程，可以求解系統(tǒng)潮流分 布 。則網(wǎng)絡(luò)共具有 10 個(gè)節(jié)點(diǎn)和 9條支路。（幾條支路斷幾次） 6. 利用 DDRTS 軟件，進(jìn)行繪制系統(tǒng)圖進(jìn)行上述各種情況潮流的分析，并進(jìn)行結(jié)果的比較。 2. 輸入各支路數(shù)據(jù)，各節(jié)點(diǎn)數(shù)據(jù)利用給定的程序進(jìn)行在變電所在某一負(fù)荷情況下的潮流計(jì)算，并對(duì)計(jì)算結(jié)果進(jìn)行分析。 線路長為 80km 。 （ 2） 牛頓 — 拉夫遜法計(jì)算程序框圖 啟動(dòng) 輸入原始數(shù)據(jù) 形成節(jié)點(diǎn)導(dǎo)納矩陣 分解各節(jié)點(diǎn)初始電壓的實(shí)部和虛部 迭代次數(shù) K=0 求 PQ 節(jié)點(diǎn)的 , ，求 PV 節(jié)點(diǎn)的 ， 置節(jié)點(diǎn)號(hào) i=0 雅克比矩陣是否形成， in 求得雅克比矩陣各元素增大節(jié)點(diǎn)號(hào) i=i+1 把雅克比矩陣單位化 求解修正方程，得 ， 求解最大修正量 ， 是否收斂 回帶各電壓新值， K=K+1 計(jì)算輸出電壓大小及相角，節(jié)電功率及支路損耗 停止 否 是 否 是 三、課程設(shè)計(jì)任務(wù) 數(shù)據(jù) ： 兩個(gè)發(fā)電廠分別通過變壓器和輸電線路與四個(gè)變電所相連。這一步在收斂到希望的值以前重復(fù)進(jìn)行，一般要反復(fù)計(jì)算滿足 ? ? ????? ?????? 1112121111 ,m a x nnnnnnnn XXXXXX ? ? 為預(yù)先規(guī)定的小正數(shù)， 1?nnX 是第 n次迭代 nX 的近似值。 )( )( )()(nnXf Xf???? 比較兩式，可以看出牛頓 — 拉夫遜法的休整量和 )(nX 的誤差的一次項(xiàng)