1、<p><b>　　摘要：</b></p><p>　　潮流計(jì)算是電力系統(tǒng)最基本最常用的計(jì)算。根據(jù)系統(tǒng)給定的運(yùn)行條件，網(wǎng)絡(luò)接線及元件參數(shù)，通過(guò)潮流計(jì)算可以確定各母線的電壓幅值和相角，各元件流過(guò)的功率，整個(gè)系統(tǒng)的功率損耗。潮流計(jì)算是實(shí)現(xiàn)電力系統(tǒng)安全經(jīng)濟(jì)發(fā)供電的必要手段和重要工作環(huán)節(jié)。因此，潮流計(jì)算在電力系統(tǒng)的規(guī)劃計(jì)算，生產(chǎn)運(yùn)行，調(diào)度管理及科學(xué)計(jì)算中都有著廣泛的應(yīng)用。</p>

2、;<p>　　首先，畫(huà)出系統(tǒng)的等效電路圖，在計(jì)算出各元件參數(shù)的基礎(chǔ)上，應(yīng)用牛頓—拉夫遜Newton-Raphson法以及MATLAB軟件進(jìn)行計(jì)算對(duì)給定系統(tǒng)圖進(jìn)行了四種不同負(fù)荷下的潮流計(jì)算，經(jīng)過(guò)調(diào)節(jié)均得到符合電壓限制及功率限制的潮流分布。</p><p>　　其次，輪流斷開(kāi)環(huán)網(wǎng)的三條支路，在新的系統(tǒng)結(jié)構(gòu)下進(jìn)行次潮流計(jì)算，結(jié)果亦均滿足潮流分布要求。牛頓—拉夫遜Newton-Raphson法具有較好的收斂

3、性，上訴計(jì)算過(guò)程經(jīng)過(guò)四到五次迭代后均能收斂。</p><p>　　電力系統(tǒng)穩(wěn)態(tài)分析包括潮流計(jì)算（或潮流分析)和靜態(tài)安全分析。潮流計(jì)算針對(duì)電力革統(tǒng)各正常運(yùn)行方式，而靜態(tài)安全分析則要研究各種運(yùn)行方式下個(gè)別系統(tǒng)元件退出運(yùn)行后系統(tǒng)的狀況。其目的是校驗(yàn)系統(tǒng)是否能安全運(yùn)行，即是否有過(guò)負(fù)荷的元件或電壓過(guò)低的母線等。原則上講，靜態(tài)安全分析也可U用潮流計(jì)算來(lái)代替。但是一般靜態(tài)安全分析需要校驗(yàn)的狀態(tài)數(shù)非常多，用嚴(yán)格的潮流計(jì)算來(lái)分析這

4、些狀態(tài)往往計(jì)算量過(guò)大，因此不得不尋求一些特殊的算法以滿足要求。</p><p>　　牛頓法是數(shù)學(xué)中解決非線性方程式的典型方法，有較好的收斂性。解決電力系統(tǒng)潮流計(jì)算問(wèn)題是以導(dǎo)納距陣為基礎(chǔ)的，因此，只要在迭代過(guò)程中盡可能保持方程式系數(shù)距陣的稀疏性，就可以大大提高牛頓法潮流程序的放率。自從20 世紀(jì)60 年代中期利用了最佳順序消去法以后，牛頓法在收斂性、內(nèi)存要求、速度方面都超過(guò)了阻抗法，成為直到目前仍在廣泛采用的優(yōu)秀方

5、法。</p><p>　　關(guān)鍵詞：牛頓-拉夫遜法 MATLAB 潮流計(jì)算</p><p><b>　　目 錄</b></p><p>　　課程設(shè)計(jì)相關(guān)要求3</p><p><b>　　課程設(shè)計(jì)題目3</b></p><p><b>　　課程設(shè)計(jì)的內(nèi)

6、容3</b></p><p>　　課程設(shè)計(jì)的要求與數(shù)據(jù)3</p><p>　　課程設(shè)計(jì)應(yīng)完成的工作3</p><p>　　潮流計(jì)算設(shè)計(jì)題目4</p><p><b>　　潮流計(jì)算題目4</b></p><p>　　對(duì)課題的分析及求解思路5</p><p&

7、gt;<b>　　算法分析5</b></p><p><b>　　潮流計(jì)算算法5</b></p><p>　　關(guān)于電力系統(tǒng)潮流計(jì)算手工計(jì)算7</p><p><b>　　節(jié)點(diǎn)導(dǎo)納矩陣7</b></p><p><b>　　簡(jiǎn)化雅可比矩陣8</b>

8、</p><p><b>　　修正、迭代9</b></p><p><b>　　程序設(shè)計(jì)10</b></p><p>　　潮流計(jì)算流程圖10</p><p>　　潮流計(jì)算源程序11</p><p><b>　　運(yùn)行計(jì)算結(jié)果16</b></

9、p><p><b>　　參考文獻(xiàn)17</b></p><p>　　一、課程設(shè)計(jì)相關(guān)要求</p><p><b>　　1、課程設(shè)計(jì)題目：</b></p><p>　　電力系統(tǒng)潮流和穩(wěn)定計(jì)算</p><p><b>　　課程設(shè)計(jì)的內(nèi)容</b></p>

10、;<p>　　掌握潮流計(jì)算的基本原理</p><p>　　根據(jù)電力系統(tǒng)網(wǎng)絡(luò)推導(dǎo)電力網(wǎng)絡(luò)數(shù)學(xué)模型，寫(xiě)出節(jié)點(diǎn)導(dǎo)納矩陣</p><p>　　賦予各節(jié)點(diǎn)電壓變量初值后，求解不平衡量</p><p><b>　　形成雅克比矩陣</b></p><p>　　求解修正量后，重新修改初值，從2開(kāi)始重新循環(huán)計(jì)算</p&

11、gt;<p>　　求解的電壓變量達(dá)到所要求的精度時(shí)，再計(jì)算各支路功率分布、功率損耗和平衡節(jié)點(diǎn)功率</p><p><b>　　上機(jī)編程調(diào)試</b></p><p>　　計(jì)算分析給定系統(tǒng)潮流，并與手工計(jì)算結(jié)果作比較分析</p><p>　　書(shū)寫(xiě)課程設(shè)計(jì)說(shuō)明書(shū)，并打印紙質(zhì)版上交，準(zhǔn)備計(jì)算機(jī)演示答辯</p><p&g

12、t;　　課程設(shè)計(jì)的要求與數(shù)據(jù)</p><p>　　運(yùn)用某種語(yǔ)言，給出程序源代碼和程序說(shuō)明</p><p>　　給定題目的輸入輸出文件</p><p><b>　　給出程序的計(jì)算過(guò)程</b></p><p>　　以某一系統(tǒng)為例，給出程序的手算過(guò)程</p><p>　　課程設(shè)計(jì)應(yīng)完成的工作</p

13、><p>　　按照格式規(guī)范，獨(dú)立完成課程設(shè)計(jì)說(shuō)明書(shū)的撰寫(xiě)</p><p>　　完成潮流計(jì)算程序的編寫(xiě)和調(diào)試</p><p>　　完成潮流計(jì)算的手算過(guò)程</p><p>　　二、潮流計(jì)算設(shè)計(jì)題目</p><p><b>　　1、潮流計(jì)算題目</b></p><p>　　圖2-1

14、電力系統(tǒng)接線圖</p><p>　　2、對(duì)課題的分析及求解思路</p><p>　　此電力系統(tǒng)是一個(gè)5節(jié)點(diǎn)，4支路的電力網(wǎng)絡(luò)。其中包含3個(gè)PQ節(jié)點(diǎn)，一個(gè)PV節(jié)點(diǎn)，和一個(gè)平衡節(jié)點(diǎn)。綜合比較牛頓拉夫遜法（直角坐標(biāo)、極坐標(biāo)）、PQ分解法等多種求解方法的特點(diǎn)，最后確定采用牛頓拉夫遜法（極坐標(biāo)）。因?yàn)榇朔椒ㄋ杞獾姆匠探M最少。</p><p><b>　　三、算法分

15、析</b></p><p><b>　　1、潮流計(jì)算算法</b></p><p>　　本題采用了題目要求的牛頓－拉夫遜潮流計(jì)算的方法。</p><p>　　牛頓-拉夫遜法潮流計(jì)算的公式。把牛頓法用于潮流計(jì)算，采用直角坐標(biāo)形式表示的如式（1-3）所示的形式。其中電壓和支路導(dǎo)納可表示為：</p><p><

16、b>　　（1-2）</b></p><p>　　將上述表示式（1-2）代入（1-1）式的右端，展開(kāi)并分出實(shí)部和虛部，便得：</p><p><b>　?。?-3）</b></p><p>　　按照以上的分類(lèi)，PQ節(jié)點(diǎn)的輸出有功功率和無(wú)功功率是給定的，則第i節(jié)點(diǎn)的給定功率設(shè)為和（稱(chēng)為注入功率）。</p><p

17、>　　假定系統(tǒng)中的第1、2、…、m節(jié)點(diǎn)為PQ節(jié)點(diǎn)，對(duì)其中每一個(gè)節(jié)點(diǎn)的N-R法表達(dá)式</p><p>　　F(x)=0[如、、]形式有些下列方程：</p><p><b>　?。?-4）</b></p><p>　　=（1、2、…、m）</p><p>　　PV節(jié)點(diǎn)的有功功率和節(jié)點(diǎn)電壓幅值是給定的。假定系統(tǒng)中的第m+

18、1、m+2、…、n-1節(jié)點(diǎn)為PV節(jié)點(diǎn)，則對(duì)其中每一PV節(jié)點(diǎn)可以列寫(xiě)方程：</p><p><b>　?。?-5）</b></p><p>　　=（m+1、m+2、…、n-1）</p><p>　　(6)形成雅可比矩陣。N-R法的思想是；本例；對(duì)F(x)求偏導(dǎo)的式（1-6）、式（1-7），即式（1-4）、式（1-5）中的、、是多維變量的函數(shù)，對(duì)多

19、維變量求偏導(dǎo)（、、、、、、、…），并以矩陣的形式表達(dá)稱(chēng)為雅可比矩陣。</p><p>　　當(dāng)j=i時(shí)，對(duì)角元素為</p><p><b>　?。?-6）</b></p><p>　　當(dāng)時(shí),矩陣非對(duì)角元素為：</p><p><b>　　（1-7）</b></p><p>　　

20、由上式不難看出，雅可比矩陣有以下特點(diǎn)。</p><p>　　雅可比矩陣中的諸元素都是節(jié)點(diǎn)電壓的函數(shù)，因此在迭代過(guò)程中，它們將隨著節(jié)點(diǎn)電壓的變化而不斷的變化。</p><p>　　雅可比矩陣具有結(jié)構(gòu)對(duì)稱(chēng)性，數(shù)據(jù)不對(duì)稱(chēng)。如非對(duì)角，，。</p><p>　　由式（1-7）可以看出，當(dāng)導(dǎo)納矩陣中非對(duì)角元素為零時(shí)，。雅可比矩陣中相應(yīng)的元素也為零，即矩陣是非常稀疏的。因此，修正

21、方程的求解同樣可以應(yīng)用稀疏矩陣的求解技巧。正是由于這一點(diǎn)才使N-R法獲得廣泛的應(yīng)用。</p><p>　　2、關(guān)于電力系統(tǒng)潮流計(jì)算手工計(jì)算</p><p><b>　　節(jié)點(diǎn)導(dǎo)納矩陣</b></p><p><b>　　求得節(jié)點(diǎn)導(dǎo)納矩陣Y</b></p><p><b>　　=</b&

22、gt;</p><p>　　各節(jié)點(diǎn)的導(dǎo)納值如下：</p><p>　　; ；</p><p>　??； ;</p><p>　　; ;</p><p>　?。?

23、 ;</p><p>　　； ;</p><p>　　; ;</p><p>　　; ;</p><p>　　; ;</p><

24、p>　　; ;</p><p>　　; ;</p><p><b>　　;</b></p><p><b>　　;</b></p><p><b>　　;</b&

25、gt;</p><p><b>　　;</b></p><p><b>　　.</b></p><p><b>　　簡(jiǎn)化雅可比矩陣</b></p><p>　　形成有功迭代和無(wú)功迭代的簡(jiǎn)化雅可比矩陣B/和B//</p><p><b>　　B/

26、= </b></p><p><b>　　B//=</b></p><p>　　將B/ 和B//進(jìn)行三角分解：</p><p><b>　　修正、迭代</b></p><p>　　給定PQ節(jié)點(diǎn)初值和各節(jié)點(diǎn)電壓相角初值</p><p>　　V1=1.05∠0。

27、 ，V2(0)=V3(0)=1.0，V4=1.1</p><p>　　δ2(0)=δ3(0)=0， δ4(0)=0</p><p>　　1 作第一次有功迭代，按公式計(jì)算節(jié)點(diǎn)有功功率不平衡量</p><p>　　△P2(0)=-0.55-（-0.024037）=-0.525963</p><p>　　△P3(0)=-0.30-（-0.02269

28、5）=-0.277305</p><p>　　△P4(0)=0.500000</p><p>　　△P1(0)/V1(0)=0.454545 △P2(0)/ V2(0)=-0.525963</p><p>　　△P3(0)/V3(0)=-0.277309</p><p>　　2 做第一次無(wú)功迭代，按公式計(jì)算無(wú)功功率不平衡量，計(jì)算時(shí)電壓

29、相角最新的 修正值。</p><p>　　△Q2(0)=-0.13-(-0.001550)=-0.039594</p><p>　　△Q3(0)=-0.18-(-0.14406)=-0.039588</p><p>　　△Q2(0)/ V2(0)=-0.131553</p><p>　　△Q3(0)/V3(0)=-0.039588&

30、lt;/p><p>　　解修正方程式，可得各節(jié)點(diǎn)電壓幅值的修正量為</p><p>　　△V3(0)）=-0.014855</p><p><b>　　于是有：</b></p><p>　　V2(1) = V2(0)+△V2(1)=0.964776</p><p>　　V3(1) = V3(0)+△V

31、3(1)=0.985145</p><p>　　到這里為止，第一輪有功迭代和無(wú)功迭代便做完了。</p><p>　　3 按公式計(jì)算平衡節(jié)點(diǎn)功率，得：</p><p>　　P1+jQ1=0.367885+j0.264696</p><p>　　經(jīng)過(guò)四輪迭代，節(jié)點(diǎn)不平衡功率也下降到10-5以下，迭代到此結(jié)束。</p><p&g

32、t;<b>　　五、 程序設(shè)計(jì)</b></p><p>　　1、 潮流計(jì)算流程圖</p><p>　　圖5-1 潮流計(jì)算流程圖</p><p>　　2、 潮流計(jì)算源程序</p><p>　　據(jù)課題題目，本程序把節(jié)點(diǎn)1設(shè)為平衡節(jié)點(diǎn)，節(jié)點(diǎn)2、3、4為PQ節(jié)點(diǎn)，節(jié)點(diǎn)5為PV節(jié)點(diǎn)。</p><p>　　G

33、(1,1)=10.834;</p><p>　　B(1,1)=-32.500;</p><p>　　G(1,2)=-1.667;</p><p>　　B(1,2)=5.000;</p><p>　　G(1,3)=-1.667;</p><p>　　B(1,3)=5.000;</p><p>　　

34、G(1,4)=-2.500;</p><p>　　B(1,4)=7.500;</p><p>　　G(1,5)=-5.000;</p><p>　　B(1,5)=15.000;</p><p>　　G(2,1)=-1.667;</p><p>　　B(2,1)=5.000;</p><p>　　

35、G(2,2)=12.917;</p><p>　　B(2,2)=-38.750;</p><p>　　G(2,3)=-10.000;</p><p>　　B(2,3)=30.000;</p><p><b>　　G(2,4)=0;</b></p><p><b>　　B(2,4)=0;&

36、lt;/b></p><p>　　G(2,5)=-1.250;</p><p>　　B(2,5)=3.750;</p><p>　　G(3,1)=-1.667;</p><p>　　B(3,1)=5.000;</p><p>　　G(3,2)=-10.000;</p><p>　　B(3,

37、2)=30.000;</p><p>　　G(3,3)=12.917;</p><p>　　B(3,3)=-38.750;</p><p>　　G(3,4)=-1.250;</p><p>　　B(3,4)=3.750;</p><p><b>　　G(3,5)=0;</b></p>

38、<p><b>　　B(3,5)=0;</b></p><p>　　G(4,1)=-2.500;</p><p>　　B(4,1)=7.500;</p><p><b>　　G(4,2)=0;</b></p><p><b>　　B(4,2)=0;</b></

39、p><p>　　G(4,3)=-1.250;</p><p>　　B(4,3)=3.750;</p><p>　　G(4,4)=3.750;</p><p>　　B(4,4)=-11.250;</p><p><b>　　G(4,5)=0;</b></p><p><b&

40、gt;　　B(4,5)=0;</b></p><p>　　G(5,1)=-5.000;</p><p>　　B(5,1)=15.000;</p><p>　　G(5,2)=-1.250;</p><p>　　B(5,2)=3.750;</p><p><b>　　G(5,3)=0;</b&g

41、t;</p><p><b>　　B(5,3)=0;</b></p><p><b>　　G(5,4)=0;</b></p><p><b>　　B(5,4)=0;</b></p><p>　　G(5,5)=6.250;</p><p>　　B(5,5)

42、=-18.750;</p><p><b>　　Y=G+j*B;</b></p><p>　　delt(1)=0;</p><p>　　delt(2)=0;</p><p>　　delt(3)=0;</p><p>　　delt(4)=0;</p><p><b>

43、;　　u(1)=1.0;</b></p><p><b>　　u(2)=1.0;</b></p><p><b>　　u(3)=1.0;</b></p><p><b>　　u(4)=1.0;</b></p><p>　　p(1)=0.20;</p>

44、<p>　　q(1)=0.20;</p><p>　　p(2)=-0.45;</p><p>　　q(2)=-0.15;</p><p>　　p(3)=-0.40;</p><p>　　q(3)=-0.05;</p><p>　　p(4)=-0.60;</p><p>　　q(4)=-

45、0.10;</p><p><b>　　k=0;</b></p><p>　　precision=1;</p><p>　　N1=4; %the N1 is the amount of the PQ bus</p><p>　　while precision>0.00001</p><p&g

46、t;　　delt(5)=0;</p><p>　　u(5)=1.06;</p><p>　　for m=1:N1</p><p>　　for n=1:N1+1</p><p>　　pt(n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));</p&g

47、t;<p>　　qt(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));</p><p><b>　　end</b></p><p>　　pp(m)=p(m)-sum(pt);</p><p>　　qq(m)=q(m)-sum(qt);&l

48、t;/p><p><b>　　end</b></p><p>　　for m=1:N1</p><p>　　for n=1:N1+1</p><p>　　h0(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));</p>

49、<p>　　n0(n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));</p><p>　　j0(n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));</p><p>　　L0(n)=-u(m)*u

50、(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));</p><p><b>　　end</b></p><p>　　H(m,m)=sum(h0)-u(m)^2*(G(m,m)*sin(delt(m)-delt(m))-B(m,m)*cos(delt(m)-delt(m</p>&l

51、t;p><b>　　)));</b></p><p>　　N(m,m)=sum(n0)-2*u(m)^2*G(m,m)+u(m)^2*(G(m,m)*cos(delt(m)-delt(m))+B(m,m)*sin(delt(m)-delt(m)));</p><p>　　J(m,m)=sum(j0)+u(m)^2*(G(m,m)*cos(delt(m)-delt

52、(m))+B(m,m)*sin(delt(m)-delt(m)));</p><p>　　L(m,m)=sum(L0)+2*u(m)^2*B(m,m)+u(m)^2*(G(m,m)*sin(delt(m)-delt(m))-B(m,m)*cos(delt(m)-delt(m)));</p><p><b>　　end</b></p><p>　

53、　for m=1:N1</p><p>　　JJ(2*m-1,2*m-1)=H(m,m);</p><p>　　JJ(2*m-1,2*m)=N(m,m);</p><p>　　JJ(2*m,2*m-1)=J(m,m);</p><p>　　JJ(2*m,2*m)=L(m,m);</p><p><b>　　e

54、nd</b></p><p>　　for m=1:N1</p><p>　　for n=1:N1</p><p><b>　　if m==n</b></p><p><b>　　else</b></p><p>　　H(m,n)=-u(m)*u(n)*(G(m,n

55、)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));</p><p>　　J(m,n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));</p><p>　　N(m,n)=-J(m,n);</p><p>　　L(m,n)=H(

56、m,n);</p><p>　　JJ(2*m-1,2*n-1)=H(m,n);</p><p>　　JJ(2*m-1,2*n)=N(m,n);</p><p>　　JJ(2*m,2*n-1)=J(m,n);</p><p>　　JJ(2*m,2*n)=L(m,n);</p><p><b>　　end<

57、/b></p><p><b>　　end</b></p><p><b>　　end</b></p><p>　　for m=1:N1</p><p>　　PP(2*m-1)=pp(m);</p><p>　　PP(2*m)=qq(m);</p><

58、;p><b>　　end</b></p><p>　　uu=-inv(JJ)*PP';</p><p>　　precision=max(abs(uu));</p><p>　　for n=1:N1</p><p>　　delt(n)=delt(n)+uu(2*n-1);</p><p&g

59、t;　　u(n)=u(n)+uu(2*n);</p><p><b>　　end</b></p><p><b>　　k=k+1;</b></p><p><b>　　end</b></p><p>　　K=k-1,delt,u'</p><p>

60、;　　%the following program is used to calculate the S5 and S(m,n)</p><p>　　for n=1:N1+1</p><p>　　U(n)=u(n)*(cos(delt(n))+j*sin(delt(n)));</p><p><b>　　end</b></p>&

61、lt;p>　　for m=1:N1+1</p><p>　　I(m)=Y(5,m)*U(m);</p><p><b>　　end</b></p><p>　　S5=U(5)*sum(conj(I))</p><p>　　for m=1:N1+1</p><p>　　for n=1:N1+1

62、</p><p>　　S(m,n)=U(m)*(conj(U(m))-conj(U(n)))*conj(-Y(m,n));</p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　Y</b></p>&l

63、t;p><b>　　S</b></p><p><b>　　3、運(yùn)行計(jì)算結(jié)果</b></p><p><b>　　K =4；</b></p><p>　　delt =[ -0.0461 -0.0839 -0.0896 -0.1044 0]；</p><p&g

64、t;　　U =[ 1.0365 1.0087 1.0073 1.0016 1.0600]’</p><p>　　S5 =1.2982 + 0.2445i</p><p><b>　　Y =</b></p><p>　　10.8340 -32.5000i -1.6670 + 5.0000i -1.6670 + 5.0000i -2.

65、5000 + 7.5000i -5.0000 +15.0000i</p><p>　　-1.6670 + 5.0000i 12.9170 -38.7500i -10.0000 +30.0000i 0 -1.2500 + 3.7500i</p><p>　　-1.6670 + 5.0000i -10.0000 +30.0000i 12.9170 -3

66、8.7500i -1.2500 + 3.7500i 0 </p><p>　　-2.5000 + 7.5000i 0 -1.2500 + 3.7500i 3.7500 -11.2500i 0 </p><p>　　-5.0000 +15.0000i -1.2500 + 3.7500i

67、 0 0 6.2500 -18.7500i</p><p><b>　　S =</b></p><p>　　0 0.2469 + 0.0815i 0.2793 + 0.0806i 0.5489 + 0.1333i -0.8751 - 0.0954i</p&g

68、t;<p>　　-0.2431 - 0.0701i 0 0.1891 - 0.0121i 0 -0.3960 - 0.0677i</p><p>　　-0.2746 - 0.0664i -0.1887 + 0.0132i 0 0.0633 + 0.0033i 0

69、 </p><p>　　-0.5370 - 0.0977i 0 -0.0630 - 0.0023i 0 0 </p><p>　　0.8895 + 0.1387i 0.4087 + 0.1058i 0 0

70、 0 </p><p><b>　　參考文獻(xiàn)</b></p><p>　　[1]于永源，楊綺雯. 電力系統(tǒng)分析[M]. 北京：中國(guó)電力出版社，2007</p><p>　　[2]孟祥萍.電力系統(tǒng)分析[M].北京：高等教育出版社，2005</p><p>　　[3]何仰贊等.電力系統(tǒng)分析[M]. 武