電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計(jì)---兩機(jī)五節(jié)點(diǎn)網(wǎng)絡(luò)潮流計(jì)算-牛拉法

1、<p>　　電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計(jì)</p><p>　　題 目：兩機(jī)五節(jié)點(diǎn)網(wǎng)絡(luò)潮流計(jì)算-牛拉法</p><p>　　前言 內(nèi)蒙古科技大學(xué)課程設(shè)計(jì)任務(wù)書(shū)………………………………………2</p><p>　　電力系統(tǒng)潮流分布簡(jiǎn)述…………………………………………………5</p><p>　　潮流計(jì)算的發(fā)展趨勢(shì)………………

2、……………………5</p><p>　　潮流計(jì)算的發(fā)展史………………………………………5</p><p>　　潮流計(jì)算的意義…………………………………………5</p><p>　　設(shè)計(jì)程序…………………………………………………………………7</p><p>　　潮流計(jì)算的基本步驟…………………………………………7</p><

3、;p>　　P-Q法潮流計(jì)算流程圖………………………………………9</p><p>　　MATLAB程序 ………………………………………10</p><p>　　潮流計(jì)算程序運(yùn)行結(jié)果示……………………………………14</p><p>　　總結(jié)……………………………………………………………………………</p><p>　　參

4、考……………………………………………………………………………</p><p>　　內(nèi)蒙古科技大學(xué)課程設(shè)計(jì)任務(wù)書(shū)</p><p><b>　　系統(tǒng)接線圖</b></p><p>　　其中節(jié)點(diǎn)1為平衡節(jié)點(diǎn)，節(jié)點(diǎn)2、3、4、5為PQ節(jié)點(diǎn)。</p><p>　　第一章 電力系統(tǒng)潮流計(jì)算簡(jiǎn)述</p><

5、;p>　　第一節(jié) 潮流計(jì)算的意義</p><p>　　(1)在電網(wǎng)規(guī)劃階段,通過(guò)潮流計(jì)算,合理規(guī)劃電源容量及接入點(diǎn),合理規(guī)劃網(wǎng)架,選擇無(wú)功補(bǔ)償方案,滿足規(guī)劃水平的大、小方式下潮流交換控制、調(diào)峰、調(diào)相、調(diào)壓的要求。 </p><p>　　(2)在編制年運(yùn)行方式時(shí),在預(yù)計(jì)負(fù)荷增長(zhǎng)及新設(shè)備投運(yùn)基礎(chǔ)上,選擇典型方式進(jìn)行潮流計(jì)算,發(fā)現(xiàn)電網(wǎng)中薄弱環(huán)節(jié),供調(diào)度員日常調(diào)度控制參考,并對(duì)規(guī)劃、基建部門(mén)

6、提出改進(jìn)網(wǎng)架結(jié)構(gòu),加快基建進(jìn)度的建議。 </p><p>　　(3)正常檢修及特殊運(yùn)行方式下的潮流計(jì)算,用于日運(yùn)行方式的編制,指導(dǎo)發(fā)電廠開(kāi)機(jī)方式,有功、無(wú)功調(diào)整方案及負(fù)荷調(diào)整方案,滿足線路、變壓器熱穩(wěn)定要求及電壓質(zhì)量要求。 </p><p>　　(4)預(yù)想事故、設(shè)備退出運(yùn)行對(duì)靜態(tài)安全的影響分析及作出預(yù)想的運(yùn)行方式調(diào)整方案。 </p><p>　　總結(jié)為在電力系統(tǒng)運(yùn)行

7、方式和規(guī)劃方案的研究中，都需要進(jìn)行潮流計(jì)算以比較運(yùn)行方式或規(guī)劃供電方案的可行性、可靠性和經(jīng)濟(jì)性。同時(shí)，為了實(shí)時(shí)監(jiān)控電力系統(tǒng)的運(yùn)行狀態(tài)，也需要進(jìn)行大量而快速的潮流計(jì)算。因此，潮流計(jì)算是電力系統(tǒng)中應(yīng)用最廣泛、最基本和最重要的一種電氣運(yùn)算。在系統(tǒng)規(guī)劃設(shè)計(jì)和安排系統(tǒng)的運(yùn)行方式時(shí)，采用離線潮流計(jì)算；在電力系統(tǒng)運(yùn)行狀態(tài)的實(shí)時(shí)監(jiān)控中，則采用在線潮流計(jì)算。 </p><p>　　第二節(jié) 潮流計(jì)算的發(fā)展史</p>&

8、lt;p>　　利用電子計(jì)算機(jī)進(jìn)行潮流計(jì)算從20世紀(jì)50年代中期就已經(jīng)開(kāi)始。此后，潮流計(jì)算曾采用了各種不同的方法，這些方法的發(fā)展主要是圍繞著對(duì)潮流計(jì)算的一些基本要求進(jìn)行的。對(duì)潮流計(jì)算的要求可以歸納為下面幾點(diǎn)： </p><p>　?。?）算法的可靠性或收斂性 </p><p>　?。?）計(jì)算速度和內(nèi)存占用量 </p><p>　　（3）計(jì)算的方便性和靈活性 &l

9、t;/p><p>　　電力系統(tǒng)潮流計(jì)算屬于穩(wěn)態(tài)分析范疇，不涉及系統(tǒng)元件的動(dòng)態(tài)特性和過(guò)渡過(guò)程。因此其數(shù)學(xué)模型不包含微分方程，是一組高階非線性方程。非線性代數(shù)方程組的解法離不開(kāi)迭代，因此，潮流計(jì)算方法首先要求它是能可靠的收斂，并給出正確答案。隨著電力系統(tǒng)規(guī)模的不斷擴(kuò)大，潮流問(wèn)題的方程式階數(shù)越來(lái)越高，目前已達(dá)到幾千階甚至上萬(wàn)階，對(duì)這樣規(guī)模的方程式并不是采用任何數(shù)學(xué)方法都能保證給出正確答案的。這種情況促使電力系統(tǒng)的研究人員不

10、斷尋求新的更可靠的計(jì)算方法。 </p><p>　　在用數(shù)字計(jì)算機(jī)求解電力系統(tǒng)潮流問(wèn)題的開(kāi)始階段，人們普遍采用以節(jié)點(diǎn)導(dǎo)納矩陣為基礎(chǔ)的高斯-賽德?tīng)柕?。這個(gè)方法的原理比較簡(jiǎn)單，要求的數(shù)字計(jì)算機(jī)的內(nèi)存量也比較小，適應(yīng)當(dāng)時(shí)的電子數(shù)字計(jì)算機(jī)制作水平和電力系統(tǒng)理論水平，于是電力系統(tǒng)計(jì)算人員轉(zhuǎn)向以阻抗矩陣為主的逐次代入法。 </p><p>　　20世紀(jì)60年代初，數(shù)字計(jì)算機(jī)已經(jīng)發(fā)展到第二代，計(jì)算機(jī)

11、的內(nèi)存和計(jì)算速度發(fā)生了很大的飛躍，從而為阻抗法的采用創(chuàng)造了條件。阻抗矩陣是滿矩陣，阻抗法要求計(jì)算機(jī)儲(chǔ)存表征系統(tǒng)接線和參數(shù)的阻抗矩陣。這就需要較大的內(nèi)存量。而且阻抗法每迭代一次都要求順次取阻抗矩陣中的每一個(gè)元素進(jìn)行計(jì)算，因此，每次迭代的計(jì)算量很大。 </p><p>　　阻抗法改善了電力系統(tǒng)潮流計(jì)算問(wèn)題的收斂性，解決了導(dǎo)納法無(wú)法解決的一些系統(tǒng)的潮流計(jì)算，在當(dāng)時(shí)獲得了廣泛的應(yīng)用，曾為我國(guó)電力系統(tǒng)設(shè)計(jì)、運(yùn)行和研究作出了

12、很大的貢獻(xiàn)。但是，阻抗法的主要缺點(diǎn)就是占用計(jì)算機(jī)的內(nèi)存很大，每次迭代的計(jì)算量很大。當(dāng)系統(tǒng)不斷擴(kuò)大時(shí)，這些缺點(diǎn)就更加突出。為了克服阻抗法在內(nèi)存和速度方面的缺點(diǎn)，后來(lái)發(fā)展了以阻抗矩陣為基礎(chǔ)的分塊阻抗法。這個(gè)方法把一個(gè)大系統(tǒng)分割為幾個(gè)小的地區(qū)系統(tǒng)，在計(jì)算機(jī)內(nèi)只需存儲(chǔ)各個(gè)地區(qū)系統(tǒng)的阻抗矩陣及它們之間的聯(lián)絡(luò)線的阻抗，這樣不僅大幅度的節(jié)省了內(nèi)存容量，同時(shí)也提高了計(jì)算速度。 </p><p>　　克服阻抗法缺點(diǎn)的另一途徑是采用

13、牛頓-拉夫遜法（以下簡(jiǎn)稱牛頓法）。牛頓法是數(shù)學(xué)中求解非線性方程式的典型方法，有較好的收斂性。解決電力系統(tǒng)潮流計(jì)算問(wèn)題是以導(dǎo)納矩陣為基礎(chǔ)的，因此，只要在迭代過(guò)程中盡可能保持方程式系數(shù)矩陣的稀疏性，就可以大大提高牛頓潮流程序的計(jì)算效率。自從20世紀(jì)60年代中期采用了最佳順序消去法以后，牛頓法在收斂性、內(nèi)存要求、計(jì)算速度方面都超過(guò)了阻抗法，成為直到目前仍被廣泛采用的方法。 </p><p>　　在牛頓法的基礎(chǔ)上，根據(jù)電

14、力系統(tǒng)的特點(diǎn)，抓住主要矛盾，對(duì)純數(shù)學(xué)的牛頓法進(jìn)行了改造，得到了P-Q分解法。P-Q分解法在計(jì)算速度方面有顯著的提高，迅速得到了推廣。 </p><p>　　牛頓法的特點(diǎn)是將非線性方程線性化。20世紀(jì)70年代后期，有人提出采用更精確的模型，即將泰勒級(jí)數(shù)的高階項(xiàng)也包括進(jìn)來(lái)，希望以此提高算法的性能，這便產(chǎn)生了保留非線性的潮流算法。另外，為了解決病態(tài)潮流計(jì)算，出現(xiàn)了將潮流計(jì)算表示為一個(gè)無(wú)約束非線性規(guī)劃問(wèn)題的模型，即非線性

15、規(guī)劃潮流算法。 </p><p>　　近20多年來(lái)，潮流算法的研究仍然非?；钴S，但是大多數(shù)研究都是圍繞改進(jìn)牛頓法和P-Q分解法進(jìn)行的。此外，隨著人工智能理論的發(fā)展，遺傳算法、人工神經(jīng)網(wǎng)絡(luò)、模糊算法也逐漸被引入潮流計(jì)算。但是，到目前為止這些新的模型和算法還不能取代牛頓法和P-Q分解法的地位。由于電力系統(tǒng)規(guī)模的不斷擴(kuò)大，對(duì)計(jì)算速度的要求不斷提高，計(jì)算機(jī)的并行計(jì)算技術(shù)也將在潮流計(jì)算中得到廣泛的應(yīng)用，成為重要的研究領(lǐng)域。

16、 </p><p>　　第三節(jié) 潮流計(jì)算的發(fā)展趨勢(shì)</p><p>　　通過(guò)幾十年的發(fā)展，潮流算法日趨成熟。近幾年，對(duì)潮流算法的研究仍然是如何改善傳統(tǒng)的潮流算法，即高斯-塞德?tīng)柗ā⑴ｎD法和快速解耦法。牛頓法，由于其在求解非線性潮流方程時(shí)采用的是逐次線性化的方法，為了進(jìn)一步提高算法的收斂性和計(jì)算速度，人們考慮采用將泰勒級(jí)數(shù)的高階項(xiàng)或非線性項(xiàng)也考慮進(jìn)來(lái)，于是產(chǎn)生了二階潮流算法。后來(lái)又提出了根據(jù)

17、直角坐標(biāo)形式的潮流方程是一個(gè)二次代數(shù)方程的特點(diǎn)，提出了采用直角坐標(biāo)的保留非線性快速潮流算法。 </p><p>　　對(duì)于一些病態(tài)系統(tǒng)，應(yīng)用非線性潮流計(jì)算方法往往會(huì)造成計(jì)算過(guò)程的振蕩或者不收斂，從數(shù)學(xué)上講，非線性的潮流計(jì)算方程組本來(lái)就是無(wú)解的。這樣，人們提出來(lái)了將潮流方程構(gòu)造成一個(gè)函數(shù)，求此函數(shù)的最小值問(wèn)題，稱之為非線性規(guī)劃潮流的計(jì)算方法。優(yōu)點(diǎn)是原理上保證了計(jì)算過(guò)程永遠(yuǎn)不會(huì)發(fā)散。如果將數(shù)學(xué)規(guī)劃原理和牛頓潮流算法

18、有機(jī)結(jié)合一起就是最優(yōu)乘子法。另外，為了優(yōu)化系統(tǒng)的運(yùn)行，從所有以上的可行潮流解中挑選出滿足一定指標(biāo)要求的一個(gè)最佳方案就是最優(yōu)潮流問(wèn)題。最優(yōu)潮流是一種同時(shí)考慮經(jīng)濟(jì)性和安全性的電力網(wǎng)絡(luò)分析優(yōu)化問(wèn)題。OPF 在電力系統(tǒng)的安全運(yùn)行、經(jīng)濟(jì)調(diào)度、可靠性分析、能量管理以及電力定價(jià)等方面得到了廣泛的應(yīng)用。</p><p><b>　　程序設(shè)計(jì)</b></p><p>　　第一節(jié) 潮流計(jì)

19、算的基本步驟</p><p>　　1.形成節(jié)點(diǎn)導(dǎo)納矩陣YB</p><p>　　由上圖可得相應(yīng)的節(jié)點(diǎn)導(dǎo)納矩陣</p><p>　　由于節(jié)點(diǎn)1為平衡節(jié)點(diǎn)，其他的節(jié)點(diǎn)均為PQ節(jié)點(diǎn)，系數(shù)矩陣，階數(shù)相同。又應(yīng)對(duì)該等值網(wǎng)絡(luò)，不存在除去與有功功率和電壓相位或無(wú)功功率和電壓大小關(guān)系較小的因素的可能性，這兩個(gè)矩陣，完全相同。他們就由導(dǎo)納矩陣的虛部</p><p&

20、gt;　　部分中除第一行第一列外的各個(gè)元素所組成.</p><p>　　2．計(jì)算各節(jié)點(diǎn)有功功率不平衡量</p><p>　　取=1.06，=0；====1.0;= ===0,</p><p>　　計(jì)算各節(jié)點(diǎn)有功功率不平衡量</p><p><b>　　=-</b></p><p>　　各節(jié)點(diǎn)電壓的

21、相位角(弧度)</p><p><b>　　由矩陣方程式</b></p><p><b>　　-()=</b></p><p>　　3．計(jì)算各節(jié)點(diǎn)無(wú)功功率不平衡量</p><p>　　按下式計(jì)算各節(jié)點(diǎn)無(wú)功功率不平衡量</p><p><b>　　=- <

22、/b></p><p>　　由矩陣方程式計(jì)算各節(jié)點(diǎn)電壓大小</p><p><b>　　-=</b></p><p>　　4．計(jì)算平衡節(jié)點(diǎn)功率和線路功率 </p><p><b>　　==()</b></p><p>　　第二節(jié) P-Q分解法的程序流程圖：</p

23、><p><b>　　MATLAB程序</b></p><p><b>　　clc</b></p><p><b>　　clear</b></p><p>　　disp('節(jié)點(diǎn)總數(shù)為：');</p><p><b>　　N=5<

24、;/b></p><p>　　disp('平衡節(jié)點(diǎn)為：');</p><p><b>　　1</b></p><p>　　disp('PQ節(jié)點(diǎn)為：');</p><p>　　JD=[2,3,4,5]</p><p>　　Y=[6.25-18.75i -5

25、+15i -1.25+3.75i 0 0;</p><p>　　-5+15i 10.834-32.5i -1.667+5i -1.667+5i -2.5+7.5i;</p><p>　　-1.25+3.75i -1.667+5i 12.917-38.75i -10+30i 0;</p&

26、gt;<p>　　0 -1.667+5i -10+30i 12.917-38.75i -1.25+3.75i;</p><p>　　0 -2.5+7.5i 0 -1.25+3.75i 3.75-11.25i]</p><p>　　G=real(Y);B=imag(Y);</p><p&

27、gt;　　B1=B(2:5, 2:5);</p><p>　　B2=B1; </p><p>　　b1=inv(B1);</p><p>　　b2=inv(B2);</p><p>　　dt(2)=0; </p><p><b>　　dt(

28、3)=0; </b></p><p><b>　　dt(4)=0;</b></p><p><b>　　dt(5)=0;</b></p><p><b>　　u(2)=1.0;</b></p><p><b>　　u(3)=1.0;</b>&l

29、t;/p><p><b>　　u(4)=1.0;</b></p><p><b>　　u(5)=1.0;</b></p><p>　　p(2)=0.20;</p><p>　　p(3)=-0.45;</p><p>　　p(4)=-0.40;</p><p&g

30、t;　　p(5)=-0.60;</p><p>　　q(2)=0.20;</p><p>　　q(3)=-0.15;</p><p>　　q(4)=-0.05;</p><p>　　q(5)=-0.10;</p><p>　　k=0;wucha=1;</p><p><b>　　N1=

31、4;</b></p><p>　　while wucha(1)>0.00001;</p><p>　　dt(1)=0;u(1)=1.06;</p><p><b>　　for m=2:5</b></p><p><b>　　for n=1:5</b></p><

32、p>　　pt(n)=u(m)*u(n)*(G(m,n)*cos(dt(m)-dt(n))+B(m,n)*sin(dt(m)-dt(n))); </p><p><b>　　end</b></p><p>　　disp('各節(jié)點(diǎn)有功率不平衡量為：') </p><p>　　dp(m)=p(m)-sum(pt)</p&

33、gt;<p><b>　　end</b></p><p><b>　　for m=2:5</b></p><p>　　u1=diag(u);</p><p>　　u2=inv(u1(2:5,2:5));</p><p>　　u3=flipud(rot90(u(2:5)));</p

34、><p>　　u4=u3(1:4,:);</p><p>　　dp1=flipud(rot90(dp(2:5)));</p><p>　　dp2=dp1(1:4,:);</p><p>　　dtt=(-b1*dp2/u4)*u2;</p><p>　　dtt=dtt(:,1:1);</p><p>

35、　　dtt=flipud(rot90(dtt))</p><p>　　disp('各節(jié)點(diǎn)電壓的相位角dt為：') </p><p>　　dt(m)=dt(m)+dtt(m-1)</p><p><b>　　end</b></p><p><b>　　for m=2:5</b><

36、/p><p><b>　　for n=1:5</b></p><p>　　qt(n)=u(m)*u(n)*(G(m,n)*sin(dt(m)-dt(n))-B(m,n)*cos(dt(m)-dt(n)));</p><p><b>　　end </b></p><p>　　disp('各節(jié)點(diǎn)無(wú)功

37、率不平衡量為：')</p><p>　　qq(m)=q(m)-sum(qt)</p><p><b>　　end</b></p><p><b>　　for m=2:5</b></p><p>　　qq1=flipud(rot90(qq(2:5)));</p><p>

38、;　　qq2=qq1(1:4,:);</p><p>　　ut=-b2*qq2/u4;</p><p>　　disp( '各節(jié)點(diǎn)電壓的修正量：') </p><p><b>　　ut</b></p><p>　　disp( '各節(jié)點(diǎn)電壓的大小為：') </p><p

39、>　　u(m)=u(m)+ut(m-1) </p><p><b>　　end</b></p><p><b>　　for n=1:5</b></p><p>　　U(n)=u(n)*(cos(dt(n))+j*sin(dt(n)));</p><p><b>　　end</b

40、></p><p><b>　　for n=1:5</b></p><p>　　I(n)=Y(1,n)*U(n);</p><p><b>　　end</b></p><p>　　disp('各節(jié)點(diǎn)功率s為：')</p><p>　　S1=U(1)*su

41、m(conj(I))</p><p><b>　　for m=1:5</b></p><p><b>　　for n=1:5</b></p><p>　　S(m,n)=U(m)*(conj(U(m))-conj(U(n)))*conj(-Y(m,n));</p><p><b>　　end

42、</b></p><p><b>　　end</b></p><p>　　disp('各節(jié)點(diǎn)功率sij為：')</p><p><b>　　S</b></p><p>　　wucha=max(abs(ut))</p><p><b>　　

43、k=k+1</b></p><p><b>　　End</b></p><p>　　第四節(jié) 程序運(yùn)行結(jié)果顯示</p><p><b>　　節(jié)點(diǎn)總數(shù)為：</b></p><p><b>　　N = 5</b></p><p><b>

44、;　　平衡節(jié)點(diǎn)為：</b></p><p><b>　　ans = 1</b></p><p><b>　　PQ節(jié)點(diǎn)為： </b></p><p>　　JD = 2 3 4 5</p><p>　　Y =6.2500 -18.7500i -5.0000 +15.

45、0000i -1.2500 + 3.7500i 0 0 </p><p>　　-5.0000 +15.0000i 10.8340 -32.5000i -1.6670 + 5.0000i -1.6670 + 5.0000i -2.5000 + 7.5000i</p><p>　　-1.2500 + 3.7500i

46、-1.6670 + 5.0000i 12.9170 -38.7500i -10.0000 +30.0000i 0 </p><p>　　0 -1.6670 + 5.0000i -10.0000 +30.0000i 12.9170 -38.7500i -1.2500 + 3.7500i</p><p>　　0

47、-2.5000 + 7.5000i 0 -1.2500 + 3.7500i 3.7500 -11.2500i</p><p>　　各節(jié)點(diǎn)有功率不平衡量為：</p><p>　　dp =0 0.5000</p><p>　　各節(jié)點(diǎn)有功率不平衡量為：</p><p>　　dp = 0 0.500

48、0 -0.3750</p><p>　　各節(jié)點(diǎn)有功率不平衡量為：</p><p>　　dp = 0 0.5000 -0.3750 -0.4000</p><p>　　各節(jié)點(diǎn)有功率不平衡量為：</p><p>　　dp = 0 0.5000 -0.3750 -0.4000 -0.6000</p>

49、<p>　　dtt =-0.0370 -0.0855 -0.0918 -0.1086</p><p>　　各節(jié)點(diǎn)電壓的相位角dt為：</p><p>　　dt =0 -0.0370 0 0 0</p><p>　　dtt =-0.0370 -0.0855 -0.0918 -0.1

50、086</p><p>　　各節(jié)點(diǎn)電壓的相位角dt為：</p><p>　　dt = 0 -0.0370 -0.0855 0 0</p><p>　　dtt =-0.0370 -0.0855 -0.0918 -0.1086</p><p>　　各節(jié)點(diǎn)電壓的相位角dt為：</p>

51、<p>　　dt = 0 -0.0370 -0.0855 -0.0918 0</p><p>　　dtt =-0.0370 -0.0855 -0.0918 -0.1086</p><p>　　各節(jié)點(diǎn)電壓的相位角dt為：</p><p>　　dt = 0 -0.0370 -0.0855 -0.0918

52、 -0.1086</p><p>　　各節(jié)點(diǎn)無(wú)功率不平衡量為：</p><p>　　qq = 0 1.2119</p><p>　　各節(jié)點(diǎn)無(wú)功率不平衡量為：</p><p>　　qq = 0 1.2119 -0.0773</p><p>　　各節(jié)點(diǎn)無(wú)功率不平衡量為：</p><p&

53、gt;　　qq = 0 1.2119 -0.0773 -0.1919</p><p>　　各節(jié)點(diǎn)無(wú)功率不平衡量為：</p><p>　　qq =0 1.2119 -0.0773 -0.1919 -0.3196</p><p>　　各節(jié)點(diǎn)電壓的修正量：</p><p><b>　　ut =</b&

54、gt;</p><p>　　0.0395 0 0 0</p><p>　　0.0081 0 0 0</p><p>　　0.0064 0 0 0</p><p>　　0.0001 0

55、 0 0</p><p>　　各節(jié)點(diǎn)電壓的大小為：</p><p>　　u =1.0600 1.0395 1.0000 1.0000 1.0000</p><p>　　各節(jié)點(diǎn)電壓的修正量：</p><p><b>　　ut =</b></p><p

56、>　　0.0395 0 0 0</p><p>　　0.0081 0 0 0</p><p>　　0.0064 0 0 0</p><p>　　0.0001 0 0 0&

57、lt;/p><p>　　各節(jié)點(diǎn)電壓的大小為：</p><p>　　u = 1.0600 1.0395 1.0081 1.0000 1.0000</p><p>　　各節(jié)點(diǎn)電壓的修正量：</p><p>　　ut =0.0395 0 0 0</p><p&g

58、t;　　0.0081 0 0 0</p><p>　　0.0064 0 0 0</p><p>　　0.0001 0 0 0</p><p>　　各節(jié)點(diǎn)電壓的大小為：</p><p>　　u = 1.

59、0600 1.0395 1.0081 1.0064 1.0000</p><p>　　各節(jié)點(diǎn)電壓的修正量：</p><p>　　ut =0.0395 0 0 0</p><p>　　0.0081 0 0 0</p><p>

60、　　0.0064 0 0 0</p><p>　　0.0001 0 0 0</p><p>　　各節(jié)點(diǎn)電壓的大小為：</p><p>　　u =1.0600 1.0395 1.0081 1.0064 1.0001</p><p&

61、gt;<b>　　各節(jié)點(diǎn)功率s為：</b></p><p>　　S1 = 1.1389 + 0.2403i</p><p>　　各節(jié)點(diǎn)功率sij為：</p><p>　　S = 0 0.7229 + 0.1332i 0.4160 + 0.1071i 0 0

62、 </p><p>　　-0.7133 - 0.1044i 0 0.3110 + 0.0850i 0.3469 + 0.0845i 0.6671 + 0.1416i</p><p>　　-0.4029 - 0.0676i -0.3052 - 0.0677i 0 0.2083 - 0.0128i

63、 0 </p><p>　　0 -0.3398 - 0.0633i -0.2078 + 0.0141i 0 0.0714 + 0.0033i</p><p>　　0 -0.6499 - 0.0900i 0 -0.0710 - 0.0021i

64、 0 </p><p><b>　　k = 1</b></p><p>　　各節(jié)點(diǎn)有功率不平衡量為：</p><p>　　dp = 0 -0.4117 -0.3750 -0.4000 -0.6000</p><p>　　各節(jié)點(diǎn)有功率不平衡量為：</p&g

65、t;<p>　　dp =0 -0.4117 0.0498 -0.4000 -0.6000</p><p>　　各節(jié)點(diǎn)有功率不平衡量為：</p><p>　　dp = 0 -0.4117 0.0498 0.0762 -0.6000</p><p>　　各節(jié)點(diǎn)有功率不平衡量為：</p><p>

66、;　　dp = 0 -0.4117 0.0498 0.0762 0.1209</p><p>　　dtt = -0.0105 0.0012 0.0017 0.0035</p><p>　　各節(jié)點(diǎn)電壓的相位角dt為：</p><p>　　dt = 0 -0.0474 -0.0855 -0.0918

67、 -0.1086</p><p>　　dtt = -0.0105 0.0012 0.0017 0.0035</p><p>　　各節(jié)點(diǎn)電壓的相位角dt為：</p><p>　　dt = 0 -0.0474 -0.0843 -0.0918 -0.1086</p><p>　　dtt = -0.

68、0105 0.0012 0.0017 0.0035</p><p>　　各節(jié)點(diǎn)電壓的相位角dt為：</p><p>　　dt = 0 -0.0474 -0.0843 -0.0901 -0.1086</p><p>　　dtt =-0.0105 0.0012 0.0017 0.0035</p><

69、;p>　　各節(jié)點(diǎn)電壓的相位角dt為：</p><p>　　dt = 0 -0.0474 -0.0843 -0.0901 -0.1050</p><p>　　各節(jié)點(diǎn)無(wú)功率不平衡量為：</p><p>　　qq = 0 -0.1367 -0.0773 -0.1919 -0.3196</p><p>　　

70、各節(jié)點(diǎn)無(wú)功率不平衡量為：</p><p>　　qq = 0 -0.1367 0.0175 -0.1919 -0.3196</p><p>　　各節(jié)點(diǎn)無(wú)功率不平衡量為：</p><p>　　qq = 0 -0.1367 0.0175 0.0237 -0.3196</p><p>　　各節(jié)點(diǎn)無(wú)功率不平衡量為

71、：</p><p>　　qq =0 -0.1367 0.0175 0.0237 0.0378</p><p>　　各節(jié)點(diǎn)電壓的修正量：</p><p>　　ut = -0.0038 0 0 0</p><p>　　0.0002 0 0

72、 0</p><p>　　0.0004 0 0 0</p><p>　　0.0009 0 0 0</p><p>　　各節(jié)點(diǎn)電壓的大小為：</p><p>　　u = 1.0600 1.0358 1.0081 1.0064

73、 1.0001</p><p>　　各節(jié)點(diǎn)電壓的修正量：</p><p>　　ut =-0.0038 0 0 0</p><p>　　0.0002 0 0 0</p><p>　　0.0004 0 0

74、 0</p><p>　　0.0009 0 0 0</p><p>　　各節(jié)點(diǎn)電壓的大小為：</p><p>　　u =1.0600 1.0358 1.0083 1.0064 1.0001</p><p>　　各節(jié)點(diǎn)電壓的修正量：</p><p&

75、gt;　　ut = -0.0038 0 0 0</p><p>　　0.0002 0 0 0</p><p>　　0.0004 0 0 0</p><p>　　0.0009 0 0

76、 0</p><p>　　各節(jié)點(diǎn)電壓的大小為：</p><p>　　u = 1.0600 1.0358 1.0083 1.0068 1.0001</p><p>　　各節(jié)點(diǎn)電壓的修正量：</p><p>　　ut =-0.0038 0 0 0</p>&

77、lt;p>　　0.0002 0 0 0</p><p>　　0.0004 0 0 0</p><p>　　0.0009 0 0 0</p><p>　　各節(jié)點(diǎn)電壓的大小為：</p><p>　　

78、u =1.0600 1.0358 1.0083 1.0068 1.0009</p><p><b>　　各節(jié)點(diǎn)功率s為：</b></p><p>　　S1 =1.3260 + 0.2508i</p><p>　　各節(jié)點(diǎn)功率sij為：</p><p>　　S = 0 0.9

79、151 + 0.1435i 0.4109 + 0.1072i 0 0 </p><p>　　-0.8998 - 0.0977i 0 0.2415 + 0.0816i 0.2741 + 0.0809i 0.5423 + 0.1344i</p><p>　　-0.3981 - 0

80、.0687i -0.2378 - 0.0707i 0 0.1903 - 0.0114i 0 </p><p>　　0 -0.2695 - 0.0672i -0.1899 + 0.0125i 0 0.0639 + 0.0036i</p><p>　　0

81、 -0.5306 - 0.0995i 0 -0.0636 - 0.0026i 0 </p><p><b>　　k =2</b></p><p>　　各節(jié)點(diǎn)有功率不平衡量為：</p><p>　　dp = 0 0.0419 0.0498 0.

82、0762 0.1209</p><p>　　各節(jié)點(diǎn)有功率不平衡量為：</p><p>　　dp = 0 0.0419 -0.0044 0.0762 0.1209</p><p>　　各節(jié)點(diǎn)有功率不平衡量為：</p><p>　　dp = 0 0.0419 -0.0044 -0.0044 0.12

83、09</p><p>　　各節(jié)點(diǎn)有功率不平衡量為：</p><p>　　dp = 0 0.0419 -0.0044 -0.0044 -0.0058</p><p>　　dtt =0.0016 0.0006 0.0006 0.0008</p><p>　　各節(jié)點(diǎn)電壓的相位角dt為：</p>

84、<p>　　dt = 0 -0.0459 -0.0843 -0.0901 -0.1050</p><p>　　dtt =0.0016 0.0006 0.0006 0.0008</p><p>　　各節(jié)點(diǎn)電壓的相位角dt為：</p><p>　　dt = 0 -0.0459 -0.0838 -0.0901

85、-0.1050</p><p>　　dtt =0.0016 0.0006 0.0006 0.0008</p><p>　　各節(jié)點(diǎn)電壓的相位角dt為：</p><p>　　dt = 0 -0.0459 -0.0838 -0.0895 -0.1050</p><p>　　dtt =0.0016 0.000

86、6 0.0006 0.0008</p><p>　　各節(jié)點(diǎn)電壓的相位角dt為：</p><p>　　dt = 0 -0.0459 -0.0838 -0.0895 -0.1043</p><p>　　各節(jié)點(diǎn)無(wú)功率不平衡量為：</p><p>　　qq = 0 0.0153 0.0175 0.02

87、37 0.0378</p><p>　　各節(jié)點(diǎn)無(wú)功率不平衡量為：</p><p>　　qq = 0 0.0153 -0.0005 0.0237 0.0378</p><p>　　各節(jié)點(diǎn)無(wú)功率不平衡量為：</p><p>　　qq = 0 0.0153 -0.0005 -0.0006 0.03

88、78</p><p>　　各節(jié)點(diǎn)無(wú)功率不平衡量為：</p><p>　　qq = 0 0.0153 -0.0005 -0.0006 -0.0001</p><p>　　各節(jié)點(diǎn)電壓的修正量：</p><p><b>　　ut =</b></p><p>　　1.0e-003 *

89、</p><p>　　0.7776 0 0 0</p><p>　　0.5118 0 0 0</p><p>　　0.5491 0 0 0</p><p>　　0.6921 0

90、 0 0</p><p>　　各節(jié)點(diǎn)電壓的大小為：</p><p>　　u = 1.0600 1.0366 1.0083 1.0068 1.0009</p><p>　　各節(jié)點(diǎn)電壓的修正量：</p><p>　　ut =1.0e-003 *</p><p>　　0.7

91、776 0 0 0</p><p>　　0.5118 0 0 0</p><p>　　0.5491 0 0 0</p><p>　　0.6921 0 0 0</p>

92、<p>　　各節(jié)點(diǎn)電壓的大小為：</p><p>　　u = 1.0600 1.0366 1.0088 1.0068 1.0009</p><p>　　各節(jié)點(diǎn)電壓的修正量：</p><p>　　ut = 1.0e-003 *</p><p>　　0.7776 0 0

93、 0</p><p>　　0.5118 0 0 0</p><p>　　0.5491 0 0 0</p><p>　　0.6921 0 0 0</p><p>　　各節(jié)點(diǎn)電壓的大小為：&

94、lt;/p><p>　　u = 1.0600 1.0366 1.0088 1.0073 1.0009</p><p>　　各節(jié)點(diǎn)電壓的修正量：</p><p>　　ut = 1.0e-003 *</p><p>　　0.7776 0 0 0</p><

95、;p>　　0.5118 0 0 0</p><p>　　0.5491 0 0 0</p><p>　　0.6921 0 0 0</p><p>　　各節(jié)點(diǎn)電壓的大小為：</p><p>　　u

96、= 1.0600 1.0366 1.0088 1.0073 1.0016</p><p><b>　　各節(jié)點(diǎn)功率s為：</b></p><p>　　S1 =1.2936 + 0.2441i</p><p>　　各節(jié)點(diǎn)功率sij為：</p><p>　　S = 0 0.

97、8855 + 0.1383i 0.4081 + 0.1057i 0 0 </p><p>　　-0.8712 - 0.0955i 0 0.2475 + 0.0815i 0.2799 + 0.0806i 0.5495 + 0.1333i</p><p>　　-0.3954 -

98、0.0678i -0.2437 - 0.0701i 0 0.1889 - 0.0122i 0 </p><p>　　0 -0.2751 - 0.0664i -0.1885 + 0.0132i 0 0.0632 + 0.0033i</p><p>　　0

99、 -0.5376 - 0.0976i 0 -0.0629 - 0.0023i 0 </p><p><b>　　k = 3</b></p><p>　　各節(jié)點(diǎn)有功率不平衡量為：</p><p>　　dp = 0 -0.0057 -0.0044 -0.

100、0044 -0.0058</p><p>　　各節(jié)點(diǎn)有功率不平衡量為：</p><p>　　dp = 0 -0.0057 0.0003 -0.0044 -0.0058</p><p>　　各節(jié)點(diǎn)有功率不平衡量為：</p><p>　　dp = 0 -0.0057 0.0003 0.0004 -0.0

101、058</p><p>　　各節(jié)點(diǎn)有功率不平衡量為：</p><p>　　dp = 0 -0.0057 0.0003 0.0004 0.0005</p><p>　　dtt = 1.0e-003 *</p><p>　　-0.2442 -0.1346 -0.1413 -0.1674</p>

102、<p>　　各節(jié)點(diǎn)電壓的相位角dt為：</p><p>　　dt = 0 -0.0461 -0.0838 -0.0895 -0.1043</p><p>　　dtt =1.0e-003 *</p><p>　　-0.2442 -0.1346 -0.1413 -0.1674</p><p>　　各節(jié)點(diǎn)電壓的

103、相位角dt為：</p><p>　　dt = 0 -0.0461 -0.0839 -0.0895 -0.1043</p><p>　　dtt = 1.0e-003 *</p><p>　　-0.2442 -0.1346 -0.1413 -0.1674</p><p>　　各節(jié)點(diǎn)電壓的相位角dt為：</p&g

104、t;<p>　　dt = 0 -0.0461 -0.0839 -0.0896 -0.1043</p><p>　　dtt = 1.0e-003 *</p><p>　　-0.2442 -0.1346 -0.1413 -0.1674</p><p>　　各節(jié)點(diǎn)電壓的相位角dt為：</p><p>　　d

105、t =0 -0.0461 -0.0839 -0.0896 -0.1044</p><p>　　各節(jié)點(diǎn)無(wú)功率不平衡量為：</p><p>　　qq = 0 -0.0020 -0.0005 -0.0006 -0.0001</p><p>　　各節(jié)點(diǎn)無(wú)功率不平衡量為：</p><p>　　qq = 0 -0.0

106、020 0.0001 -0.0006 -0.0001</p><p>　　各節(jié)點(diǎn)無(wú)功率不平衡量為：</p><p>　　qq = 0 -0.0020 0.0001 0.0001 -0.0001</p><p>　　各節(jié)點(diǎn)無(wú)功率不平衡量為：</p><p>　　qq = 0 -0.0020 0.000

107、1 0.0001 0.0001</p><p>　　各節(jié)點(diǎn)電壓的修正量：</p><p><b>　　ut =</b></p><p>　　1.0e-004 *</p><p>　　-0.9362 0 0 0</p><p>　　-0.5

108、617 0 0 0</p><p>　　-0.6002 0 0 0</p><p>　　-0.7119 0 0 0</p><p>　　各節(jié)點(diǎn)電壓的大小為：</p><p>　　u = 1.0600

109、 1.0365 1.0088 1.0073 1.0016</p><p>　　各節(jié)點(diǎn)電壓的修正量：</p><p>　　ut =1.0e-004 *</p><p>　　-0.9362 0 0 0</p><p>　　-0.5617 0 0

110、 0</p><p>　　-0.6002 0 0 0</p><p>　　-0.7119 0 0 0</p><p>　　各節(jié)點(diǎn)電壓的大小為：</p><p>　　u =1.0600 1.0365 1.0087 1.

111、0073 1.0016</p><p>　　各節(jié)點(diǎn)電壓的修正量：</p><p>　　ut =1.0e-004 *</p><p>　　-0.9362 0 0 0</p><p>　　-0.5617 0 0 0</p><

112、p>　　-0.6002 0 0 0</p><p>　　-0.7119 0 0 0</p><p>　　各節(jié)點(diǎn)電壓的大小為：</p><p>　　u =1.0600 1.0365 1.0087 1.0072 1.0016</p>

113、<p>　　各節(jié)點(diǎn)電壓的修正量：</p><p>　　ut =1.0e-004 *</p><p>　　-0.9362 0 0 0</p><p>　　-0.5617 0 0 0</p><p>　　-0.6002 0

114、 0 0</p><p>　　-0.7119 0 0 0</p><p>　　各節(jié)點(diǎn)電壓的大小為：</p><p>　　u =1.0600 1.0365 1.0087 1.0072 1.0015</p><p><b>　　各節(jié)點(diǎn)功

115、率s為：</b></p><p>　　S1 =1.2987 + 0.2445i</p><p>　　各節(jié)點(diǎn)功率sij為：</p><p>　　S = 0 0.8900 + 0.1387i 0.4087 + 0.1058i 0 0 </p>&

116、lt;p>　　-0.8755 - 0.0954i 0 0.2468 + 0.0815i 0.2792 + 0.0806i 0.5488 + 0.1333i</p><p>　　-0.3960 - 0.0677i -0.2430 - 0.0701i 0 0.1891 - 0.0121i 0 &

117、lt;/p><p>　　0 -0.2745 - 0.0665i -0.1887 + 0.0132i 0 0.0633 + 0.0033i</p><p>　　0 -0.5369 - 0.0977i 0 -0.0630 - 0.0023i 0 &

118、lt;/p><p><b>　　k =4</b></p><p>　　各節(jié)點(diǎn)有功率不平衡量為：</p><p>　　dp = 1.0e-003 *</p><p>　　0 0.7371 0.2618 0.4425 0.5138</p><p>　　各節(jié)點(diǎn)有功率不平衡量為：<

119、;/p><p>　　dp =1.0e-003 *</p><p>　　0 0.7371 -0.0482 0.4425 0.5138</p><p>　　各節(jié)點(diǎn)有功率不平衡量為：</p><p>　　dp = 1.0e-003 *</p><p>　　0 0.7371 -0.0482 -0

120、.0639 0.5138</p><p>　　各節(jié)點(diǎn)有功率不平衡量為：</p><p>　　dp =1.0e-003 *</p><p>　　0 0.7371 -0.0482 -0.0639 -0.1022</p><p>　　dtt =1.0e-004 *</p><p>　　0.2903

121、 0.1366 0.1431 0.1567</p><p>　　各節(jié)點(diǎn)電壓的相位角dt為：</p><p>　　dt = 0 -0.0461 -0.0839 -0.0896 -0.1044</p><p>　　dtt =1.0e-004 *</p><p>　　0.2903 0.1366 0.14

122、31 0.1567</p><p>　　各節(jié)點(diǎn)電壓的相位角dt為：</p><p>　　dt = 0 -0.0461 -0.0839 -0.0896 -0.1044</p><p>　　dtt = 1.0e-004 *</p><p>　　0.2903 0.1366 0.1431 0.1567</