復雜網(wǎng)絡N-R法潮流分析計算

1、編號 課 程 設 計（ 2012級本科）題 目： 復雜網(wǎng)絡N-R法潮流分析與計算的設計 系（部）院： 物理與機電工程學院 專 業(yè)： 電氣工程及其自動化 作者姓名： 指導教師： 職稱： 副教授 完成日期： 2015 年 6 月 30 日 二一五 年 六 月河西學院本科生課程設計任務書設 計 題 目復雜網(wǎng)絡N-R法潮流分析與計算的設計作 者 姓 名學院、專業(yè)、年級物電學院電氣工程專業(yè) 12 級指導教師姓名、職稱 副教授任務下達日期2015年5月20日一、設計資料1.系統(tǒng)圖的確定選擇六節(jié)點、環(huán)網(wǎng)、兩電源和多引出的電力系統(tǒng)，簡化電力系統(tǒng)圖如圖1所示，等值阻抗圖如圖2所示。運用以直角坐標表示的牛頓-拉夫

2、遜計算如圖1所示系統(tǒng)中的潮流分布。計算精度要求各節(jié)點電壓的誤差或修正量不大于。 圖1 電力系統(tǒng)圖圖2 電力系統(tǒng)等值阻抗圖2.各節(jié)點的初值及阻抗參數(shù)該系統(tǒng)中，節(jié)點為平衡節(jié)點，保持=1.05+j0為定值，節(jié)點為PV節(jié)點，其他四個節(jié)點都是PQ節(jié)點。給定的注入電壓標幺值、線路阻抗標幺值、線路阻抗標幺值、輸出功率標幺值和變壓器變比標幺值如圖2所示的注釋。表1 各節(jié)點電壓標幺值參數(shù)UUUUUU1.051.001.001.001.001.05表2 線路、變壓器阻抗標幺值線路TL2L3L4L5L6阻抗j0.030.06+j0.0250.04+j0.250.08+j0.300.1+j0.35j0.015表3 節(jié)

3、點輸出功率節(jié)點功率2+j11.8+j0.401.6+j0.83.7+j1.35注：各PQ節(jié)點的電壓取1是為了方便計算和最后驗證程序的正確性。二、設計的基本要求3.1設計及計算說明書（1）說明書要求書寫整齊，條理分明，表達正確、語言正確。（2）計算書內(nèi)容：為各設計內(nèi)容最終成果、確定提供依據(jù)進行的技術分析、論證和定量計算，如。（3）計算書要求：計算無誤，分析論證過程簡單明了，各設計內(nèi)容列表匯總。3.2圖紙（1）繪制分析所需的必要圖紙 （2）圖紙要求：用標準符號繪制，布置均勻，設備符號大小合適，清晰美觀。三、論文(設計)進度安排階段論文（設計）各階段名稱起止日期1熟悉設計任務書、設計題目及設計背景資

4、料5.205.252查閱有關資料5.265.273閱讀設計要求必讀的參考資料5.285.294書寫設計說明書5.306.155小組答辯質(zhì)疑6.216.226上交設計成果6.30四、需收集和閱讀的資料及參考文獻（指導教師指定）1: 陳珩.電力系統(tǒng)穩(wěn)態(tài)分析（第三版）M，北京，中國電力出版社，20072：何仰贊.溫增銀.電力系統(tǒng)分析第三版M，武漢,華中科技大學出版社，20023：陳悅.電氣工程畢業(yè)設計指南電力系統(tǒng)分冊M，北京，中國水利水電出版社，20084：5：教 研 室 意 見 負責人簽名： 年 月 日 目錄1 牛頓拉夫遜法概述11.1牛頓-拉夫遜法基本原理11.2牛頓-拉夫遜法潮流求解過程22手

5、算潮流計算62.1用圖中數(shù)據(jù)和等值網(wǎng)絡形成節(jié)點導納矩陣62.2設各節(jié)點電壓初始值為：72.3用公式72.4求取雅可比矩陣72.5求修正量矩陣82.6計算修正各節(jié)點電壓83計算機算法潮流計算93.1牛頓拉夫遜法的程序框圖93.2結果顯示9總結19附件20參考文獻251 牛頓拉夫遜法概述電力系統(tǒng)潮流計算是電力系統(tǒng)分析中的一種最基本的計算，是對復雜電力系統(tǒng)正常和故障條件下穩(wěn)態(tài)運行狀態(tài)的計算。潮流計算的目標是求取電力系統(tǒng)在給定運行狀態(tài)的計算。即節(jié)點電壓和功率分布，用以檢查系統(tǒng)各元件是否過負荷。各點電壓是否滿足要求，功率的分布和分配是否合理以及功率損耗等。對現(xiàn)有電力系統(tǒng)的運行和擴建，對新的電力系統(tǒng)進行規(guī)

6、劃設計以及對電力系統(tǒng)進行靜態(tài)和暫態(tài)穩(wěn)定分析都是以潮流計算為基礎。潮流計算結果可用如電力系統(tǒng)穩(wěn)態(tài)研究，安全估計或最優(yōu)潮流等對潮流計算的模型和方法有直接影響。實際電力系統(tǒng)的潮流技術那主要采用牛頓-拉夫遜法。1.1牛頓-拉夫遜法基本原理牛頓-拉夫遜法(簡稱NR法)在數(shù)學上是求解非線性代數(shù)方程式的有效方法。其要點是把非線性方程式的求解過程變成反復地對相應的線性方程式進行求解的過程。即通常所稱的逐次線性化過程。對于非線性代數(shù)方程組：fx=0即 fix1,x2,xn=0 i=1,2,n (1-1-1)在待求量x的某一個初始估計值附近，將上式展開成泰勒級數(shù)并略去二階及以上的高階項，得到如下的經(jīng)線性化的方程組

7、：fx0+f'x0x0=0 (1-1-2)上式稱之為牛頓法的修正方程式。由此可以求得第一次迭代的修正量 x0=-f'x0-1fx0 (1-1-3)將和相加，得到變量的第一次改進值。接著就從出發(fā)，重復上述計算過程。因此從一定的初值出發(fā)，應用牛頓法求解的迭代格式為： f'xkxk=-fxk (1-1-4) xk+1=xk+xk (1-1-5)上兩式中：是函數(shù)對于變量x的一階偏導數(shù)矩陣，即雅可比矩陣Jk為迭代次數(shù)。有上式可見，牛頓法的核心便是反復形式并求解修正方程式。牛頓法當初始估計值和方程的精確解足夠接近時，收斂速度非?？欤哂衅椒绞諗刻匦?。牛頓潮流算法突出的優(yōu)點是收斂速度

8、快，若選擇到一個較好的初值，算法將具有平方收斂特性，一般迭代45次便可以收斂到一個非常精確的解。而且其迭代次數(shù)與所計算網(wǎng)絡的規(guī)?；緹o關。牛頓法也具有良好的收斂可靠性，對于對以節(jié)點導納矩陣為基礎的高斯法呈病態(tài)的系統(tǒng)，牛頓法也能可靠收斂。牛頓法所需的內(nèi)存量及每次迭代所需時間均較高斯法多。牛頓法的可靠收斂取決于有一個良好的啟動初值。如果初值選擇不當，算法有可能根本不收斂或收斂到一個無法運行的節(jié)點上。對于正常運行的系統(tǒng)，各節(jié)點電壓一般均在額定值附近，偏移不會太大，并且各節(jié)點間的相位角差也不大，所以對各節(jié)點可以采用統(tǒng)一的電壓初值(也稱為平直電壓)，如假定：Ui0=1 i0=0 或 ei0=1 fi0=

9、0 i=q1,2,n,is (1-1-6)這樣一般能得到滿意的結果。但若系統(tǒng)因無功緊張或其它原因?qū)е码妷嘿|(zhì)量很差或有重載線路而節(jié)點間角差很大時，仍用上述初始電壓就有可能出現(xiàn)問題。解決這個問題的辦法可以用高斯法迭代12次，以此迭代結果作為牛頓法的初值。也可以先用直流法潮流求解一次以求得一個較好的角度初值，然后轉(zhuǎn)入牛頓法迭代。1.2牛頓-拉夫遜法潮流求解過程以下討論的是用直角坐標形式的牛頓拉夫遜法潮流的求解過程。當采用直角坐標時，潮流問題的待求量為各節(jié)點電壓的實部和虛部兩個分量e1,f1,e2,f2,en,fn由于平衡節(jié)點的電壓向量是給定的，因此待求兩共需要2(n-1)個方程式。事實上，除了平衡節(jié)

10、點的功率方程式在迭代過程中沒有約束作用以外，其余每個節(jié)點都可以列出兩個方程式。對PQ節(jié)點來說Pis和Qis是給定的，因而可以寫出Pi=Pis-eijiGijej-Bijfj-fjjiGijfj+Bijej=0Qi=Qis-fijiGijej-Bijfj+ejjiGijfj+Bijej=0 （1-2-1）對PV節(jié)點來說，給定量是Pis和Qis，因此可以列出 Pi=Pis-eijiGijej-Bijfj-fijiGijfj+Bijej=0Vi2=Vis2-ei2+fi2=0 （1-2-2）求解過程大致可以分為以下步驟：（1）形成節(jié)點導納矩陣（2）將各節(jié)點電壓設初值U，（3）將節(jié)點初值代入相關求式，

11、求出修正方程式的常數(shù)項向量（4）將節(jié)點電壓初值代入求式，求出雅可比矩陣元素（5）求解修正方程，求修正向量（6）求取節(jié)點電壓的新值（7）檢查是否收斂，如不收斂，則以各節(jié)點電壓的新值作為初值自第3步重新開始進行狹義次迭代，否則轉(zhuǎn)入下一步（8）計算支路功率分布，PV節(jié)點無功功率和平衡節(jié)點注入功率。以直角坐標系形式表示迭代推算式 采用直角坐標時,節(jié)點電壓相量及復數(shù)導納可表示為:Vi=ei+jfi Yij=Gij+jBij (1-2-3) 將以上二關系式代入上式中,展開并分開實部和虛部;假定系統(tǒng)中的第1,2,m號為PQ節(jié)點, m+1,m+2,n-1為PV節(jié)點,根據(jù)節(jié)點性質(zhì)的不同,得到如下迭代推算式: 對

12、于PQ節(jié)點Pi=Pi-eij=1nGijej-Bijfj-fjj=1nGijfj+BijejQi=Qi-fij=1nGijej-Bijfj+eij=1nGijfj+Bijej (1-2-4)i=1,2,m對于PV節(jié)點 Pi=Pi-eij=1nGijej-Bijfj-fij=1nGijfj+BijejVi2=Vi2-ei2+fi2 (1-2-5)i=m+1,m+2,n-1對于平衡節(jié)點 平衡節(jié)點只設一個,電壓為已知,不參見迭代,其電壓為: Vn=en+jfn (1-2-6)修正方程式(2-3-5)和(2-3-6)兩組迭代式共包括2(n-1)個方程.選定電壓初值及變量修正量符號之后代入式(2-3-5

13、)和(2-3-6),并將其按泰勒級數(shù)展開,略去ei,fi二次方程及以后各項,得到修正方程如下:W=-JU (1-2-7)W= U= (1-2-8) 雅可比矩陣各元素的算式式(1-2-8)中, 雅可比矩陣中的各元素可通過對式(1-2-4)和(1-2-5)進行偏導而求得.當時, 雅比矩陣中非對角元素為 Piej=-Qifj=-Gijei+BijfiPifj=Qiej=Bijei-GijfiU2ej=U2fj=0 (1-2-9)當時,雅可比矩陣中對角元素為: (1-2-10)由式(3-2-9和(3-2-10)看出,雅可比矩陣的特點:陣中各元素是節(jié)點電壓的函數(shù),在迭代過程中,這些元素隨著節(jié)點電壓的變化

14、而變化;導納矩陣中的某些非對角元素為零時,雅可比矩陣中對應的元素也是為零.若,則 必有Jij=0;雅可比矩陣不是對稱矩陣;(i=q1,2,n,is)雅可比矩陣各元素的表示如下:Hij=Piej=-Gijei+Bijfi (ji)-jiGijej-Bijfj-Giiei-Biifi(j=i)Nij=Pifj=Bijei-Gijfi (ji)-jiGijfj+Bijej+Biiei-Giifi(j=i)Mij=Qiej=Bijei-Gijfi (ji)jiGijfi+Bijej+Biiei-Giifi(j=i)Lij=Qifj=Gijei+Bijfi (ji)-jiGijej-Bijfj+Giie

15、i+Biifi(j=i)Rij=Ui2ej=0 (ji)-2ei(j=i) Sij=Ui2fj=0 (ji)-2fi(j=i)2手算潮流計算2.1用圖中數(shù)據(jù)和等值網(wǎng)絡形成節(jié)點導納矩陣節(jié)點導納矩陣YB對角線上的元素為：Y11=1KYT1+K-1KYT1=-j33.3333Y22=1KYT1+1-KK2YT1+Y0+y23+y25=1.53174-j37.4166Y33=Y0+y23+y34=1.7376-j6.3942Y44=Y0+y34+y45=1.5846-j5.2535Y55=1KYT2+1-KK2YT2+Y0+y45+y25=1.3787-j66.5103Y66=1KYT2+K-1KYT

16、2=-j66.6667非對角線上的元素為：Y12=Y21=1KYT1=j31.7460Y13=Y31=0;Y14=Y41=0;Y15=Y51=0;Y16=Y61=0Y23=Y32=-10.06+j0.25=-0.9077+j3.7822Y25=Y52=-10.04+j0.25=-0.6240+j3.9002Y24=Y42=0;Y26=Y62=0Y34=Y43=-10.08+j0.3=-0.8299+j3.1120Y46=Y64=0Y56=Y65=1KYT2=j63.4921所以導納矩陣為Z= 2.2設各節(jié)點電壓初始值為：U=e+fe1=1.05 f1=0e2=1f2=0e3=1f3=0e4=1

17、f4=0e5=1f5=0e6=1.05 f6=02.3用公式Pi=Pi-j=1j=neiGijej-Bijfj+fiGijfj+BijejQi=Qi-j=1j=nfiGijfj-Bijfj-eiGijfj+BijejUi2=Ui2-ei2+fi2對PQ和PV節(jié)點求取Pi,Qi,Ui2得2.4求取雅可比矩陣 2.5求修正量矩陣J-1=2.6計算修正各節(jié)點電壓e11=1.05000 e21=1.09470 e31= 0.93882 e41= 0.86841 e51= 1.09625 e61=1.05000 f11=0 f21=-0.12300 f31=-0.60441 f41=-0.63777 f

18、51=-0.15857 f61=-0.087753計算機算法潮流計算3.1牛頓拉夫遜法的程序框圖3.2結果顯示：Y = 0 -33.3333i 0 +31.7460i 0 0 0 0 0 +31.7460i 1.5317 -37.4166i -0.9077 + 3.7822i 0 -0.6240 + 3.9002i 0 0 -0.9077 + 3.7822i 1.7376 - 6.3942i -0.8299 + 3.1120i 0 0 0 0 -0.8299 + 3.1120i 1.5846 - 5.2535i -0.7547 + 2.6415i 0 0 -0.6240 + 3.9002i 0

19、 -0.7547 + 2.6415i 1.3787 -66.5103i 0 +63.4921i 0 0 0 0 0 +63.4921i 0 -66.6667i次數(shù)time= 1雅可比矩陣JJ= 0.0312 29.5393 -1.5317 -37.4166 0.9077 3.7822 0 0 0.6240 3.9002 33.9528 0.0312 -37.4166 1.5317 3.7822 -0.9077 0 0 3.9002 -0.6240 0 0 0.9077 3.2822 -1.7376 -6.3942 0.8299 3.1120 0 0 0 0 4.2821 -0.9077 -6.

20、3942 1.7376 3.1120 -0.8299 0 0 0 0 0 0 0.8676 2.4800 -1.5846 -5.2535 0.7547 2.6415 0 0 0 0 3.7441 -0.7921 -5.2535 1.5846 2.6415 -0.7547 0 0 0.6240 3.9002 0 0 0.6858 -0.7310 -1.3787 -66.5103 0 0 3.9002 -0.6240 0 0 6.0140 -0.8237 -66.5103 1.3787 0 0 0 0 0 0 0 0 0 70.0000 0 0 0 0 0 0 0 0 -2.1000 0E =0.

21、8137 0.8366 0.6713 0.6636 1.0500 1.0500F = -0.1896 -0.2825 -0.7053 -0.6191 -0.0714 0U =0.8137 - 0.1896i 0.8366 - 0.2825i 0.6713 - 0.7053i 0.6636 - 0.6191i 1.0500 - 0.0714i 1.0500 dU = 0.1863 0.1896 0.1634 0.2825 0.3287 0.7053 0.3364 0.6191 00.0714PQ = -1.9688 1.2067 -1.8000 0.1000 -1.5623 -0.1679 -3

22、.7689 2.0725 5.0000 0precision = 5次數(shù)time= 2雅可比矩陣JJ= 7.6062 24.4159 -8.3403 -30.1569 1.4557 2.9056 0 0 1.2472 3.0554 27.2498 -4.4313 -30.1569 8.3403 2.9056 -1.4557 0 0 3.0554 -1.2472 0 0 3.4865 2.4260 -3.2599 -4.8584 1.5734 2.3690 0 0 0 0 3.3894 -0.1690 -4.8584 3.2599 2.3690 -1.5734 0 0 0 0 0 0 3.918

23、7 0.4693 -4.7692 -2.4089 2.3697 1.2409 0 0 0 0 2.5381 -1.5855 -2.4089 4.7692 1.2409 -2.3697 0 0 2.8286 2.2020 0 0 3.7251 -1.1059 -42.0899 -43.2856 0 0 2.2020 -2.8286 0 0 3.6775 -0.5472 -43.2856 42.0899 0 0 0 0 0 0 0 0 0 70.0000 0 0 0 0 0 0 0 0 -2.1000 -0.1429E = 1.1261 1.1333 0.8698 0.7324 1.0476 1.

24、0500F = -0.1055 -0.2062 -0.5336 -0.4376 -0.0714 0U = 1.1261 - 0.1055i 1.1333 - 0.2062i 0.8698 - 0.5336i 0.7324 - 0.4376i 1.0476 - 0.0714i 1.0500 dU = -0.3124 -0.0841 -0.2967 -0.0762 -0.1985 -0.1717 -0.0688 -0.1815 0.0024 0PQ =-0.4396 -0.1479 -0.2762 -0.4656 -0.0873 -0.9285 -1.1649 -0.6965 0 -0.0051p

25、recision = 1.1649次數(shù)time= 3雅可比矩陣JJ= 5.1274 34.8149 -5.6738 -41.9740 1.4214 4.1634 0 0 1.1143 4.3262 36.6848 -1.5735 -41.9740 5.6738 4.1634 -1.4214 0 0 4.3262 -1.1143 0 0 3.2041 3.4717 -3.2879 -6.8880 1.5823 3.3556 0 0 0 0 4.7263 -0.4133 -6.8880 3.2879 3.3556 -1.5823 0 0 0 0 0 0 3.1837 0.8345 -4.1815

26、-3.7238 2.0659 1.8948 0 0 0 0 3.6933 -1.5811 -3.7238 4.1815 1.8948 -2.0659 0 0 2.1636 2.5835 0 0 4.3148 -2.1031 -30.1120 -48.1109 0 0 2.5835 -2.1636 0 0 5.3120 0.8976 -48.1109 30.1120 0 0 0 0 0 0 0 0 0 70.0000 0 0 0 0 0 0 0 0 -2.0951 -0.1429E = 1.1221 1.1336 0.8716 0.7159 1.0476 1.0500F =-0.0117 -0.

27、1248 -0.5167 -0.4513 -0.0714 0U =1.1221 - 0.0117i 1.1336 - 0.1248i 0.8716 - 0.5167i 0.7159 - 0.4513i 1.0476 - 0.0714i 1.0500 dU = 0.0040 -0.0939 -0.0004 -0.0814 -0.0019 -0.0169 0.0165 0.0137 0.0000 -0.0000PQ = 0.0997 -0.1346 -0.0892 0.0231 -0.1403 0.0157 -0.1689 0.2752 -0.0000 -0.0000precision = 0.2

28、752次數(shù)time= 4雅可比矩陣JJ= 2.1458 34.7123 -2.1557 -41.9669 1.0627 4.2333 0 0 0.7458 4.3690 36.5315 1.4043 -41.9669 2.1557 4.2333 -1.0627 0 0 4.3690 -0.7458 0 0 3.0372 3.6423 -2.7679 -7.0318 1.3292 3.4243 0 0 0 0 4.7064 0.0350 -7.0318 2.7679 3.4243 -1.3292 0 0 0 0 0 0 3.2852 0.7975 -4.0959 -3.7604 2.0228 1

29、.9125 0 0 0 0 3.7700 -1.3777 -3.7604 4.0959 1.9125 -2.0228 0 0 2.2069 2.5105 0 0 4.6076 -2.0746 -31.0036 -46.9932 0 0 2.5105 -2.2069 0 0 5.1756 1.1428 -46.9932 31.0036 0 0 0 0 0 0 0 0 0 70.0000 0 0 0 0 0 0 0 0 -2.0951 -0.1429E =1.1201 1.1321 0.8731 0.7174 1.0476 1.0500F = -0.0113 -0.1244 -0.5167 -0.

30、4519 -0.0714 0U =1.1201 - 0.0113i 1.1321 - 0.1244i 0.8731 - 0.5167i 0.7174 - 0.4519i 1.0476 - 0.0714i 1.0500 dU = 0.0020 -0.0004 0.0016 -0.0004 -0.0014 -0.0001 -0.0015 0.0006 0.0000 0.0000PQ =0.0024 -0.0001 0.0078 0.0114 -0.0007 0.0027 -0.0056 -0.0023 0.0000 -0.0000precision = 0.0114次數(shù)time= 5雅可比矩陣JJ

31、= 2.1354 34.6468 -2.1387 -41.8916 1.0595 4.2260 0 0 0.7430 4.3614 36.4683 1.4174 -41.8916 2.1387 4.2260 -1.0595 0 0 4.3614 -0.7430 0 0 3.0309 3.6468 -2.7627 -7.0223 1.3267 3.4197 0 0 0 0 4.6905 0.0345 -7.0223 2.7627 3.4197 -1.3267 0 0 0 0 0 0 3.2881 0.8064 -4.0978 -3.7680 2.0237 1.9163 0 0 0 0 3.770

32、1 -1.3767 -3.7680 4.0978 1.9163 -2.0237 0 0 2.2100 2.5161 0 0 4.6104 -2.0690 -31.0428 -47.0929 0 0 2.5161 -2.2100 0 0 5.1771 1.1403 -47.0929 31.0428 0 0 0 0 0 0 0 0 0 70.0000 0 0 0 0 0 0 0 0 -2.0951 -0.1429E = 1.1201 1.1320 0.8731 0.7174 1.0476 1.0500F =-0.0113 -0.1244 -0.5167 -0.4519 -0.0714 0U = 1

33、.1201 - 0.0113i 1.1320 - 0.1244i 0.8731 - 0.5167i 0.7174 - 0.4519i 1.0476 - 0.0714i 1.0500 dU = 1.0e-005 * 0.1966 -0.0481 0.1552 -0.0453 -0.1877 -0.1040 -0.1272 -0.1305 -0.0000 0PQ =1.0e-004 * -0.0319 -0.0180 0.0939 0.1495 0.0312 -0.0620 -0.0088 -0.0317 0 0.0000precision = 1.4954e-005次數(shù)time= 6雅可比矩陣J

34、J= 2.1354 34.6467 -2.1387 -41.8915 1.0595 4.2260 0 0 0.7430 4.3614 36.4682 1.4175 -41.8915 2.1387 4.2260 -1.0595 0 0 4.3614 -0.7430 0 0 3.0309 3.6468 -2.7627 -7.0223 1.3267 3.4197 0 0 0 0 4.6904 0.0345 -7.0223 2.7627 3.4197 -1.3267 0 0 0 0 0 0 3.2881 0.8064 -4.0978 -3.7681 2.0237 1.9163 0 0 0 0 3.77

35、01 -1.3767 -3.7681 4.0978 1.9163 -2.0237 0 0 2.2100 2.5161 0 0 4.6104 -2.0690 -31.0427 -47.0930 0 0 2.5161 -2.2100 0 0 5.1771 1.1403 -47.0930 31.0427 0 0 0 0 0 0 0 0 0 70.0000 0 0 0 0 0 0 0 0 -2.0951 -0.1429E =1.1201 1.1320 0.8731 0.7174 1.0476 1.0500F = -0.0113 -0.1244 -0.5167 -0.4519 -0.0714 0U =1

36、.1201 - 0.0113i 1.1320 - 0.1244i 0.8731 - 0.5167i 0.7174 - 0.4519i 1.0476 - 0.0714i 1.0500 dU = 1.0e-011 * 0.1748 0.3054 0.1691 0.2967 -0.0946 0.1358 0.0568 0.0980 -0.0000 0PQ = 1.0e-010 *-0.1365 -0.0187 0.1313 0.1907 -0.0803 -0.0356 0.1180 0.0176 0 0.0000precision =1.9074e-011S =-19.5876 +27.9344i線

37、路流動功率S12 =-0.0000 + 2.5240iS25 = -0.0000 + 2.5240iS23 = 4.7133 + 1.7798iS34 = 3.4698 + 1.3835iS45 = 0.3937 + 1.3666iS56 = -28.1644 +18.3564i功率損耗S =-19.5876 +27.9344i總結近兩周的課程設計時間是緊迫的，在此過程首先讓我對所學知識進行重新回顧，明白了使電力系統(tǒng)穩(wěn)定的運行，必須經(jīng)過精密的設計和運算。在進行此次課程設計的過程中，加深了我對潮流計算的認識，尤其是對牛頓拉夫遜潮流計算的求解思路有了比較透徹的理解。通過自己動手查資料，分析原理，繪

38、制電路乃至撰寫文檔感到了設計的難度。綜合檢驗和鞏固自己學到的知識。而且在此次課程設計中，我發(fā)現(xiàn)了自己的基礎知識不足。在以后的學習中應注意基礎知識的掌握。此次設計中采用Matlab進行潮流計算，提高了準確性，使用方便，已得到成功應用。在此，衷心的感謝老師對我們的悉心指導，謝謝老師!19附件計算機算法程序：y11=-33.33333j;y12=31.74603j;y13=0;y14=0;y15=0;y16=0;y21=31.74603j;y22=1.53174-37.41662j;y23=-0.90772+3.78215j;y24=0;y25=-0.62402+3.90016j;y26=0;y31

39、=0;y32=-0.90772+3.78215j;y33=1.7376-6.39418j;y34=-0.8299+3.11203j;y35=0; y36=0;y41=0;y42=0;y43=-0.82988+3.11203j;y44=1.5846-5.25354j;y45=-0.75472+2.64151j;y46=0;y51=0;y52=-0.62402+3.90016j;y53=0;y54=-0.75472+2.64151j;y55=1.37874-66.5103j;y56=63.49206j;y61=0;y62=0;y63=0;y64=0;y65=63.4921j;y66=-66.66667j;Y=y11 y12 y13 y14 y15 y16; y21 y22 y23 y24 y25 y26;y31 y32 y33 y34 y35 y36; y41 y42 y43 y44 y45 y46;y51 y52 y53 y54 y55 y56;y61 y62 y63 y64 y65 y66W=zeros(6); for i=1:6 %導納矩陣第一行放到最后一行，方便下面