n-r潮流計算分析電力系統(tǒng)課程設(shè)計

1、<p><b>　　課程設(shè)計任務(wù)書</b></p><p><b>　　摘 要</b></p><p>　　電力系統(tǒng)穩(wěn)態(tài)分析包括潮流計算和靜態(tài)安全分析。本文主要運用的事潮流計算，潮流計算是電力網(wǎng)絡(luò)設(shè)計與運行中最基本的運算，對電力網(wǎng)絡(luò)的各種設(shè)計方案及各種運行方式進行潮流計算，可以得到各種電網(wǎng)各節(jié)點的電壓，并求得網(wǎng)絡(luò)的潮流及網(wǎng)絡(luò)中的各元件

2、的電力損耗，進而求得電能損耗。本位就是運用潮流計算具體分析，并有MATLAB仿真。 </p><p>　　關(guān)鍵詞： 電力系統(tǒng) 潮流計算 MATLAB仿真</p><p><b>　　目錄</b></p><p><b>　　1.052</b></p><p><b>　　摘

3、 要5</b></p><p><b>　　目錄6</b></p><p>　　第1章 電力系統(tǒng)的基本概念6</p><p>　　第二章 潮流計算7</p><p>　　2.1潮流計算概述與發(fā)展8</p><p>　　2.2復(fù)雜電力系統(tǒng)潮流計算8</p>&

4、lt;p>　　2.3潮流計算的方法及優(yōu)、缺點9</p><p>　　2.4 MATLAB概述9</p><p>　　2.6牛頓-拉夫遜法原理10</p><p>　　2.7牛頓-拉夫遜法解決潮流計算問題12</p><p>　　2.8計算機潮流計算的步驟13</p><p>　　2.9計算機程序的實現(xiàn)

5、14</p><p>　　程序運行結(jié)果：16</p><p>　　第四章 學(xué)習(xí)心得29</p><p><b>　　附錄：30</b></p><p>　　第1章 電力系統(tǒng)的基本概念</p><p>　　電力系統(tǒng)：發(fā)電機把機械能轉(zhuǎn)化為電能，電能經(jīng)變壓器和電力線路輸送并分配到用戶，在那里經(jīng)電

6、動機、電爐和電燈等設(shè)備又將電能轉(zhuǎn)化為機械能、熱能和光能等。這些生產(chǎn)、變換、輸送、分配、消費電能的發(fā)電機、變壓器、變換器、電力線路及各種用電設(shè)備等聯(lián)系在一起組成的統(tǒng)一整體稱為電力系統(tǒng)。電力網(wǎng)：電力系統(tǒng)中除發(fā)電機和用電設(shè)備外的部分。動力系統(tǒng)：電力系統(tǒng)和“動力部分”的總和。</p><p>　　電力系統(tǒng)潮流計算是研究電力系統(tǒng)穩(wěn)態(tài)運行情況的一種計算，它根據(jù)給定的運行條件及系統(tǒng)接線情況確定整個電力系統(tǒng)各部分的運行狀態(tài)：各母

7、線的電壓，各元件中流過的功率，系統(tǒng)的功率損耗等等。在電力系統(tǒng)規(guī)劃的設(shè)計和現(xiàn)有電力系統(tǒng)運行方式的研究中，都需要利用潮流計算來定量地分析比較供電方案或運行方式的合理性。可靠性和經(jīng)濟性。此外，電力系統(tǒng)潮流計算也是計算系統(tǒng)動態(tài)穩(wěn)定和靜態(tài)穩(wěn)定的基礎(chǔ)。所以潮流計算是研究電力系統(tǒng)的一種很重要和很基礎(chǔ)的計算。</p><p>　　隨著科學(xué)技術(shù)的發(fā)展，電力系統(tǒng)變得越來越復(fù)雜，電氣工程師掌握一種好的能對電力系統(tǒng)進行仿真的軟件是學(xué)習(xí)和

8、研究的需要。文章簡要介紹了MATLAB發(fā)展歷史、組成和強大的功能，并用簡單例子分別就編程和仿真兩方面分析了MATIAB軟件在電力系統(tǒng)研究中的具體應(yīng)用。采取等效電路法，能對特殊、復(fù)雜地電力系統(tǒng)進行高效仿真研究，因此，掌握編程和仿真是學(xué)好MATLAB的基礎(chǔ)。與眾多專門的電力系統(tǒng)仿真軟件相比，MATLAB軟件具有易學(xué)、功能強大和開放性好，是電力系統(tǒng)仿真研究的有力工具。</p><p><b>　　第2章 潮流

9、計算</b></p><p>　　2.1潮流計算概述與發(fā)展</p><p>　　電力系統(tǒng)潮流計算也分為離線計算和在線計算兩種，前者主要用于系統(tǒng)規(guī)劃設(shè)計和安排系統(tǒng)的運行方式，后者則用于正在運行系統(tǒng)的經(jīng)常監(jiān)視及實時控制。</p><p>　　利用電子數(shù)字計算機進行電力系統(tǒng)潮流計算從50年代中期就已經(jīng)開始。在這20年內(nèi)，潮流計算曾采用了各種不同的方法，這些方法

10、的發(fā)展主要圍繞著對潮流計算的一些基本要求進行的。</p><p>　　牛頓-拉夫遜法作為一種實用的，有競爭力的電力系統(tǒng)潮流計算方法，是在應(yīng)用了稀疏矩陣技巧和高斯消去法求修正方程后。牛頓-拉夫遜法是求解非線性代數(shù)方程有效的迭代計算。</p><p>　　2.2復(fù)雜電力系統(tǒng)潮流計算</p><p>　　電力系統(tǒng)潮流計算是對復(fù)雜電力系統(tǒng)正常和故障條件下穩(wěn)態(tài)運行狀態(tài)的計算。

11、潮流計算的目標(biāo)是求取電力系統(tǒng)在給定運行方式下的節(jié)點電壓和功率分布，用以檢查系統(tǒng)各元件是否過負(fù)荷、各點電壓是否滿足要求、功率的分布和分配是否合理以及功率損耗等。對現(xiàn)有電力系統(tǒng)的運行和擴建，對新的電力系統(tǒng)進行規(guī)劃設(shè)計都是以潮流計算為基礎(chǔ)。</p><p>　　潮流計算結(jié)果的用途，例如用于電力系統(tǒng)穩(wěn)定研究、安全估計或最優(yōu)潮流等也對潮流計算的模型和方法有直接影響。</p><p><b>

12、;　　節(jié)點類型：</b></p><p>　　1）PV節(jié)點：柱入有功功率P為給定值，電壓也保持在給定數(shù)值。</p><p>　　2）PQ節(jié)點：諸如有功功率和無功功率是給定的。</p><p>　　3）平衡節(jié)點：用來平衡全電網(wǎng)的功率。選一容量足夠大的發(fā)電機擔(dān)任平衡全電網(wǎng)功率的職責(zé)。平衡節(jié)點的電壓大小與相位是給定的，通常以它的相角為參考量，即取其電壓相角為0

13、。一個獨立的電力網(wǎng)中只設(shè)一個平衡點。</p><p><b>　　基本步驟：</b></p><p>　　1）形成節(jié)點導(dǎo)納矩陣；</p><p>　　2）將各節(jié)點電壓設(shè)初值U；</p><p>　　3）將節(jié)點初值代入相關(guān)求式，求出修正方程式的常數(shù)項向量；</p><p>　　4）將節(jié)點電壓初值代入

14、求式，求出雅可比矩陣元素；</p><p>　　5）求解修正方程，求修正向量；</p><p>　　6）求取節(jié)點電壓的新值；</p><p>　　7）檢查是否收斂，如不收斂，則以各節(jié)點電壓的新值作為初值自第3步重新開始進行狹義次迭代，否則轉(zhuǎn)入下一步；</p><p>　　8）計算支路功率分布，PV節(jié)點無功功率和平衡節(jié)點柱入功率。</p&

15、gt;<p>　　2.3潮流計算的方法及優(yōu)、缺點</p><p>　　潮流計算法有，簡化梯度法、二次規(guī)劃法、牛頓-拉夫遜法(Newton-Raphson)等。簡化梯度法是采用梯度法進行搜索，用罰函數(shù)處理違約的不等式約束。該方法程序編制簡便，所需存儲量小，對初始點無特殊要求，曾獲得普遍重視，成為第一種有效的優(yōu)化潮流方法。簡化梯度法的缺點：迭代過程中，尤其是在接近最優(yōu)點附近會出現(xiàn)鋸齒現(xiàn)象，收斂性較差，收

16、斂速度很慢；每次迭代都要重新計算潮流，計算量很大，耗時較多。二次規(guī)劃法是二階的方法，解決最優(yōu)潮流問題收斂精度較好，能很好地解決耦合的最優(yōu)潮流問題，但缺點是計算Lagrange函數(shù)的二階偏導(dǎo)數(shù)，計算量大、計算復(fù)雜。</p><p>　　2.4 MATLAB概述</p><p>　　目前電子計算機已廣泛應(yīng)用于電力系統(tǒng)的分析計算，潮流計算是其基本應(yīng)用之一?，F(xiàn)有很多潮流計算方法。對潮流計算方法有五

17、方面的要求：</p><p><b>　?。?）計算速度快；</b></p><p><b>　?。?）內(nèi)存需要少；</b></p><p>　　（3）計算結(jié)果有良好的可靠性和可信性；</p><p>　　（4）適應(yīng)性好，亦即能處理變壓器變比調(diào)整、系統(tǒng)元件的不同描述和與其它程序配合的能力強；<

18、/p><p><b>　?。?）簡單。</b></p><p>　　MATLAB是一種交互式、面向?qū)ο蟮某绦蛟O(shè)計語言，廣泛應(yīng)用于工業(yè)界與學(xué)術(shù)界，主要用于矩陣運算，同時在數(shù)值分析、自動控制模擬、數(shù)字信號處理、動態(tài)分析、繪圖等方面也具有強大的功能。</p><p>　　MATLAB程序設(shè)計語言結(jié)構(gòu)完整，且具有優(yōu)良的移植性，它的基本數(shù)據(jù)元素是不需要定義的

19、數(shù)組。它可以高效率地解決工業(yè)計算問題，特別是關(guān)于矩陣和矢量的計算。MATLAB與C語言和FORTRAN語言相比更容易被掌握。通過M語言，可以用類似數(shù)學(xué)公式的方式來編寫算法，大大降低了程序所需的難度并節(jié)省了時間，從而可把主要的精力集中在算法的構(gòu)思而不是編程上。</p><p>　　2.6牛頓-拉夫遜法原理</p><p>　　假設(shè)有n個聯(lián)立的非線性代數(shù)方程：</p><p

20、>　　假設(shè)以給出各變量的初值，，……，，令其分別為個變量的修正量，使?jié)M足以上方程，所以：</p><p>　　將上式中的n個多元函數(shù)在初始值附近分別展開成泰勒級數(shù)，并略去含有，，……，的二次及以上階次的各項，便得：</p><p><b>　　方程可寫成：</b></p><p>　　以上方程是對于修正量，，……，的線性方程組，稱為牛頓法

21、的修正方程，可解出，，……，。對初始近似解進行修正：</p><p>　?。╥=1,2,……，n）</p><p>　　反復(fù)迭代，在進行k+1次迭代時，從求解修正方程式：</p><p>　　得到修正量,,……,,對各量進行修正 （i=1,2,……，n）迭代過程一直進行到滿足收斂判據(jù)</p><p>　　2.7牛頓-拉夫遜法解決潮流計算問

22、題</p><p>　　節(jié)點總數(shù)為n；PQ節(jié)點有m，；PV節(jié)點有n-m-1，平衡節(jié)點有1個，節(jié)點編號按照先PQ節(jié)點，再PV節(jié)點，最后平衡節(jié)點的順序進行編號，即:1，2，…，m為PQ節(jié)點；m+1，m+2，…，n-1為PV節(jié)點；n為平衡節(jié)點。可形成結(jié)點導(dǎo)納矩陣。導(dǎo)納矩陣元素可表示為，本文中節(jié)點電壓以直角坐標(biāo)形式表示，即。由此下列公式可求出Pi，Qi</p><p>　　假設(shè)系統(tǒng)中的第1，2，…

23、，m號節(jié)點為PQ節(jié)點，第i個節(jié)點的給定功率為和，對該節(jié)點可列方程：</p><p>　　假設(shè)系統(tǒng)中的第m+1，m+2，…，n-1號節(jié)點為PV節(jié)點，則對其中每一個節(jié)點可列方程:</p><p>　　第n號節(jié)點為平衡節(jié)點，其電壓為是給定的，故不參加迭代。</p><p>　　修正方程可寫成分塊矩陣的形式：</p><p>　　通過反復(fù)求解修正方程

24、，解出各節(jié)點的未知量，再通過收斂判據(jù)判定是否已為真值。從而求得PQ節(jié)點的電壓V及相角δ的真值，PV節(jié)點的Q、δ真值，平衡節(jié)點的P、Q真值，以上即為牛頓-拉夫遜迭代法的潮流計算過程，其優(yōu)點為計算精確，運行速度快。其中的各個環(huán)節(jié)都可通過MATLAB程序來實現(xiàn)。</p><p>　　2.8計算機潮流計算的步驟</p><p>　　（1）對電力網(wǎng)絡(luò)的所有參數(shù)設(shè)初值,包括電壓、相角、有功、無功等。&

25、lt;/p><p>　?。?）處理非標(biāo)準(zhǔn)變比支路,使其變成標(biāo)準(zhǔn)變比為1的變壓器支路。</p><p>　?。?）形成節(jié)點導(dǎo)納矩陣Y。</p><p>　?。?）計算有功功率的不平衡量ΔPi,從而求出。</p><p>　?。?）根據(jù)節(jié)點的類型形成J。</p><p>　　（6）解修正方程式,求各節(jié)點的電壓的變化量Δei(i

26、=1,2,3...n,i≠s)</p><p>　?。?）求各節(jié)點相角的新值ei=ei+Δei (i=1,2,3... n,i≠s)</p><p>　　（8）計算無功功率的不平衡量ΔQi,從而求出 (i=1,2,3...n,i≠s)</p><p>　?。?）解修正方程式,求各節(jié)點的電壓大小的變化量 (i=1,2,3...,n,i≠s)。</p>&

27、lt;p>　　（10）求各節(jié)點的電壓大小的新值 (i=1,2,3...,n,i≠s)。</p><p>　　（11）運用個節(jié)點的電壓的新值自第四步開始下一次迭代。</p><p>　　計算平衡節(jié)點的功率和線路功率。其中平衡節(jié)點的功率的計算公式為</p><p><b>　　線路上的功率為:</b></p><p>

28、　　從而線路上的損耗的功率為:</p><p>　　2.9計算機程序的實現(xiàn)</p><p>　　節(jié)點導(dǎo)納矩陣是方陣,其階數(shù)等于網(wǎng)絡(luò)中出參考節(jié)點外的節(jié)點數(shù)n。</p><p>　　節(jié)點導(dǎo)納矩陣是稀疏矩陣,其各行非零非對角元數(shù)就等于該行相對應(yīng)節(jié)點所連接的不接地支路數(shù)。</p><p>　　節(jié)點導(dǎo)納矩陣的對角元就等于各該節(jié)點所連接導(dǎo)納的總和。<

29、;/p><p>　　節(jié)點導(dǎo)納矩陣的非對角元等于連接節(jié)點i,j支路導(dǎo)納的負(fù)值。</p><p>　　點導(dǎo)納矩陣一般是對稱矩陣。</p><p>　　對于支路中有非標(biāo)準(zhǔn)變比變壓器的支路來說,利用下面的公式來計算它的導(dǎo)納。變壓器的變比,假如已知非標(biāo)準(zhǔn)變比支路i,j上的阻抗(以下沒有特殊說明所有的參數(shù)都用標(biāo)幺值)為,則線路導(dǎo)納為,線路上的對地半導(dǎo)納為。</p>&

30、lt;p><b>　　J的形成</b></p><p>　　Y是由最終形成的導(dǎo)納矩陣的虛部組成的,但是pv節(jié)點以及平衡節(jié)點不參加Q～V迭代,因此Y中不包含與這些節(jié)點有關(guān)的元素。</p><p><b>　　迭代條件和約束方程</b></p><p>　　迭代條件就是如果ΔQ≤時就停止迭代。</p>&l

31、t;p>　　對節(jié)點的約束條件分為三類:即對節(jié)點注入功率的約束、對節(jié)點電壓大小的約束和對相角的約束。其中對節(jié)點注入功率的約束,主要是對電源注入功率的約束條件不能滿足時,將威脅到發(fā)電機的安全運行。對電壓大小的約束不能滿足時,將影響電能的質(zhì)量,嚴(yán)重時將影響系統(tǒng)運行的穩(wěn)定性。對相對相角的約束條件不能滿足時,也將危及系統(tǒng)運行的穩(wěn)定性。</p><p>　　圖2.2 程序流程圖</p><p>

32、　　第3章 潮流計算及運行結(jié)果</p><p>　　3.0 手算潮流計算</p><p>　　用圖1和圖2的數(shù)據(jù)和等值網(wǎng)絡(luò)形成節(jié)點導(dǎo)納矩陣</p><p><b>　　1. 節(jié)點導(dǎo)納矩陣</b></p><p><b>　　程序運行結(jié)果：</b></p><p><b

33、>　　導(dǎo)納矩陣 Y=</b></p><p>　　0 -30.2343i 0 +31.7460i 0 0 0 0 </p><p>　　0 +31.7460i 14.8252 -42.6506i -14.2012 + 5.9172i 0

34、 -0.6240 + 3.9002i 0 0 -14.2012+5.9172i 15.0311-8.5292i -0.8299+3.1120i 0 0 </p><p>　　0 0 -0.8299 + 3.1120i 1.5846 - 5.5035i -0.7547 + 2.641

35、5i 0 </p><p>　　0 -0.6240 + 3.9002i 0 -0.7547 + 2.6415i 1.3787 -72.9583i 0 +63.4921i</p><p>　　0 0 0 0 0 +63.4921i

36、 0 -60.4686i</p><p><b>　　B1 =</b></p><p>　　1.0000 2.0000 0 + 0.0300i 0 1.0500 1.0000 1.0000 </p><p>　　2.0000 3.0000 0.0600 + 0.025

37、0i 0 + 0.5000i 1.0000 0 0 </p><p>　　2.0000 5.0000 0.0400 + 0.2500i 0 + 0.5000i 1.0000 0 0 </p><p>　　3.0000 4.0000 0.0800 + 0.3000i 0 + 0.5000i

38、 1.0000 0 0 </p><p>　　4.0000 5.0000 0.1000 + 0.3500i 0 1.0000 0 0 </p><p>　　6.0000 5.0000 0 + 0.0150i 0 1.0500 1.0000 1

39、.0000 </p><p><b>　　B2 =</b></p><p>　　0 0 1.2000 0 1.0000 </p><p>　　0 2.1000 + 1.0000i 1.0000

40、 0 2.0000 </p><p>　　0 1.8000 + 0.4000i 1.0000 0 2.0000 </p><p>　　0 1.6000 + 0.8000i 1.0000 0 2.000

41、0 </p><p>　　0 3.7000 + 1.3000i 1.0000 0 2.0000 </p><p>　　0 -5.0000 1.2000 0 3.0000 </p>&l

42、t;p>　　功率方程第(1)次差值：</p><p>　　0 0 4.2619 -2.1000 0.1000 -1.8000 -0.5500 -1.6000 8.4738 -3.7000 0 5.0000</p><p>　　形成的第(1)次Jacobi矩陣：</p><p>　　Co

43、lumns 1 through 12 </p><p>　　0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0</p><p>

44、　　0 0 -37.3887 14.8252 5.9172 -14.2012 0 0 3.9002 -0.6240 0 0</p><p>　　0 0 -14.8252 -47.9126 14.2012 5.9172 0 0 0.6240 3.9002 0 0</p><p>　　0 0 5.

45、9172 -14.2012 -8.0292 15.0311 3.1120 -0.8299 0 0 0 0 0 0 14.2012 5.9172 -15.0311 -9.0292 0.8299 3.1120 0 0 0 0</p><p>　　0 0 0 0 3.1120 -0.8299 -5.2535

46、 1.5846 2.6415 -0.7547 0 0 0 0 0 0 0.8299 3.1120 -1.5846 -5.7535 0.7547 2.6415 0 0 0 0 3.9002 -0.6240 0 0 2.6415 -0.7547 -63.1845 1.3787 63.4921 0 0 0 .6

47、240 3.9002 0 0 0.7547 2.6415 -1.3787 -82.7321 0 63.4921</p><p>　　0 0 0 0 0 0 0 0 0 0 -2.4000 0</p><p>　　0 0 0 0 0 0

48、 0 0 0 76.1905 0 -63.4921</p><p>　　Column 13 </p><p><b>　　0</b></p><p><b>　　0</b></p><p><b>　　4.2619</b><

49、/p><p><b>　　-2.1000</b></p><p><b>　　0.1000</b></p><p><b>　　-1.8000</b></p><p><b>　　-0.5500</b></p><p><b>

50、;　　-1.6000</b></p><p><b>　　8.4738</b></p><p><b>　　-3.7000</b></p><p><b>　　0</b></p><p><b>　　5.0000</b></p>

51、<p>　　Jacobi矩陣第(1)次消去運算</p><p>　　Columns 1 through 12 </p><p>　　0 0 0 0 0 0 0 0 0 0 0 0</p>&l

52、t;p>　　0 0 0 0 0 0 0 0 0 0 0 0</p><p>　　0 0 1.0000 -0.3965 -0.1583 0.3798 0 0 -0.

53、1043 0.0167 0 0</p><p>　　0 0 0 1.0000 -0.2204 -0.2147 0 0 0.0171 -0.0771 0 0</p><p>　　0 0 0

54、 0 1.0000 -1.0549 -0.3206 0.0855 -0.0845 0.1044 0 0</p><p>　　0 0 0 0 0 1.0000 0.1079 -0.1753 -0.0184 -0.0757 0 0<

55、/p><p>　　0 0 0 0 0 0 1.0000 -0.3868 -0.6525 0.1977 0 0</p><p>　　0 0 0 0 0 0 0 1.

56、0000 0.0419 -0.5520 0 0</p><p>　　0 0 0 0 0 0 0 0 1.0000 -0.0135 -1.0440 0</p><p>　　0 0

57、 0 0 0 0 0 0 0 1.0000 0.0106 -0.7856</p><p>　　0 0 0 0 0 0 0 0 0 0 1.0000

58、 0</p><p>　　0 0 0 0 0 0 0 0 0 0 0 1.0000</p><p>　　Column 13 </p><p><b>　　0</b><

59、/p><p><b>　　0</b></p><p><b>　　-0.1140</b></p><p><b>　　0.0705</b></p><p><b>　　-0.1659</b></p><p><b>　　0.1

60、184</b></p><p><b>　　0.0717</b></p><p><b>　　0.3118</b></p><p><b>　　-0.1439</b></p><p><b>　　0.0636</b></p>&l

61、t;p><b>　　0</b></p><p><b>　　-0.0431</b></p><p>　　Jacobi矩陣第(1)次回代運算</p><p>　　0 0 -0.1490 0.1103 -0.0024 0.1657 0.1015 0.3342 -0.14

62、35 0.0297 0 -0.0431</p><p><b>　　各個節(jié)點電壓模</b></p><p>　　0 1.1543 1.0160 0.9586 1.1439 1.2008</p><p>　　功率方程第(2)次差值：</p><p>　　0

63、 0 -0.9196 -0.3298 -0.2847 -0.0127 -0.5652 0.0384 -1.5185 0.1960 -0.0019 -0.3928</p><p>　　形成的第(2)次Jacobi矩陣。</p><p>　　Column 13 </p><p><b>　　0</b>

64、</p><p><b>　　0</b></p><p><b>　　-0.9196</b></p><p><b>　　-0.3298</b></p><p><b>　　-0.2847</b></p><p><b>

65、　　-0.0127</b></p><p><b>　　-0.5652</b></p><p><b>　　0.0384</b></p><p><b>　　-1.5185</b></p><p><b>　　0.1960</b></p&

66、gt;<p><b>　　-0.0019</b></p><p><b>　　-0.3928</b></p><p>　　Jacobi矩陣第(2)次消去運算。</p><p>　　Column 13 </p><p><b>　　0</b></p>

67、<p><b>　　0</b></p><p><b>　　0.0195</b></p><p><b>　　-0.0011</b></p><p><b>　　0.0425</b></p><p><b>　　-0.0053</

68、b></p><p><b>　　0.2265</b></p><p><b>　　-0.0488</b></p><p><b>　　0.0271</b></p><p><b>　　-0.0013</b></p><p>

69、<b>　　0.0008</b></p><p><b>　　0.0487</b></p><p>　　Jacobi矩陣第(2)次回代運算</p><p>　　0 0 0.0441 0.0155 0.0700 -0.0210 0.1790 -0.0275 0.0288

70、 0.0422 -0.0010 0.0487</p><p><b>　　各個節(jié)點電壓模</b></p><p>　　0 1.1120 0.9435 0.7821 1.1170 1.2010</p><p>　　功率方程第(3)次差值：</p><p>　　0

71、0 -0.0417 0.0017 -0.0185 -0.0126 -0.1340 -0.0217 -0.0367 0.0646 -0.0024 -0.0916</p><p>　　形成的第(3)次Jacobi矩陣。</p><p>　　Column 13 </p><p><b>　　0</b><

72、;/p><p><b>　　0</b></p><p><b>　　-0.0417</b></p><p><b>　　0.0017</b></p><p><b>　　-0.0185</b></p><p><b>　　-0

73、.0126</b></p><p><b>　　-0.1340</b></p><p><b>　　-0.0217</b></p><p><b>　　-0.0367</b></p><p><b>　　0.0646</b></p>

74、<p><b>　　-0.0024</b></p><p><b>　　-0.0916</b></p><p>　　Jacobi矩陣第(3)次消去運算。</p><p>　　Column 13 </p><p><b>　　0</b></p><

75、;p><b>　　0</b></p><p><b>　　0.0009</b></p><p><b>　　-0.0004</b></p><p><b>　　0.0034</b></p><p><b>　　-0.0004</b&g

76、t;</p><p><b>　　0.1046</b></p><p><b>　　-0.0114</b></p><p><b>　　0.0036</b></p><p><b>　　-0.0002</b></p><p><

77、;b>　　0.0010</b></p><p><b>　　0.0181</b></p><p>　　Jacobi矩陣第(3)次回代運算</p><p>　　0 0 0.0078 0.0043 0.0191 -0.0063 0.0768 -0.0055 0.0056

78、0.0156 0.0011 0.0181</p><p><b>　　各個節(jié)點電壓模</b></p><p>　　0 1.1048 0.9237 0.7098 1.1126 1.2001</p><p>　　功率方程第(4)次差值：</p><p>　　0 0

79、 -0.0004 0.0002 -0.0002 -0.0012 -0.0276 -0.0058 0.0006 0.0016 -0.0003 -0.0053</p><p>　　形成的第(4)次Jacobi矩陣</p><p>　　Column 13 </p><p><b>　　0</b></p&

80、gt;<p><b>　　0</b></p><p><b>　　-0.0004</b></p><p><b>　　0.0002</b></p><p><b>　　-0.0002</b></p><p><b>　　-0.001

81、2</b></p><p><b>　　-0.0276</b></p><p><b>　　-0.0058</b></p><p><b>　　0.0006</b></p><p><b>　　0.0016</b></p><

82、;p><b>　　-0.0003</b></p><p><b>　　-0.0053</b></p><p>　　Jacobi矩陣第(4)次消去運算</p><p>　　Column 13 </p><p><b>　　0</b></p><p>

83、<b>　　0</b></p><p><b>　　0.0000</b></p><p><b>　　-0.0000</b></p><p><b>　　0.0000</b></p><p><b>　　0.0000</b></

84、p><p><b>　　0.0450</b></p><p><b>　　-0.0033</b></p><p><b>　　0.0010</b></p><p><b>　　0.0002</b></p><p><b>　　

85、0.0001</b></p><p><b>　　0.0053</b></p><p>　　Jacobi矩陣第(4)次回代運算</p><p>　　0 0 0.0021 0.0014 0.0056 -0.0020 0.0259 -0.0018 0.0016 0.0048

86、 0.0002 0.0053</p><p><b>　　各個節(jié)點電壓模</b></p><p>　　0 1.1029 0.9178 0.6857 1.1114 1.2000</p><p>　　功率方程第(5)次差值：</p><p>　　0 0 -0.0000

87、 0.0000 0.0000 -0.0001 -0.0032 -0.0007 -0.0000 0.0001 -0.0000 -0.0005</p><p>　　形成的第(5)次Jacobi矩陣</p><p>　　Column 13 </p><p><b>　　0</b></p><

88、;p><b>　　0</b></p><p><b>　　-0.0000</b></p><p><b>　　0.0000</b></p><p><b>　　0.0000</b></p><p><b>　　-0.0001</b&g

89、t;</p><p><b>　　-0.0032</b></p><p><b>　　-0.0007</b></p><p><b>　　-0.0000</b></p><p><b>　　0.0001</b></p><p>&l

90、t;b>　　-0.0000</b></p><p><b>　　-0.0005</b></p><p>　　Jacobi矩陣第(5)次消去運算</p><p>　　Column 13 </p><p><b>　　0</b></p><p><b>

91、;　　0</b></p><p><b>　　0.0000</b></p><p><b>　　-0.0000</b></p><p><b>　　-0.0000</b></p><p><b>　　0.0000</b></p>

92、<p><b>　　0.0087</b></p><p><b>　　-0.0005</b></p><p><b>　　0.0002</b></p><p><b>　　0.0000</b></p><p><b>　　0.0000&

93、lt;/b></p><p><b>　　0.0008</b></p><p>　　Jacobi矩陣第(5)次回代運算</p><p>　　0 0 0.0003 0.0002 0.0008 -0.0003 0.0039 -0.0003 0.0002 0.0007 0.000

94、0 0.0008</p><p><b>　　各個節(jié)點電壓模</b></p><p>　　0 1.1026 0.9169 0.6820 1.1113 1.2000</p><p>　　功率方程第(6)次差值：</p><p>　　1.0e-004 *</p><p&

95、gt;　　0 0 -0.0026 0.0007 0.0032 -0.0182 -0.7419 -0.1647 -0.0016 0.0184 -0.0063 -0.1036</p><p>　　形成的第(6)次Jacobi矩陣</p><p>　　Column 13 </p><p><b>　　

96、0</b></p><p><b>　　0</b></p><p><b>　　-0.0000</b></p><p><b>　　0.0000</b></p><p><b>　　0.0000</b></p><p>

97、<b>　　-0.0000</b></p><p><b>　　-0.0001</b></p><p><b>　　-0.0000</b></p><p><b>　　-0.0000</b></p><p><b>　　0.0000</b&

98、gt;</p><p><b>　　-0.0000</b></p><p><b>　　-0.0000</b></p><p>　　Jacobi矩陣第(6)次消去運算</p><p>　　Column 13 </p><p><b>　　0</b><

99、;/p><p><b>　　0</b></p><p><b>　　0.0000</b></p><p><b>　　-0.0000</b></p><p><b>　　-0.0000</b></p><p><b>　　0.

100、0000</b></p><p><b>　　0.0002</b></p><p><b>　　-0.0000</b></p><p><b>　　0.0000</b></p><p><b>　　0.0000</b></p>&

101、lt;p><b>　　0.0000</b></p><p><b>　　0.0000</b></p><p>　　Jacobi矩陣第(6)次回代運算</p><p>　　1.0e-004 * 0 0 0.0764 0.0496 0.2048 -0.0728 0.9503

102、 -0.0641 0.0564 0.1753 0.0074 0.1923</p><p><b>　　各個節(jié)點電壓模</b></p><p>　　0 1.1026 0.9169 0.6819 1.1113 1.2000</p><p>　　功率方程第(7)次差值：</p>&l

103、t;p>　　1.0e-007 * 0 0 -0.0015 0.0004 0.0019 -0.0106 -0.4343 -0.0965 -0.0010 0.0108 -0.0037 -0.0606</p><p><b>　　迭代次數(shù)：</b></p><p><b>　　6</b&g

104、t;</p><p>　　沒有達到精度要求的個數(shù)：</p><p>　　9 10 10 10 8 3 0</p><p>　　各節(jié)點的電壓復(fù)數(shù)值E為(節(jié)點號從小到大排列)：</p><p>　　1.2000 1.0947 - 0.1316i 0.9068 - 0.1361i

105、 0.6128 - 0.2991i 1.1073 - 0.0932i 1.1996 - 0.0298i</p><p>　　各節(jié)點的電壓模值大小V為(節(jié)點號從小到大排列)：</p><p>　　1.2000 1.1026 0.9169 0.6819 1.1113 1.2000</p><p>　　各節(jié)點的電壓相角sida為(

106、節(jié)點號從小到大排列)：</p><p>　　0 -6.8549 -8.5347 -26.0187 -4.8097 -1.4241</p><p>　　各節(jié)點的功率S為(節(jié)點號從小到大排列)：</p><p>　　5.0135 + 1.8332i -2.1000 - 1.0000i -1.8000 - 0.4000i -1.6000 - 0.8

107、000i -3.7000 - 1.3000i 5.0000 + 2.5555i</p><p>　　----------------------------------------------------</p><p>　　各條支路的首端功率Si為(順序同您輸入B1時一致)：</p><p>　　S(1,2)=5.0135+1.8332i----------

108、-------------------------------------------</p><p>　　S(2,3)=3.0895+0.48941i----------------------------------------------------</p><p>　　S(2,5)=-0.176-0.31075i-----------------------------------

109、------------------</p><p>　　S(3,4)=0.78738+0.39432i-----------------------------------------------------</p><p>　　S(4,5)=-0.90639-0.43087i----------------------------------------------------<

110、/p><p>　　S(6,5)=5+2.5555i-----------------------------------------------------</p><p>　　各條支路的末端功率Sj為(順序同您輸入B1時一致)：</p><p>　　S(2,1)=-5.0135-1.1787i--------------------------------------

111、--------------</p><p>　　S(3,2)=-2.5874-0.79432i-----------------------------------------------------</p><p>　　S(5,2)=0.17703-0.29553i-----------------------------------------------------</p&g

112、t;<p>　　S(4,3)=-0.69361-0.36913i-----------------------------------------------------</p><p>　　S(5,4)=1.123+1.1889i-----------------------------------------------------</p><p>　　S(5,6)=-5

113、-2.1934i-----------------------------------------------------</p><p>　　各條支路的功率損耗DS為(順序同您輸入B1時一致)：</p><p>　　DS(1,2)=-8.8818e-016+0.65451i----------------------------------------------------</p

114、><p>　　DS(2,3)=0.50212-0.3049i-----------------------------------------------------</p><p>　　DS(2,5)=0.0010207-0.60628i-----------------------------------------------------</p><p>　　DS

115、(3,4)=0.093766+0.025185i-----------------------------------------------------</p><p>　　DS(4,5)=0.21659+0.75805i-----------------------------------------------------</p><p>　　DS(6,5)=0+0.36211i--

116、---------------------------------------------------</p><p><b>　　ZF =</b></p><p>　　DS(6,5)=0+0.36211i</p><p><b>　　第4章 學(xué)習(xí)心得</b></p><p>　　近一周的課程設(shè)計到

117、現(xiàn)在就要告一段落，通過課程設(shè)計，使我對這學(xué)期所學(xué)的知識又重新回顧了一遍。課程設(shè)計不同于以前的實驗，它是一個系統(tǒng)的工程，雖然以書本理論為核心，卻涉及了眾多實用性領(lǐng)域和經(jīng)驗化常識，待到真正要靠自己動手去查找資料、分析原理、繪制電路、選取器件乃至撰寫文檔時，才會感到設(shè)計的難度，不會的東西其實還很多很多，在以后的工作中還需要認(rèn)真的學(xué)習(xí)和實踐。時光飛逝，我們的學(xué)習(xí)到了最后一個環(huán)節(jié)，也是一個很重要的環(huán)節(jié)—畢業(yè)論文設(shè)計。因為我們可以通過畢業(yè)設(shè)計來進一

118、步綜合檢驗和鞏固自己學(xué)到的知識。</p><p><b>　　附錄A：</b></p><p><b>　　clear;</b></p><p>　　isb=1; %input('請輸入平衡母線節(jié)點號：isb=');</p><p>　　pr=0.00001; %in

119、put('請輸入誤差精度：pr=');%---------------------------------------------------</p><p>　　n=6;%input('請輸入節(jié)點數(shù)：n=');</p><p>　　nl=6;%input('請輸入支路數(shù)：nl=');</p><p>　　B1=[1

120、2 0+0.03i 0 1.05 1 1;</p><p>　　2 3 0.06+0.025i 0+0.5i 1 0 0;</p><p>　　2 5 0.04+0.25i 0+0.5i 1 0 0;</p><p>　　3 4 0.08+0.30i 0+0.5i

121、 1 0 0; </p><p>　　4 5 0.1+0.35i 0 1 0 0;</p><p>　　6 5 0+0.015i 0 1.05 1 1]</p><p>　　B2=[0 0 1.20 0 1;</p><

122、p>　　0 2.1+1i 1.00 0 2;</p><p>　　0 1.8+0.40i 1.00 0 2;</p><p>　　0 1.6+0.8i 1.00 0 2;</p><p>　　0 3.7+1.3i 1.00 0 2;</p><

123、p>　　0 -5+0i 1.20 0 3] </p><p>　　Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl);</p><p>　　for i=1:nl%從1到n1（總支路數(shù)）</p><p>　　if

124、 B1(i,7)==1 %-----------如果是變壓器支路--------</p><p>　　if B1(i,6)==0%左節(jié)點(首端)處于1側(cè)</p><p>　　p=B1(i,1);q=B1(i,2);</p><p>　　else %左節(jié)點(首端)處于K側(cè)</

125、p><p>　　p=B1(i,2);q=B1(i,1);</p><p><b>　　end</b></p><p>　　Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5));%非對角元</p><p>　　Y(q,p)=Y(p,q); %非對角元</p>

126、;<p>　　Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2);%對角元K側(cè)</p><p>　　Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4);%對角元1側(cè)+勵磁導(dǎo)納 </p><p>　　else %------------否則為線路支路--------------------</p><p