1、3減我孕畢業(yè)設(shè)計（論文）題,目配電網(wǎng)潮流計算與程序設(shè)計摘要 TOC o 1-5 h z 一配電網(wǎng)潮流概述5 HYPERLINK l bookmark83 o Current Document 1.1配電網(wǎng)潮流計算的目的與意義5 HYPERLINK l bookmark86 o Current Document 1.2潮流計算方法概述5 HYPERLINK l bookmark89 o Current Document 牛頓拉夫遜法6 HYPERLINK l bookmark107 o Current Document 快速解耦法6 HYPERLINK l bookmark117 o Curre

2、nt Document 回路阻抗法9前推回代法11 HYPERLINK l bookmark133 o Current Document 1.3本文工作11 HYPERLINK l bookmark136 o Current Document 二.配電網(wǎng)網(wǎng)絡(luò)模型11 HYPERLINK l bookmark139 o Current Document 2.1元件模型11電力線路的數(shù)學(xué)模型11 HYPERLINK l bookmark145 o Current Document 變壓器的等值電路13 HYPERLINK l bookmark169 o Current Document 2.2網(wǎng)絡(luò)

3、模型15 HYPERLINK l bookmark172 o Current Document 三：基于mat lab的配電網(wǎng)潮流計算算法16 HYPERLINK l bookmark175 o Current Document 3.1配電網(wǎng)潮流計算算法原理16 HYPERLINK l bookmark329 o Current Document 3.2 matlab 的概述19 HYPERLINK l bookmark332 o Current Document 3.3程序設(shè)計21牛頓一拉夫遜法潮流求解過程21牛頓一拉夫遜法的程序框圖25 HYPERLINK l bookmark381 o

4、Current Document 四：算例27 HYPERLINK l bookmark384 o Current Document 參考文獻(xiàn)28致謝29配電網(wǎng)潮流計算與程序設(shè)計學(xué) 生：石昊晨指導(dǎo)教師：劉會家（三峽大學(xué) 國際文化交流學(xué)院）摘 要：本文首先分析了配電網(wǎng)的特點及對算法的要求，然后建立配電網(wǎng)潮流計 算模型。針對配電網(wǎng)潮流計算的現(xiàn)狀進(jìn)行了全面分析，深入討論了目前各方法的 特點，并從收斂性及其他性能指標(biāo)進(jìn)行了比較分析；詳細(xì)研究用的比較廣泛的牛 頓拉夫遜法，并以廣度優(yōu)先順序搜索策略作為理論基礎(chǔ)。針對某地區(qū)配電網(wǎng) 的具體情況，選取IOKV的配電網(wǎng)子系統(tǒng)進(jìn)行潮流計算。利用MATLAB 2009

5、a進(jìn) 行了基于牛頓一一拉夫遜法的配電網(wǎng)的潮流計算程序。由計算結(jié)果可知，該算法 具有一定的優(yōu)越性，軟件的開發(fā)具有一定的實用性。關(guān)鍵詞：電力系統(tǒng)，配電網(wǎng)潮流，牛頓一一拉夫遜法，MATLAB程序設(shè)計Abstract: In this paper, ungrounded system, the characteristics of non-zero sequence path, a three-phase decoupled power flow calculation method. This method ignores the influence of zero sequence compone

6、nts, making the three-phase asymmetrical load caused by phase coupling decoupling to be achieved by the phase flow calculation. The algorithm flow algorithm to the existing distribution network in the three-phase node voltage equation 3n-order decomposition of the node voltage equation of three n-or

7、der, so no matter what kind of algorithm can greatly save memory and computation for the distribution network to achieve by phase analysis provides a good way. In this paper, a system of 36 nodes to verify the results show that the method can fully into account the impact of unbalanced three-phase l

8、oads, a better computational speed and accuracy.Keywords: power systems, phase decoupling, power flow, back/forward sweep algorithm一.電力系統(tǒng)潮流概述1.1配電網(wǎng)潮流計算的目的與意義電力系統(tǒng)潮流計算是研究電力系統(tǒng)穩(wěn)態(tài)運行情況的一種基本電氣計算，是電力系統(tǒng)規(guī)劃和運 營中不可缺少的一個重要組成部分?？梢哉f，潮流計算是電力系統(tǒng)分析中最基本、最重要的 計算，是電力系統(tǒng)安全、經(jīng)濟(jì)分析和實時控制與調(diào)度的基礎(chǔ)。潮流計算的任務(wù)是根據(jù)給定的 運行條件和網(wǎng)絡(luò)結(jié)構(gòu)確定整個系統(tǒng)的運行狀態(tài)

9、，如各母線上的電壓（幅值及相角）、網(wǎng)絡(luò)中 的功率分布及功率損耗等。即潮流計算是對電力系統(tǒng)的功率分布和電壓分布的計算其具體 任務(wù)就是編制系統(tǒng)的調(diào)度計劃和電氣設(shè)備檢修計劃，確定電力系統(tǒng)中變壓器分接頭位置和 系統(tǒng)中樞點與電壓控制點的電壓曲線，進(jìn)行事故運行方式的分析，為電力系統(tǒng)短路和穩(wěn)定 計算提供數(shù)據(jù)，為繼電保護(hù)及自動裝置整定與電力系統(tǒng)設(shè)計和規(guī)劃提供依據(jù)等。潮流計算的 目的是對現(xiàn)有電力系統(tǒng)的正常運行狀態(tài)進(jìn)行分析，以提示必要的改進(jìn)措施，同時為新建系統(tǒng) 或擴(kuò)建系統(tǒng)的有關(guān)分析、計算打下基礎(chǔ)。配電網(wǎng)潮流計算是配電網(wǎng)經(jīng)濟(jì)運行、系統(tǒng)分析等的重要基礎(chǔ)，但由于配電網(wǎng)與輸電網(wǎng)有 著明顯的差異:配電網(wǎng)具有環(huán)形結(jié)構(gòu)，而通常

10、以開環(huán)方式運行。通常呈輻射狀，支路比值較大， 分支線較多；配電線路中的R/X比值偏大使輸電網(wǎng)中常用的潮流計算算法如傳統(tǒng)的牛頓法 和快速分解法在應(yīng)用于配電網(wǎng)潮流計算時容易形成病態(tài)而無法收斂，因此，研究適合于配電 網(wǎng)的潮流算法也是至關(guān)重要的。目前，輸電系統(tǒng)潮流計算方法已較為成熟，而且獲得了廣泛的實際應(yīng)用。但隨著電力系 統(tǒng)規(guī)模的不斷擴(kuò)大，潮流方程的階數(shù)越來越高，對這種規(guī)模的方程并不是采用任何數(shù)學(xué)方法 都能保證給出正確答案的，因此，這也成為促使電力系統(tǒng)研究人員不斷尋求新的、更可靠的 潮流計算方法的動力。隨著現(xiàn)代電力系統(tǒng)大系統(tǒng)、強非線性與多元件的特點日益突出，其計算量與計算復(fù)雜度急劇 增加。舊的計算機(jī)軟

11、件在處理潮流計算時，其速度已無法滿足大電網(wǎng)模擬和實時控制的仿真 要求，而高效的潮流問題的相關(guān)軟件的研究已成為大規(guī)模電力系統(tǒng)仿真計算的關(guān)鍵。1.2潮流計算方法概述與輸電網(wǎng)相比，配電網(wǎng)的網(wǎng)絡(luò)結(jié)構(gòu)有著明顯的差異：配電網(wǎng)的網(wǎng)絡(luò)呈現(xiàn)輻射 狀，在正常運行是開環(huán)的，只有在倒換負(fù)荷或發(fā)生故障時才有可能出現(xiàn)短時環(huán)網(wǎng) 運行或多電源運行的情況；配電線路的總長度較輸電網(wǎng)絡(luò)要長且分支較多，配電 線的線徑比輸電線細(xì)，導(dǎo)致配電網(wǎng)的R/ X較大，無法滿足氣 fi(x )(3-3)將Ax(0)和x(0)相加，得到變量的第一次改進(jìn)值x(i)。接著就從x(i)出發(fā)，重復(fù)上述計算過程。因此從一定的初值x(0)出發(fā)，應(yīng)用牛頓法求解的迭

12、代格式為： HYPERLINK l bookmark101 o Current Document f( xk )A xk( = f( x ()(3-4)x(k + 1 = x k(書 Ax k (3-5)上兩式中：f(x)是函數(shù)f (x)對于變量x的一階偏導(dǎo)數(shù)矩陣，即雅可比矩陣J; k為迭代次數(shù)。有上式可見，牛頓法的核心便是反復(fù)形式并求解修正方程式。牛頓法當(dāng)初始估計值雙0）和方程的精確解足夠接近時，收斂速度非常快，具有平方收斂特性。牛頓潮流算法突出的優(yōu)點是收斂速度快，若選擇到一個較好的初值，算法將 具有平方收斂特性，一般迭代45次便可以收斂到一個非常精確的解。而且其迭 代次數(shù)與所計算網(wǎng)絡(luò)的規(guī)模

13、基本無關(guān)。牛頓法也具有良好的收斂可靠性，對于對 以節(jié)點導(dǎo)納矩陣為基礎(chǔ)的高斯法呈病態(tài)的系統(tǒng)，牛頓法也能可靠收斂。牛頓法所 需的內(nèi)存量及每次迭代所需時間均較高斯法多。牛頓法的可靠收斂取決于有一個良好的啟動初值。如果初值選擇不當(dāng)，算法 有可能根本不收斂或收斂到一個無法運行的節(jié)點上。對于正常運行的系統(tǒng)，各節(jié) 點電壓一般均在額定值附近，偏移不會太大，并且各節(jié)點間的相位角差也不大， 所以對各節(jié)點可以采用統(tǒng)一的電壓初值（也稱為平直電壓），如假定：U（0）= 10（ 0）= 0 或 e（ 0）= 1 f（。匚 0 （i = q1, 2 , n,誅;s （3-6）iiii這樣一般能得到滿意的結(jié)果。但若系統(tǒng)因無功

14、緊張或其它原因?qū)е码妷嘿|(zhì)量 很差或有重載線路而節(jié)點間角差很大時，仍用上述初始電壓就有可能出現(xiàn)問題。 解決這個問題的辦法可以用高斯法迭代12次，以此迭代結(jié)果作為牛頓法的初 值。也可以先用直流法潮流求解一次以求得一個較好的角度初值，然后轉(zhuǎn)入牛頓 法迭代?？焖俳怦罘榱烁倪M(jìn)牛頓法在內(nèi)容占用量及計算速度方面的不足，早在1974年有人提 出的快速解耦法（對稱P-Q分解法）是較成功的一種算法；它是密切結(jié)合高壓電 力系統(tǒng)固有特點，對牛頓法改進(jìn)后得到的一種方法。原理是根據(jù)系統(tǒng)有功決定于 電壓相角的變化，而無功主要決定于電壓模值的變化這一特性，并進(jìn)行合理假設(shè):（1）線路兩端的相角差不大，且|氣 |Bj ,即認(rèn)為

15、cos0tj = 1 ；G sin 0 B ；（2）與節(jié)點無功功率對應(yīng)的導(dǎo)納Q /遠(yuǎn)小于節(jié)點的自導(dǎo)納b，即Q Q / AP / U = B&Q / U = BAU,(215)式中：B，、b是由節(jié)點導(dǎo)納矩陣的虛部構(gòu)成的常數(shù)對稱矩陣，可有XB、BX等這種方法具有簡單、快速、內(nèi)存節(jié)省且收斂可靠的優(yōu)點，是廣泛應(yīng)用于高壓網(wǎng)在 線處理計算的方法。該方法存在的問題是R/X比值敏感，用于配電網(wǎng)可能迭代 次數(shù)過多或不收斂。針對這一問題，提出了一種改進(jìn)的快速解耦法。該方法的特點是，它根據(jù)配 電網(wǎng)的輻射型特點，從一種新概念上構(gòu)造出潮流方程，即前一節(jié)點的電壓電流用 含后一節(jié)點的電壓和電流的關(guān)系式表示，即Wk-1kk(

16、216)其中，wk為第k個節(jié)點的電壓和對應(yīng)的支路電流矩陣，gk為前后兩個節(jié)點的關(guān)系方程。根據(jù)邊界條件 七=0,U= U，可建立潮流方程如下所示:(217)其中，U0(氣)為按(2-16)從末端遞推到始端形成的以末端電壓氣為變量的方程，f (U)的雅可比矩陣可以表示為從饋線末端到始端所有支路雅可比矩陣的乘積， 即其中J (U)=afaUaUaUawaUaUaUk-1a其中J (U)=afaUaUaUawaUaUaUk-1aIaUaUk-1k+1aIaIk+1agn-1 aU(218)(219)這樣，一方面可以減少方程的數(shù)目，使之等于支路數(shù)；另一方面能夠充分利 用配電網(wǎng)的輻射型結(jié)構(gòu)導(dǎo)致的數(shù)值特性，

17、將雅可比矩陣簡化為一個三角矩陣，使 其求解的實質(zhì)變?yōu)橐环N前推回推算法，從而簡化了運算，并極大提高了其收斂性 能。文獻(xiàn)6通過以下假設(shè)將式(2-19)中的G化為單位矩陣：k節(jié)點K電壓氣的微小變化，將引起前一節(jié)點Uki幾乎相同的變化，因此 左上角項的所有元素近似為1；電流Ik+i的微小變化對氣i影響很小，因此右上角項的所有元素近似為0；Uk的微小變化對Ik影響很小，因此左下角項的所有元素近似為0；(4)電流L的微小變化時，將引起Ik幾乎相同的變化，因此左下角項的所有 元素近似為1?；芈纷杩狗ㄔ谝话汶娏ο到y(tǒng)(發(fā)、輸電網(wǎng)絡(luò))中，各節(jié)點和大地間有發(fā)電機(jī)、負(fù)荷、線路 電容等對地支路，節(jié)點和節(jié)點間也有輸電線路

18、和變壓器支路，使得系統(tǒng)的節(jié)點方 程式數(shù)小于回路方程式數(shù)。因而，一般電力系統(tǒng)的分析計算采用節(jié)點電壓方程為 宜。但對于低電壓配電網(wǎng)絡(luò)，由于一般不計配電線路對地充電電容的影響，并忽 略變壓器的對地導(dǎo)納，網(wǎng)絡(luò)中樹支數(shù)將總大于連支數(shù)，因而適合采用回路電流方 程進(jìn)行分析。因此提出了一種基于回路方程的潮流算法，并稱之為直接解:Direct SolutionMethod)。由于它基于回路阻抗方程，稱之為回路阻抗法。該方法將各節(jié) 點的負(fù)荷用恒定阻抗表示，從饋線節(jié)點到每一個負(fù)荷節(jié)點形成一條回路，以回路 電流為變量，根據(jù)基爾霍夫電壓定律，可列出回路電流方程式組：= Z I + Z I + + Z I1(1,1) 1

19、(1,2) 2 n n:= Z I + Z I + + Z In(n , 1 ) 1 n ( , 2 ) 2 n n n ( , )(220)式中，*為根節(jié)點電壓，匕為第i條回路上的回路電流(等于負(fù)荷節(jié)點i的 負(fù)荷電流)，Z”為第i條回路的自阻抗(等于節(jié)點i與根節(jié)點s之間的支路阻抗 和，加上節(jié)點i的負(fù)荷阻抗)，Z為第i條回路和第j條回路的互阻抗(等于節(jié) 點i與節(jié)點j到根節(jié)點s的共同支路阻抗和)。設(shè)負(fù)荷節(jié)點數(shù)為L，則回路阻抗矩陣Z是一個LXL維的不含零元素的方陣。采用LU分解方法對方程式(2-20)進(jìn)行分解，可求出回路電流，也就得到各 個負(fù)荷節(jié)點的負(fù)荷電流。然后可求出各條支路上的電壓降，進(jìn)而可求

20、得各節(jié)點的 電壓和負(fù)荷節(jié)點的功率，反復(fù)迭代，直到求得的負(fù)荷節(jié)點功率與給定負(fù)荷的差值 滿足一定的精度要求為止。在回路阻抗陣中有許多相同的元素，實際上只有網(wǎng)絡(luò) 支路數(shù)目個不同元素。但是在一般的編號方式下，這些不同的元素交叉混雜，無 規(guī)律性可言。為了減少占用計算機(jī)的存儲容量，文獻(xiàn)8采用了一種特別的節(jié)點 和支路編號方案，在這種編號方案下，回路阻抗矩陣Z和它三角分解得到的上 三角矩陣U中的元素能夠有規(guī)律地排列，即許多相同的元素集中排列在一起， 因而可以借用“稀疏存儲”技術(shù)，只存儲其中不同的元素，只是這種編號方案太 復(fù)雜而不易實現(xiàn)。在求U矩陣的元素時，文獻(xiàn)6也通過采用一些求解技巧，提 高了計算速度。但這些

21、技巧不適用于在U矩陣中占很大比例的對角元素和同一 行與它緊相鄰的元素，因而限制了求解速度的提高。特別地，回路阻抗法處理網(wǎng) 孔的能力較強，它對增加一條環(huán)路后的處理方法比較簡單：假定連接節(jié)點i和i (i ; B = I %S /100U2T 0N T 0 NN(22)式中：rt 一變壓器高低壓繞組的總電阻(Q )；XT 一變壓器高低壓繞組的總電抗(Q )；G 一變壓器的電導(dǎo)(S)；TB 變壓器的電納(S)；TPk 一變壓器的短路損耗(kW)；Sn 一變壓器的額定容量(MVA)；U n 一變壓器的額定電壓(kV)；U%一變壓器的短路電壓百分值；10% 一變壓器的空載電流百分值；R|r-j |GrjB

22、rUN 圖2-4雙繞組變壓器的T型等值電路 n型等值電路也就是等值變壓器模型：| y 21 |2_yiO Wy20圖2-5雙繞組變壓器n型等值電路不論采用有名制或標(biāo)幺制，凡涉及多電壓級網(wǎng)絡(luò)的計算，都必須將網(wǎng)絡(luò)中所 有參數(shù)和變量歸算至同一電壓級。這是因為r形或t形等值電路做變壓器模型 時，這些等值電路模型并不能體現(xiàn)變壓器實際具有的電壓變換功能。但是等值變 壓器模型則具有這種電壓變換功能，它也是運用計算機(jī)進(jìn)行電力系統(tǒng)分析時采用 的變壓器模型，雖然運用這種模型時并不排斥手算。既然這種模型體現(xiàn)電壓變換， 在多電壓等級網(wǎng)絡(luò)計算中采用這種變壓器模型后，就可以不必進(jìn)行參數(shù)和變量的 歸算，這正是這種變壓器模型

23、的主要特點之一。以下，即介紹這種變壓器模型。首先，從一個未作電壓歸算的簡單網(wǎng)絡(luò)入手。設(shè)圖2-6、圖2-7中變壓器 的導(dǎo)納或勵磁支路和線路的導(dǎo)納支路都可略去；設(shè)變壓器兩側(cè)線路的阻抗都未經(jīng) 歸算，即分別為高低壓側(cè)或I、II側(cè)線路的實際阻抗，變壓器本身的阻抗歸在低 壓側(cè)；設(shè)變壓器的變比為k，其值為高、低壓繞組電壓之比。K:12I 廠*正圖2-6變壓器模型(1)圖2-7變壓器模型(2)顯然,在這些假設(shè)條件下,如在變壓器阻抗Z/側(cè)串聯(lián)一變比為K的理想 變壓器如圖2-8：凸U1. K1己2占圖2-8變壓器模型(3)其效果就如同將變壓器及其低壓側(cè)線路的阻抗都?xì)w算至高壓側(cè)，或?qū)⒏邏簜?cè) 線路的阻抗歸算至低壓側(cè)，

24、從而實際上獲得將所有參數(shù)和變量都?xì)w算到同一側(cè)的 等值網(wǎng)絡(luò)，只要變壓器的變比取的是實際變比，這一等值網(wǎng)絡(luò)無疑是嚴(yán)格的。因 此很容易知道圖2-5中的參數(shù)：七2 =七1 = I上.犬y10 = (1-k)/Z, - K2七。=(k -1)/Z K (23) 附帶指出，可以證明，變壓器不僅有改變電壓大小而且有移相功能時，其變 比k將為復(fù)數(shù)，這時，仍將得到上面所示的y、y、y、y，但其中y、y 102021121221不相等，無源電路的互易特性不復(fù)存在，不能用n形等值電路表示這種變壓器模 型，雖然這樣不影響運用這種模型進(jìn)行計算。2.2網(wǎng)絡(luò)模型有名制：所有參數(shù)和變量都以有名單位，如Q、S、kV(V)、kA

25、(A)、MVA(VA) 等表示。標(biāo)幺制：所有參數(shù)和變量都以與他們同名基準(zhǔn)值相對的標(biāo)幺值表示，因此都 沒有單位。對多電壓級網(wǎng)絡(luò)，變壓器模型：采用等值變壓器模型時，所有參數(shù)和變量可 不進(jìn)行歸算；采用有名制或標(biāo)幺制取決于習(xí)慣。在我國，電力工程界使用標(biāo)幺值 已有多年；但在國外，有名制的使用也很普遍。至于變壓器模型的使用范圍，則 涇渭分明。手算時，都是用r形或T型等值電路模型；計算機(jī)計算時，都是用 等值變壓器或n型等值電路模型。此外，在制定電力網(wǎng)絡(luò)等值電路模型時，有時還同時作某些簡化，常見的有： 線路的電導(dǎo)通常都被略去；變壓器的電導(dǎo)有時以具有定值的有功功率損耗的 形式出現(xiàn)在電路中；100km以下架空線路

26、的電納被略去；100300km架空線 路或變壓器的電納有時以具有定值的容性或感性無功功率損耗的形式出現(xiàn)在電 路中。有時，整個元件，甚至部分系統(tǒng)都可能不包括在等值電路中。例如，將某些發(fā)電廠的高壓母線看作為可維持給定電壓、輸出給定功率的等值電源時，這些 發(fā)電廠內(nèi)部的元件就不再包括在等值電路中。三：基于matlab的配電網(wǎng)潮流計算算法3.1配電網(wǎng)潮流計算算法原理牛頓一拉夫遜迭代法是常用的解非線性方程組的方法，也是當(dāng)前廣泛采用的 計算潮流的方法，其標(biāo)準(zhǔn)模式如下。設(shè)有非線性方程組f （x ,x ,x ）= y112 n 1f （x ,x ,x ）= y式（3.1）212 n式（3.1）f （x ,x ,

27、x ）= yn 12 n n其近似解為xi（其近似解為xi（0）, x2）,-,平的關(guān)式應(yīng)該成立：f C（0）+ Ax , x（ 0 + Ax，x（+Ax ）= y1122 nn、f x（o）+ Ax , x（ d + Ax ； x（+Ax）= y1122 nn設(shè)近似解與精確解分別相差,從則如下式（3.2）f (G e B f ) + G e + B f (j = i)jij j 司 Ji ，I jei(3-19)5A20 (j/i)-2e (j = i) jIc 6A72 f 0 (S)S =r = S .可2f. (j =，)(3-20)牛頓一拉夫遜法的程序框圖四：算例原始資料及主要參數(shù):

28、輻射配電網(wǎng)的線路參數(shù)和負(fù)載參數(shù)支路起點終點電阻（歐姆）電抗（歐姆）:W無功(KVA)1120.00580.002910013182230.03080.0157907503340.02280.01151204794450.02380.0121606605560.05110.0441606306670.01170.038620010497780.10680.077120014258890.06430.0462160110699100.06510.0462606021010110.01230.0041452981111120.02340.0077608951212130.09160.07216096

29、11313140.03380.04451205991414150.03690.0328606211515160.04660.0340609471616170.08040. 1074609361717180.04570.035890371182190.01020.0998909701919200.09390.08469010992020210.02550.0298906582121220.04420.052590431223230.02820.0192907762323240.05600.044242010532424250.05590.0437420673原始資料及主要參數(shù):22256260.

30、01270.0065604952626270.01770.0090607722727280.06610.0583609632828290.05020.04371206842929300.03170.01611007323030310.06080.06011509093131320.01940.02262104083232330.02130.0331100261參考文獻(xiàn)1王淳.一種實用的輻射網(wǎng)潮流算法.繼電器.2008年12期2王明崗.基于前推后代法解決環(huán)網(wǎng)的配電網(wǎng)潮流計算.南京工程學(xué)院學(xué)報（自然科學(xué)版）. 2009年01期3余加喜，郭志忠.一種面向?qū)ο蟮妮椛錉钆潆娋W(wǎng)潮流計算方法.哈爾濱理工大學(xué)

31、學(xué)報.2007 年06期4謝開貴，周平._配電網(wǎng)絡(luò)潮流計算的遞推算法.電力系統(tǒng)自動化.2007年04期5張榮，王秀和.改進(jìn)的帶二階項配電網(wǎng)快速潮流算法.電工技術(shù)學(xué)報.2007年07期6劉耀年，豈小梅.基于回路阻抗法的配電網(wǎng)潮流計算.繼電器.2007年08期7盧本初.陳紅坤.輻射配電網(wǎng)潮流計算實用方法.福建電力與電工.2007年04期8馬瑾.徐建源.輻射狀配電網(wǎng)的潮流計算.沈陽工業(yè)大學(xué)學(xué)報.2008年01期9咎加強.駱平.配電網(wǎng)的疊加法潮流計算.現(xiàn)代電力.2009年06期10顏偉.劉方.輻射型網(wǎng)絡(luò)潮流的分層前推回代算法.中國電機(jī)工程學(xué)報.2009年08期11汪衛(wèi)華.楊期余.配電網(wǎng)潮流計算中簡單編

32、號方法.電力自動化設(shè)備.2010年01期12曹亮.孔峰.一種配電網(wǎng)的實用潮流算法.繼電器.2012年05期13戴雯霞.吳捷.基于支路電流的配網(wǎng)潮流前推后代法.繼電器.2012年05期14許夢黔.趙武智.基于圖形化界面配電網(wǎng)潮流計算的開發(fā).電力標(biāo)準(zhǔn)化與計量.2010年 03期15蘇宏田.王秀玲.二叉樹理論在配電網(wǎng)潮流計算中的應(yīng)用.內(nèi)蒙古電力技術(shù).2010年02 期16 Ulas Eminoglu and M. Hakan Hocaoglu A new power flow method for radial distribution systems including voltage depen

33、dent load models Electric power system 2005.5衷心感謝導(dǎo)師劉會家老師在學(xué)術(shù)上的悉心指導(dǎo)、深深教誨，在生活上的熱情 關(guān)懷，人生路上給予我的指引。本人大學(xué)學(xué)習(xí)期間在劉老師教導(dǎo)下求學(xué)，做人， 劉老師淵博的學(xué)識、敏銳的思維、民主而嚴(yán)謹(jǐn)?shù)淖黠L(fēng)使學(xué)生受益匪淺、永生難忘; 追求真理、獻(xiàn)身科學(xué)、學(xué)而不厭、誨人不倦的崇高品質(zhì)對學(xué)生將是永遠(yuǎn)的鞭策， 他不但教會了我許多理論上的知識、實踐中的經(jīng)驗，同時我還從他那里學(xué)到了人 生的哲理、指引我今后的人生方向。感謝我的父母多年來對我的支持和鼓勵，比起他們的付出，我所做的一切是 那么的微不足道。最后向鼓勵和幫助的各位老師、同學(xué)致以

34、最誠摯的謝意和最衷心的祝福。附錄：源程序n=input(請輸入節(jié)點數(shù)：n=);n1=input(,請輸入支路數(shù)：n1=)；isb=input(,請輸入平衡母線節(jié)點號：isb=)；pr=input(,請輸入誤差精度：pr=)；B1二input(請輸入由支路參數(shù)形成的矩陣：B1=)；B2=input(,請輸入各節(jié)點參數(shù)形成的矩陣：B2=)；X=input(請輸入由節(jié)點參數(shù)形成的矩陣：X=)；Y=zeros(n)；e=zeros(1,n)；f=zeros(1,n)；V=seros(1,n)；O=zeros(1,n)；S1=zeros(n1)；for i=1:nif X(i,2)=0;P=X(i,1)

35、;Y(p,p)=1./X(i,2);endendfor i=1:n1if B1(i,6)=0p=B1(i,1);q=B1(i,2);else p=B1(i,2);q=B1(i,1);endY(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5);Y(p,q)=Y(p,q);Y(p,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)2)+B1(i,4)./2;Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;end%求導(dǎo)納矩陣G=real(Y);B=imag(Y);for i=1:ne(i)=real(B2(i,3);f(i)=imag(B2(i,3);V

36、(i)=B2(i,4);endfor i=1:nS(i)=B2(i,1)-B2(i,2);B(i,i)=B(i,i)+B2(i,5);endP=rea(S);Q=imag(S);ICT1=0;IT2=1;NO=2*n;N=NO+1;a=0;while IT2=0IT2=0;a=a+1;for i=1:n;if i=isbC(i) =0;D(i) =0;for j1=1:nC(i)=C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1);D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1);endP1=C(i)*e(i)+f(i)*D(i);Q1=f(i)*C(i

37、)-D(i)*e(i);%求P,QV2=e(i)2+f(i)2;if B2(i,6)=3DP=P(i)-P1;DQ=Q(i)-Q1;for j1=1:nif j1=isb&j1=iX1=-G(i,j1)*e(i)-B(i,j1)*f(i);X2=B(i,j1)*e(i)-G(i,j1)*f(i);X3=X2;X4=-X1;p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;m=p+1;J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X4;J(m,q)=X2;endendelseDPHPEP1;DVHV(i)2V2;for jlHlmif jl Hisbe

38、l HiX1PGPjl)*e(i)Bpjl)*f(i)； X2HBPjl)*e(i)Gpjl)*f(i)；X5H0;X6H0;- qH2* j T1 二pq) HX5 二pN) HDV - nlHP+1: JFq) Hxl 二FN) hdp - qHq+1 二pq) HX6 二Fq) Hx elseif jlHHi中jlHisbX1HC(i)Gpi)*e(i)Bpi)*f (i)： X2HDE+BPi)*e(i)Gpi)*f (i)： x5H2*e(i)：X6H2M (i)：- qH2* j T1 二pq) HX5 二pN) HDV WHP+1: JFq) Hxl 二(目 N) hdp - q

39、Hq+1 二pq) HX6 二(目 q) HX2 -endendendendend K沛回氏咨帛for kH3N0kIHk+rNIHN;for k2HklNlJFk2TJFk2)/JFk)：endJFkTl;if k3k4Hkl:for k3H3k432for k2=k1:N1J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endendfor k3=k1:N0for k2=k1:N1J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endendendfor k=3:2:N0-1L=(k+1)./2;e(L

40、)=e(L)-J(k,N);k1=k+1;f(L)=f(L)-J(k1,N);endfor k=3:N0DET=abs(J(k,N);if DET=prIT2=IT2+1endendICT2(a)=IT2ICT1=ICT1+1;for k=1:ndy(k)=sqrt(e(k)”2+f(k)”2);endfor i=1:nDy(k)=sqrt(e(k)”2+f(k)”2);endfor i=1:nDy(ICT1,i)=dy(i);endend%用高斯消去法解“ w=-J*V”disp(迭代次數(shù))；disp(ICT1);disp(沒有達(dá)到精度要求的個數(shù))；disp(ICT2);for k=1:nV

41、(k)二sqrt(e(k)”2+f(k)”2);0(k)=atan(f(k)./e(k)*180./pi;endE=e+f*j;disp(各節(jié)點的實際電壓標(biāo)么值E為(節(jié)點號從小到大的排列)：);disp(E);disp(各節(jié)點的電壓大小V為(節(jié)點號從小到大的排列)：)； disp(V);disp(各節(jié)點的電壓相角0為(節(jié)點號從小到大的排列)：)； disp(O);for p=1:nC(p)=0;for q=1:nC(p)=C(p)+conj(Y(p,q)*conj(E(q);endS(p)=E(p)*C(p);enddisp(各節(jié)點的功率S為(節(jié)點號從小到大排列)：)；disp(S);disp

42、(各條支路的首端功率Si為(順序同您輸AB1時一樣)：)；for i=1:n1if B1(i,6)=0p=B1(i,1);q=B1(i,2);else p=B1(i,2);q=B1(i,1);endSi(p,q)=E(p)*(conj(E(p)*conj(B1(i,4)./2) + (conj(E(p)*B1(i,5)-conj(E(q)*conj(1./(B1(i,3)*B1(i,5);disp(Si(p.q);enddisp(各條支路的末端功率Sj為(順序同您的輸入B1時一樣)：)；for i=1:n1if B1(i,6)=0p=B1(i,1);q=B1(i,2);else p=B1(i,

43、2);q=B1(i,1);endSj(q,p)=E(q)*(conj(E(q)*conj(B1(i,4)./2) + (xonj(E(q)./B1(i,5)-con j(E(p)*xonj(1./(B1(i,3)*B1(i,5);disp(Sj(q,p);enddisp(各條支路的功率損耗DS為(順序同您輸入B1時一樣)：:)；for i=1:n1if B1(i,6)=0p=B1(i,1);q=B1(i,2);else p=B1(i,2);q=B1(i,1);endDS(i)=Si(p,q)+Sj(q,p);disp(DS(i);endfor i=1:ICT1Cs(i)=i;enddisp(以

44、下是每次迭代后各節(jié)點的電壓值(如圖所示)； plot(Cs,Dy),xlabel(迭代次數(shù))，ylabel(電壓)，title(電壓迭代次數(shù)曲線 );畢業(yè)設(shè)計（論文）原創(chuàng)性聲明和使用授權(quán)說明原創(chuàng)性聲明本人鄭重承諾：所呈交的畢業(yè)設(shè)計（論文），是我個人在指導(dǎo)教 師的指導(dǎo)下進(jìn)行的研究工作及取得的成果。盡我所知，除文中特別加 以標(biāo)注和致謝的地方外，不包含其他人或組織已經(jīng)發(fā)表或公布過的研 究成果，也不包含我為獲得 及其它教育機(jī)構(gòu)的學(xué)位或?qū)W歷而使用過的材料。對本研究提供過幫助和做出過貢獻(xiàn)的個人或集體， 均已在文中作了明確的說明并表示了謝意。作者簽名： 日 期：指導(dǎo)教師簽名：日 期：使用授權(quán)說明本人完全了解

45、 大學(xué)關(guān)于收集、保存、使用畢業(yè)設(shè)計（論文）的規(guī)定，即：按照學(xué)校要求提交畢業(yè)設(shè)計（論文）的印刷本和電 子版本；學(xué)校有權(quán)保存畢業(yè)設(shè)計（論文）的印刷本和電子版，并提供 目錄檢索與閱覽服務(wù)；學(xué)校可以采用影印、縮印、數(shù)字化或其它復(fù)制 手段保存論文；在不以贏利為目的前提下，學(xué)校可以公布論文的部分 或全部內(nèi)容。作者簽名： 日 期：學(xué)位論文原創(chuàng)性聲明本人鄭重聲明：所呈交的論文是本人在導(dǎo)師的指導(dǎo)下獨立進(jìn)行研 究所取得的研究成果。除了文中特別加以標(biāo)注引用的內(nèi)容外，本論文 不包含任何其他個人或集體已經(jīng)發(fā)表或撰寫的成果作品。對本文的研 究做出重要貢獻(xiàn)的個人和集體，均已在文中以明確方式標(biāo)明。本人完 全意識到本聲明的法律

46、后果由本人承擔(dān)。作者簽名：日期： 年 月 日本學(xué)位論文作者完全了解學(xué)校有關(guān)保留、使用學(xué)位論文的規(guī)定， 同意學(xué)校保留并向國家有關(guān)部門或機(jī)構(gòu)送交論文的復(fù)印件和電子版， 允許論文被查閱和借閱。本人授權(quán) 大學(xué)可以將本學(xué)位論文的全部或部分內(nèi)容編入有關(guān)數(shù)據(jù)庫進(jìn)行檢索，可以采用影印、縮 印或掃描等復(fù)制手段保存和匯編本學(xué)位論文。涉密論文按學(xué)校規(guī)定處理。作者簽名：日期：年月日導(dǎo)師簽名：日期：年月日獨創(chuàng)聲明本人鄭重聲明：所呈交的畢業(yè)設(shè)計（論文），是本人在指導(dǎo)老師的 指導(dǎo)下，獨立進(jìn)行研究工作所取得的成果，成果不存在知識產(chǎn)權(quán)爭議。 盡我所知，除文中已經(jīng)注明引用的內(nèi)容外，本設(shè)計（論文）不含任何 其他個人或集體已經(jīng)發(fā)表或

47、撰寫過的作品成果。對本文的研究做出重 要貢獻(xiàn)的個人和集體均已在文中以明確方式標(biāo)明。本聲明的法律后果由本人承擔(dān)。作者簽名:二00年九月二十日畢業(yè)設(shè)計（論文）使用授權(quán)聲明本人完全了解*學(xué)院關(guān)于收集、保存、使用畢業(yè)設(shè)計（論文）的 規(guī)定。本人愿意按照學(xué)校要求提交學(xué)位論文的印刷本和電子版，同意學(xué) 校保存學(xué)位論文的印刷本和電子版，或采用影印、數(shù)字化或其它復(fù)制 手段保存設(shè)計（論文）；同意學(xué)校在不以營利為目的的前提下，建立 目錄檢索與閱覽服務(wù)系統(tǒng)，公布設(shè)計（論文）的部分或全部內(nèi)容，允 許他人依法合理使用。（保密論文在解密后遵守此規(guī)定）作者簽名:二00年九月二十日基本要求：寫畢業(yè)論文主要目的是培養(yǎng)學(xué)生綜合運用所

48、學(xué)知識 和技能，理論聯(lián)系實際，獨立分析，解決實際問題的能力，使學(xué)生得 到從事本專業(yè)工作和進(jìn)行相關(guān)的基本訓(xùn)練。畢業(yè)論文應(yīng)反映出作者能 夠準(zhǔn)確地掌握所學(xué)的專業(yè)基礎(chǔ)知識，基本學(xué)會綜合運用所學(xué)知識進(jìn)行 科學(xué)研究的方法，對所研究的題目有一定的心得體會，論文題目的范 圍不宜過寬，一般選擇本學(xué)科某一重要問題的一個側(cè)面。畢業(yè)論文的基本教學(xué)要求是：1、培養(yǎng)學(xué)生綜合運用、鞏固與擴(kuò)展所學(xué)的基礎(chǔ)理論和專業(yè)知識， 培養(yǎng)學(xué)生獨立分析、解決實際問題能力、培養(yǎng)學(xué)生處理數(shù)據(jù)和信息的 能力。2、培養(yǎng)學(xué)生正確的理論聯(lián)系實際的工作作風(fēng)，嚴(yán)肅認(rèn)真的科 學(xué)態(tài)度。3、培養(yǎng)學(xué)生進(jìn)行社會調(diào)查研究；文獻(xiàn)資料收集、閱讀和整 理、使用；提出論點、綜

49、合論證、總結(jié)寫作等基本技能。畢業(yè)論文是畢業(yè)生總結(jié)性的獨立作業(yè)，是學(xué)生運用在校學(xué)習(xí)的基 本知識和基礎(chǔ)理論，去分析、解決一兩個實際問題的實踐鍛煉過程， 也是學(xué)生在校學(xué)習(xí)期間學(xué)習(xí)成果的綜合性總結(jié)，是整個教學(xué)活動中不 可缺少的重要環(huán)節(jié)。撰寫畢業(yè)論文對于培養(yǎng)學(xué)生初步的科學(xué)研究能 力，提高其綜合運用所學(xué)知識分析問題、解決問題能力有著重要意義。 畢業(yè)論文在進(jìn)行編寫的過程中，需要經(jīng)過開題報告、論文編寫、論文 上交評定、論文答辯以及論文評分五個過程，其中開題報告是論文進(jìn) 行的最重要的一個過程，也是論文能否進(jìn)行的一個重要指標(biāo)。撰寫意義：1.撰寫畢業(yè)論文是檢驗學(xué)生在校學(xué)習(xí)成果的重要措 施，也是提高教學(xué)質(zhì)量的重要環(huán)節(jié)

50、。大學(xué)生在畢業(yè)前都必須完成畢業(yè) 論文的撰寫任務(wù)。申請學(xué)位必須提交相應(yīng)的學(xué)位論文，經(jīng)答辯通過后， 方可取得學(xué)位?？梢赃@么說，畢業(yè)論文是結(jié)束大學(xué)學(xué)習(xí)生活走向社會 的一個中介和橋梁。畢業(yè)論文是大學(xué)生才華的第一次顯露，是向祖國 和人民所交的一份有份量的答卷，是投身社會主義現(xiàn)代化建設(shè)事業(yè)的 報到書。一篇畢業(yè)論文雖然不能全面地反映出一個人的才華，也不一 定能對社會直接帶來巨大的效益，對專業(yè)產(chǎn)生開拓性的影響。但是， 實踐證明，撰寫畢業(yè)論文是提高教學(xué)質(zhì)量的重要環(huán)節(jié)，是保證出好人 才的重要措施。通過撰寫畢業(yè)論文，提高寫作水平是干部隊伍“四化”建設(shè)的需 要。黨中央要求，為了適應(yīng)現(xiàn)代化建設(shè)的需要，領(lǐng)導(dǎo)班子成員應(yīng)當(dāng)逐

51、 步實現(xiàn)“革命化、年輕化、知識化、專業(yè)化”。這個“四化”的要求，也 包含了對干部寫作能力和寫作水平的要求。提高大學(xué)生的寫作水平是社會主義物質(zhì)文明和精神文明建設(shè) 的需要。在新的歷史時期，無論是提高全族的科學(xué)文化水平，掌握現(xiàn) 代科技知識和科學(xué)管理方法，還是培養(yǎng)社會主義新人，都要求我們的 干部具有較高的寫作能力。在經(jīng)濟(jì)建設(shè)中，作為領(lǐng)導(dǎo)人員和機(jī)關(guān)的辦 事人員，要寫指示、通知、總結(jié)、調(diào)查報告等應(yīng)用文；要寫說明書、 廣告、解說詞等說明文；還要寫科學(xué)論文、經(jīng)濟(jì)評論等議論文。在當(dāng) 今信息社會中，信息對于加快經(jīng)濟(jì)發(fā)展速度，取得良好的經(jīng)濟(jì)效益發(fā) 揮著愈來愈大的作用。寫作是以語言文字為信號，是傳達(dá)信息的方式。 信息

52、的來源、信息的收集、信息的儲存、整理、傳播等等都離不開寫 作。論文種類：畢業(yè)論文是學(xué)術(shù)論文的一種形式，為了進(jìn)一步探討和 掌握畢業(yè)論文的寫作規(guī)律和特點，需要對畢業(yè)論文進(jìn)行分類。由于畢 業(yè)論文本身的內(nèi)容和性質(zhì)不同，研究領(lǐng)域、對象、方法、表現(xiàn)方式不 同，因此，畢業(yè)論文就有不同的分類方法。按內(nèi)容性質(zhì)和研究方法的不同可以把畢業(yè)論文分為理論性論文、 實驗性論文、描述性論文和設(shè)計性論文。后三種論文主要是理工科大 學(xué)生可以選擇的論文形式，這里不作介紹。文科大學(xué)生一般寫的是理 論性論文。理論性論文具體又可分成兩種：一種是以純粹的抽象理論 為研究對象，研究方法是嚴(yán)密的理論推導(dǎo)和數(shù)學(xué)運算，有的也涉及實 驗與觀測，用

53、以驗證論點的正確性。另一種是以對客觀事物和現(xiàn)象的 調(diào)查、考察所得觀測資料以及有關(guān)文獻(xiàn)資料數(shù)據(jù)為研究對象，研究方 法是對有關(guān)資料進(jìn)行分析、綜合、概括、抽象，通過歸納、演繹、類 比，提出某種新的理論和新的見解。按議論的性質(zhì)不同可以把畢業(yè)論文分為立論文和駁論文。立論性 的畢業(yè)論文是指從正面闡述論證自己的觀點和主張。一篇論文側(cè)重于 以立論為主，就屬于立論性論文。立論文要求論點鮮明，論據(jù)充分， 論證嚴(yán)密，以理和事實服人。駁論性畢業(yè)論文是指通過反駁別人的論 點來樹立自己的論點和主張。如果畢業(yè)論文側(cè)重于以駁論為主，批駁 某些錯誤的觀點、見解、理論，就屬于駁論性畢業(yè)論文。駁論文除按 立論文對論點、論據(jù)、論證的

54、要求以外，還要求針鋒相對，據(jù)理力爭。 按研究問題的大小不同可以把畢業(yè)論文分為宏觀論文和微觀論文。凡 屆國家全局性、帶有普遍性并對局部工作有一定指導(dǎo)意義的論文，稱 為宏觀論文。它研究的面比較寬廣，具有較大范圍的影響。反之，研 究局部性、具體問題的論文，是微觀論文。它對具體工作有指導(dǎo)意義， 影響的面窄一些。另外還有一種綜合型的分類方法，即把畢業(yè)論文分為專題型、論 辯型、綜述型和綜合型四大類：1.專題型論文。這是分析前人研究成果的基礎(chǔ)上，以直接論述 的形式發(fā)表見解，從正面提出某學(xué)科中某一學(xué)術(shù)問題的一種論文。如 本書第十二章例文中的淺析領(lǐng)導(dǎo)者突出工作重點的方法與藝術(shù)一 文，從正面論述了突出重點的工作方

55、法的意義、方法和原則，它表明 了作者對突出工作重點方法的肯定和理解。2.論辯型論文。這是針 對他人在某學(xué)科中某一學(xué)術(shù)問題的見解，憑借充分的論據(jù)，著重揭露 其不足或錯誤之處，通過論辯形式來發(fā)表見解的一種論文。3.綜述 型論文。這是在歸納、總結(jié)前人或今人對某學(xué)科中某一學(xué)術(shù)問題已有 研究成果的基礎(chǔ)上，加以介紹或評論，從而發(fā)表自己見解的一種論文。綜合型論文。這是一種將綜述型和論辯型兩種形式有機(jī)結(jié)合起來 寫成的一種論文。如關(guān)于中國民族關(guān)系史上的幾個問題一文既介 紹了研究民族關(guān)系史的現(xiàn)狀，又提出了幾個值得研究的問題。因此， 它是一篇綜合型的論文。寫作步驟：畢業(yè)論文是高等教育自學(xué)考試本科專業(yè)應(yīng)考者完成本 科階段學(xué)業(yè)的最后一個環(huán)節(jié)，它是應(yīng)考者的總結(jié) 性獨立作業(yè)，目的 在于總結(jié)學(xué)習(xí)專業(yè)的成果，培養(yǎng)綜合運用所學(xué)知識解決實際問題 的 能力。從文體而言，它也是對某一專業(yè)領(lǐng)域的現(xiàn)實問題或理論問題 進(jìn)行 科學(xué)研究探索的具有一定意義的論說文。完成畢業(yè)論文的撰 寫