基于Matlab電力系統(tǒng)潮流計算

II目錄引言 11緒論 11.1電力系統(tǒng)敘述 11.2潮流計算簡介 21.3潮流計算的意義及其發(fā)展 22電力系統(tǒng)潮流計算基本原理 32.1電力網(wǎng)絡(luò)的基本方程式 32.1.2自導(dǎo)納和互導(dǎo)納的確定方法 42.1.3節(jié)點導(dǎo)納矩陣的性質(zhì)及意義 52.1.4非標(biāo)準(zhǔn)變比變壓器等值電路 72.2潮流計算的數(shù)學(xué)模型 82.2.1潮流計算的節(jié)點類型 82.2.2潮流計算基本方程 92.3潮流計算方法 102.3.1牛頓——拉夫遜法 102.3.2高斯——賽德爾法 102.3.3PQ分解法 123牛頓拉夫遜潮流計算理論分析 133.1概述 133.2牛頓法基本原理 133.3牛頓法潮流計算方程 173.3.1節(jié)點功率方程 173.3.2修正方程 184基于matlab的實例仿真 214.1潮流計算原始資料參數(shù) 214.2參數(shù)計算及等值電路的繪制 224.2.1節(jié)點設(shè)置及分類 224.2.3支路參數(shù)計算并求解 224.3求解方法 244.4牛頓–拉夫遜潮流計算法的求解過程 254.4.1牛頓–拉夫遜潮流計算法的計算框圖 254.4.2牛頓法計算潮流的步驟如下 254.4.3利用已知MATLAB程序求解，并修改相應(yīng)程序變量 264.4.4變電所負(fù)荷為題目所給數(shù)據(jù)進(jìn)行求解 264.4.5修改程序 274.5運行matlab程序輸出結(jié)果 284.6matlab仿真結(jié)果分析 295總結(jié) 31致謝 32參考文獻(xiàn) 32PAGE34基于Matlab電力系統(tǒng)潮流計算摘要:電力系統(tǒng)潮流計算是研究穩(wěn)態(tài)運情況一種重要分析計算，它根據(jù)給定的運行條件及系統(tǒng)接線情況確定整個電力系統(tǒng)各部分的運行狀態(tài)。在電力系統(tǒng)規(guī)劃設(shè)計和現(xiàn)有電力系統(tǒng)運行方式的研究中，都需要利用潮流計算來定量地分析比較供電方案或運行方式的合理性，可靠性和經(jīng)濟(jì)性。MATLAB使用方便，有著其他高級語言無法比擬的強(qiáng)大的矩陣處理功能。這樣使MATLAB成為電力系統(tǒng)潮流計算的首選計算機(jī)語言。牛頓-拉夫遜法是電力系統(tǒng)潮流計算的常用算法之一，它收斂性好，迭代次數(shù)少。介紹了電力系統(tǒng)潮流計算機(jī)輔分析的基本知識及潮流計算牛頓-拉普遜法,最后介紹了利用matlab進(jìn)行實例仿真。關(guān)鍵詞：電力系統(tǒng)潮流計算；牛頓-拉普遜法；matlab1緒論1.1電力系統(tǒng)敘述電力工業(yè)發(fā)展初期，電能是直接在用戶附近的發(fā)電站(或稱發(fā)電廠)中生產(chǎn)的，各發(fā)電站孤立運行。隨著工農(nóng)業(yè)生產(chǎn)和城市的發(fā)展，電能的需要量迅速增加，而熱能資源(如煤田)和水能資源豐富的地區(qū)又往往遠(yuǎn)離用電比較集中的城市和工礦區(qū)，為了解決這個矛盾，就需要在動力資源豐富的地區(qū)建立大型發(fā)電站，然后將電能遠(yuǎn)距離輸送給電力用戶。同時，為了提高供電可靠性以及資源利用的綜合經(jīng)濟(jì)性，又把許多分散的各種形式的發(fā)電站，通過送電線路和變電所聯(lián)系起來。這種由發(fā)電機(jī)、升壓和降壓變電所，送電線路以及用電設(shè)備有機(jī)連接起來的整體，即稱為電力系統(tǒng)。電力系統(tǒng)加上發(fā)電機(jī)的原動機(jī)(如汽輪機(jī)、水輪機(jī))，原動機(jī)的力能部分(如熱力鍋爐、水庫、原子能電站的反應(yīng)堆)、供熱和用熱設(shè)備，則稱為動力系統(tǒng)。現(xiàn)代電力系統(tǒng)提出了“靈活交流輸電與新型直流輸電”的概念。靈活交流輸電技術(shù)是指運用固態(tài)電子器件與現(xiàn)代自動控制技術(shù)對交流電網(wǎng)的電壓、相位角、阻抗、功率以及電路的通斷進(jìn)行實時閉環(huán)控制，從而提高高壓輸電線路的輸送能力和電力系統(tǒng)的穩(wěn)定水平。新型直流輸電技術(shù)是指應(yīng)用現(xiàn)電力電子技術(shù)的最新成果，改善和簡化變流站的造價等。1.2潮流計算簡介電力系統(tǒng)潮流計算是研究電力系統(tǒng)穩(wěn)態(tài)運行情況的一種計算，它根據(jù)給定的運行條件及系統(tǒng)接線情況確定整個電力系統(tǒng)各部分的運行狀態(tài)：各母線的電壓，各元件中流過的功率，系統(tǒng)的功率損耗等等。在電力系統(tǒng)規(guī)劃的設(shè)計和現(xiàn)有電力系統(tǒng)運行方式的研究中，都需要利用潮流計算來定量地分析比較供電方案或運行方式的合理性?？煽啃院徒?jīng)濟(jì)性。此外，電力系統(tǒng)潮流計算也是計算系統(tǒng)動態(tài)穩(wěn)定和靜態(tài)穩(wěn)定的基礎(chǔ)。所以潮流計算是研究電力系統(tǒng)的一種很重要和基礎(chǔ)的計算。電力系統(tǒng)潮流計算也分為離線計算和在線計算兩種，前者主要用于系統(tǒng)規(guī)劃設(shè)計和安排系統(tǒng)的運行方式，后者則用于正在運行系統(tǒng)的經(jīng)常監(jiān)視及實時控制。利用電子數(shù)字計算機(jī)進(jìn)行電力系統(tǒng)潮流計算從50年代中期就已經(jīng)開始。在這20年內(nèi)，潮流計算曾采用了各種不同的方法，這些方法的發(fā)展主要圍繞著對潮流計算的一些基本要求進(jìn)行的。對潮流計算的要求可以歸納為下面幾點：（1）計算方法的可靠性或收斂性；（2）對計算機(jī)內(nèi)存量的要求；（3）計算速度；（4）計算的方便性和靈活性。電力系統(tǒng)潮流計算問題在數(shù)學(xué)上是一組多元非線性方程式求解問題，其解法都離不開迭代。因此，對潮流計算方法，首先要求它能可靠地收斂，并給出正確答案。由于電力系統(tǒng)結(jié)構(gòu)及參數(shù)的一些特點，并且隨著電力系統(tǒng)不斷擴(kuò)大，潮流計算的方程式階數(shù)也越來越高，對這樣的方程式并不是任何數(shù)學(xué)方法都能保證給出正確答案的。這種情況成為促使電力系統(tǒng)計算人員不斷尋求新的更可靠方法的重要因素。1.3潮流計算的意義及其發(fā)展作為電力系統(tǒng)分析中最基本的計算，潮流計算能夠在復(fù)雜電力系統(tǒng)故障和正常條件下對穩(wěn)態(tài)運行狀態(tài)進(jìn)行分析計算。計算給定運行狀態(tài)是其目標(biāo)。潮流計算的基礎(chǔ)是對電力系統(tǒng)的規(guī)劃，對電力系統(tǒng)靜、暫態(tài)穩(wěn)定分析。潮流計算的結(jié)果可于穩(wěn)態(tài)研究，安全估計或著最優(yōu)潮流等。這些對其方法和模型有直接影響，牛-拉發(fā)是潮流計算主要方法。在用數(shù)字計算機(jī)解統(tǒng)潮流問題的開始階段，通常采取以節(jié)點導(dǎo)納矩陣為基礎(chǔ)的逐次代入法。此方法原理簡單，占用計算機(jī)內(nèi)存較小，適應(yīng)50年代電子計算機(jī)制造水平和當(dāng)時電力系統(tǒng)理論水平。它的缺點是收斂性比較差，系統(tǒng)規(guī)模擴(kuò)大，迭代次數(shù)就會驟然上升，計算過程常常會出現(xiàn)不收斂現(xiàn)象。為了解決這個困擾，我們運用阻抗法。它改善了計算過程中收斂性問題，解決了導(dǎo)納法無法求解的一些系統(tǒng)的潮流計算，在60年代應(yīng)用廣泛。占用計算機(jī)內(nèi)存大是它的主要不足，隨著迭代次數(shù)變大，系統(tǒng)擴(kuò)大時，這些缺點就更加明顯。60年代中期以阻抗矩陣為基礎(chǔ)的分塊阻抗法很好的克服了它的缺點。這中方法把大系統(tǒng)分割成幾個小系統(tǒng)，這樣大大減少了計算機(jī)內(nèi)存，同時也提高了計算速度。此外還有一種方法也可以很好的解決上面出現(xiàn)的各種問題。牛—拉發(fā)是導(dǎo)納矩陣為基礎(chǔ)，因此只要我們在計算過程中保持矩陣的稀疏性，就能提高這種方法的2電力系統(tǒng)潮流計算基本原理2.1電力網(wǎng)絡(luò)的基本方程式電力網(wǎng)絡(luò)可以用結(jié)點方程式或回路方程式表示出來。在結(jié)點方程式中表示網(wǎng)絡(luò)狀態(tài)的變量是各節(jié)點的電壓，在回路方程式中是各回路中的回路電流。一般若給出網(wǎng)絡(luò)的支路數(shù)b，結(jié)點數(shù)n，則回路方程式數(shù)m為m=b-n+1結(jié)點方程式數(shù)為=n-1因此，回路方程式數(shù)比結(jié)點方程式數(shù)多d=m-=b-2n+2在一般電力系統(tǒng)中，各結(jié)點(母線)和大地間有發(fā)電機(jī)、負(fù)荷、線路電容等對地支路，還有結(jié)點和結(jié)點之間也有輸電線路和變壓器之路，一般b>2n，用結(jié)點方程式表示比用回路方程式表示方程式數(shù)目要少。而且如以下所示，用結(jié)點方程式表示容易建立直觀的方程式，輸電線的連接狀態(tài)等變化時也很容易變更網(wǎng)絡(luò)方程式。把電力系統(tǒng)的發(fā)電機(jī)端子和負(fù)荷端子（同步調(diào)相機(jī)等的端子也作為發(fā)電機(jī)端來處理）抽出來，剩下的輸電線路及其它輸電系統(tǒng)概括為網(wǎng)絡(luò)Ｎet表示。在發(fā)電機(jī)結(jié)點和負(fù)荷結(jié)點上標(biāo)出任意順序的記號：1,2,…,I,…,n.在輸電系統(tǒng)Net的內(nèi)部不包含電源，并且各節(jié)點和大地間連接的線路對地電容、電力電容器等都作為負(fù)荷來處理。令端子1,2……,n的對地電壓分別為，由各端子流向輸電系統(tǒng)Net的電流相應(yīng)為，則此網(wǎng)絡(luò)方程組可以表示為(2-1)(2-1)式可以簡單寫成(I=1,2,…,n)(2-2)或者寫成I＝Y(jié)V(2-3)其中(2-4)(2-4)的Y稱為節(jié)點導(dǎo)納矩陣。因輸電系統(tǒng)Net只是由無源元件構(gòu)成的，而導(dǎo)納矩陣是對稱矩陣，于是有以下關(guān)系（2-5）電壓V和電流Ｉ的關(guān)系用式(2-1)～(2-5)表示時稱為節(jié)點導(dǎo)納方程式。這里電壓Ｖ用電流Ｉ的方程式表示時，則(2-3)式化為V＝ZI(2-6)其中(2-6)式稱為結(jié)點阻抗方程式，當(dāng)然，阻抗矩陣也是對稱矩陣。2.1.2自導(dǎo)納和互導(dǎo)納的確定方法電力網(wǎng)絡(luò)的節(jié)點電壓方程：(2-7)式(2-7)為節(jié)點注入電流列向量，注入電流有正有負(fù)，注入網(wǎng)絡(luò)的電流為正，流出網(wǎng)絡(luò)的電流為負(fù)。根據(jù)這一規(guī)定，電源節(jié)點的注入電流為正，負(fù)荷節(jié)點為負(fù)。既無電源又無負(fù)荷的聯(lián)絡(luò)節(jié)點為零，帶有地方負(fù)荷的電源節(jié)點為二者代數(shù)之和。式(2-7)為節(jié)點電壓列向量，由于節(jié)點電壓是對稱于參考節(jié)點而言的，因而需先選定參考節(jié)點。在電力系統(tǒng)中一般以地為參考節(jié)點。如整個網(wǎng)絡(luò)無接地支路，則需要選定某一節(jié)點為參考。設(shè)網(wǎng)絡(luò)中節(jié)點數(shù)為（不含參考節(jié)點），則，均為n*n列向量。為n*n階節(jié)點導(dǎo)納矩陣。節(jié)電導(dǎo)納矩陣的節(jié)點電壓方程：，展開為：(2-8)是一個n*n階節(jié)點導(dǎo)納矩陣，其階數(shù)就等于網(wǎng)絡(luò)中除參考節(jié)點外的節(jié)點數(shù)。節(jié)點導(dǎo)納矩陣的對角元素(i=1，2，n)成為自導(dǎo)納。自導(dǎo)納數(shù)值上就等于在i節(jié)點施加的電壓，其余節(jié)點接地時，經(jīng)過i節(jié)點注入的電流，因此，它可以定義為：（2-9）節(jié)點i的自導(dǎo)納數(shù)值上就等于與節(jié)點直接連接的所有支路導(dǎo)納的總和。節(jié)點導(dǎo)納矩陣的非對角元素(j=1，2，…，n;i=1，2，…。，n;j=i)稱互導(dǎo)納，由此可得互導(dǎo)納數(shù)值上就等于在節(jié)點i施加單位電壓，其他節(jié)點全部接地時，經(jīng)節(jié)點j注入網(wǎng)絡(luò)的電流，因此可定義為：(2-10)節(jié)點j，i之間的互導(dǎo)納數(shù)值上就等于連接節(jié)點j，i支路到導(dǎo)納的負(fù)值。顯然，恒等于?；?dǎo)納的這些性質(zhì)決定了節(jié)點導(dǎo)納矩陣是一個對稱稀疏矩陣。而且，由于每個節(jié)點所連接的支路數(shù)總有一個限度，隨著網(wǎng)絡(luò)中節(jié)點數(shù)的增加非零元素相對愈來愈少，節(jié)點導(dǎo)納矩陣的稀疏度，即零元素數(shù)與總元素的比值就愈來愈高。2.1.3節(jié)點導(dǎo)納矩陣的性質(zhì)及意義節(jié)點導(dǎo)納矩陣的性質(zhì)：（1）為對稱矩陣，=。如網(wǎng)絡(luò)中含有源元件，如移相變壓器，則對稱性不再成立。（2）對無接地支路的節(jié)點，其所在行列的元素之和均為零，即。對于有接地支路的節(jié)點，其所在行列的元素之和等于該點接地支路的導(dǎo)納。利用這一性質(zhì)，可以檢驗所形成節(jié)點導(dǎo)納矩陣的正確性。（3）具有強(qiáng)對角性：對角元素的值不小于同一行或同一列中任一元素。（4）為稀疏矩陣，因節(jié)點i，j之間無支路直接相連時=0，這種情況在實際電力系統(tǒng)中非常普遍。矩陣的稀疏性用稀疏度表示，其定義為矩陣中的零元素與全部元素之比，即，式中Z為中的零元素。S隨節(jié)點數(shù)n的增加而增加：n=50，S可達(dá)92%；n=100，S可達(dá)90%；n=500，S可達(dá)99%，充分利用節(jié)點導(dǎo)納矩陣的稀疏性可節(jié)省計算機(jī)內(nèi)存，加快計算速度，這種技巧稱為稀疏技術(shù)。節(jié)點導(dǎo)納矩陣的意義：是n*n階方陣，其對角元素(i=1，2，n)稱為自導(dǎo)納，非對角元素(i，j=1，2，n，)稱為互導(dǎo)納。將節(jié)點電壓方程展開為：（2-11）可見(2-12)表明，自導(dǎo)納在數(shù)值上等于僅在節(jié)點i施加單位電壓而其余節(jié)點電壓均為零（即其余節(jié)點全部接地）時，經(jīng)節(jié)點i注入網(wǎng)絡(luò)的電流。其顯然等于與節(jié)點i直接相連的所有支路的導(dǎo)納之和。同時可見。表明，互導(dǎo)納在數(shù)值上等于僅在節(jié)點j施加單位電壓而其余節(jié)點電壓均為零時，經(jīng)節(jié)點i注入網(wǎng)絡(luò)的電流，其顯然等于()即=。為支路的導(dǎo)納，負(fù)號表示該電流流出網(wǎng)絡(luò)。如節(jié)點ij之間無支路直接相連，則該電流為0，從而=0。注意字母幾種不寫法的不同意義：粗體黑字表示導(dǎo)納矩陣，大寫字母代矩陣中的第i行第j列元素，即節(jié)點i和節(jié)點j之間的互導(dǎo)納。小寫字母i，j支路的導(dǎo)納等于支路阻抗的倒數(shù)數(shù)，。2.1.4非標(biāo)準(zhǔn)變比變壓器等值電路變壓器型等值電路更便于計算機(jī)反復(fù)計算,更適宜于復(fù)雜網(wǎng)絡(luò)的潮流計算.雙繞組變壓器可用阻抗與一個理想變壓器串聯(lián)的電路表示.理想變壓器只是一個參數(shù),那就是變比。現(xiàn)在變壓器阻抗按實際變比歸算到低壓側(cè)為例,推導(dǎo)出變壓器型等值電路。圖2-2雙繞組變壓器原理圖圖2-3變壓器阻抗歸算到低壓側(cè)等值模型流入和流出理想變壓器的功率相等(2-13)式(2-13)中,是理想變壓器的變比,和分別為變壓器高,低繞組的實際電壓.從圖2-3直接可得：（2-14）從而可得:（2-15）式（2-14）中,又因節(jié)點電流方程應(yīng)具有如下形式:（2-16）將式（2-14）與（2-15）比較，得：，；，。因此可得各支路導(dǎo)納為：（2-17）由此可得用導(dǎo)納表示的變壓器型等值電路：圖2-4變壓器型等值電路2.2潮流計算的數(shù)學(xué)模型2.2.1潮流計算的節(jié)點類型根據(jù)電力系統(tǒng)中各節(jié)點性質(zhì)的不同，很自然地把節(jié)點分成三類：PQ節(jié)點對于這類節(jié)點，等值負(fù)荷功率PLi、QLi和等值電源功率PGi、QGi是給定的，從而注入功率Pi、Qi也是給定的，待求的則是節(jié)點電壓的大小Ui和相位角δi。這類節(jié)點稱為PQ節(jié)點PU節(jié)點對于這類節(jié)點，等值負(fù)荷和等值電源的有功功率PLi、QGi，從而注入有功功率Pi是給定。等值負(fù)荷的無功功率QLi和節(jié)點電壓Ui的大小事給定。待求的則是等值電源的無功功率QGi，從而注入無功功率Qi和節(jié)點電壓相位角δi。這類節(jié)點稱為PU節(jié)點平衡節(jié)點潮流計算時，一般都只設(shè)一個平衡節(jié)點。對這個節(jié)點，等值負(fù)荷功率PLS、QLS是給定的，節(jié)點電壓大小和相位角Us、δs.也是給定的。待求的則是等值電源功率PGS、QGS，從而注入功率Ps、Qs。擔(dān)負(fù)調(diào)整系統(tǒng)頻率任務(wù)的電廠母線往往被選作平衡節(jié)點。2.2.2潮流計算基本方程采用導(dǎo)納矩陣時，節(jié)點注入電流和節(jié)點電壓構(gòu)成如式(2-7)所示線性方程組可展開如下形式：(2-18)由于實際電網(wǎng)中測量的節(jié)點注入量一般不是電流而是功率，因此必須將式中的注入電流用節(jié)點注入功率來表示。節(jié)點功率與節(jié)點電流之間的關(guān)系為：(2-19)式中，因此用導(dǎo)納矩陣時，PQ節(jié)點可以表示為把這個關(guān)系代入式中，得（2-20）式（2-20）就是電力系統(tǒng)潮流計算的數(shù)學(xué)模型潮流方程。它具有如下特點：1：它是一組代數(shù)方程，因而表征的是電力系統(tǒng)的穩(wěn)定運行特性。2：它是一組非線性方程，因而只能用迭代方法求其數(shù)值解。3：由于方程中的電壓和導(dǎo)納既可以表為直角坐標(biāo)，又可表為極坐標(biāo)，因而潮流方程有多種表達(dá)形式極坐標(biāo)形式，直角坐標(biāo)形式和混合坐標(biāo)形式。（1）取，，得到潮流方程的極坐標(biāo)形式：(2-21)(2)取，，得到潮流方程的直角坐標(biāo)形式：(2-22)(3)取，得到潮流方程的混合坐標(biāo)形式：(2-23)不同坐標(biāo)形式的潮流方程適用于不同的迭代解法。例如：利用牛頓拉夫遜迭代法求解，以直角坐標(biāo)和混合坐標(biāo)形式的潮流方程為方便；而P-Q解耦法是在混合坐標(biāo)形式的基礎(chǔ)上發(fā)展而成，故當(dāng)然采用混合坐標(biāo)形式。4:它是一組n個復(fù)數(shù)方程，因而實數(shù)方程數(shù)為2n個但方程中共含4n個變量：P，Q，U和，i=1，2，，n，故必須先指定2n個變量才能求解。2.3潮流計算方法2.3.1牛頓——拉夫遜法牛頓法是數(shù)學(xué)中求解非線性方程式的典型方法，它是通過泰勒級數(shù)展開，忽略二階以上高階項，原理是逐次將非線性方程組線性，在多次形成和求解修正方程，直至滿足要求，具體的內(nèi)容參照第三章。2.3.2高斯——賽德爾法高斯-塞德爾法原理比較簡單，主要以節(jié)點導(dǎo)納矩陣為基礎(chǔ)。下面簡單介紹下其原理和潮流計算過程。（1）高斯-塞德爾法的基本原理設(shè)有n個聯(lián)立的非線性方程（2-24）解此方程組可得（2-25）若已經(jīng)求得各變量的第k此迭代值，則第（k+1）次迭代值為（2-26）只要給定變量的初值就可以按式（2-10）迭代計算，一直進(jìn)行到所有變量都滿足收斂條件：即可。（2）高斯-塞德爾潮流計算過程假設(shè)有n個節(jié)點的電力系統(tǒng)，沒有PV節(jié)點，平衡節(jié)點編號為s，功率方程可寫成下列復(fù)數(shù)方程式：（2-27）對每一個PQ節(jié)點都可列出一個方程式，因而有n-1個方程式。在這些方程式中，注入功率和都是給定的，平衡節(jié)點電壓也是已知的，因而只有n-1個節(jié)點的電壓為未知量，從而有可能求得唯一解。將上式寫成高斯-塞德爾法的迭代形式(2-28)如系統(tǒng)內(nèi)存在PV節(jié)點，假設(shè)節(jié)點p為PV節(jié)點，設(shè)定的節(jié)點電壓為Up0。假定高斯-塞德爾迭代法已完成第k次迭代，接著要做第k+1次迭代前，先按下式求出節(jié)點p的注入無功功率：(2-29)然后代入下式，求出p點電壓(2-30)在迭代過程中，按上式求得的節(jié)點p的電壓大小不一定等于設(shè)定的節(jié)點電壓Up0，所有在下一次的迭代中，應(yīng)以設(shè)定的Up0對電壓進(jìn)行修正，但其相角仍保持上式所求得的值，使得(2-31)如果所求得PV節(jié)點的無功功率越限，則無功功率在限，該P(yáng)V節(jié)點轉(zhuǎn)化為PQ節(jié)點。歸納起來，高斯-塞德爾迭代法計算潮流的步驟為：1.設(shè)定各節(jié)點電壓的初值，并給定迭代誤差判據(jù)；2.對每一個PQ節(jié)點，以前一次迭代的節(jié)點電壓值代入功率迭代方程式求出新值；3．對于PV節(jié)點，求出其無功功率，并判斷是否越限，如越限則將PV節(jié)點轉(zhuǎn)化為PQ節(jié)點；4.判別各節(jié)點電壓前后二次迭代值相量差的模是否小于給定誤差，如不小于，則回到第2步，繼續(xù)進(jìn)行計算，否則轉(zhuǎn)到第5步；5.根據(jù)功率方程求出平衡節(jié)點注入功率；6求支路功率分布和支路功率損耗。2.3.3PQ分解法PQ分解法是牛頓法的一種簡化方法，它利用了電力系統(tǒng)特有的運行特性，改進(jìn)和提高了運行速度。由牛頓法的修正方程進(jìn)行展開可得：(2-32)根據(jù)電力系統(tǒng)的運行特性進(jìn)行簡化：考到電力系統(tǒng)中有功功率分布主要受節(jié)點電壓相角的影響，無功功率分布主要受節(jié)點電壓幅值的影響，所以可以近似的忽略電壓幅值變化對有功功率和電壓相位變化對無功功率分布的影響，即(2-33)根據(jù)電力系統(tǒng)的正常運行條件作如假設(shè)：電壓相位角在電力系統(tǒng)運行時的變化不大（不超過10~20度）；通常架空線路的電抗遠(yuǎn)遠(yuǎn)大于電阻節(jié)點無功功率相應(yīng)的導(dǎo)納Q/U*U遠(yuǎn)小于該節(jié)點的自導(dǎo)納的虛部。用算式表示如下：（2-34）由上面假設(shè)可得雅克比矩陣的表達(dá)式：（2-35）修正式為(2-36)3牛頓拉夫遜潮流計算理論分析3.1概述牛頓法收斂性好，迭代次數(shù)少，在潮流計算方法中得到廣泛的應(yīng)用，目前為止還沒有更好的方法能夠完全取代它。牛頓拉夫遜法（下面簡稱牛頓法）是數(shù)學(xué)中求解非線性方程的典型方法，能快速求出其他方法求不出或者難以求出的解。本章將主要針對牛頓法的理論進(jìn)行具體介紹。3.2牛頓法基本原理牛頓-拉夫遜法是解非線性方程式的有效方法。牛頓拉夫遜法潮流計算是目前最為廣泛、效果最好的一種潮流計算方法。這種把非線性方程式的求解過程變成反復(fù)對相應(yīng)的線性方程式的求解過程，即逐次線性化過程，這就是牛頓法的核心。我們以如下非線性方程式的求解過程為例來說明：（3-1）設(shè)為該方程式的初值。而真正解x在它的近旁：（3-2）式中：為初始值的修正量。如果求得，則由式（3-2）就可以得到真正解x。為此將式（3-3）按泰勒級數(shù)展開（3-4）當(dāng)我們選擇的初始值比較好，即很小時，式（3-4）中包含的和更高階次項可以略去不計。因此，式（3-4）可以簡化為（3-5）這是對于變量的形式方程式，用它可以求出修正量。由于式（3-5）是式（3-4）的簡化結(jié)果，所以由式（3-5）解出后，還不能得到方程式（3-1）的真正解。實際上，用對修正后得到的：（3-6）現(xiàn)在如果再以作為初值，解式（3-5就能得到更趨近真正解的：（3-7）這樣反復(fù)下去，就構(gòu)成了不斷求解非線性方程式的逐次線性化過程。第t次迭代時的參數(shù)方程為（3-8）（3-9）上式左端可以看成是近似解引起的誤差，當(dāng)時，就滿足了原方程式（3-1），因而就成為該方程的解。式中是函數(shù)在點的一次導(dǎo)數(shù)，也就是曲線在點的斜率，如圖（3-1）所示，修正量則是由點的切線與橫軸的交點來確定，由圖（3-1）可以直觀的看出牛頓法的求解過程。圖3-1牛頓-拉夫遜法幾何解釋現(xiàn)在把牛頓法推廣到多變量非線性方程組的情況。設(shè)有變量的非線性聯(lián)立方程組：（3-10）給定各變量初值，假設(shè)為其修正量，并使其滿足（3-11）對以上n個方程式分別按泰勒級數(shù)展開，當(dāng)忽略所組成的二次項和高次項時，可以得到（3-12）式中：為函數(shù)對自變量的偏導(dǎo)數(shù)在點（）處的值。把上式寫成矩陣形式：（3-13）這是變量的線性方程組，稱為牛頓法的修正方程，通過它可以解出，并可以進(jìn)一步求得（3-14）式中向真正解逼近了一步，如果再以它們作為初值重復(fù)解式（3-13）修正方程式，等到更接近真解的，如此迭代下去，并按式（3-14）進(jìn)行修正，直到滿足收斂要求為止并停止迭代計算，這就構(gòu)成了牛頓法的迭代過程。一般第t次迭代式的修正方程為（3-15）上式可以簡寫為(3-16)其中，其中的為第t次迭代時的雅克比矩陣；同理可以得到第t次迭代時的修正量：（3-17）同樣，也可以寫出類似（3-14）的算式（3-18）這樣反復(fù)交替的解式（3-16）及式（3-18）就可以使逐步趨近方程式的真正解。當(dāng)滿足人為收斂條件時，即或（3-19）迭代結(jié)束，式中為預(yù)先給定的小正數(shù)。3.3牛頓法潮流計算方程3.3.1節(jié)點功率方程電力系統(tǒng)的負(fù)荷習(xí)慣用功率表示，對于有n個節(jié)點的電力系統(tǒng)，系統(tǒng)中各節(jié)點注入電流與注入功率以標(biāo)幺值表示的關(guān)系為i=1，2，……，n（3-20）式中表示其共軛復(fù)數(shù)。將此關(guān)系式代入節(jié)點電壓方程的通式，可得到以節(jié)點注入功率表示的節(jié)點電壓方程:(3-21)上述的方程式，通常稱為功率方程。根據(jù)方程中的節(jié)點電壓向量表示的不同，可以得到不同形式的功率方程。若節(jié)點電壓向量以直角坐標(biāo)表示，即以復(fù)數(shù)平面上實軸與虛軸上的投影表示可寫成（3-22）其共軛值為（3-23）導(dǎo)納表示為（3-24）把這兩關(guān)系式代回式（3-21）的功率方程中，展開后再將功率方程的實部和虛部分別寫成有功、無功功率分離的節(jié)點方功率方程：（3-25）式中：i=1，2，……，n為各節(jié)點的編號。若節(jié)點電壓以極坐標(biāo)表示，則或?qū)懗桑?-26）將其同導(dǎo)納的復(fù)數(shù)表達(dá)式一起代入式（3-21）的功率方程，進(jìn)整理可以得到（3-27）式中：——i與j節(jié)點電壓的相角差。由式（3-25）和（3-27）給出的功率方程表示方法避免了復(fù)數(shù)運算，因此，在潮流計算中普遍采用。3.3.2修正方程采用牛頓法計算潮流時，需要對功率方程進(jìn)行修改。下面將根據(jù)在不同坐標(biāo)內(nèi)的修改進(jìn)行討論：（1）在直角坐標(biāo)系內(nèi)時，由PQ節(jié)點功率方程（3-25）可知：節(jié)點i的注入功率是各點電壓的函數(shù)，設(shè)節(jié)點的電壓已知，代入式（3-25），可以求出節(jié)點i的有功及無功功率，它們與給定的PQ節(jié)點的注入功率的差值應(yīng)滿足以下方程（3-28）對于PV節(jié)點，已知節(jié)點的注入有功功率及節(jié)點電壓大小，記作，其節(jié)點的有功功率應(yīng)滿方程：（3-29）對于平衡節(jié)點，因為其電壓給定，故不需要迭代求解。通過以上分析可見，式（3-28）和式（3-29）共2（n-1）個方程，待求量共2（n-1）個。將上述2（n-1）個方程按泰勒級數(shù)展開，并略去修正量的高次方項后得到修正方程如下：（3-30）其中雅克比矩陣的各元素可以對式（3-28）和式（3-29）求偏導(dǎo)數(shù)獲得。對于非對角元素（）有（3-31）對于對角元素（有（3-32）由上述表達(dá)式可以看到，雅克比矩陣具有以下特點：各元素是各節(jié)點電壓的函數(shù)，迭代過程中每迭代一次各節(jié)點電壓都要變化，因而各元素每次也變化；雅克比矩陣不具有對稱性；互導(dǎo)納，與之對應(yīng)的非對角元素亦為零，此外因非對角元素，故雅克比矩陣是稀疏矩。當(dāng)在極坐標(biāo)系內(nèi)時，由功率方程（3-27）可知節(jié)點i的注入功率是各節(jié)點電壓幅值和相角的函數(shù)。代入式（3-27）可以求出節(jié)點i的有功功率和無功功率，它們與給定的PQ節(jié)點的注入功率的差值滿足下面方程：（3-33）式中：——i與j節(jié)點電壓的相角差。在有n個節(jié)點的系統(tǒng)中，假定第號節(jié)點為PQ節(jié)點，第m+1~n-1號節(jié)點為PV節(jié)點，第n號節(jié)點為平衡節(jié)點。和是給定的，PV節(jié)點的電壓幅值也是給定的，因此，只剩下n-1個節(jié)點的電壓相角和m個節(jié)點的電壓幅值是未知量。由（3-33）可知一共包含了n-1+m方程式，正好同未知量的數(shù)目相等，而直角坐標(biāo)形式的方程少了n-1-m個。由方程（3-33）可以寫出修正方程(3-34)式中（3-35）其中：H是階方陣，其元素為；N是階矩陣，其元素為；K是階矩陣，其元素為；L是階矩陣，其元素為。對式（3-33）求偏導(dǎo)數(shù)，可得雅克比矩陣元素的表達(dá)式如下：非對角元素（）（3-36）對角元素（）（3-37）4基于matlab的實例仿真4.1潮流計算原始資料參數(shù)系統(tǒng)地理接線如圖1所示線長為：1－火廠：70km 2－火廠：120km4－火廠：80km 2－3：100km 3－水廠：60km 3－火廠：200km均為雙回線，電壓為220KV，單回線路參數(shù)為：r1=0.085Ω/km x1=0.32Ω/km b1=3.5*10-6s/km圖1系統(tǒng)地理接線圖四個變電所的負(fù)荷cosФ=0.85,年初負(fù)荷分配如下：①70MW②90MW ③60MW ④80MW10KV35KV 10KV 35KV變壓器配置：P=60MW時，配二臺50MVA變壓器；每臺PK=254KW，uk%=14.5 P=70MW時，配二臺63MVA變壓器；每臺PK=245KW，uk%=12.5 P=80MW時，配二臺63MVA變壓器；每臺PK=245KW，uk%=12.5 P=90MW時，配二臺75MVA變壓器；每臺PK=414KW，uk%=16.7要求：①、變電所的負(fù)荷為年初數(shù)值，年末最大運行方式數(shù)值應(yīng)比年初數(shù)值高10％年中最大運行方式數(shù)值應(yīng)比年初數(shù)值低3％各時期最小運行方式數(shù)值應(yīng)比當(dāng)時的最大運行方式數(shù)值低30％。②、火電廠作為平衡節(jié)點，水電廠作為PV節(jié)點。③、各節(jié)點電壓應(yīng)在UN－1.05UN之間。④、水廠夏季滿發(fā)P=120MW，－40<Q<75Mvar，冬季最大調(diào)峰功率P=60MW，－40<Q<75Mvar。=5\*GB3⑤、夏季最小運行方式。系統(tǒng)的網(wǎng)絡(luò)接線圖如下圖所示：圖2系統(tǒng)網(wǎng)絡(luò)接線圖4.2參數(shù)計算及等值電路的繪制4.2.1節(jié)點設(shè)置及分類根據(jù)電力系統(tǒng)圖將火電廠母線設(shè)為節(jié)點節(jié)點1，水電廠節(jié)點設(shè)為節(jié)點8，將變電所1、2、3、4的高壓側(cè)分別設(shè)為節(jié)點2、5、6、9、低壓側(cè)為3、4、7、10。并且，將節(jié)點1設(shè)為平衡節(jié)點，將節(jié)點8設(shè)為PV節(jié)點，其余節(jié)點設(shè)為PQ節(jié)點，即：節(jié)點分類標(biāo)號平衡節(jié)點：1（電壓初始設(shè)置為231V）PV節(jié)點：2（電壓初始設(shè)置為231V，由于總負(fù)荷為200MW電廠1和電廠2裝機(jī)容量比約為3:2，所以設(shè)置電廠2的發(fā)電功率為120MW）PQ節(jié)點：2、3、4、5、6、7、9、10（電壓設(shè)置為220V）。電壓基準(zhǔn)值為220V，視在功率基準(zhǔn)值為100WVA。4.2.3支路參數(shù)計算并求解設(shè)定電壓基準(zhǔn)值為220V，視在功率基準(zhǔn)值為100WVA，根據(jù)題目原始資料，計算發(fā)電廠、變壓器及線路的參數(shù)有名值。參數(shù)計算公式如下：1.變壓器阻抗：2.變壓器電抗：3.線路電阻：4.電抗：5.線路電納：B=b·L6.變電所負(fù)荷分別為：變電所1=47.53+j29.456變電所2=61.11+j37.873變電所3=40.74+j25.248變電所4=54.32+j33.665將參數(shù)整理，雙回路支路阻抗除以2，對地電納乘以2，見下表：表1各支路等值參數(shù)首端號 末端號支路的阻抗（R+jx）/支路的對地電納B/122.975+j11.2j0.00049231.494+j48.01601510.2+j38.4j0.00042541.718+j53.8850568.5+j32J0.00035682.55+j9.6j0.00042672.459+j70.1801617+j64J0.0007193.4+j12.8j0.000569101.494+j48.01607.計算變壓器分接頭：0.95-1.05變壓器有5個抽頭，電壓調(diào)節(jié)范圍為2*2.5%，對應(yīng)的分接頭開始時設(shè)變壓器高壓側(cè)接主接頭,降壓變壓器5個分接頭時的非標(biāo)準(zhǔn)變比以備調(diào)壓時選用。4.3求解方法畫出系統(tǒng)的等效電路圖，在計算出各元件參數(shù)的基礎(chǔ)上，應(yīng)用牛頓—拉夫遜Newton-Raphson法以及MATLAB軟件進(jìn)行計算對給定系統(tǒng)圖進(jìn)行潮流計算，經(jīng)過調(diào)節(jié)均得到符合電壓限制及功率限制的潮流分布。采用標(biāo)么值進(jìn)行計算,根據(jù)題中最高電壓等級為220KV,所以選擇基準(zhǔn)電壓及功率分別為：據(jù)此可以計算各部分標(biāo)幺值。繪制等效電路時，雙回線阻抗減半，容抗加倍；并聯(lián)變壓器阻抗變?yōu)橐话?，無容抗。繪制的等效電路如下：圖3等效電路圖4.4牛頓–拉夫遜潮流計算法的求解過程4.4.1牛頓–拉夫遜潮流計算法的計算框圖根據(jù)給出的程序畫出的潮流計算框圖如下圖：輸入原始數(shù)據(jù)輸入原始數(shù)據(jù)形成節(jié)點導(dǎo)納陣給定電壓初值e（0）、f（0）k=0根據(jù)公式計算用公式計算雅克比矩陣各元素高斯法解修正方程，求修正節(jié)點電壓是計算平衡節(jié)點功率及全部線路功率輸出k=k+1否圖4程序流程圖4.4.2牛頓法計算潮流的步驟如下⑴給這各節(jié)點電壓初始值；⑵將以上電壓初始值代入公式，求修正方程的常數(shù)項向量；⑶將電壓初始值在帶入公式，求出修正方程中系數(shù)矩陣的各元素。⑷解修正方程式；⑸修正各節(jié)點電壓，；⑹將，在帶入方程式，求出；⑺檢驗是否收斂，即⑻如果收斂，迭代到此結(jié)束，進(jìn)一步計算各線路潮流和平衡節(jié)點功率，并打印輸出結(jié)果。如果不收斂，轉(zhuǎn)回（2）進(jìn)行下次迭代計算，直到收斂為止。4.4.3利用已知MATLAB程序求解，并修改相應(yīng)程序變量需要輸入的數(shù)據(jù)（1）n—節(jié)點數(shù)、n1—線路數(shù)、isb—平衡母線節(jié)點號、pr—誤差精度，如果要輸入則輸入eps即可。（2）請輸入由支路參數(shù)形成的矩陣B1。矩陣B1的每行是由如下元素構(gòu)成的：某支路的首端號；某支路的末端號；支路的阻抗；支路的對地容抗；支路的變比K；支路首端處于K側(cè)為1，1側(cè)為0。（3）請輸入各節(jié)點參數(shù)形成的矩陣B2。矩陣B2的每行是由下列參數(shù)構(gòu)成的：該節(jié)點所接發(fā)電機(jī)的功率；該節(jié)點的負(fù)荷功率；節(jié)點電壓的初始值；PV節(jié)點電壓V的給定值；節(jié)點所接的無功補(bǔ)償裝置的容量；節(jié)點分類標(biāo)號節(jié)點分類標(biāo)號4.4.4變電所負(fù)荷為題目所給數(shù)據(jù)進(jìn)行求解1.形成B1、B2矩陣：根據(jù)所求參數(shù)，以及B1矩陣的含義，列寫B(tài)1矩陣如下：B1=為了平衡兩發(fā)電廠發(fā)出的電量，令發(fā)水電廠的功率為120MW，為了減小線路上的損耗，令發(fā)電機(jī)的電壓為額定電壓的1.05倍。并且，根據(jù)前面敘述的節(jié)點分類，形成B2矩陣如下：B2=4.4.5修改程序修改程序：由于變壓器變比發(fā)生改變，兩側(cè)標(biāo)幺值同時發(fā)生變化，就要對標(biāo)幺值進(jìn)行折算：在程序中添加程序段如下：Vv=V./[11K10K911K811K7];plot(Vv);………[2202201138.52202201122022038.5].*VvVVV=[2202201138.52202201122022038.5].*V該程序段先把高壓側(cè)標(biāo)么值換算為低壓側(cè)標(biāo)么值，再把所有電壓的標(biāo)么值轉(zhuǎn)化為有名值，輸出到變量表中。4.5運行matlab程序輸出結(jié)果表1未調(diào)整時變壓器分接頭均選為1參數(shù)節(jié)點電壓節(jié)點功率支路支路功率123186.0775+15.5794i1-295.3985+39.57222i2229.4166-0-0i2-395.1443+61.9441i311.126-47.53-29.456i1-513.9470+38.6380i437.9755-61.11-37.873i5-4122.4048+81.338i5227.3954-0-0i5-6-109.2134-10.0545i6230.1655-0-0i6-8-239+23.8450i711.0681-40.74-25.248i6-781.5953+53.7859i8231120-21.6392i1-6-46.9966+18.3199i9228.88840-0i1-9109.1168+46.6099i1038.6594-54.32-33.665i9-10108.7650+71.3483i調(diào)整合適后的結(jié)果如下：表2變壓器分接頭均選為0.975參數(shù)節(jié)點電壓節(jié)點功率支路支路功率123186.0652+15.3868i1-295.3985+39.57222i2229.4166-0-0i2-395.1443+61.9441i311.126-47.53-29.456i1-513.9470+38.6380i440.0571-61.11-37.873i5-4122.4048+81.338i5227.42430-0i5-6-109.2134-10.0545i6230.1715-0-0i6-8-239+23.8450i711.0684-40.74-25.248i6-781.5953+53.7859i8231120-21.7828i1-6-46.9966+18.3199i9228.88840-0i1-9109.1168+46.6099i1038.6594-54.32-33.665i9-10108.7650+71.3483i當(dāng)分接頭為0.95時，輸出結(jié)果如下：表3變壓器分接頭均選為0.95參數(shù)節(jié)點電壓節(jié)點功率支路支路功率123186.0534+15.2013i1-295.3985+39.57222i2229.4166-0-0i2-395.1443+61.9441i311.126-47.53-29.456i1-513.9470+38.6380i442.3041-61.11-37.873i5-4122.4048+81.338i5227.4522-0-0i5-6-109.2134-10.0545i6230.17730-0i6-8-239+23.8450i711.0687-40.74-25.248i6-781.5953+53.7859i8231120-21.9211i1-6-46.9966+18.3199i9228.88840-0i1-9109.1168+46.6099i1038.6594-54.32-33.665i9-10108.7650+71.3483i4.6matlab仿真結(jié)果分析為比較電壓的變化，制作電壓表格如下：表4電壓比較表電壓節(jié)點號未調(diào)節(jié)調(diào)節(jié)適當(dāng)后過調(diào)節(jié)12312312312229.4166229.4166229.4166311.12611.12611.126437.975540.057142.30415227.3954227.4243227.45226230.1655230.1715230.1773711.068111.068411.068782312312319228.8884228.8884228.88841038.659438.659438.6594由表可知未調(diào)節(jié)時，電壓超出允許范圍，調(diào)節(jié)適當(dāng)后，電壓在適當(dāng)范圍內(nèi)，過調(diào)節(jié)時電壓有超出給定范圍。支路損耗對比表如下：表5支路損耗比較表支路損耗未調(diào)節(jié)調(diào)節(jié)適當(dāng)后過調(diào)節(jié)1-20.1499-25.4039i0.1499-25.4039i0.1499-25.4039i2-30.0943+3.0321i0.0943+3.0321i0.0943-3.0321i1-50.0845-21.7464i0.0832-21.7542i0.0819-21.7616i5-40.1955+5.9146i0.1848+5.592i0.1745+5.2811i5-60.489-16.479i0.4886-16.4832i0.4882-16.4872i6-80.6933-19.7206i0.6935-19.7206i0.6936-19.7206i6-70.1153+3.2901i0.1153+3.2899i0.1153+3.2897i1-60.2039-36.4503i0.2038-36.4517i0.2037-36.453i1-90.2267-28.7566i0.2267-28.7566i0.2267-28.7566i9-100.125+4.0183i0.125+4.0183i0.125+4.0183i總和2.3775-132.3i2.3652-132.64i2.3534-132.96i對比以上參數(shù)可知，分接頭越低，對應(yīng)節(jié)點電壓越高，同時損耗越小表1形成的圖像如下圖5變壓器分接頭未調(diào)節(jié)為1時仿真圖形從圖5可以直觀的看出4號節(jié)點電壓不在1-1.05范圍之內(nèi)，所以節(jié)點4不滿足題目要求。表2形成的圖像如下圖6變壓器分接頭調(diào)節(jié)為0.975時仿真圖形從圖中可以看出各節(jié)點電壓都在UN－1.05UN之間，因此滿足題目要求。表3形成的圖像如下圖7變壓器分接頭調(diào)節(jié)為0.95時仿真圖形綜上所述在分接頭為0.95和1時均超出允許范圍，在分接頭為0.975時滿足電壓在1-1.05UN電壓在允許的范圍內(nèi)，滿足題目要求的條件。另外,在調(diào)節(jié)電壓時，如果單獨調(diào)節(jié)分接頭不能輸出滿足條件時，可以同時調(diào)節(jié)發(fā)電機(jī)的電壓，以使輸出參數(shù)滿足要求的條件。5總結(jié)潮流計算是電力系統(tǒng)的最基本、最常用的計算。根據(jù)系統(tǒng)給定的運行條件、網(wǎng)絡(luò)接線及元件參數(shù)，通過潮流計算可以確定各母線的電壓（相角及幅值）、各元件中流過的功率、整個系統(tǒng)的功率損耗等。潮流計算是實現(xiàn)電力系統(tǒng)安全經(jīng)濟(jì)發(fā)展的必要手段和重要環(huán)節(jié)。牛頓—拉夫遜法是一種經(jīng)典的數(shù)學(xué)方法，用來解決矩陣的計算問題，在數(shù)值分析課程中我們也接觸過這種方法，這次在電力系統(tǒng)的潮流計算中得到了應(yīng)用，更讓我覺得知識靈活運用的重要性。MATLAB是一套高性能的數(shù)學(xué)計算軟件，它集數(shù)值分析、矩陣計算、信號處理和圖形顯示于一身，構(gòu)成了一個方便的界面和友好的用戶環(huán)境，其強(qiáng)大的擴(kuò)展功能為各個領(lǐng)域的應(yīng)用提供了方便。在潮流計算的大量數(shù)值計算過程中更顯示出其優(yōu)勢，是我們應(yīng)該掌握的一門基本技術(shù)。在這次課程設(shè)計中，我不僅學(xué)會了matlab的使用方法，還學(xué)會了很多的實用的matlab指令。在這次畢業(yè)設(shè)計中我還將平時學(xué)的調(diào)壓方法運用在電力系統(tǒng)的調(diào)節(jié)，真正做到了學(xué)以致用。實際的情況遠(yuǎn)比我們計算的情況復(fù)雜，這讓我深刻了解了潮流計算的重要性。精準(zhǔn)的潮流計算不僅可以節(jié)省能源并充分發(fā)揮能源的作用，還可以使電網(wǎng)處于穩(wěn)定且平衡的狀態(tài)，這對電網(wǎng)的安全運行起到關(guān)鍵性的作用。所以我認(rèn)為學(xué)好電力系統(tǒng)的關(guān)鍵在于學(xué)好潮流計算。致謝畢業(yè)設(shè)計是培養(yǎng)學(xué)生綜合運用所學(xué)知識，發(fā)現(xiàn)，提出，分析和解決實際問題,鍛煉實踐能力的重要環(huán)節(jié),是對學(xué)生實際工作能力的具體訓(xùn)練和考察過程。回顧此次畢業(yè)設(shè)計，從選題到定稿，從理論到實踐，在近幾個月的日子里，雖然工作辛苦，但是我掌握了很多新的東西，這都是在四年的理論學(xué)習(xí)中沒有接觸過的，不僅鞏固了以前所學(xué)過的知識，而且拓展了本專業(yè)的知識面。通過這次畢業(yè)設(shè)計我懂得了理論與實際相結(jié)合的重要性，只有把所學(xué)的理論知識與實踐相結(jié)合起來，從實驗中得出結(jié)論，驗證理論，應(yīng)用理論，才能提高自己的實際動手能力和獨立思考的能力。在設(shè)計的過程中不斷遇到各種問題，正是這些問題使我發(fā)現(xiàn)了自己的不足之處，解決每一個問題的過程使我收獲頗豐，也督促我還要繼續(xù)深入學(xué)習(xí)。這次畢業(yè)設(shè)計的順利完成離不開劉玉娟老師的悉心指導(dǎo)，劉老師在設(shè)計方向以及很多細(xì)節(jié)方面都給予了我很大幫助，劉老師對待工作一絲不茍的嚴(yán)謹(jǐn)工作作風(fēng)深深影響了我，特在此表示感謝！同時，對給過我?guī)椭乃型瑢W(xué)和各位指導(dǎo)老師再次表示忠心的感謝！參考文獻(xiàn)[1]于永源,楊綺雯，電力系統(tǒng)分析（第三版）［M].中國電力出版社，2007[2]張勇軍，任震等，基于災(zāi)變遺傳算法的無功規(guī)劃優(yōu)化，〔J〕電力系統(tǒng)自動化，2002年第23期。[3]王守相，劉玉田電力系統(tǒng)潮流計算研究現(xiàn)狀--《山東電力技術(shù)》1996年05期[4]何仰贊，溫增銀.電力系統(tǒng)分析（上冊）第三版［M].湖北：華中科技大學(xué)出版社，2002[5]劉同娟.MATLAB在電路分析中的應(yīng)用.電氣電子教學(xué)學(xué)報.2002[6]西安交通大學(xué)等.電力系統(tǒng)計算[M].北