1、摘 要潮流計算是電力系統的一項重要分析功能，是進行故障計算，繼電保護整定，安全分析的必要工具。是對復雜電力系統正常和故障條件下穩態運行狀態的計算。電力系統潮流計算的結果是電力系統穩定計算和故障分析的基礎。潮流計算是電力系統的各種計算的基礎，同時它又是研究電力系統的一項重要分析功能，在電力系統規劃設計和現有電力系統運行方式的研究中，都需要利用電力系統潮流計算來定量的比較供電方案和現有電力系統運行方式的合理性，可靠性和經濟性。實際電力系統的潮流技術那主要采用牛頓-拉夫遜法。傳統的潮流計算程序缺乏圖形用戶界面，結果顯示不直觀，難于與其他分析功能集成。網絡原始數據輸入工作量大且易于出錯。隨著計算機技術

2、的飛速發展，MICROSOFT WINDOWS操作系統早已被大家所熟悉，其友好的圖形用戶界面已成為PC機的標準，而DOS操作系統下的應用程序因其界面不夠友好，開發具有WINDOWS風格界面的電力系統分析軟件已成為當前的主流趨勢。另外，傳統的程序設計方法是結構化程序設計方法，該方法基于功能分解，把整個軟件工程看作是一個個對象的組合，由于對某個特定問題域來說，該對象組成基本不變，因此，這種基于對象分解方法設計的軟件結構上比較穩定，易于維護和擴充。本文介紹了圖形化潮流計算軟件的開發設計思想和總體結構，闡述了該軟件所具備的功能和特點。結合電力系統的特點，軟件采用 MATLAB語言運行于WINDOWS操

3、作系統的圖形化潮流計算軟件。本系統的主要特點是操作簡單，圖形界面直觀，運行穩定。計算準確。計算中，算法做了一些改進，提高了計算速度，各個類的有效封裝又使程序具有很好的模塊性。可維護性和可重用性。關鍵詞：電力系統潮流計算；牛頓拉夫遜法潮流計算； MATLABAbstractPower flow calculation is an important function of the electric power system is analyzed, fault calculation, relay protection, safety analysis tools necessary. Is t

4、he calculation of the steady state operation of normal and fault conditions of complex power system under the. The results of power flow calculation is the basis of computation and fault analysis of power system stability. Power flow calculation is the basis of all kinds of power system calculation,

5、 and it is also an important function of power system analysis, the research on the operating mode of power system planning and design of the existing power system, need to use power flow calculation to quantitative comparison of the rationality of power supply scheme and operation mode of the exist

6、ing power system, reliability and economy. The actual power flow technology that mainly uses the Newton-Raphson method.The traditional flow calculation program lacks a graphical user interface, the display is not intuitive, it is difficult to integrate with other analysis function. The network input

7、 data and heavy workload and error prone. With the rapid development of computer technology, Microsoft windows operating system has long been familiar, its friendly graphical user interface has become PC standard, and the application of the DOS operating system because of its interface is not friend

8、ly enough, the power system analysis software development with windows style interface has become the main trend. In addition, the traditional design method is a structured program design method, this method is based on the function decomposition, the software engineering as a combination of objects

9、, due to a particular problem domain, the composition of the object is essentially the same, therefore, based on the software structure of the object decomposition method on the design of stable, easy to maintain and expansion.The development of design idea and overall structure of this paper introd

10、uces the graphical power flow calculation software, expounds the function and features of the software. According to the characteristics of the power system, the graphical trend MATLAB language operating system running on a windows calculation software. The main feature of this system is simple, int

11、uitive graphical interface, accurate calculation of stable operation. In the calculation, the algorithm made some improvements, improve the calculation speed, and the effective encapsulation of class program module has good maintainability and reusability.Keywords: power flow calculation; Newton Rap

12、hson power flow calculation; MATLAB目錄引言1第1章電力系統潮流計算概述21.1 電力系統敘述21.2 潮流計算簡介21.3 潮流計算的意義及其發展3第2章潮流計算的數學模型42.1 導納矩陣的原理及計算方法42.2 潮流計算的基本方程72.3 電力系統節點分類92.4 潮流計算的約束條件10第3章牛頓拉夫遜法概述123.1 牛頓拉夫遜法基本原理123.2 牛頓-拉夫遜法潮流求解過程133.3 牛頓拉夫遜法的程序框圖17第4章 MATLAB概述184.1 MATLAB簡介184.2 MATLAB應用在潮流計算中的優勢184.3 矩陣的運算19第5章 潮流計算主

13、界面設計實現215.1 主界面介紹215.2 數據初始化215.3 潮流計算225.4 數據處理235.5 數據傳遞的問題235.6 例：某電網接線圖及給定的參數245.7運算結果24第6章 系統潮流計算的前沿算法及發展前景316.1 保留非線性算法316.2 最優潮流分析法316.3 OPF分析法32結論與展望33致謝34參考文獻35附錄36附錄A:基于MATLAB的牛頓拉夫遜法潮流計算程序清單36附錄B:外文文獻及譯文39附錄C:參考文獻的題錄摘要45插圖清單圖2-1雙繞組變壓器原理圖6圖2-2變壓器等值電路7圖2-3潮流計算用的電網結構圖8圖2-4潮流計算等值網絡8圖3-1牛頓拉夫遜法的

14、程序框圖17圖5-1某電網接線圖24引言潮流計算是在給定電力系統網絡結構、參數和決定系統運行狀態的邊界條件的情況下確定系統穩態運行狀態的一種基本方法，是電力系統規劃和運營中不可缺少的一個重要組成部分。可以說，它是電力系統分析中最基本、最重要的計算，是系統安全、經濟分析和實時控制與調度的基礎。是電力系統研究人員長期研究的一個課題。MATLAB自1980年問世以來，它的強大的矩陣處理功能給電力系統的分析、計算帶來許多方便。在處理潮流計算時，其計算機軟件的速度已無法滿足大電網模擬和實時控制的仿真要求，而高效的潮流問題相關軟件的研究已成為大規模電力系統仿真計算的關鍵。隨著計算機技術的不斷發展和成熟，對

15、MATLAB潮流計算的研究為快速、詳細地解決大電網的計算問題開辟了新思路。電力系統潮流計算是電力系統分析中的一種最基本的計算，是對復雜電力系統正常和故障條件下穩態運行狀態的計算。潮流計算的目標是求取電力系統在給定運行狀態的計算。即節點電壓和功率分布，用以檢查系統各元件是否過負荷。各點電壓是否滿足要求，功率的分布和分配是否合理以及功率損耗等。對現有電力系統的運行和擴建，對新的電力系統進行規劃設計以及對電力系統進行靜態和暫態穩定分析都是以潮流計算為基礎。潮流計算結果可用如電力系統穩態研究，安全估計或最優潮流等對潮流計算的模型和方法有直接影響。在用數字計算機解電力系統潮流問題的開始階段，普遍采取以節

16、點導納矩陣為基礎的逐次代入法。這個方法的原理比較簡單，要求的數字計算機內存量比較下，適應50年代電子計算機制造水平和當時電力系統理論水平。但它的收斂性較差，當系統規模變大時，迭代次數急劇上升，在計算中往往出現迭代不收斂的情況。這就迫使電力系統計算人員轉向以阻抗矩陣為基礎的逐次代入法。阻抗法改善了系統潮流計算問題的收斂性，解決了導納法無法求解的一些系統的潮流計算，在60年代獲得了廣泛的應用。阻抗法的主要缺點是占用計算機內存大，每次迭代的計算量大。當系統不斷擴大時，這些缺點就更加突出。為了克服阻抗法在內存和速度方面的缺點，60年代中期發展了以阻抗矩陣為基礎的分塊阻抗法。這個方法把一個大系統分割為幾

17、個小的地區系統，在計算機內只需要存儲各個地區系統的阻抗矩陣及它們之間聯絡線的阻抗，這樣不僅大幅度地節省了內存容量，同時也提高了計算速度。克服阻抗法缺點的另一途徑是采用牛頓-拉夫遜法。這是數學中解決非線性方程式的典型方法，有較好的收斂性。在解決電力系統潮流計算問題時，是以導納矩陣為基礎的，因此，只要我們能在迭代過程中盡可能保持方程式系數矩陣的稀疏性，就可以大大提高牛頓法潮流程序的效率。自從60年代中期，在牛頓法中利用了最佳順序消去法以后，牛頓法在收斂性。內存要求。速度方面都超過了阻抗法，成為60年代末期以后廣泛采用的優秀方法。 第1章 電力系統潮流計算概述1.1 電力系統敘述電力工業發展初期，電

18、能是直接在用戶附近的發電站(或稱發電廠)中生產的，各發電站孤立運行。隨著工農業生產和城市的發展，電能的需要量迅速增加，而熱能資源(如煤田)和水能資源豐富的地區又往往遠離用電比較集中的城市和工礦區，為了解決這個矛盾，就需要在動力資源豐富的地區建立大型發電站，然后將電能遠距離輸送給電力用戶。同時，為了提高供電可靠性以及資源利用的綜合經濟性，又把許多分散的各種形式的發電站，通過送電線路和變電所聯系起來。這種由發電機、升壓和降壓變電所，送電線路以及用電設備有機連接起來的整體，即稱為電力系統。電力系統加上發電機的原動機(如汽輪機、水輪機)，原動機的力能部分(如熱力鍋爐、水庫、原子能電站的反應堆)、供熱和

19、用熱設備，則稱為動力系統。現代電力系統提出了“靈活交流輸電與新型直流輸電”的概念。靈活交流輸電技術是指運用固態電子器件與現代自動控制技術對交流電網的電壓、相位角、阻抗、功率以及電路的通斷進行實時閉環控制，從而提高高壓輸電線路的輸送能力和電力系統的穩定水平。新型直流輸電技術是指應用現電力電子技術的最新成果，改善和簡化變流站的造價等。運行方式管理中，潮流是確定電網運行方式的基本出發點；在規劃領域，需要進行潮流分析驗證規劃方案的合理性；在實時運行環境，調度員潮流提供了電網在預想操作情況下電網的潮流分布以校驗運行可靠性。在電力系統調度運行的多個領域都涉及到電網潮流計算。潮流是確定電力網絡運行狀態的基本

20、因素，潮流問題是研究電力系統穩態問題的基礎和前提。1.2 潮流計算簡介電力系統潮流計算是研究電力系統穩態運行情況的一種計算，它根據給定的運行條件及系統接線情況確定整個電力系統各部分的運行狀態：各母線的電壓，各元件中流過的功率，系統的功率損耗等等。在電力系統規劃的設計和現有電力系統運行方式的研究中，都需要利用潮流計算來定量地分析比較供電方案或運行方式的合理性。可靠性和經濟性。此外，電力系統潮流計算也是計算系統動態穩定和靜態穩定的基礎。所以潮流計算是研究電力系統的一種很重要和基礎的計算。電力系統潮流計算也分為離線計算和在線計算兩種，前者主要用于系統規劃設計和安排系統的運行方式，后者則用于正在運行系

21、統的經常監視及實時控制。利用電子數字計算機進行電力系統潮流計算從50年代中期就已經開始。在這20年內，潮流計算曾采用了各種不同的方法，這些方法的發展主要圍繞著對潮流計算的一些基本要求進行的。對潮流計算的要求可以歸納為下面幾點：(1)計算方法的可靠性或收斂性。(2)對計算機內存量的要求。(3)計算速度。(4)計算的方便性和靈活性。電力系統潮流計算問題在數學上是一組多元非線性方程式求解問題，其解法都離不開迭代。因此，對潮流計算方法，首先要求它能可靠地收斂，并給出正確答案。由于電力系統結構及參數的一些特點，并且隨著電力系統不斷擴大，潮流計算的方程式階數也越來越高，對這樣的方程式并不是任何數學方法都能

22、保證給出正確答案的。這種情況成為促使電力系統計算人員不斷尋求新的更可靠方法的重要因素。1.3 潮流計算的意義及其發展電力系統潮流計算是電力系統分析中的一種最基本的計算，是對復雜電力系統正常和故障條件下穩態運行狀態的計算。潮流計算的目標是求取電力系統在給定運行狀態的計算。即節點電壓和功率分布，用以檢查系統各元件是否過負荷。各點電壓是否滿足要求，功率的分布和分配是否合理以及功率損耗等。對現有電力系統的運行和擴建，對新的電力系統進行規劃設計以及對電力系統進行靜態和暫態穩定分析都是以潮流計算為基礎。潮流計算結果可用如電力系統穩態研究，安全估計或最優潮流等對潮流計算的模型和方法有直接影響。實際電力系統的

23、潮流技術那主要采用牛頓-拉夫遜法。在運行方式管理中，潮流是確定電網運行方式的基本出發點；在規劃領域，需要進行潮流分析驗證規劃方案的合理性；在實時運行環境，調度員潮流提供了多個在預想操作情況下電網的潮流分布以校驗運行可靠性。在電力系統調度運行的多個領域都涉及到電網潮流計算。潮流是確定電力網絡運行狀態的基本因素，潮流問題是研究電力系統穩態問題的基礎和前提。在用數字計算機解電力系統潮流問題的開始階段，普遍采取以節點導納矩陣為基礎的逐次代入法。這個方法的原理比較簡單，要求的數字計算機內存量比較下，適應50年代電子計算機制造水平和當時電力系統理論水平。但它的收斂性較差，當系統規模變大時，迭代次數急劇上升

24、，在計算中往往出現迭代不收斂的情況。這就迫使電力系統計算人員轉向以阻抗矩陣為基礎的逐次代入法。阻抗法改善了系統潮流計算問題的收斂性，解決了導納法無法求解的一些系統的潮流計算，在60年代獲得了廣泛的應用。阻抗法的主要缺點是占用計算機內存大，每次迭代的計算量大。當系統不斷擴大時，這些缺點就更加突出。為了克服阻抗法在內存和速度方面的缺點，60年代中期發展了以阻抗矩陣為基礎的分塊阻抗法。這個方法把一個大系統分割為幾個小的地區系統，在計算機內只需要存儲各個地區系統的阻抗矩陣及它們之間聯絡線的阻抗，這樣不僅大幅度地節省了內存容量，同時也提高了計算速度。克服阻抗法缺點的另一途徑是采用牛頓-拉夫遜法。這是數學

25、中解決非線性方程式的典型方法，有較好的收斂性。在解決電力系統潮流計算問題時，是以導納矩陣為基礎的，因此，只要我們能在迭代過程中盡可能保持方程式系數矩陣的稀疏性，就可以大大提高牛頓法潮流程序的效率。自從60年代中期，在牛頓法中利用了最佳順序消去法以后，牛頓法在收斂性。內存要求。速度方面都超過了阻抗法，成為60年代末期以后廣泛采用的優秀方法。第2章 潮流計算的數學模型2.1 導納矩陣的原理及計算方法 自導納和互導納的確定方法電力網絡的節點電壓方程： (2-1)為節點注入電流列向量，注入電流有正有負，注入網絡的電流為正，流出網絡的電流為負。根據這一規定，電源節點的注入電流為正，負荷節點為負。既無電源

26、又無負荷的聯絡節點為零，帶有地方負荷的電源節點為二者代數之和。為節點電壓列向量，由于節點電壓是對稱于參考節點而言的，因而需先選定參考節點。在電力系統中一般以地為參考節點。如整個網絡無接地支路，則需要選定某一節點為參考。設網絡中節點數為（不含參考節點）,則，均為n*n列向量。為n*n階節點導納矩陣。節電導納矩陣的節點電壓方程： 展開為： (2-2)是一個n*n階節點導納矩陣，其階數就等于網絡中除參考節點外的節點數。節點導納矩陣的對角元素 (i=1，2，n)成為自導納。自導納數值上就等于在i節點施加單位電壓，其他節點全部接地時，經節點i注入網絡的電流，因此，它可以定義為： (2-3)節點i的自導納

27、數值上就等于與節點直接連接的所有支路導納的總和。節點導納矩陣的非對角元素 (j=1，2，n;i=1，2，。，n;j=i)稱互導納，由此可得互導納數值上就等于在節點i施加單位電壓，其他節點全部接地時，經節點j注入網絡的電流，因此可定義為： (2-4)節點j，i之間的互導納數值上就等于連接節點j，i支路到導納的負值。顯然，恒等于。互導納的這些性質決定了節點導納矩陣是一個對稱稀疏矩陣。而且，由于每個節點所連接的支路數總有一個限度，隨著網絡中節點數的增加非零元素相對愈來愈少，節點導納矩陣的稀疏度，即零元素數與總元素的比值就愈來愈高。 節點導納矩陣的性質及意義節點導納矩陣的性質：(1)為對稱矩陣，=。如

28、網絡中含有源元件，如移相變壓器，則對稱性不再成立。(2)對無接地支路的節點，其所在行列的元素之和均為零，即：對于有接地支路的節點，其所在行列的元素之和等于該點接地支路的導納。利用這一性質，可以檢驗所形成節點導納矩陣的正確性。(3)具有強對角性：對角元素的值不小于同一行或同一列中任一元素。(4)為稀疏矩陣，因節點i ，j 之間無支路直接相連時=0，這種情況在實際電力系統中非常普遍。矩陣的稀疏性用稀疏度表示，其定義為矩陣中的零元素與全部元素之比，即 ， 式中 為中的零元素。 隨節點數的增加而增加：=50，可達92%；=100，可達90%；=500，可達99%，充分利用節點導納矩陣的稀疏性可節省計算

29、機內存，加快計算速度，這種技巧稱為稀疏技術。節點導納矩陣的意義：是n×n階方陣，其對角元素 (i=1，2，-n)稱為自導納，非對角元素(i，j=1，2，n，)稱為互導納。將節點電壓方程展開為 (2-5)表明，自導納在數值上等于僅在節點i施加單位電壓而其余節點電壓均為零（即其余節點全部接地）時，經節點i注入網絡的電流。其顯然等于與節點i直接相連的所有支路的導納之和。同時可見。表明，互導納在數值上等于僅在節點j施加單位電壓而其余節點電壓均為零時，經節點i注入網絡的電流，其顯然等于()即。為支路的導納，負號表示該電流流出網絡。如節點ij之間無支路直接相連，則該電流為0，從而=0。注意字母幾

30、種不寫法的不同意義：粗體黑字表示導納矩陣，大寫字母代矩陣中的第i行第j列元素，即節點i和節點j之間的互導納。小寫字母i，j支路的導納等于支路阻抗的倒數數。根據定義直接求取節點導納矩陣時，注意以下幾點：(1)節點導納矩陣是方陣，其階數就等于網絡中除去參考節點外的節點數。參考節點一般取大地，編號為零。(2)節點導納矩陣是稀疏矩陣，其各行非零非對角元素就等于與該行相對應節點所連接的不接地支路數。(3)節點導納矩陣的對角元素就等于各該節點所連接導納的總和。因此，與沒有接地支路的節點對應的行或列中，對角元素為非對角元素之和的負值。(4)節點導納矩陣的非對角元素等于連接節點i，j支路導納的負值。因此，一般

31、情況下，節點導納矩陣的對角元素往往大于非對角元素的負值。(5)節點導納矩陣一般是對稱矩陣，這是網絡的互易特性所決定的。從而，一般只要求求取這個矩陣的上三角或下三角部分。 非標準變比變壓器等值電路變壓器型等值電路更便于計算機反復計算,更適宜于復雜網絡的潮流計算.雙繞組變壓器可用阻抗與一個理想變壓器串聯的電路表示.理想變壓器只是一個參數,那就是變比。現在變壓器阻抗按實際變比歸算到低壓側為例,推導出變壓器型等值電路.圖2-1雙繞組變壓器原理圖流入和流出理想變壓器的功率相等： (2-6)式中,是理想變壓器的變比,和 分別為變壓器高,低繞組的實際電壓.從圖b直接可得： （2-7）從而可得: （2-8）式

32、中,又因節點電流方程應具有如下形式: （2-9）將式（1-8）與（1-9）比較，得：因此可得各支路導納為: （2-10）由此可得用導納表示的變壓器型等值電路:圖2-2變壓器型等值電路2.2 潮流計算的基本方程在潮流問題中，任何復雜的電力系統都可以歸納為以下元件（參數）組成：（1）發電機（注入電流或功率）（2）負荷（注入負的電流或功率）（3）輸電線支路（電阻，電抗）（4）變壓器支路（電阻，電抗，變比）（5）母線上的對地支路（阻抗和導納）（6）線路上的對地支路（一般為線路充電點容導納）集中了以上各類型的元件的簡單網絡如圖圖2-3潮流計算用的電網結構圖圖2-4潮流計算等值網絡采用導納矩陣時，節點注入

33、電流和節點電壓構成以下線性方程組:,其中：可展開如下形式: (2-12)由于實際電網中測量的節點注入量一般不是電流而是功率，因此必須將式中的注入電流用節點注入功率來表示。節點功率與節點電流之間的關系為: (2-13)式中，因此用導納矩陣時，節點可以表示為把這個關系代入式中得： （2-14）式（3-4 ）就是電力系統潮流計算的數學模型-潮流方程。它具有如下特點：(1)它是一組代數方程，因而表征的是電力系統的穩定運行特性。(2)它是一組非線性方程，因而只能用迭代方法求其數值解。(3)由于方程中的電壓和導納既可以表為直角坐標，又可表為極坐標，因而潮流方程有多種表達形式-極坐標形式，直角坐標形式和混合

34、坐標形式。取 ，得到潮流方程的極坐標形式： (2-15)取 ， ，得到潮流方程的直角坐標形式： (2-16)取，得到潮流方程的混合坐標形式： (2-17)不同坐標形式的潮流方程適用于不同的迭代解法。例如：利用牛頓-拉夫遜迭代法求解，以直角坐標和混合坐標形式的潮流方程為方便；而P-Q解耦法是在混合坐標形式的基礎上發展而成，故當然采用混合坐標形式。(4)它是一組n個復數方程，因而實數方程數為2n個但方程中共含4n個變量：P，Q，U和，i=1，2，n，故必須先指定2n個變量才能求解。2.3 電力系統節點分類用一般的電路理論求解網絡方程，目的是給出電壓源(或電流源)研究網絡內的電流(或電壓)分布，作為

35、基礎的方程式，一般用線性代數方程式表示。然而在電力系統中，給出發電機或負荷連接母線上電壓或電流(都是向量)的情況是很少的，一般是給出發電機母線上發電機的有功功率(P)和母線電壓的幅值(U)，給出負荷母線上負荷消耗的有功功率(P)和無功功率(Q)。主要目的是由這些已知量去求電力系統內的各種電氣量。所以，根據電力系統中各節點性質的不同，很自然地把節點分成三類：(1)PQ節點對這一類點，事先給定的是節點功率(P，Q)，待求的未知量是節點電壓向量(U，)，所以叫PQ節點。通常變電所母線都是PQ節點，當某些發電機的輸出功率P。Q給定時，也作為PQ節點。PQ節點上的發電機稱之為PQ機(或PQ給定型發電機)

36、。在潮流計算中，系統大部分節點屬于PQ節點。(2)PU節點這類節點給出的參數是該節點的有功功率P及電壓幅值U，待求量為該節點的無功功率Q及電壓向量的相角。這類節點在運行中往往要有一定可調節的無功電源。用以維持給定的電壓值。通常選擇有一定無功功率儲備的發電機母線或者變電所有無功補償設備的母線做PU節點處理。PU節點上的發電機稱為PU機(或PU給定型發電機)。(3)平衡節點在潮流計算中，這類節點一般只設一個。對該節點，給定其電壓值，并在計算中取該節點電壓向量的方向作為參考軸，相當于給定該點電壓向量的角度為零。也就是說，對平衡節點給定的運行參數是U和，因此有城為U節點，而待求量是該節點的P。Q，整個

37、系統的功率平衡由這一節點承擔。關于平衡節點的選擇，一般選擇系統中擔任調頻調壓的某一發電廠(或發電機)，有時也可能按其他原則選擇，例如，為提高計算的收斂性。可以選擇出線數多或者靠近電網中心的發電廠母線作平衡節點。以上三類節點4個運行參數P。Q。U。中，已知量都是兩個，待求量也是兩個，只是類型不同而已。2.4 潮流計算的約束條件電力系統運行必須滿足一定技術和經濟上的要求。這些要求夠成了潮流問題中某些變量的約束條件，常用的約束條件如下：(1)節點電壓應滿足： (2-18)(2)從保證電能質量和供電安全的要求來看，電力系統的所有電氣設備都必須運行在額定電壓附近。PU節點電壓幅值必須按上述條件給定。因此

38、，這一約束條件對PQ節點而言(3)節點的有功功率和無功功率應滿足： (2-19)PQ節點的有功功率和無功功率，以及PU節點的有功功率，在給定是就必須滿足上述條件，因此，對平衡節點的P和Q以及PU節點的Q應按上述條件進行檢驗。(4)節點之間電壓的相位差應滿足： (2-30)為了保證系統運行的穩定性，要求某些輸電線路兩端的電壓相位不超過一定的數值。這一約束的主要意義就在于此。因此，潮流計算可以歸結為求解一組非線性方程組，并使其解答滿足一定的約束條件。常用的方法是迭代法和牛頓法，在計算過程中，或得出結果之后用約束條件進行檢驗。如果不能滿足要求，則應修改某些變量的給定值，甚至修改系統的運行方式，重新進

39、行計算。第3章 牛頓拉夫遜法概述3.1 牛頓拉夫遜法基本原理電力系統潮流計算是電力系統分析中的一種最基本的計算，是對復雜電力系統正常和故障條件下穩態運行狀態的計算。潮流計算的目標是求取電力系統在給定運行狀態的計算。即節點電壓和功率分布，用以檢查系統各元件是否過負荷。各點電壓是否滿足要求，功率的分布和分配是否合理以及功率損耗等。對現有電力系統的運行和擴建，對新的電力系統進行規劃設計以及對電力系統進行靜態和暫態穩定分析都是以潮流計算為基礎。潮流計算結果可用如電力系統穩態研究，安全估計或最優潮流等對潮流計算的模型和方法有直接影響。實際電力系統的潮流技術那主要采用牛頓-拉夫遜法。牛頓-拉夫遜法(簡稱牛

40、頓法)在數學上是求解非線性代數方程式的有效方法。其要點是把非線性方程式的求解過程變成反復地對相應的線性方程式進行求解的過程。即通常所稱的逐次線性化過程。對于非線性代數方程組： 即 (3-1)在待求量x的某一個初始估計值附近，將上式展開成泰勒級數并略去二階及以上的高階項，得到如下的經線性化的方程組： (3-2)上式稱之為牛頓法的修正方程式。由此可以求得第一次迭代的修正量： (3-3)將和相加，得到變量的第一次改進值。接著就從出發，重復上述計算過程。因此從一定的初值出發，應用牛頓法求解的迭代格式為： (3-4) (3-5)上兩式中：是函數對于變量x的一階偏導數矩陣，即雅可比矩陣J;k為迭代次數有上

41、式可見，牛頓法的核心便是反復形式并求解修正方程式。牛頓法當初始估計值和方程的精確解足夠接近時，收斂速度非常快，具有平方收斂特性。牛頓潮流算法突出的優點是收斂速度快，若選擇到一個較好的初值，算法將具有平方收斂特性，一般迭代45次便可以收斂到一個非常精確的解。而且其迭代次數與所計算網絡的規模基本無關。牛頓法也具有良好的收斂可靠性，對于對以節點導納矩陣為基礎的高斯法呈病態的系統，牛頓法也能可靠收斂。牛頓法所需的內存量及每次迭代所需時間均較高斯法多。牛頓法的可靠收斂取決于有一個良好的啟動初值。如果初值選擇不當，算法有可能根本不收斂或收斂到一個無法運行的節點上。對于正常運行的系統，各節點電壓一般均在額定

42、值附近，偏移不會太大，并且各節點間的相位角差也不大，所以對各節點可以采用統一的電壓初值(也稱為平直電壓)，如假定： 或 (3-6)這樣一般能得到滿意的結果。但若系統因無功緊張或其它原因導致電壓質量很差或有重載線路而節點間角差很大時，仍用上述初始電壓就有可能出現問題。解決這個問題的辦法可以用高斯法迭代12次，以此迭代結果作為牛頓法的初值。也可以先用直流法潮流求解一次以求得一個較好的角度初值，然后轉入牛頓法迭代。3.2 牛頓-拉夫遜法潮流求解過程以下討論的是用直角坐標形式的牛頓拉夫遜法潮流的求解過程。當采用直角坐標時，潮流問題的待求量為各節點電壓的實部和虛部兩個分量由于平衡節點的電壓向量是給定的，

43、因此待求兩共需要2(n-1)個方程式。事實上，除了平衡節點的功率方程式在迭代過程中沒有約束作用以外，其余每個節點都可以列出兩個方程式。對PQ節點來說，和是給定的，因而可以寫出：（3-7）對PV節點來說，給定量是和，因此可以列出：（3-8）求解過程大致可以分為以下步驟：(1)形成節點導納矩陣。(2)將各節點電壓設初值U。(3)將節點初值代入相關求式，求出修正方程式的常數項向量。(4)將節點電壓初值代入求式，求出雅可比矩陣元素。(5)求解修正方程，求修正向量。(6)求取節點電壓的新值。(7)檢查是否收斂，如不收斂，則以各節點電壓的新值作為初值自第3步重新開始進行狹義次迭代，否則轉入下一步。(8)計

44、算支路功率分布，PV節點無功功率和平衡節點柱入功率。以直角坐標系形式表示：迭代推算式采用直角坐標時,節點電壓相量及復數導納可表示為: (3-9)將以上二關系式代入上式中,展開并分開實部和虛部;假定系統中的第1,2,m號為PQ節點,第m+1,m+2,n-1為PV節點,根據節點性質的不同,得到如下迭代推算式:對于PQ節點： (3-10)對于PV節點：(3-11)對于平衡節點：平衡節點只設一個,電壓為已知,不參見迭代,其電壓為: (3-12)修正方程式(2-3-5)和(2-3-6)兩組迭代式工包括2(n-1)個方程.選定電壓初值及變量修正量號之后代入式(2-3-5)和(2-3-6),并將其按泰勒級數

45、展開,略去二次方程及以各項,得到修正方程如下: (3-13)(3-14)雅可比矩陣各元素的算式：式(3-2-8)中, 雅可比矩陣中的各元素可通過對式(3-2-4)和(3-2-5)進行偏導而求得。當時，雅可比矩陣中非對角元素為： (3-15)當時,雅可比矩陣中對角元素為: (3-16)由式(3-2-9和(3-2-10)看出,雅可比矩陣的特點:(1)陣中各元素是節點電壓的函數,在迭代過程中,這些元素隨著節點電壓的變化而變化。(2)導納矩陣中的某些非對角元素為零時,雅可比矩陣中對應的元素也是為零。若,則必有。(3)可比矩陣不是對稱矩陣; 雅可比矩陣各元素的表示如下:3.3 牛頓拉夫遜法的程序框圖以對

46、于PU節點，計算。對于PQ節點，計算形成導納矩陣輸入原始數據啟動給定電壓初值置按系統的潮流分布計算平衡節點的功率及線路功率是否YN計算雅可比矩陣各元素輸出求修正方程式，求以用，求修正節點電壓 圖3-1 牛頓拉夫遜法程序框圖 第4章 MATLAB概述4.1 MATLAB簡介目前電子計算機已廣泛應用于電力系統的分析計算，潮流計算是其基本應用軟件之一。現有很多潮流計算方法。對潮流計算方法有五方面的要求：(1)計算速度快。(2)內存需要少。(3)計算結果有良好的可靠性和可信性。(4)適應性好，亦即能處理變壓器變比調整、系統元件的不同描述和與其它程序配合的能力強。(5)簡單。MATLAB是一種交互式、面

47、向對象的程序設計語言，廣泛應用于工業界與學術界，主要用于矩陣運算，同時在數值分析、自動控制模擬、數字信號處理、動態分析、繪圖等方面也具有強大的功能。MATLAB程序設計語言結構完整，且具有優良的移植性，它的基本數據元素是不需要定義的數組。它可以高效率地解決工業計算問題，特別是關于矩陣和矢量的計算。MATLAB與C語言和FORTRAN語言相比更容易被掌握。通過M語言，可以用類似數學公式的方式來編寫算法，大大降低了程序所需的難度并節省了時間，從而可把主要的精力集中在算法的構思而不是編程上。另外，MATLAB提供了一種特殊的工具：工具箱（TOOLBOXES）.這些工具箱主要包括：信號處理（SIGNA

48、L PROCESSING）、控制系統（CONTROL SYSTEMS）、神經網絡（NEURAL NETWORKS）、模糊邏輯(FUZZY LOGIC)、小波(WAVELETS)和模擬（SIMULATION）等等。不同領域、不同層次的用戶通過相應工具的學習和應用，可以方便地進行計算、分析及設計工作。MATLAB設計中，原始數據的填寫格式是很關鍵的一個環節，它與程序使用的方便性和靈活性有著直接的關系。原始數據輸入格式的設計，主要應從使用的角度出發，原則是簡單明了，便于修改。4.2 MATLAB應用在潮流計算中的優勢MATLAB由一系列工具組成。這些工具方便用戶使用MATLAB的函數和文件，其中許多

49、工具采用的是圖形用戶界面。包括MATLAB桌面和命令窗口、歷史命令窗口、編輯器和調試器、路徑搜索和用于用戶瀏覽幫助、工作空間、文件的瀏覽器。隨著MATLAB的商業化以及軟件本身的不斷升級，MATLAB的用戶界面也越來越精致，更加接近Windows的標準界面，人機交互性更強，操作更簡單。而且新版本的MATLAB提供了完整的聯機查詢、幫助系統，極大的方便了用戶的使用。簡單的編程環境提供了比較完備的調試系統，程序不必經過編譯就可以直接運行，而且能夠及時地報告出現的錯誤及進行出錯原因分析。MATLAB是一個高級的矩陣/陣列語言，它包含控制語句、函數、數據結構、輸入和輸出和面向對象編程特點。用戶可以在命

50、令窗口中將輸入語句與執行命令同步，也可以先編寫好一個較大的復雜的應用程序（M文件）后再一起運行。新版本的MATLAB語言是基于最為流行的C+語言基礎上的，因此語法特征與C+語言極為相似，而且更加簡單，更加符合科技人員對數學表達式的書寫格式。使之更利于非計算機專業的科技人員使用。MATLAB是一個包含大量計算算法的集合。其擁有600多個工程中要用到的數學運算函數，可以方便的實現用戶所需的各種計算功能。函數中所使用的算法都是科研和工程計算中的最新研究成果，而前經過了各種優化和容錯處理。在通常情況下，可以用它來代替底層編程語言，如C和C+ 。在計算要求相同的情況下，使用MATLAB的編程工作量會大大

51、減少。基于MATLAB的電力系統潮流計算使計算機在計算、分析、研究復雜的電力系統潮流分布問題上又前進了一步。不管采用什么算法,所有的潮流計算都是基于矩陣的迭代運算。而MATLAB 語言正是以處理矩陣見長, 實踐證明，MATLAB 語言在電力系統潮流計算仿真研究中的應用是可行的,而且由于其強大的矩陣處理功能,完全可以應用于電力系統的其它分析計算中；用MATLAB語言編程效率高, 程序調試十分方便,可大大縮減軟件開發周期，如果像控制界一樣開發出電力系統自己的專用工具箱,將系統分析用的一些基本計算以函數的形式直接調用,那么更高層次的系統軟件也可以很容易地實現。4.3 矩陣的運算矩陣是MATLAB數據

52、存儲的基本單元，而矩陣的運算是MATLAB語言的核心，在MATLAB語言系統中幾乎一切運算均是以對矩陣的操作為基礎的。矩陣的基本數學運算包括矩陣的四則運算、與常數的運算、逆運算、行列式運算、秩運算、特征值運算等基本函數運算，這里進行簡單介紹。(1)四則運算：矩陣的加、減、乘運算符分別為“+，*”，用法與數字運算幾乎相同，但計算時要滿足其數學要求在MATLAB中矩陣的除法有兩種形式：左除“”和右除“/”。在傳統的MATLAB算法中，右除是先計算矩陣的逆再相乘，而左除則不需要計算逆矩陣直接進行除運算。通常右除要快一點，但左除可避免被除矩陣的奇異性所帶來的麻煩。在MATLAB6中兩者的區別不太大。(

53、2)與常數的運算：常數與矩陣的運算即是同該矩陣的每一元素進行運算。但需注意進行數除時，常數通常只能做除數。(3)基本函數運算：矩陣的函數運算是矩陣運算中最實用的部分，常用的主要有以下幾個：det(a) 求矩陣a的行列式eig(a) 求矩陣a的特征值inv(a)或a (-1) 求矩陣a的逆矩陣rank(a) 求矩陣a的秩trace(a) 求矩陣a的跡（對角線元素之和）我們在進行工程計算時常常遇到矩陣對應元素之間的運算。這種運算不同于前面講的數學運算，為有所區別，我們稱之為數組運算。(4)基本數學運算：數組的加、減與矩陣的加、減運算完全相同。而乘除法運算有相當大的區別，數組的乘除法是指兩同維數組對

54、應元素之間的乘除法，它們的運算符為“.*”和“./”或“.”。前面講過常數與矩陣的除法運算中常數只能做除數。在數組運算中有了“對應關系”的規定，數組與常數之間的除法運算沒有任何限制。另外，矩陣的數組運算中還有冪運算（運算符為. ）、指數運算（exp）、對數運算（log）、和開方運算（sqrt）等。有了“對應元素”的規定，數組的運算實質上就是針對數組內部的每個元素進行的。矩陣的冪運算與數組的冪運算有很大的區別。(5)邏輯關系運算：邏輯運算是MATLAB中數組運算所特有的一種運算形式，也是幾乎所有的高級語言普遍適用的一種運算。第5章 潮流計算主界面設計實現5.1 主界面介紹主界面是進行潮流計算的的

55、操作界面，能完成數據的初始化、進行潮流計算、并顯示潮流計算的相關結果以供用戶分析。具體來講就是點擊第一個彈出框選擇潮流計算數據組，并在主界面上顯示所選擇數據組的電路接線圖，單擊計算按鈕就能完成潮流計算并將潮流計算結果保存下來，單擊第二個彈出框可以查看具體的結果，同時還將畫出潮流計算電壓迭代曲線。5.2 數據初始化本設計將數據保存在Excel表格里，在主界面上單擊彈出框選擇數據組從而從Excel表格里調入數據完成數據初始化。本設計定義了7個數據供主程序調用，他們分別是節點數n，支路數n1，平衡節點標號isb，精度值pr，節點參數矩陣B1，支路參數矩陣B2和節點號及其對地阻抗形成的舉證X。從Excel表格中調入數據采用Xlsread函數，這樣Excel表格中的數據就以