電力系統(tǒng)穩(wěn)態(tài)潮流計(jì)算課程設(shè)計(jì)

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

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

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

4、p>　　電力系統(tǒng)的出現(xiàn)，使高效、無(wú)污染、使用方便、易于調(diào)控的電能得到廣泛應(yīng)用，推動(dòng)了社會(huì)生產(chǎn)各個(gè)領(lǐng)域的變化，開創(chuàng)了電力時(shí)代，發(fā)生了第二次技術(shù)革命。電力系統(tǒng)的規(guī)模和技術(shù)水準(zhǔn)已成為一個(gè)國(guó)家經(jīng)濟(jì)發(fā)展水平的標(biāo)志之一。</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)

5、是否能安全運(yùn)行，即是否有過負(fù)荷的元件或電壓過低的母線等。原則上講，靜態(tài)安全分析也可U用潮流計(jì)算來(lái)代替。但是一般靜態(tài)安全分析需要校驗(yàn)的狀態(tài)數(shù)非常多，用嚴(yán)格的潮流計(jì)算來(lái)分析這些狀態(tài)往往計(jì)算量過大，因此不得不尋求一些特殊的算法以滿足要求。</p><p>　　牛頓法是數(shù)學(xué)中解決非線性方程式的典型方法，有較好的收斂性。解決電力系統(tǒng)潮流計(jì)算問題是以導(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è)計(jì)目的與要求4</p><

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

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

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

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

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

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

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

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

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

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

17、<p>　　2.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>　　第三章 潮流計(jì)算算法及手工計(jì)算</p><

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

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

20、</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+1、m+2、…、n-1節(jié)點(diǎn)為PV節(jié)點(diǎn)，則對(duì)其中每一PV節(jié)點(diǎn)可以列寫方程：</p&g

21、t;<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ì)多維變量求偏導(dǎo)（、、、、、、、…），并以矩陣的形式表達(dá)稱為雅可比矩陣。</p>

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

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

24、t;/p><p>　　3.2 關(guān)于電力系統(tǒng)潮流計(jì)算手工計(jì)算</p><p>　　3.2.1 節(jié)點(diǎn)導(dǎo)納矩陣</p><p><b>　　求得節(jié)點(diǎn)導(dǎo)納矩陣Y</b></p><p><b>　　=</b></p><p>　　各節(jié)點(diǎn)的導(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 簡(jiǎn)化雅可比矩陣</p><p>　　形成有功迭代和無(wú)功迭代的簡(jiǎn)化雅可比矩陣B/和B//</p><p><b>　　B/= </b></p><p><b>　　B//=</b><

29、/p><p>　　將B/ 和B//進(jìn)行三角分解：</p><p>　　3.2.3 修正、迭代</p><p>　　給定PQ節(jié)點(diǎn)初值和各節(jié)點(diǎn)電壓相角初值</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 作第一次有功迭代，按公式計(jì)算節(jié)點(diǎn)有功功率不平衡量</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做第一次無(wú)功迭代，按公式計(jì)算無(wú)功功率不平衡量，計(jì)算時(shí)電壓相角最新的修正值。</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é)點(diǎn)電壓幅值的修正量為</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>　　到這里為止，第一輪有功迭代和無(wú)功迭代便做完了。</p><p>　　3

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

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

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

37、上。</p><p>　　另外，MATLAB提供了一種特殊的工具：工具箱（TOOLBOXES）.這些工具箱主要包括：信號(hào)處理（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)用，可以方便地進(jìn)行計(jì)算

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

39、的每一元素進(jìn)行運(yùn)算。但需注意進(jìn)行數(shù)除時(shí)，常數(shù)通常只能做除數(shù)。</p><p>　　基本函數(shù)運(yùn)算中，矩陣的函數(shù)運(yùn)算是矩陣運(yùn)算中最實(shí)用的部分，常用的主要有以下幾個(gè)：</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的跡（對(duì)角線元素之和）</p><p>　　我們?cè)谶M(jìn)行工程計(jì)算時(shí)常常遇

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

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

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

44、題題目，本程序把節(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(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 運(yùn)行計(jì)算結(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è)計(jì)首先讓我明白了要使電力系統(tǒng)運(yùn)行的穩(wěn)定，必須經(jīng)過精密的設(shè)計(jì)和計(jì)算。在進(jìn)行課題設(shè)計(jì)的過程中，加深了我對(duì)潮流計(jì)算的認(rèn)識(shí)，尤其是對(duì)牛頓拉夫遜潮流計(jì)算的求解思路有了比較透徹的理

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

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