1、<p><b>　　課程設(shè)計(論文)</b></p><p>　　題 目 名 稱 電力系統(tǒng)潮流計算 </p><p>　　課 程 名 稱 電力系統(tǒng)穩(wěn)態(tài)分析 </p><p>　　學(xué) 生 姓 名 </p&g

2、t;<p>　　學(xué) 號 </p><p>　　系 、專 業(yè) 電氣工程系電氣工程及其自動化—電力方向 </p><p>　　指 導(dǎo) 教 師 </p><p><b>　　2009年1月6日</b></p><p

3、><b>　　前 言</b></p><p>　　在如今的社會，電力已經(jīng)成為人們必不可少的需求，而建立結(jié)構(gòu)合理的大型電力系統(tǒng)不僅便于電能生產(chǎn)與消費的集中管理、統(tǒng)一調(diào)度和分配，減少總裝機容量，節(jié)省動力設(shè)施投資，且有利于地區(qū)能源資源的合理開發(fā)利用，更大限度地滿足地區(qū)國民經(jīng)濟日益增長的用電需要。電力系統(tǒng)建設(shè)往往是國家及地區(qū)國民經(jīng)濟發(fā)展規(guī)劃的重要組成部分。</p><

4、p>　　電力系統(tǒng)的出現(xiàn)，使高效、無污染、使用方便、易于調(diào)控的電能得到廣泛應(yīng)用，推動了社會生產(chǎn)各個領(lǐng)域的變化，開創(chuàng)了電力時代，發(fā)生了第二次技術(shù)革命。電力系統(tǒng)的規(guī)模和技術(shù)水準(zhǔn)已成為一個國家經(jīng)濟發(fā)展水平的標(biāo)志之一。</p><p>　　電力系統(tǒng)穩(wěn)態(tài)分析包括潮流計算（或潮流分析)和靜態(tài)安全分析。潮流計算針對電力革統(tǒng)各正常運行方式，而靜態(tài)安全分析則要研究各種運行方式下個別系統(tǒng)元件退出運行后系統(tǒng)的狀況。其目的是校驗系統(tǒng)

5、是否能安全運行，即是否有過負荷的元件或電壓過低的母線等。原則上講，靜態(tài)安全分析也可U用潮流計算來代替。但是一般靜態(tài)安全分析需要校驗的狀態(tài)數(shù)非常多，用嚴格的潮流計算來分析這些狀態(tài)往往計算量過大，因此不得不尋求一些特殊的算法以滿足要求。</p><p>　　牛頓法是數(shù)學(xué)中解決非線性方程式的典型方法，有較好的收斂性。解決電力系統(tǒng)潮流計算問題是以導(dǎo)納距陣為基礎(chǔ)的，因此，只要在迭代過程中盡可能保持方程式系數(shù)距陣的稀疏性，就

6、可以大大提高牛頓法潮流程序的放率。自從20 世紀(jì)60 年代中期利用了最佳順序消去法以后，牛頓法在收斂性、內(nèi)存要求、速度方面都超過了阻抗法，成為直到目前仍在廣泛采用的優(yōu)秀方法。</p><p><b>　　目 錄</b></p><p>　　第一章 系統(tǒng)概述4</p><p>　　1.1 設(shè)計目的與要求4</p><

7、p>　　1.1.1 設(shè)計目的4</p><p>　　1.1.2 設(shè)計要求4</p><p>　　1.2 設(shè)計題目4</p><p>　　1.3 設(shè)計內(nèi)容4</p><p>　　第二章 潮流計算設(shè)計題目5</p><p>　　2.1 潮流計算題目5</p><p>　　2.2 對

8、課題的分析及求解思路6</p><p>　　第三章 潮流計算算法及手工計算6</p><p>　　3.1 潮流計算算法6</p><p>　　3.2 關(guān)于電力系統(tǒng)潮流計算手工計算8</p><p>　　3.2.1 節(jié)點導(dǎo)納矩陣8</p><p>　　3.2.2 簡化雅可比矩陣9</p>&

9、lt;p>　　3.2.3修正、迭代10</p><p>　　第四章 Matlab概述10</p><p>　　4.1 Matlab簡介10</p><p>　　4.2 矩陣的運算11</p><p>　　4.2.1 四則運算12</p><p>　　4.2.2 與常數(shù)的運算12</p>

10、<p>　　4.2.3 基本數(shù)學(xué)運算12</p><p>　　4.2.4 邏輯關(guān)系運算12</p><p>　　第五章 潮流計算流程圖及源程序13</p><p>　　5.1 潮流計算流程圖13</p><p>　　5.2 潮流計算源程序14</p><p>　　5.3 運行計算結(jié)果19<

11、/p><p><b>　　總 結(jié)20</b></p><p><b>　　參考文獻21</b></p><p><b>　　第一章 系統(tǒng)概述</b></p><p>　　1.1 設(shè)計目的與要求</p><p>　　1.1.1 設(shè)計目的</p>

12、;<p>　　1. 掌握電力系統(tǒng)潮流計算的基本原理；</p><p>　　2. 掌握并能熟練運用一門計算機語言（MATLAB語言或FORTRAN或C語言或C++語言）；</p><p>　　3. 采用計算機語言對潮流計算進行計算機編程計算。 </p><p>　　1.1.2 設(shè)計要求</p><p>　　1. 程序源代碼； &l

13、t;/p><p>　　2. 給定題目的輸入，輸出文件； </p><p><b>　　3. 程序說明； </b></p><p>　　4. 給定系統(tǒng)的程序計算過程； </p><p>　　5. 給定系統(tǒng)的手算過程（至少迭代2次）。 </p><p><b>　　1.2 設(shè)計題目</b&

14、gt;</p><p>　　電力系統(tǒng)潮流計算（牛頓-拉夫遜法、P-Q 分解法）</p><p><b>　　1.3 設(shè)計內(nèi)容</b></p><p>　　1.根據(jù)電力系統(tǒng)網(wǎng)絡(luò)推導(dǎo)電力網(wǎng)絡(luò)數(shù)學(xué)模型，寫出節(jié)點導(dǎo)納矩陣； </p><p>　　2.賦予各節(jié)點電壓變量（直角坐標(biāo)系形式）初值后，求解不平衡量；</p>

15、<p>　　3.形成雅可比矩陣； </p><p>　　4.求解修正量后，重新修改初值，從2開始重新循環(huán)計算； </p><p>　　5.求解的電壓變量達到所要求的精度時，再計算各支路功率分布、功率損耗和平衡節(jié)</p><p><b>　　點功率; </b></p><p>　　6.上機編程調(diào)試；連調(diào)； &l

16、t;/p><p>　　7.計算分析給定系統(tǒng)潮流分析并與手工計算結(jié)果作比較分析。 </p><p>　　8.準(zhǔn)備計算機演示答辯，書寫該課程設(shè)計說明書（必須計算機打?。?lt;/p><p>　　第二章 潮流計算設(shè)計題目</p><p>　　2.1 潮流計算題目</p><p>　　圖2-1 電力系統(tǒng)接線圖</p>

17、<p>　　2.2 對課題的分析及求解思路</p><p>　　此電力系統(tǒng)是一個5節(jié)點，4支路的電力網(wǎng)絡(luò)。其中包含3個PQ節(jié)點，一個PV節(jié)點，和一個平衡節(jié)點。綜合比較牛頓拉夫遜法（直角坐標(biāo)、極坐標(biāo)）、PQ分解法等多種求解方法的特點，最后確定采用牛頓拉夫遜法（極坐標(biāo)）。因為此方法所需解的方程組最少。</p><p>　　第三章 潮流計算算法及手工計算</p><

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

19、　將上述表示式（1-2）代入（1-1）式的右端，展開并分出實部和虛部，便得：</p><p><b>　?。?-3）</b></p><p>　　按照以上的分類，PQ節(jié)點的輸出有功功率和無功功率是給定的，則第i節(jié)點的給定功率設(shè)為和（稱為注入功率）。</p><p>　　假定系統(tǒng)中的第1、2、…、m節(jié)點為PQ節(jié)點，對其中每一個節(jié)點的N-R法表達式

20、</p><p>　　F(x)=0[如、、]形式有些下列方程：</p><p><b>　?。?-4）</b></p><p>　　=（1、2、…、m）</p><p>　　PV節(jié)點的有功功率和節(jié)點電壓幅值是給定的。假定系統(tǒng)中的第m+1、m+2、…、n-1節(jié)點為PV節(jié)點，則對其中每一PV節(jié)點可以列寫方程：</p&g

21、t;<p><b>　　（1-5）</b></p><p>　　=（m+1、m+2、…、n-1）</p><p>　　(6)形成雅可比矩陣。N-R法的思想是；本例；對F(x)求偏導(dǎo)的式（1-6）、式（1-7），即式（1-4）、式（1-5）中的、、是多維變量的函數(shù)，對多維變量求偏導(dǎo)（、、、、、、、…），并以矩陣的形式表達稱為雅可比矩陣。</p>

22、<p>　　當(dāng)j=i時，對角元素為</p><p><b>　?。?-6）</b></p><p>　　當(dāng)時,矩陣非對角元素為：</p><p><b>　?。?-7）</b></p><p>　　由上式不難看出，雅可比矩陣有以下特點。</p><p>　　雅可

23、比矩陣中的諸元素都是節(jié)點電壓的函數(shù)，因此在迭代過程中，它們將隨著節(jié)點電壓的變化而不斷的變化。</p><p>　　雅可比矩陣具有結(jié)構(gòu)對稱性，數(shù)據(jù)不對稱。如非對角，，。</p><p>　　由式（1-7）可以看出，當(dāng)導(dǎo)納矩陣中非對角元素為零時，。雅可比矩陣中相應(yīng)的元素也為零，即矩陣是非常稀疏的。因此，修正方程的求解同樣可以應(yīng)用稀疏矩陣的求解技巧。正是由于這一點才使N-R法獲得廣泛的應(yīng)用。&l

24、t;/p><p>　　3.2 關(guān)于電力系統(tǒng)潮流計算手工計算</p><p>　　3.2.1 節(jié)點導(dǎo)納矩陣</p><p><b>　　求得節(jié)點導(dǎo)納矩陣Y</b></p><p><b>　　=</b></p><p>　　各節(jié)點的導(dǎo)納值如下：</p><p

25、>　　; ；</p><p>　?。?;</p><p>　　; ;</p><p>　??； ;</p><p>　??；

26、 ;</p><p>　　; ;</p><p>　　; ;</p><p>　　; ;</p><p>　　; ;</p><p&g

27、t;　　; ;</p><p><b>　　;</b></p><p><b>　　;</b></p><p><b>　　;</b></p><p><b>　　;</b></p&

28、gt;<p><b>　　.</b></p><p>　　3.2.2 簡化雅可比矩陣</p><p>　　形成有功迭代和無功迭代的簡化雅可比矩陣B/和B//</p><p><b>　　B/= </b></p><p><b>　　B//=</b><

29、/p><p>　　將B/ 和B//進行三角分解：</p><p>　　3.2.3 修正、迭代</p><p>　　給定PQ節(jié)點初值和各節(jié)點電壓相角初值</p><p>　　V1=1.05∠0。 ，V2(0)=V3(0)=1.0，V4=1.1</p><p>　　δ2(0)=δ3(0)=0， δ4(0)=0</p&g

30、t;<p>　　1 作第一次有功迭代，按公式計算節(jié)點有功功率不平衡量</p><p>　　△P2(0)=-0.55-（-0.024037）=-0.525963</p><p>　　△P3(0)=-0.30-（-0.022695）=-0.277305</p><p>　　△P4(0)=0.500000</p><p>　　△P1(

31、0)/V1(0)=0.454545 △P2(0)/ V2(0)=-0.525963</p><p>　　△P3(0)/V3(0)=-0.277309</p><p>　　2做第一次無功迭代，按公式計算無功功率不平衡量，計算時電壓相角最新的修正值。</p><p>　　△Q2(0)=-0.13-(-0.001550)=-0.039594</p>

32、<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</p><p>　　解修正方程式，可得各節(jié)點電壓幅值的修正量為</p><p>　　△V3(0)）=-0.014855&

33、lt;/p><p><b>　　于是有：</b></p><p>　　V2(1) = V2(0)+△V2(1)=0.964776</p><p>　　V3(1) = V3(0)+△V3(1)=0.985145</p><p>　　到這里為止，第一輪有功迭代和無功迭代便做完了。</p><p>　　3

34、按公式計算平衡節(jié)點功率，得：</p><p>　　P1+jQ1=0.367885+j0.264696</p><p>　　經(jīng)過四輪迭代，節(jié)點不平衡功率也下降到10-5以下，迭代到此結(jié)束。</p><p>　　第四章 Matlab概述</p><p>　　4.1 Matlab簡介</p><p>　　目前電子計算機已廣

35、泛應(yīng)用于電力系統(tǒng)的分析計算，潮流計算是其基本應(yīng)用軟件之一。現(xiàn)有很多潮流計算方法。對潮流計算方法有五方面的要求：（1）計算速度快（2）內(nèi)存需要少（3）計算結(jié)果有良好的可靠性和可信性（4）適應(yīng)性好，亦即能處理變壓器變比調(diào)整、系統(tǒng)元件的不同描述和與其它程序配合的能力強（5）簡單。</p><p>　　MATLAB是一種交互式、面向?qū)ο蟮某绦蛟O(shè)計語言，廣泛應(yīng)用于工業(yè)界與學(xué)術(shù)界，主要用于矩陣運算，同時在數(shù)值分析、自動控制模

36、擬、數(shù)字信號處理、動態(tài)分析、繪圖等方面也具有強大的功能。</p><p>　　MATLAB程序設(shè)計語言結(jié)構(gòu)完整，且具有優(yōu)良的移植性，它的基本數(shù)據(jù)元素是不需要定義的數(shù)組。它可以高效率地解決工業(yè)計算問題，特別是關(guān)于矩陣和矢量的計算。MATLAB與C語言和FORTRAN語言相比更容易被掌握。通過M語言，可以用類似數(shù)學(xué)公式的方式來編寫算法，大大降低了程序所需的難度并節(jié)省了時間，從而可把主要的精力集中在算法的構(gòu)思而不是編程

37、上。</p><p>　　另外，MATLAB提供了一種特殊的工具：工具箱（TOOLBOXES）.這些工具箱主要包括：信號處理（SIGNAL PROCESSING）、控制系統(tǒng)（CONTROL SYSTEMS）、神經(jīng)網(wǎng)絡(luò)（NEURAL NETWORKS）、模糊邏輯(FUZZY LOGIC)、小波(WAVELETS)和模擬（SIMULATION）等等。不同領(lǐng)域、不同層次的用戶通過相應(yīng)工具的學(xué)習(xí)和應(yīng)用，可以方便地進行計算

38、、分析及設(shè)計工作。</p><p>　　MATLAB設(shè)計中，原始數(shù)據(jù)的填寫格式是很關(guān)鍵的一個環(huán)節(jié)，它與程序使用的方便性和靈活性有著直接的關(guān)系。</p><p>　　原始數(shù)據(jù)輸入格式的設(shè)計，主要應(yīng)從使用的角度出發(fā)，原則是簡單明了，便于修改。</p><p>　　4.2.2 與常數(shù)的運算 </p><p>　　常數(shù)與矩陣的運算即是同該矩陣

39、的每一元素進行運算。但需注意進行數(shù)除時，常數(shù)通常只能做除數(shù)。</p><p>　　基本函數(shù)運算中，矩陣的函數(shù)運算是矩陣運算中最實用的部分，常用的主要有以下幾個：</p><p>　　det(a) 求矩陣a的行列式</p><p>　　eig(a) 求矩陣a的特征值</p

40、><p>　　inv(a)或a ^ (-1) 求矩陣a的逆矩陣</p><p>　　rank(a) 求矩陣a的秩</p><p>　　trace(a) 求矩陣a的跡（對角線元素之和）</p><p>　　我們在進行工程計算時常常遇

41、到矩陣對應(yīng)元素之間的運算。這種運算不同于前面講的數(shù)學(xué)運算，為有所區(qū)別，我們稱之為數(shù)組運算。</p><p>　　4.2.3 基本數(shù)學(xué)運算</p><p>　　數(shù)組的加、減與矩陣的加、減運算完全相同。而乘除法運算有相當(dāng)大的區(qū)別，數(shù)組的乘除法是指兩同維數(shù)組對應(yīng)元素之間的乘除法，它們的運算符為“.*”和“./”或“.\”。前面講過常數(shù)與矩陣的除法運算中常數(shù)只能做除數(shù)。在數(shù)組運算中有了“對應(yīng)關(guān)系”

42、的規(guī)定，數(shù)組與常數(shù)之間的除法運算沒有任何限制。</p><p>　　另外，矩陣的數(shù)組運算中還有冪運算（運算符為 .^ ）、指數(shù)運算（exp）、對數(shù)運算（log）、和開方運算（sqrt）等。有了“對應(yīng)元素”的規(guī)定，數(shù)組的運算實質(zhì)上就是針對數(shù)組內(nèi)部的每個元素進行的。矩陣的冪運算與數(shù)組的冪運算有很大的區(qū)別。</p><p>　　4.2.4 邏輯關(guān)系運算 </p><p

43、>　　邏輯運算是MATLAB中數(shù)組運算所特有的一種運算形式，也是幾乎所有的高級語言普遍適用的一種運算。</p><p>　　第五章 潮流計算流程圖及源程序</p><p>　　5.1 潮流計算流程圖</p><p>　　圖5-1 潮流計算流程圖</p><p>　　5.2 潮流計算源程序</p><p>　　據(jù)課

44、題題目，本程序把節(jié)點1設(shè)為平衡節(jié)點，節(jié)點2、3、4為PQ節(jié)點，節(jié)點5為PV節(jié)點。</p><p>　　G(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

45、.667;</p><p>　　B(1,3)=5.000;</p><p>　　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

46、.667;</p><p>　　B(2,1)=5.000;</p><p>　　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>

47、;　　G(2,4)=0;</b></p><p><b>　　B(2,4)=0;</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;&

48、lt;/p><p>　　G(3,2)=-10.000;</p><p>　　B(3,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.75

49、0;</p><p><b>　　G(3,5)=0;</b></p><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;&l

50、t;/b></p><p><b>　　B(4,2)=0;</b></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>&l

51、t;p><b>　　G(4,5)=0;</b></p><p><b>　　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>

52、　　B(5,2)=3.750;</p><p><b>　　G(5,3)=0;</b></p><p><b>　　B(5,3)=0;</b></p><p><b>　　G(5,4)=0;</b></p><p><b>　　B(5,4)=0;</b>&

53、lt;/p><p>　　G(5,5)=6.250;</p><p>　　B(5,5)=-18.750;</p><p><b>　　Y=G+j*B;</b></p><p>　　delt(1)=0;</p><p>　　delt(2)=0;</p><p>　　delt(3)=

54、0;</p><p>　　delt(4)=0;</p><p><b>　　u(1)=1.0;</b></p><p><b>　　u(2)=1.0;</b></p><p><b>　　u(3)=1.0;</b></p><p><b>　　

55、u(4)=1.0;</b></p><p>　　p(1)=0.20;</p><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

56、;</p><p>　　p(4)=-0.60;</p><p>　　q(4)=-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<

57、/p><p>　　while precision>0.00001</p><p>　　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

58、(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));</p><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>　

59、　pp(m)=p(m)-sum(pt);</p><p>　　qq(m)=q(m)-sum(qt);</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

60、)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));</p><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

61、(m,n)*sin(delt(m)-delt(n)));</p><p>　　L0(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>　　H(m,m)=sum(h0)-u(m)^2*(G(m

62、,m)*sin(delt(m)-delt(m))-B(m,m)*cos(delt(m)-delt(m</p><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)));</

63、p><p>　　J(m,m)=sum(j0)+u(m)^2*(G(m,m)*cos(delt(m)-delt(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)));<

64、;/p><p><b>　　end</b></p><p>　　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>

65、<p>　　JJ(2*m,2*m)=L(m,m);</p><p><b>　　end</b></p><p>　　for m=1:N1</p><p>　　for n=1:N1</p><p><b>　　if m==n</b></p><p><b>

66、;　　else</b></p><p>　　H(m,n)=-u(m)*u(n)*(G(m,n)*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

67、><p>　　N(m,n)=-J(m,n);</p><p>　　L(m,n)=H(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

68、>　　JJ(2*m,2*n)=L(m,n);</p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　end</b></p><p>　　for m=1:N1</p><p>　　PP(2

69、*m-1)=pp(m);</p><p>　　PP(2*m)=qq(m);</p><p><b>　　end</b></p><p>　　uu=-inv(JJ)*PP';</p><p>　　precision=max(abs(uu));</p><p>　　for n=1:N1<

70、/p><p>　　delt(n)=delt(n)+uu(2*n-1);</p><p>　　u(n)=u(n)+uu(2*n);</p><p><b>　　end</b></p><p><b>　　k=k+1;</b></p><p><b>　　end</b

71、></p><p>　　K=k-1,delt,u'</p><p>　　%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

72、(n)));</p><p><b>　　end</b></p><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

73、><p>　　for m=1:N1+1</p><p>　　for n=1:N1+1</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<

74、;/b></p><p><b>　　Y</b></p><p><b>　　S</b></p><p>　　5.3 運行計算結(jié)果</p><p><b>　　K =4；</b></p><p>　　delt =[ -0.0461 -0.083

75、9 -0.0896 -0.1044 0]；</p><p>　　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.5

76、000i -1.6670 + 5.0000i -1.6670 + 5.0000i -2.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>　　-

77、1.6670 + 5.0000i -10.0000 +30.0000i 12.9170 -38.7500i -1.2500 + 3.7500i 0 </p><p>　　-2.5000 + 7.5000i 0 -1.2500 + 3.7500i 3.7500 -11.2500i 0 </p>

78、<p>　　-5.0000 +15.0000i -1.2500 + 3.7500i 0 0 6.2500 -18.7500i</p><p><b>　　S =</b></p><p>　　0 0.2469 + 0.0815i 0.2793 + 0.0806

79、i 0.5489 + 0.1333i -0.8751 - 0.0954i</p><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

80、 0 0.0633 + 0.0033i 0 </p><p>　　-0.5370 - 0.0977i 0 -0.0630 - 0.0023i 0 0 </p><p>　　0.8895 + 0.1387i 0.4087 + 0.

81、1058i 0 0 0 </p><p><b>　　總 結(jié)</b></p><p>　　此次課程設(shè)計首先讓我明白了要使電力系統(tǒng)運行的穩(wěn)定，必須經(jīng)過精密的設(shè)計和計算。在進行課題設(shè)計的過程中，加深了我對潮流計算的認識，尤其是對牛頓拉夫遜潮流計算的求解思路有了比較透徹的理

82、解。同時由于求解過程中用到求節(jié)點導(dǎo)鈉矩陣，求矩陣的逆等等，又使我對以前所學(xué)的知識有了一次很好的溫習(xí)。同時也看到了研究性學(xué)習(xí)的效果，從研究中去學(xué)習(xí)，理論結(jié)合實際，將理論運用到實際，同時在實踐中發(fā)現(xiàn)問題，然后解決問題，</p><p>　　而且在此次課程設(shè)計中，我發(fā)現(xiàn)了自己的基礎(chǔ)知識有很多的不足。這些基礎(chǔ)的缺乏給我的設(shè)計計劃造成了不小的障礙。在這個過程中，我明白了，只要用心去做，認真去做，持之以恒，就會有新的發(fā)現(xiàn)，有

83、意外的收獲。</p><p>　　此次課程設(shè)計的完成，要感謝老師的幫助以及同學(xué)的協(xié)力合作，團結(jié)才能出成績。 </p><p><b>　　參考文獻</b></p><p>　　[1] 何仰贊等.電力系統(tǒng)分析[M]. 武漢：華中理工大學(xué)出版社，2002.3 </p><p>　　[2] 西安交通大學(xué)等.電力系統(tǒng)計算[M].

84、北京：水利電力出版社，1993.12 </p><p>　　[3] 陳 衍.電力系統(tǒng)穩(wěn)態(tài)分析[M].北京：水利電力出版社，2004.1 </p><p>　　[4] 李光琦.電力系統(tǒng)暫態(tài)分析[M].北京： 水利電力出版社，2002.5 </p><p>　　[5] 于永源，楊綺雯. 電力系統(tǒng)分析（第二版）[M]. 北京：中國電力出版社，2004.3</p>