復雜網絡NR法潮流分析與計算的設計畢業(yè)設計

1、摘 要電力系統(tǒng)穩(wěn)態(tài)分析包括潮流計算和靜態(tài)安全分析。本文主要運用的事潮流計算，潮流計算是電力網絡設計與運行中最基本的運算，對電力網絡的各種設計方案及各種運行方式進行潮流計算，可以得到各種電網各節(jié)點的電壓，并求得網絡的潮流及網絡中的各元件的電力損耗，進而求得電能損耗。本位就是運用潮流計算具體分析，并有matlab仿真。關鍵詞： 電力系統(tǒng) 潮流計算 matlab仿真abstractelectric power system steady flow calculation and analysis of the static safety analysis. this paper, by means

2、of the calculation, flow calculation is the trend of the power network design and operation of the most basic operations of electric power network, various design scheme and the operation ways to tide computation, can get all kinds of each node of the power grid voltage and seek the trend of the net

3、work and the network of the components of the power loss, and getting electric power. the standard is to use the power flow calculation and analysis, the specific have matlab simulation.key words: power system; flow calculation; matlab simulation目錄摘 要1abstract1目錄2第一章 電力系統(tǒng)的基本概念2第二章 潮流計算32.1潮流計算概述與發(fā)展3

4、2.2復雜電力系統(tǒng)潮流計算42.3潮流計算的方法及優(yōu)、缺點42.4 matlab概述52.6牛頓-拉夫遜法原理52.7牛頓-拉夫遜法解決潮流計算問題72.8計算機潮流計算的步驟82.9計算機程序的實現9程序運行結果：11第四章 學習心得20附錄：21第一章 電力系統(tǒng)的基本概念電力系統(tǒng)：發(fā)電機把機械能轉化為電能，電能經變壓器和電力線路輸送并分配到用戶，在那里經電動機、電爐和電燈等設備又將電能轉化為機械能、熱能和光能等。這些生產、變換、輸送、分配、消費電能的發(fā)電機、變壓器、變換器、電力線路及各種用電設備等聯(lián)系在一起組成的統(tǒng)一整體稱為電力系統(tǒng)。電力網：電力系統(tǒng)中除發(fā)電機和用電設備外的部分。動力系統(tǒng)：

5、電力系統(tǒng)和“動力部分”的總和。電力系統(tǒng)潮流計算是研究電力系統(tǒng)穩(wěn)態(tài)運行情況的一種計算，它根據給定的運行條件及系統(tǒng)接線情況確定整個電力系統(tǒng)各部分的運行狀態(tài)：各母線的電壓，各元件中流過的功率，系統(tǒng)的功率損耗等等。在電力系統(tǒng)規(guī)劃的設計和現有電力系統(tǒng)運行方式的研究中，都需要利用潮流計算來定量地分析比較供電方案或運行方式的合理性?？煽啃院徒洕?。此外，電力系統(tǒng)潮流計算也是計算系統(tǒng)動態(tài)穩(wěn)定和靜態(tài)穩(wěn)定的基礎。所以潮流計算是研究電力系統(tǒng)的一種很重要和很基礎的計算。隨著科學技術的發(fā)展，電力系統(tǒng)變得越來越復雜，電氣工程師掌握一種好的能對電力系統(tǒng)進行仿真的軟件是學習和研究的需要。文章簡要介紹了matlab發(fā)展歷史、組

6、成和強大的功能，并用簡單例子分別就編程和仿真兩方面分析了matiab軟件在電力系統(tǒng)研究中的具體應用。采取等效電路法，能對特殊、復雜地電力系統(tǒng)進行高效仿真研究，因此，掌握編程和仿真是學好matlab的基礎。與眾多專門的電力系統(tǒng)仿真軟件相比，matlab軟件具有易學、功能強大和開放性好，是電力系統(tǒng)仿真研究的有力工具。第二章 潮流計算2.1潮流計算概述與發(fā)展電力系統(tǒng)潮流計算也分為離線計算和在線計算兩種，前者主要用于系統(tǒng)規(guī)劃設計和安排系統(tǒng)的運行方式，后者則用于正在運行系統(tǒng)的經常監(jiān)視及實時控制。利用電子數字計算機進行電力系統(tǒng)潮流計算從50年代中期就已經開始。在這20年內，潮流計算曾采用了各種不同的方法，

7、這些方法的發(fā)展主要圍繞著對潮流計算的一些基本要求進行的。牛頓-拉夫遜法作為一種實用的，有競爭力的電力系統(tǒng)潮流計算方法，是在應用了稀疏矩陣技巧和高斯消去法求修正方程后。牛頓-拉夫遜法是求解非線性代數方程有效的迭代計算。2.2復雜電力系統(tǒng)潮流計算電力系統(tǒng)潮流計算是對復雜電力系統(tǒng)正常和故障條件下穩(wěn)態(tài)運行狀態(tài)的計算。潮流計算的目標是求取電力系統(tǒng)在給定運行方式下的節(jié)點電壓和功率分布，用以檢查系統(tǒng)各元件是否過負荷、各點電壓是否滿足要求、功率的分布和分配是否合理以及功率損耗等。對現有電力系統(tǒng)的運行和擴建，對新的電力系統(tǒng)進行規(guī)劃設計都是以潮流計算為基礎。潮流計算結果的用途，例如用于電力系統(tǒng)穩(wěn)定研究、安全估計或

8、最優(yōu)潮流等也對潮流計算的模型和方法有直接影響。節(jié)點類型：1）pv節(jié)點：柱入有功功率p為給定值，電壓也保持在給定數值。2）pq節(jié)點：諸如有功功率和無功功率是給定的。3）平衡節(jié)點：用來平衡全電網的功率。選一容量足夠大的發(fā)電機擔任平衡全電網功率的職責。平衡節(jié)點的電壓大小與相位是給定的，通常以它的相角為參考量，即取其電壓相角為0。一個獨立的電力網中只設一個平衡點?；静襟E：1）形成節(jié)點導納矩陣；2）將各節(jié)點電壓設初值u；3）將節(jié)點初值代入相關求式，求出修正方程式的常數項向量；4）將節(jié)點電壓初值代入求式，求出雅可比矩陣元素；5）求解修正方程，求修正向量；6）求取節(jié)點電壓的新值；7）檢查是否收斂，如不收斂

9、，則以各節(jié)點電壓的新值作為初值自第3步重新開始進行狹義次迭代，否則轉入下一步；8）計算支路功率分布，pv節(jié)點無功功率和平衡節(jié)點柱入功率。2.3潮流計算的方法及優(yōu)、缺點潮流計算法有，簡化梯度法、二次規(guī)劃法、牛頓-拉夫遜法(newton-raphson)等。簡化梯度法是采用梯度法進行搜索，用罰函數處理違約的不等式約束。該方法程序編制簡便，所需存儲量小，對初始點無特殊要求，曾獲得普遍重視，成為第一種有效的優(yōu)化潮流方法。簡化梯度法的缺點：迭代過程中，尤其是在接近最優(yōu)點附近會出現鋸齒現象，收斂性較差，收斂速度很慢；每次迭代都要重新計算潮流，計算量很大，耗時較多。二次規(guī)劃法是二階的方法，解決最優(yōu)潮流問題收

10、斂精度較好，能很好地解決耦合的最優(yōu)潮流問題，但缺點是計算lagrange函數的二階偏導數，計算量大、計算復雜。2.4 matlab概述目前電子計算機已廣泛應用于電力系統(tǒng)的分析計算，潮流計算是其基本應用之一?，F有很多潮流計算方法。對潮流計算方法有五方面的要求：（1）計算速度快；（2）內存需要少；（3）計算結果有良好的可靠性和可信性；（4）適應性好，亦即能處理變壓器變比調整、系統(tǒng)元件的不同描述和與其它程序配合的能力強；（5）簡單。matlab是一種交互式、面向對象的程序設計語言，廣泛應用于工業(yè)界與學術界，主要用于矩陣運算，同時在數值分析、自動控制模擬、數字信號處理、動態(tài)分析、繪圖等方面也具有強大的

11、功能。matlab程序設計語言結構完整，且具有優(yōu)良的移植性，它的基本數據元素是不需要定義的數組。它可以高效率地解決工業(yè)計算問題，特別是關于矩陣和矢量的計算。matlab與c語言和fortran語言相比更容易被掌握。通過m語言，可以用類似數學公式的方式來編寫算法，大大降低了程序所需的難度并節(jié)省了時間，從而可把主要的精力集中在算法的構思而不是編程上。2.6牛頓-拉夫遜法原理假設有n個聯(lián)立的非線性代數方程：假設以給出各變量的初值，令其分別為個變量的修正量，使?jié)M足以上方程，所以：將上式中的n個多元函數在初始值附近分別展開成泰勒級數，并略去含有，的二次及以上階次的各項，便得：方程可寫成：以上方程是對于修

12、正量，的線性方程組，稱為牛頓法的修正方程，可解出，。對初始近似解進行修正： （i=1,2,，n）反復迭代，在進行k+1次迭代時，從求解修正方程式：得到修正量,對各量進行修正 （i=1,2,，n）迭代過程一直進行到滿足收斂判據2.7牛頓-拉夫遜法解決潮流計算問題節(jié)點總數為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é)點?？尚纬山Y點導納矩陣。導納矩陣元素可表示為，本文中節(jié)點電壓以直角坐標形式表示，即。由此下列公式可求出pi，qi假設系統(tǒng)中的第1，2，m號

13、節(jié)點為pq節(jié)點，第i個節(jié)點的給定功率為和，對該節(jié)點可列方程：假設系統(tǒng)中的第m+1，m+2，n-1號節(jié)點為pv節(jié)點，則對其中每一個節(jié)點可列方程: 第n號節(jié)點為平衡節(jié)點，其電壓為是給定的，故不參加迭代。修正方程可寫成分塊矩陣的形式：通過反復求解修正方程，解出各節(jié)點的未知量，再通過收斂判據判定是否已為真值。從而求得pq節(jié)點的電壓v及相角的真值，pv節(jié)點的q、真值，平衡節(jié)點的p、q真值，以上即為牛頓-拉夫遜迭代法的潮流計算過程，其優(yōu)點為計算精確，運行速度快。其中的各個環(huán)節(jié)都可通過matlab程序來實現。2.8計算機潮流計算的步驟（1）對電力網絡的所有參數設初值,包括電壓、相角、有功、無功等。（2）處理

14、非標準變比支路,使其變成標準變比為1的變壓器支路。（3）形成節(jié)點導納矩陣y。（4）計算有功功率的不平衡量pi,從而求出。（5）根據節(jié)點的類型形成j。（6）解修正方程式,求各節(jié)點的電壓的變化量ei(i=1,2,3.n,is)（7）求各節(jié)點相角的新值ei=ei+ei (i=1,2,3. n,is)（8）計算無功功率的不平衡量qi,從而求出 (i=1,2,3.n,is)（9）解修正方程式,求各節(jié)點的電壓大小的變化量 (i=1,2,3.,n,is)。（10）求各節(jié)點的電壓大小的新值 (i=1,2,3.,n,is)。（11）運用個節(jié)點的電壓的新值自第四步開始下一次迭代。計算平衡節(jié)點的功率和線路功率。其中

15、平衡節(jié)點的功率的計算公式為線路上的功率為:從而線路上的損耗的功率為:2.9計算機程序的實現導納矩陣的形成節(jié)點導納矩陣是方陣,其階數等于網絡中出參考節(jié)點外的節(jié)點數n。節(jié)點導納矩陣是稀疏矩陣,其各行非零非對角元數就等于該行相對應節(jié)點所連接的不接地支路數。節(jié)點導納矩陣的對角元就等于各該節(jié)點所連接導納的總和。節(jié)點導納矩陣的非對角元等于連接節(jié)點i,j支路導納的負值。點導納矩陣一般是對稱矩陣。對于支路中有非標準變比變壓器的支路來說,利用下面的公式來計算它的導納。變壓器的變比,假如已知非標準變比支路i,j上的阻抗(以下沒有特殊說明所有的參數都用標幺值)為,則線路導納為,線路上的對地半導納為。j的形成y是由最

16、終形成的導納矩陣的虛部組成的,但是pv節(jié)點以及平衡節(jié)點不參加qv迭代,因此y中不包含與這些節(jié)點有關的元素。迭代條件和約束方程迭代條件就是如果q時就停止迭代。對節(jié)點的約束條件分為三類:即對節(jié)點注入功率的約束、對節(jié)點電壓大小的約束和對相角的約束。其中對節(jié)點注入功率的約束,主要是對電源注入功率的約束條件不能滿足時,將威脅到發(fā)電機的安全運行。對電壓大小的約束不能滿足時,將影響電能的質量,嚴重時將影響系統(tǒng)運行的穩(wěn)定性。對相對相角的約束條件不能滿足時,也將危及系統(tǒng)運行的穩(wěn)定性。圖2.2 程序流程圖3.0 手算潮流計算用圖1和圖2的數據和等值網絡形成節(jié)點導納矩陣1. 節(jié)點導納矩陣程序運行結果：導納矩陣 y=

17、0 -30.2343i 0 +31.7460i 0 0 0 0 0 +31.7460i 14.8252 -42.6506i -14.2012 + 5.9172i 0 -0.6240 + 3.9002i 0 0 -14.2012+5.9172i 15.0311-8.5292i -0.8299+3.1120i 0 0 0 0 -0.8299 + 3.1120i 1.5846 - 5.5035i -0.7547 + 2.6415i 0 0 -0.6240 + 3.9002i 0 -0.7547 + 2.6415i 1.3787 -72.9583i 0 +63.4921i0 0 0 0 0 +63.4

18、921i 0 -60.4686ib1 =1.0000 2.0000 0 + 0.0300i 0 1.0500 1.0000 1.0000 2.0000 3.0000 0.0600 + 0.0250i 0 + 0.5000i 1.0000 0 0 2.0000 5.0000 0.0400 + 0.2500i 0 + 0.5000i 1.0000 0 0 3.0000 4.0000 0.0800 + 0.3000i 0 + 0.5000i 1.0000 0 0 4.0000 5.0000 0.1000 + 0.3500i 0 1.0000 0 0 6.0000 5.0000 0 + 0.0150i

19、 0 1.0500 1.0000 1.0000 b2 =0 0 1.2000 0 1.0000 0 2.1000 + 1.0000i 1.0000 0 2.0000 0 1.8000 + 0.4000i 1.0000 0 2.0000 0 1.6000 + 0.8000i 1.0000 0 2.0000 0 3.7000 + 1.3000i 1.0000 0 2.0000 0 -5.0000 1.2000 0 3.0000 功率方程第(1)次差值： 0 0 4.2619 -2.1000 0.1000 -1.8000 -0.5500 -1.6000 8.4738 -3.7000 0 5.0000

20、形成的第(1)次jacobi矩陣： columns 1 through 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 00 0 -37.3887 14.8252 5.9172 -14.2012 0 0 3.9002 -0.6240 0 00 0 -14.8252 -47.9126 14.2012 5.9172 0 0 0.6240 3.9002 0 00 0 5.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.

21、8299 3.1120 0 0 0 00 0 0 0 3.1120 -0.8299 -5.2535 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 .6240 3.9002 0 0 0.7547 2.6415 -1.3787 -82.7321 0 63.49210 0 0 0 0 0 0 0 0 0 -2.4000 00 0 0 0 0 0 0

22、0 0 76.1905 0 -63.4921 column 13 0 0 4.2619 -2.1000 0.1000 -1.8000 -0.5500 -1.6000 8.4738 -3.7000 0 5.0000jacobi矩陣第(1)次消去運算 columns 1 through 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1.0000 -0.3965 -0.1583 0.3798 0 0 -0.1043 0.0167 0 0 0 0 0 1.0000 -0.2204 -0.2147 0 0 0.0171 -0.0771 0

23、0 0 0 0 0 1.0000 -1.0549 -0.3206 0.0855 -0.0845 0.1044 0 0 0 0 0 0 0 1.0000 0.1079 -0.1753 -0.0184 -0.0757 0 0 0 0 0 0 0 0 1.0000 -0.3868 -0.6525 0.1977 0 0 0 0 0 0 0 0 0 1.0000 0.0419 -0.5520 0 0 0 0 0 0 0 0 0 0 1.0000 -0.0135 -1.0440 0 0 0 0 0 0 0 0 0 0 1.0000 0.0106 -0.7856 0 0 0 0 0 0 0 0 0 0 1.

24、0000 0 0 0 0 0 0 0 0 0 0 0 0 1.0000 column 13 0 0 -0.1140 0.0705 -0.1659 0.1184 0.0717 0.3118 -0.1439 0.0636 0 -0.0431jacobi矩陣第(1)次回代運算 0 0 -0.1490 0.1103 -0.0024 0.1657 0.1015 0.3342 -0.1435 0.0297 0 -0.0431各個節(jié)點電壓模 0 1.1543 1.0160 0.9586 1.1439 1.2008功率方程第(2)次差值： 0 0 -0.9196 -0.3298 -0.2847 -0.0127

25、 -0.5652 0.0384 -1.5185 0.1960 -0.0019 -0.3928形成的第(2)次jacobi矩陣。 column 13 0 0 -0.9196 -0.3298 -0.2847 -0.0127 -0.5652 0.0384 -1.5185 0.1960 -0.0019 -0.3928jacobi矩陣第(2)次消去運算。 column 13 0 0 0.0195 -0.0011 0.0425 -0.0053 0.2265 -0.0488 0.0271 -0.0013 0.0008 0.0487jacobi矩陣第(2)次回代運算 0 0 0.0441 0.0155 0.0

26、700 -0.0210 0.1790 -0.0275 0.0288 0.0422 -0.0010 0.0487各個節(jié)點電壓模 0 1.1120 0.9435 0.7821 1.1170 1.2010功率方程第(3)次差值： 0 0 -0.0417 0.0017 -0.0185 -0.0126 -0.1340 -0.0217 -0.0367 0.0646 -0.0024 -0.0916形成的第(3)次jacobi矩陣。 column 13 0 0 -0.0417 0.0017 -0.0185 -0.0126 -0.1340 -0.0217 -0.0367 0.0646 -0.0024 -0.09

27、16jacobi矩陣第(3)次消去運算。 column 13 0 0 0.0009 -0.0004 0.0034 -0.0004 0.1046 -0.0114 0.0036 -0.0002 0.0010 0.0181jacobi矩陣第(3)次回代運算 0 0 0.0078 0.0043 0.0191 -0.0063 0.0768 -0.0055 0.0056 0.0156 0.0011 0.0181各個節(jié)點電壓模 0 1.1048 0.9237 0.7098 1.1126 1.2001功率方程第(4)次差值： 0 0 -0.0004 0.0002 -0.0002 -0.0012 -0.0276

28、 -0.0058 0.0006 0.0016 -0.0003 -0.0053形成的第(4)次jacobi矩陣。 column 13 0 0 -0.0004 0.0002 -0.0002 -0.0012 -0.0276 -0.0058 0.0006 0.0016 -0.0003 -0.0053jacobi矩陣第(4)次消去運算。 column 13 0 0 0.0000 -0.0000 0.0000 0.0000 0.0450 -0.0033 0.0010 0.0002 0.0001 0.0053jacobi矩陣第(4)次回代運算 0 0 0.0021 0.0014 0.0056 -0.0020

29、 0.0259 -0.0018 0.0016 0.0048 0.0002 0.0053各個節(jié)點電壓模 0 1.1029 0.9178 0.6857 1.1114 1.2000功率方程第(5)次差值： 0 0 -0.0000 0.0000 0.0000 -0.0001 -0.0032 -0.0007 -0.0000 0.0001 -0.0000 -0.0005形成的第(5)次jacobi矩陣。 column 13 0 0 -0.0000 0.0000 0.0000 -0.0001 -0.0032 -0.0007 -0.0000 0.0001 -0.0000 -0.0005jacobi矩陣第(5)

30、次消去運算。 column 13 0 0 0.0000 -0.0000 -0.0000 0.0000 0.0087 -0.0005 0.0002 0.0000 0.0000 0.0008jacobi矩陣第(5)次回代運算 0 0 0.0003 0.0002 0.0008 -0.0003 0.0039 -0.0003 0.0002 0.0007 0.0000 0.0008各個節(jié)點電壓模 0 1.1026 0.9169 0.6820 1.1113 1.2000功率方程第(6)次差值： 1.0e-004 * 0 0 -0.0026 0.0007 0.0032 -0.0182 -0.7419 -0.1

31、647 -0.0016 0.0184 -0.0063 -0.1036形成的第(6)次jacobi矩陣。 column 13 0 0 -0.0000 0.0000 0.0000 -0.0000 -0.0001 -0.0000 -0.0000 0.0000 -0.0000 -0.0000jacobi矩陣第(6)次消去運算。 column 13 0 0 0.0000 -0.0000 -0.0000 0.0000 0.0002 -0.0000 0.0000 0.0000 0.0000 0.0000jacobi矩陣第(6)次回代運算 1.0e-004 * 0 0 0.0764 0.0496 0.2048

32、 -0.0728 0.9503 -0.0641 0.0564 0.1753 0.0074 0.1923各個節(jié)點電壓模 0 1.1026 0.9169 0.6819 1.1113 1.2000功率方程第(7)次差值： 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迭代次數： 6沒有達到精度要求的個數： 9 10 10 10 8 3 0各節(jié)點的電壓復數值e為(節(jié)點號從小到大排列)： 1.2000 1.0947 - 0.1316i 0.9068 - 0.1361i

33、0.6128 - 0.2991i 1.1073 - 0.0932i 1.1996 - 0.0298i各節(jié)點的電壓模值大小v為(節(jié)點號從小到大排列)： 1.2000 1.1026 0.9169 0.6819 1.1113 1.2000各節(jié)點的電壓相角sida為(節(jié)點號從小到大排列)： 0 -6.8549 -8.5347 -26.0187 -4.8097 -1.4241各節(jié)點的功率s為(節(jié)點號從小到大排列)： 5.0135 + 1.8332i -2.1000 - 1.0000i -1.8000 - 0.4000i -1.6000 - 0.8000i -3.7000 - 1.3000i 5.0000

34、 + 2.5555i-各條支路的首端功率si為(順序同您輸入b1時一致)：s(1,2)=5.0135+1.8332i-s(2,3)=3.0895+0.48941i-s(2,5)=-0.176-0.31075i-s(3,4)=0.78738+0.39432i-s(4,5)=-0.90639-0.43087i-s(6,5)=5+2.5555i-各條支路的末端功率sj為(順序同您輸入b1時一致)：s(2,1)=-5.0135-1.1787i-s(3,2)=-2.5874-0.79432i-s(5,2)=0.17703-0.29553i-s(4,3)=-0.69361-0.36913i-s(5,4)=

35、1.123+1.1889i-s(5,6)=-5-2.1934i-各條支路的功率損耗ds為(順序同您輸入b1時一致)：ds(1,2)=-8.8818e-016+0.65451i-ds(2,3)=0.50212-0.3049i-ds(2,5)=0.0010207-0.60628i-ds(3,4)=0.093766+0.025185i-ds(4,5)=0.21659+0.75805i-ds(6,5)=0+0.36211i-zf = ds(6,5)=0+0.36211i第四章 學習心得近一周的課程設計到現在就要告一段落，通過課程設計，使我對這學期所學的知識又重新回顧了一遍。課程設計不同于以前的實驗，它

36、是一個系統(tǒng)的工程，雖然以書本理論為核心，卻涉及了眾多實用性領域和經驗化常識，待到真正要靠自己動手去查找資料、分析原理、繪制電路、選取器件乃至撰寫文檔時，才會感到設計的難度，不會的東西其實還很多很多，在以后的工作中還需要認真的學習和實踐。時光飛逝，我們的學習到了最后一個環(huán)節(jié)，也是一個很重要的環(huán)節(jié)畢業(yè)論文設計。因為我們可以通過畢業(yè)設計來進一步綜合檢驗和鞏固自己學到的知識。我們本次的課程設計時間應該是比較緊迫的。所以我們也做了比較仔細和充分的準備。先是到學校圖書館和學校的網上數字圖書館進到了各種資料，在按照老師給我們的技術要求的基礎上然后再分析整理加以組織，這樣構成了自己課題的主體部份。在這過程中我

37、態(tài)度端正，積極上進，踏實認真。認為這是一次自己學習的好機會 。同時也得到老師親自指導，這是很寶貴的。附錄：clear;isb=1; %input(請輸入平衡母線節(jié)點號：isb=);pr=0.00001; %input(請輸入誤差精度：pr=);%-n=6;%input(請輸入節(jié)點數：n=);nl=6;%input(請輸入支路數：nl=);b1=1 2 0+0.03i 0 1.05 1 1; 2 3 0.06+0.025i 0+0.5i 1 0 0; 2 5 0.04+0.25i 0+0.5i 1 0 0; 3 4 0.08+0.30i 0+0.5i 1 0 0; 4 5 0.1+0.35i 0

38、 1 0 0; 6 5 0+0.015i 0 1.05 1 1b2=0 0 1.20 0 1; 0 2.1+1i 1.00 0 2; 0 1.8+0.40i 1.00 0 2; 0 1.6+0.8i 1.00 0 2; 0 3.7+1.3i 1.00 0 2; 0 -5+0i 1.20 0 3 y=zeros(n);e=zeros(1,n);f=zeros(1,n);v=zeros(1,n);sida=zeros(1,n);s1=zeros(nl);for i=1:nl%從1到n1（總支路數）if b1(i,7)=1 %-如果是變壓器支路- if b1(i,6)=0%左節(jié)點(首端)處于1側 p

39、=b1(i,1);q=b1(i,2); else %左節(jié)點(首端)處于k側 p=b1(i,2);q=b1(i,1); end y(p,q)=y(p,q)-1./(b1(i,3)*b1(i,5);%非對角元 y(q,p)=y(p,q); %非對角元 y(q,q)=y(q,q)+1./(b1(i,3)*b1(i,5)2);%對角元k側 y(p,p)=y(p,p)+1./b1(i,3)+b1(i,4);%對角元1側+勵磁導納 else %-否則為線路支路-p=b1(i,1);q=b1(i,2);y(p,q)=y(p,q)-1./b1(i,3);%非對角元 y(q,p)=y(p,q); %非對角元 y(q,q)=y(q,q)+1./b1(i,3)+b1(i,4)./2.0000;%對角元j側+線路電納的一半 y(p,p)=y(p,p)+1./b1(i,3)+b1(i,4)./2.0000;%對角元i側+線路電納的一半endenddisp(導納矩陣 y=);disp(y);%-給定各節(jié)點初始電壓及給定各節(jié)點注入功率-g=real(y);b=imag(y);%分解出導納陣的實部和虛部 for i=1:n%給定各節(jié)點初始電壓的實部和虛部 e(i)=real(b2(i,3); f(i)=imag(b2(i,3); v(i)