【正文】 變電所資料：(一) 變電所4低壓母線的電壓等級(jí)分別為：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分別配有兩臺(tái)容量為75MVA的變壓器，短路損耗414KW，短路電壓（%）=；變電所3和變電所4分別配有兩臺(tái)容量為63MVA的變壓器，短路損耗為245KW，短路電壓（%）=；輸電線路資料：發(fā)電廠和變電所之間的輸電線路的電壓等級(jí)及長(zhǎng)度標(biāo)于圖中，單位長(zhǎng)度的電阻為，單位長(zhǎng)度的電抗為，單位長(zhǎng)度的電納為。（；電壓35KV要求調(diào)整范圍在3536之間）5. 輪流斷開環(huán)網(wǎng)一回線，分析潮流的分布。4. 對(duì)不同的負(fù)荷變化，進(jìn)行潮流的調(diào)節(jié)控制，并說明調(diào)節(jié)控制的方法，并列表表示調(diào)節(jié)控制的參數(shù)變化。 牛頓—拉夫遜NewtonRaphson法是數(shù)學(xué)上解非線性方程組的有效方法，有較好的收斂性。關(guān)鍵詞 牛頓拉夫遜法（NewtonRaphson） 潮流計(jì)算與S仿真 潮流計(jì)算結(jié)果分析與調(diào)控利用MATLAB程序運(yùn)行的結(jié)果未調(diào)整時(shí)各節(jié)點(diǎn)電壓的標(biāo)么值如下： 節(jié)點(diǎn)12345678910電壓圖像如下： 由調(diào)壓范圍可知，變電所4的節(jié)點(diǎn)電壓均不滿足電壓范圍。 2 5 + 1 0。 4 5 + 1 0。6 10 + 0 1]。 0 0 1 0 0 2。 0 + 1 0 0 2。 0 0 1 0 0 2。 0 + 1 0 0 2。 1 4 + 1 0。 3 7 + 0 1。 5 9 + 0 1。 0 0 1 0 0 2。 0 + 1 0 0 2。 0 0 1 0 0 2。 0 + 1 0 0 2。 1 4 + 1 0。 3 7 + 0 1。 5 9 + 0 1。 0 0 1 0 0 2。 0 + 1 0 0 2。 0 0 1 0 0 2。 0 + 1 0 0 2。 1 4 + 1 0。 3 7 + 0 1。 5 9 + 0 1。 0 0 1 0 0 2。 0 + 1 0 0 2。 2 6 + 1 0。 4 8 + 0 1。B2=[0 0 0 1。 0 0 1 0 0 2。0 + 1 0 0 2]電壓調(diào)節(jié)方法：通過改變發(fā)電機(jī)極端電壓、變壓器變比或增加無功功率補(bǔ)償來調(diào)整電壓。 3 4 + 1 0。 5 6 + 1 0。調(diào)整前后各節(jié)點(diǎn)電壓的標(biāo)幺值如下：節(jié)點(diǎn)電壓12345678910調(diào)整前0,9716調(dià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)號(hào)從小到大排列)： + + + + 0 + 調(diào)整后的具體潮流分布各條支路的首端功率Si各條支路的末端功率Sj 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,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)椋築1=[1 3 + 1 0。 3 7 + 0 1。 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é)果：B1=[1 3 + 1 0。 3 7 + 0 1。 5 9 + 0 1。 0 0 1 0 0 2。 0 + 1 0 0 2。 2 5 + 1 0。 4 8 + 0 1。B2矩陣不變。 2 5 + 1 0。 4 8 + 0 1。B2=[0 0 0 1。 0 0 1 0 0 2。0 + 1 0 0 2]。 3 4 + 1 0。 5 6 + 1 0。電壓調(diào)節(jié)方法：通過改變發(fā)電機(jī)極端電壓、變壓器變比或增加無功功率補(bǔ)償來調(diào)整電壓。 3 4 + 1 0。 5 6 + 1 0。調(diào)整前后各節(jié)點(diǎn)的電壓標(biāo)幺值如下：節(jié)點(diǎn)電壓12345678910調(diào)整前0,99730,9998調(diào)整后各支路的功率損耗如下表：功率損耗未調(diào)整調(diào)整后（1,3）支路 （1,4）支路 (2,6) 支路 (3,4) 支路 (3,7) 支路 + +(4,5) 支路 (4,8) 支路 + + (5,6) 支路 (5,9) 支路 + + (6,10) 支路 + +各節(jié)點(diǎn)的功率S為(節(jié)點(diǎn)號(hào)從小到大排列)： + + + 調(diào)整后的具體潮流分布各條支路的首端功率Si各條支路的末端功率Sj S(1,3)=+ S(3,1)= S(1,4)=+ S(4,1)= S(2,6)=110+ 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)=圖像如下： 斷開6節(jié)點(diǎn)之間一回線時(shí)的潮流計(jì)算B1陣變?yōu)椋築1=[1 3 + 1 0。 3 7 + 0 1。 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 7 + 0 1。 5 9 + 0 1。 1 4 + 1 0。 3 7 + 0 1。6 10 + 0 1]。 1 4 + 1 0。 3 7 + 0 1。6 10 + 0 1]。 0 0 1 0 0 2。 0 + 1 0 0 2。潮流計(jì)算是實(shí)現(xiàn)電力系統(tǒng)安全經(jīng)濟(jì)發(fā)展的必要手段和重要環(huán)節(jié)。在潮流計(jì)算的大量數(shù)值計(jì)算過程中更顯示出其優(yōu)勢(shì)，是我們應(yīng)該掌握的一門基本技術(shù)。%input(39。%input(39。%input(39。%input(39。 1 4 + 1 0。 3 7 + 0 1。 5 9 + 0 1。)。 0 0 1 0 0 2。 0 + 1 0 0 2。請(qǐng)輸入各節(jié)點(diǎn)參數(shù)形成的矩陣： B2=39。f=zeros(1,n)。%－－－－－－－修改部分－－－－－－－－－－－－ym=0。您輸入的參數(shù)是標(biāo)么值？（若不是則輸入一個(gè)不為零的數(shù)值）39。)。 YB=SB./UB./UB。 B1(i,4)=B1(i,4)./YB。 disp(B1) for i=1:n B2(i,1)=B2(i,1)./SB。 B2(i,5)=B2(i,5)./SB。 disp(B2)end% % %for i=1:nl %支路數(shù) if B1(i,6)==0 %左節(jié)點(diǎn)處于低壓側(cè) p=B1(i,1)。 end Y(p,q)=Y(p,q)1./(B1(i,3)*B1(i,5))。 %對(duì)角元1側(cè) end%求導(dǎo)納矩陣disp(39。B=imag(Y)。 %PV節(jié)點(diǎn)電壓給定模值 endfor i=1:n %給定各節(jié)點(diǎn)注入功率 S(i)=B2(i,1)B2(i,2)。ICT1=0。a=0。D(i)=0。%節(jié)點(diǎn)功率P計(jì)算eiΣ(Gij*ejBij*fj)+fiΣ(Gij*fj+Bij*ej) Q1=C(i)*f(i)e(i)*D(i)。 %電壓模平方 %========= 以下針對(duì)非PV節(jié)點(diǎn)來求取功率差及Jacobi矩陣元素 ========= if B2(i,6)~=3 %非PV節(jié)點(diǎn) DP=P(i)P1。非對(duì)角元 X1=G(i,j1)*e(i)B(i,j1)*f(i)。 % X1=dP/de X4=dQ/df p=2*i1。m=p+1。 J(p,q)=X4。對(duì)角元 X1=C(i)G(i,i)*e(i)B(i,i)*f(i)。% dQ/df p=2*i1。%擴(kuò)展列△Q m=p+1。J(m,N)=DP。 % PV節(jié)點(diǎn)電壓誤差 for j1=1:n if j1~=isbamp。 % dP/df X5=0。J(p,q)=X5。J(m,N)=DP。 elseif j1==iamp。% dP/df X5=2*e(i)。J(p,q)=X5。J(m,