【正文】 910電壓實(shí)際值230230224226227220 節(jié)點(diǎn)12345678910電壓相角0結(jié)果分析：變電所1的低壓母線電壓與要求的電壓相比略微偏低，需要進(jìn)行無(wú)功補(bǔ)償，將節(jié)點(diǎn)1的無(wú)功功率降低3var；變電所3的低壓母線電壓與要求電壓相比略微偏低，需要進(jìn)行無(wú)功補(bǔ)償，調(diào)節(jié)9節(jié)點(diǎn)的無(wú)功功率，將9節(jié)點(diǎn)的無(wú)功功率降低4var，即符合條件。 0 (50+)* 220 0 0 2。 0 (70+)* 220 0 0 2。 0 0 220 0 0 2。 0 0 220 0 0 2。⑴ 4個(gè)變電所的負(fù)荷同時(shí)以2%比例增大，支路參數(shù)變化及其調(diào)節(jié)方法；支路矩陣B2的變?yōu)椋築2=[0 0 230 230 0 1。由于程序較復(fù)雜，故將其放入附錄中，具體程序和結(jié)果分別見(jiàn)附錄一，二。為了平衡兩發(fā)電廠發(fā)出的電量，令發(fā)電廠二的功率為100 MW，為了減小線路上的損耗。數(shù)據(jù)輸入格式如下：NN——節(jié)點(diǎn)數(shù)，NL——線路數(shù)、ISB——平衡母線節(jié)點(diǎn)號(hào)、PR——誤差精度；支路參數(shù)形成的矩陣B1 1）某之路的首端號(hào)P； 2）末端號(hào)Q；且PQ； 3）支路的阻抗（R+jX）； 4）支路的對(duì)地容抗； 5）支路的變比； 6）折算到哪一側(cè)的標(biāo)志 （如果支路的首端P出入高壓側(cè)輸入“1”，否則輸入“0”）。選擇電壓基準(zhǔn)值為 和功率基準(zhǔn)值，所以 。⑶ 確定完節(jié)點(diǎn)及編號(hào)后，各條支路也相應(yīng)確定了，網(wǎng)絡(luò)中總計(jì)有10條支路，我們對(duì)各支路參數(shù)進(jìn)行了計(jì)算。⑴ 系統(tǒng)的等值電路圖我們首先對(duì)給定的程序輸入部分作了簡(jiǎn)要的分析，程序開(kāi)始需要我們確定輸入節(jié)點(diǎn)數(shù)、支路數(shù)、平衡母線號(hào)、支路參數(shù)矩陣、節(jié)點(diǎn)參數(shù)矩陣。當(dāng)系統(tǒng)的電壓偏離允許值時(shí)，電力系統(tǒng)必須應(yīng)用電壓調(diào)節(jié)技術(shù)調(diào)節(jié)系統(tǒng)電壓的大小，使其維持在允許值范圍內(nèi)。電壓過(guò)高或過(guò)低，都將對(duì)人身及其用電設(shè)備產(chǎn)生重大的影響。7. 最終形成課程設(shè)計(jì)成品說(shuō)明書(shū)。（；電壓35KV要求調(diào)整范圍在3536之間）5. 輪流斷開(kāi)環(huán)網(wǎng)一回線，分析潮流的分布。3. 跟隨變電所負(fù)荷按一定比例發(fā)生變化，進(jìn)行潮流計(jì)算分析。變電所資料：(一) 變電所4低壓母線的電壓等級(jí)分別為：10KV 10KV 35KV 35KV(二) 變電所的負(fù)荷分別為：（1）50MW 50MW 70MW 60MW （2）60MW 80MW 50MW 45MW （3）50MW 40MW 55MW 70MW （4）70MW 50MW 60MW 60MW （5）70MW 60MW 50MW 45MW（6）40MW 50MW 60MW 70MW（7）70MW 40MW 50MW 60MW（8）30MW 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)度的電納為?！獭颉?. 算例電力系統(tǒng)潮流計(jì)算課程設(shè)計(jì)任務(wù)書(shū)一、 題目原始資料系統(tǒng)圖：兩個(gè)發(fā)電廠分別通過(guò)變壓器和輸電線路與四個(gè)變電所相連。)。ylabel(39。支路首端有功39。bar(Siz)。grid on。節(jié)點(diǎn)注入無(wú)功39。)。xlabel(39。subplot(2,2,3)。)。ylabel(39。節(jié)點(diǎn)號(hào)39。plot(sida)。grid on。電壓標(biāo)幺值39。)。xlabel(39。subplot(2,2,1)。)。 disp(39。 disp(ZF)。)=39。,39。DS(39。 DDS(i)=DS(i)*SB。 DS(i)=Si(p,q)+Sj(q,p)。for i=1:nl p=B1(i,1)。各條支路的功率損耗DS為(順序同您輸入B1時(shí)一致)：39。)。 disp(39。 disp(ZF)。)=39。,39。S(39。 SSj(q,p)=Sj(q,p)*SB。 Sjy(i)=Sj(q,p)。 Sjy(i)=Sj(q,p)。q=B1(i,2)。)。enddisp(39。39。 %disp(SSi(p,q))。,num2str(SSi(p,q))]。,num2str(q),39。,num2str(p),39。 ZF=[39。 end disp(Si(p,q))。 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))))。 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))))。for i=1:nl p=B1(i,1)。各條支路的首端功率Si為(順序同您輸入B1時(shí)一致)：39。)。disp(39。SS=S*SB。~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~39。disp(S)。各節(jié)點(diǎn)的功率S為(節(jié)點(diǎn)號(hào)從小到大排列)：39。 end S(p)=E(p)*C(p)。for p=1:n C(p)=0。)。disp(39。39。disp(VV)。disp(V)。各節(jié)點(diǎn)的電壓大小V為(節(jié)點(diǎn)號(hào)從小到大排列)：39。)。disp(39。EE=E*UB。)。end%=============== 計(jì)算各輸出量 ===========================disp(39。 sida(k)=atan(f(k)./e(k))*180./pi。disp(ICT2)。沒(méi)有達(dá)到精度要求的個(gè)數(shù)：39。disp(ICT1)。迭代次數(shù)：39。 ICT1=ICT1+1。 if DET=pr 電壓偏差量是否滿足要求 IT2=IT2+1。 f(L)=f(L)J(k1,N)。 e(L)=e(L)J(k,N)。%消去運(yùn)算 end J(k3,k)=0。%消去運(yùn)算 end J(k3,k)=0。 end if k==N0 break。 for k3=3:k4 用k3行從第三行開(kāi)始到當(dāng)前行前的k4行消去 for k2=k1:N1 k3行后各行下三角元素 J(k3,k2)=J(k3,k2)J(k3,k)*J(k,k2)。 非對(duì)角元規(guī)格化 end J(k,k)=1。N1=N。 J(m,q)=X2。q=q+1。 J(m,q)=X1。J(p,N)=DV。q=2*j11。 X6=2*f(i)。 dP/de X2=D(i)+B(i,i)*e(i)G(i,i)*f(i)。j1~=isb 非平衡節(jié)點(diǎn)amp。 J(m,q)=X2。q=q+1。 J(m,q)=X1。J(p,N)=DV。q=2*j11。X6=0。 dP/de X2=B(i,j1)*e(i)G(i,j1)*f(i)。j1~=i 非平衡節(jié)點(diǎn)amp。 PV節(jié)點(diǎn)有功誤差 DV=V(i)^2V2。 擴(kuò)展列△P J(m,q)=X2。J(p,q)=X4。 J(m,q)=X1。J(p,N)=DQ。q=2*j11。 dQ/de X4=C(i)+G(i,i)*e(i)+B(i,i)*f(i)。 dP/de X2=D(i)+B(i,i)*e(i)G(i,i)*f(i)。j1~=isb 非平衡節(jié)點(diǎn)amp。J(m,q)=X2。q=q+1。 J(m,q)=X1。J(p,N)=DQ。q=2*j11。 % X2=dp/df X3=dQ/de X4=X1。 % dP/de=dQ/df X2=B(i,j1)*e(i)G(i,j1)*f(i)。j1~=i 非平衡節(jié)點(diǎn)amp。 節(jié)點(diǎn)有功功率差 DQ=Q(i)Q1。 V2=e(i)^2+f(i)^2。%節(jié)點(diǎn)功率Q計(jì)算fiΣ(Gij*ejBij*fj)eiΣ(Gij*fj+Bij*ej)求P39。%Σ(Gij*fj+Bij*ej) end P1=C(i)*e(i)+f(i)*D(i)。 for j1=1:n C(i)=C(i)+G(i,j1)*e(j1)B(i,j1)*f(j1)。 for i=1:n if i~=isb 非平衡節(jié)點(diǎn) C(i)=0。while IT2~=0 IT2=0。N=N0+1。IT2=1。Q=imag(S)。 i節(jié)點(diǎn)注入功率SGSL B(i,i)=B(i,i)+B2(i,5)。 V(i)=B2(i,4)。 分解出導(dǎo)納陣的實(shí)部和虛部 for i=1:n 給定各節(jié)點(diǎn)初始電壓的實(shí)部和虛部 e(i)=real(B2(i,3))。disp(Y)G=real(Y)。導(dǎo)納矩陣 Y=39。 對(duì)角元K側(cè) Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2。 非對(duì)角元 Y(q,p)=Y(p,q)。q=B1(i,1)。q=B1(i,2)。)。 end disp(39。 B2(i,4)=B2(i,4)./UB。 B2(i,2)=B2(i,2)./SB。)。 end disp(39。 for i=1:nl B1(i,3)=B1(i,3)*YB。 BB1=B1。)。 UB=input(39。請(qǐng)輸入功率基準(zhǔn)值:SB=39。)。ym=input(39。SB=100。S1=zeros(nl)。V=zeros(1,n)。e=zeros(1,n)。)。 input(39。 0 50+ 220 0 0 2。 0 70+ 220 0 0 2。 0 0 220 0 0 2。 0 0 220 0 0 2。B2=[0 0 230 230 0 1。請(qǐng)輸入由支路參數(shù)形成的矩陣： B1=39。 6 10 + 0 1 1]。 4 8 + 0 1 1。