電力系統(tǒng)分析課程設(shè)計

1、精品資料推薦W過由3滎火厚力系統(tǒng)分析課程設(shè)計學(xué)院電氣工程學(xué)院班級學(xué)號姓名指導(dǎo)教師時間Z、.前言電力系統(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)濟性。此外，電力系統(tǒng)潮流計算也是計算系統(tǒng)動態(tài)穩(wěn)定和靜態(tài)穩(wěn)定的基礎(chǔ)。所以潮流計算是研究電力系統(tǒng)的一種很重要和很基礎(chǔ)的計算。隨著科學(xué)技術(shù)的發(fā)展，電力系統(tǒng)變得越來越復(fù)雜，電氣工程師掌握一種好的能對電力系統(tǒng)進行仿真

2、的軟件是學(xué)習(xí)和研究的需要。與眾多專門的電力系統(tǒng)仿真軟件相比，MATLA瞅件具有易學(xué)、功能強大和開放性好，是電力系統(tǒng)仿真研究的有力工具。49目錄1設(shè)計題目31.1 系統(tǒng)圖的確定31.2 各節(jié)點的初值及阻抗參數(shù)42潮流計算52.1 潮流計算概述與發(fā)展52.2 復(fù)雜電力系統(tǒng)潮流計算52.3 MATLAB概述62.4 牛頓-拉夫遜法原理72.5 牛頓-拉夫遜法解決潮流計算問題82.6 計算機潮流計算的步驟92.7 計算機潮流計算流程圖113手算潮流計算123.1 確定節(jié)點類型123.2 求節(jié)點導(dǎo)納矩陣Yb123.3 計算各節(jié)點功率的修正方程的初始初（不平衡量）134計算機算法潮流計算154.1 計算機

3、法機流計算過程154.2 計算機法機流計算結(jié)果234.3 系統(tǒng)功率分布圖25個人心得29附錄：源程序30參考文獻401設(shè)計題目1.1系統(tǒng)圖的確定選擇六節(jié)點、環(huán)網(wǎng)、兩電源和多引出的電力系統(tǒng)，簡化電力系統(tǒng)圖如圖1-1所示，等值導(dǎo)納圖如圖1-2所示。運用以直角坐標表示的牛頓-拉夫遜計算如圖1-1所示系統(tǒng)中的潮流分布。計算精度要求各節(jié)點電壓的誤差或修正量不大于105O2+j13.7+J1.31.8+j4.01.6+J0.8Hkj0.25p25Jj0.25j0.25二兩.25市0.252+j10.04+j0.251.05:11.05:13.7+j1.3J Lj 0. 03U11 . 050o圖1-2電力

4、系統(tǒng)等值導(dǎo)納圖該系統(tǒng)中，節(jié)點為平衡節(jié)點,1.2各節(jié)點的初值及阻抗參數(shù)保持U1=1.05+j0為定值，節(jié)點為PV節(jié)點,其他四個節(jié)點都是PQ節(jié)點。給定的注入電壓標幺值如表1-1、線路阻抗標幺值如表1-2、輸出功率標幺值如表1-3。表1-1各節(jié)點電壓標幺值參數(shù)U1U2U3U4U5U61.051.001.001.001.001.05表1-2線路、變壓器阻抗標幺值線路L2L3L4L5T1T2Y/2阻抗0.06+j0.250.04+j0.250.08+j0.300.1+j0.35j0.03j0.015j0.25表1-3節(jié)點輸出功率節(jié)點功率2+j11.8+j0.41.6+j0.83.7+j1.35注：各PQ

5、節(jié)點的電壓取1是為了方便計算和最后驗證程序的正確性。2潮流計算2.1 潮流計算概述與發(fā)展電力系統(tǒng)潮流計算也分為離線計算和在線計算兩種，前者主要用于系統(tǒng)規(guī)劃設(shè)計和安排系統(tǒng)的運行方式，后者則用于正在運行系統(tǒng)的經(jīng)常監(jiān)視及實時控制。利用電子數(shù)字計算機進行電力系統(tǒng)潮流計算從50年代中期就已經(jīng)開始。在這20年內(nèi)，潮流計算曾采用了各種不同的方法，這些方法的發(fā)展主要圍繞著對潮流計算的一些基本要求進行的。牛頓-拉夫遜法作為一種實用的，有競爭力的電力系統(tǒng)潮流計算方法，是在應(yīng)用了稀疏矩陣技巧和高斯消去法求修正方程后。牛頓-拉夫遜法是求解非線性代數(shù)方程有效的迭代計算。2.2 復(fù)雜電力系統(tǒng)潮流計算電力系統(tǒng)潮流計算是對復(fù)

6、雜電力系統(tǒng)正常和故障條件下穩(wěn)態(tài)運行狀態(tài)的計算。潮流計算的目標是求取電力系統(tǒng)在給定運行方式下的節(jié)點電壓和功率分布，用以檢查系統(tǒng)各元件是否過負荷、各點電壓是否滿足要求、功率的分布和分配是否合理以及功率損耗等。對現(xiàn)有電力系統(tǒng)的運行和擴建，對新的電力系統(tǒng)進行規(guī)劃設(shè)計都是以潮流計算為基礎(chǔ)。潮流計算結(jié)果的用途，例如用于電力系統(tǒng)穩(wěn)定研究、安全估計或最優(yōu)潮流等也對潮流計算的模型和方法有直接影響。節(jié)點類型：(1) PV節(jié)點：柱入有功功率P為給定值，電壓也保持在給定數(shù)值。(2) PQ節(jié)點：諸如有功功率和無功功率是給定的。（ 3）平衡節(jié)點：用來平衡全電網(wǎng)的功率。選一容量足夠大的發(fā)電機擔(dān)任平衡全電網(wǎng)功率的職責(zé)。平衡節(jié)

7、點的電壓大小與相位是給定的，通常以它的相角為參考量，即取其電壓相角為0。一個獨立的電力網(wǎng)中只設(shè)一個平衡點?；静襟E：（ 1）形成節(jié)點導(dǎo)納矩陣；（ 2）將各節(jié)點電壓設(shè)初值U；（ 3）將節(jié)點初值代入相關(guān)求式，求出修正方程式的常數(shù)項向量；（ 4）將節(jié)點電壓初值代入求式，求出雅可比矩陣元素；（ 5）求解修正方程，求修正向量；（ 6）求取節(jié)點電壓的新值；（ 7）檢查是否收斂，如不收斂，則以各節(jié)點電壓的新值作為初值自第3步重新開始進行狹義次迭代，否則轉(zhuǎn)入下一步；（ 8）計算支路功率分布，PV節(jié)點無功功率和平衡節(jié)點柱入功率。2.3 MATLAB概述目前電子計算機已廣泛應(yīng)用于電力系統(tǒng)的分析計算，潮流計算是其基

8、本應(yīng)用之一?，F(xiàn)有很多潮流計算方法。對潮流計算方法有五方面的要求：（ 1）計算速度快；（ 2）內(nèi)存需要少；（ 3）計算結(jié)果有良好的可靠性和可信性；（ 4）適應(yīng)性好，亦即能處理變壓器變比調(diào)整、系統(tǒng)元件的不同描述和與其它程序配合的能力強；（ 5）簡單。MATLAB是一種交互式、面向?qū)ο蟮某绦蛟O(shè)計語言，廣泛應(yīng)用于工業(yè)界與學(xué)術(shù)界，主要用于矩陣運算，同時在數(shù)值分析、自動控制模擬、數(shù)字信號處理、動態(tài)分析、繪圖等方面也具有強大的功能。MATLAB程序設(shè)計語言結(jié)構(gòu)完整，且具有優(yōu)良的移植性，它的基本數(shù)據(jù)元素是不需要定義的數(shù)組。它可以高效率地解決工業(yè)計算問題，特別是關(guān)于矩陣和矢量的計算。MATLAB!C語言和FOR

9、TRAN言相比更容易被掌握。通過M語言,可以用類似數(shù)學(xué)公式的方式來編寫算法，大大降低了程序所需的難度并節(jié)省了時問，從而可把主要的精力集中在算法的構(gòu)思而不是編程上。2.4 牛頓-拉夫遜法原理假設(shè)有n個聯(lián)立的非線性代數(shù)方程:fl(x1,x2xn)0f2(x1,x2xn)0Y(0),xn ,令其分別為個變量的修正fn(x1,x2xn)0(0)(0)假設(shè)以給出各變量的初值x1,x2,量，使?jié)M足以上方程，所以：9X2o1X1X7(1X/V2fn(xix?x20)x20),xn0)xn0)0(01X將上式中的n個多元函數(shù)在初始值附近分別展開成泰勒級數(shù)，并略去含有xn的二次及以上階次的各項，使得:fi(xi

10、(0),x20),xn0)fixi0xi(0)fix20x20).fix310x30)0f2(xi(0),x20),xn0)上0xi(0)f20x20).f20x3(0)0xix2x3fn(x1(0),x20),xn0)fn0xi(0)fn0x20).fn0x3(0)0xix2x3方程可寫成:(0)(0)(0)、3(x1,x2,xn)(0)(0)(0)、f2(x1,x2,xn)fn(xr,x20),xn0)f1f1f1000x1x2xnf2f2f2_*000x1x2xnfnfnfn000x1x2xn(o)Xi(0)X2x(0)xnx20)以上方程是對于修正量Xi（0）xn0)的線性方程組，稱為

11、牛頓法的修正方程，可解出x1（0）x20）,x,）。對初始近似解進行修正:xix(0)x(0)(i=1,2,，n)反復(fù)迭代，在進行k+1次迭代時，從求解修正方程式:f1(x1(k),x2k),.xnk)f2(x1(k),x2k),.xnk)fn(x1(k),x2k),.xnk)f1x1f2x1fnx1f1kx2f2kx2fnkx2f1k一xnf2k一xf,kxn得到修正量(k)x1(k)x2Jk）xn，對各量進行修正xi(k1)xi(k)xi(k)(i=1,2,）迭代過程一直進行到ma刈fi(xi(k),x2k)，xnk)i12.5牛頓-拉夫遜法解決潮流計算問題節(jié)點總數(shù)為n；PQ節(jié)點有m,;P

12、V節(jié)點有n-m-1,平衡節(jié)點有1個，節(jié)點編號按照先PQ節(jié)點，再PV節(jié)點，最后平衡節(jié)點的順序進行編號，即：1,2,，m為PQ節(jié)點；m+1,m+Z，n-1為PV節(jié)點;n為平衡節(jié)點?？尚纬山Y(jié)點導(dǎo)納矩陣。導(dǎo)納矩陣元素可表示為YjgjjBj,本文中節(jié)點電壓以直角坐標形式表示，即0eijfio由此下列公式可求出Pi，QiPiQinei(Gjejj1nfi(Gjejj1Bjfj)fi(GjfjBjej)j1nBjfj)ei(GjfjBjej)j1假設(shè)系統(tǒng)中的第1,2,，m號節(jié)點為PQ節(jié)點，第i個節(jié)點的給定功率為P和Qis，：對該節(jié)點可列方程:PinPisPiPisei(Gijejj1nBifj)QiQisQ

13、iQisfi(GijejBifj)1nfi(Gifjj1nei(Gifjj1Bijej)(i1,2,m)Bijej)假設(shè)系統(tǒng)中的第m+1，m+Z，n-1號節(jié)點為PV節(jié)點，則對其中每一個節(jié)點可列方程:PiPisPiPisVi2nei(GijejBi fj)j1Vis2Vi2Vis2(ei2nfi(Gi fjj1fi2)Bij ej ) (i m 1,m 2, ,n 1)第 n 號節(jié)點為平衡節(jié)點，其電壓為是給定的Vnjf n ，故不參加迭代。修正方程 W通過反復(fù)求解修正方程，W1J11J12J1,n1V1W2J21J22J2,n1V2W3J31J32J3,n1V3Wn1Jn1,1Jn1,2Jn1,

14、n1Vn1JV可寫成分塊矩陣的形式：解出各節(jié)點的未知量，再通過收斂判據(jù)判定是否已為真值。從而求得PQ節(jié)點的電壓V及相角6的真值，PV節(jié)點的Q6真值，平衡節(jié)點的P、Q真值，以上即為牛頓-拉夫遜迭代法的潮流計算過程，其優(yōu)點為計算精確，運行速度快。其中的各個環(huán)節(jié)都可通過MATLAB?序來實現(xiàn)。2.6計算機潮流計算的步驟1)對電力網(wǎng)絡(luò)的所有參數(shù)設(shè)初值,包括電壓、相角、有功、無功等。2)處理非標準變比支路,使其變成標準變比為1的變壓器支路。3)形成節(jié)點導(dǎo)納矩陣Y。(4)計算有功功率的不平衡量APi,從而求出Pi/Ui(i1,2,3,nJs)5)根據(jù)節(jié)點的類型形成J。(6)解修正方程式，求各節(jié)點的電壓的變

15、化量Aei(i=1,2,3.n,iws)(7)求各節(jié)點相角的新值ei=ei+Aei(i=1,2,3n,iws)(8)計算無功功率的不平衡量AQi,從而求出Qi/Ui(i=1,2,3n,iws)(9)解修正方程式,求各節(jié)點的電壓大小的變化量Ui(i=1,2,3.,n,i5) 。(10)求各節(jié)點的電壓大小的新值UiUiUi(i=1,2,3,n,iws)。(11)運用個節(jié)點的電壓的新值自第四步開始下一次迭代。計算平衡節(jié)點的功率和線路功率。其中平衡節(jié)點的功率的計算公式為in.SSUSYUiPiJQii1si線路上的功率為:.SijUiIijUiUiyi0(UiUj)yijPijJQij.SjiUjIj

16、iUjUiyj0(UjUi)yjiPjiJQji從而線路上的損耗的功率為:SijSijSjiPijJQij2.7計算機潮流計算流程圖選用牛頓-拉夫遜法進行計算機潮流計算。牛頓-拉夫遜法潮流計算流程圖如圖2-1。3手算潮流計算3.1 確定節(jié)點類型節(jié)點為平衡節(jié)點節(jié)點為PV節(jié)點，其他四個節(jié)點都是PQ節(jié)點3.2 求節(jié)點導(dǎo)納矩陣Yb由圖1可知，該系統(tǒng)以串聯(lián)支路的導(dǎo)納標幺值和對地并聯(lián)導(dǎo)納標幺值得等值電路如圖2所示?？傻脠D2的等值節(jié)點導(dǎo)納矩陣。對角線上的元素為：y33y44y55y66yiiiy22y 30K21YT1YT1j33.33y 20y30y23 y251.53-j37.42y,40y23y341

17、.74j6.39y40 y,50y34y451.58YT1'2- YT2y 20 y50K K1 K -1Yt1- Yt2j66.67KKj5.25y 25 y451.37 j66.51非對角線上的元素為:y510； y16y610yi2y2ij31.75乂3y310;y14y，10;y15y23y320.91j3.78y24y420；y26y620y25y520.62j3.90y34y430.83j3.11y35y530；y36y630y45y540.75j2.64y46y640y56y65j63.49所以節(jié)點導(dǎo)納矩陣為:Y=-j33.33j31.750000j31.751.53j3

18、7.420.91j3.7800.62j3.90000.91j3.781.74j6.390.83j3.1100000.83j3.111.58j5.250.75j2.64000.62j3.9000.75j2.641.37j6651j63.490000j6349j66.673.3 計算各節(jié)點功率的修正方程的初始值（不平衡量）設(shè)各節(jié)點電壓初始值為：?U11.05j0=L00+/0.=L00+加武=LOO+jQUT=iM+只猶=105+/）根據(jù)上述N-R法的求解過程，將各節(jié)點電壓的初始值代入（1）式和（2）式,進行潮流計算，得節(jié)點功率和節(jié)點電壓：P22;Q22.61;P31.8;Q30.10;P41.6

19、;Q40.30;P53.7;Q55.40;P65.0;U620.00;雅可比矩陣：-1.53-41.020.913.78000.623.9000-33.811.533.78-0.91003.90-0.64000.913.78-1.746.890.833.1100003.780.915.891.743.110.830000000.833.111.585.750.752.6400J=003.110.834.751.582.640.75000.623.90000.752.641.3873.21063.493.900.62002.640.7559.811.3863.490000000066.67063

20、.49000000002.10修正各節(jié)點電壓：e2(0)0.095e2(1)1.095f2(0)0.123f2(1)0.123e3(0)0.061e3(1)0.939f(30)0.604f3(1)0.6403e4(0)0.131e4(1)0.869f(40)0.638f4(1)0.638e5(0)0.096e5(1)1.096f(50)0.159f5(1)0.159e6(0)0e6(1)1f(60)0.88f6(1)0.884計算機算法潮流計算4.1計算機法潮流計算過程導(dǎo)納矩陣Y=0-30.2343i0+31.7460i0000+31.7460i14.8252-42.6506i-14.2012

21、+5.9172i0-0.62403.9002i00-0.8299+3.1120i000-0.829900-0.6240+3.9002i0+63.4921i0000-60.4686i-14.2012+5.9172i15.0311-8.5292i0+3.1120i1.5846-5.5035i-0.75470-0.7547+2.6415i1.3787-72.9583i0+63.4921i0+2.6415i功率方程第(1)次差值：Columns1through1204.2619-2.10000.1000-1.8000-0.5500-1.60008.4738-3.700005.0000形成的第(1)次J

22、acobi矩陣：Columns1through1100000000000003.90020000.624003.1120-0.829900.82993.11200-5.25351.58460-1.5846-5.75350000000000000-37.3887-0.624000-14.82523.9002005.9172000014.2012000002.6415-0.75470000.75472.6415003.900214.8252-47.9126-14.20125.917200-0.62405.917214.2012-8.0292-15.03113.11200.829900-14.201

23、25.917215.0311-9.0292-0.82993.112002.6415-0.7547-63.18451.378763.4921000.62403.90020.75472.6415-1.3787-82.732100000000000-2.400000000000076.19050Columns12through13000004.26190-2.100000.10000-1.80000-0.55000-1.600008.473863.4921-3.700000-63.49215.0000Jacobi矩陣第(1)次回代運算Columns1through120.1103-0.00240.1

24、657-0.043100-0.14900.10150.3342-0.14350.02970各個節(jié)點電壓模01.15431.01600.95861.14391.2008功率方程第(2)次差值：Columns1through12-0.2847-0.012700-0.9196-0.3298-0.56520.0384-1.51850.1960-0.0019-0.3928形成的第(2)次Jacobi矩陣：Columns1through110000000000000000000-47.156723.256800 4.4126000 1.147002.9819 -1.3474 01.3474 2.9819

25、0-3.5899 4.7795 0-1.7466 -5.2407 0-1.147000-20.21704.4126003.5786000015.2150000002.1211-1.56090001.56092.1211004.4413-47.5883-15.21503.578600-0.82942.9982-0.9415-83.48957.152972.6032000.82944.44130.94152.9982-0.3333-83.28441.88540000000000-2.400000000000-2.737576.1905-1.88545.233316.9700-5.6605-14.7

26、6232.51881.785700-16.97005.233318.1967-6.4582-1.78572.518800Columns12through1300000-0.91960-0.32980-0.28470-0.01270-0.565200.0384-1.8854-1.518572.60320.1960-0.0862-0.0019-72.6032-0.3928Jacobi矩陣第(2)次回代運算Columns1through12000.04410.01550.07000.1790-0.02750.02880.0422-0.0010.0487-0.0210各個節(jié)點電壓模01.11200.9

27、4350.78211.11701.2010功率方程第(3)次差值：Columns1through120-0.04170.0017-0.0185-0.0126-0.1340-0.0217-0.03670.0646-0.0024-0.0916形成的第(3)次Jacobi矩陣：Columns1through1100000000000004.23090001.179802.7814-1.224001.22402.78140-1.89914.35040-1.3054-5.0484000000000000-44.1917-1.179800-19.96214.2309003.4621000014.09630

28、00001.6691-1.35310001.35311.6691004.302623.5230-46.3315-14.09633.462100-0.97622.8902-1.0313-79.882410.076070.775100.97624.3026004.752616.4350-5.0872-13.43791.98451.551500-16.43504.752617.0579-6.4674-1.55151.9845001.03132.8902-3.4952-82.57374.56770000000000-2.4019000000000.353576.2524-4.5677Columns12

29、through1300000-0.041700.00170-0.01850-0.01260-0.13400-0.0217-4.5677-0.036770.77510.06460.0111-0.0024-70.7751-0.0916Jacobi矩陣第(3)次回代運算Columns1through12000.00780.00430.0768-0.00550.00560.01560.00110.0181各個節(jié)點電壓模01.10480.92370.70981.11261.20010.0191-0.0063功率方程第(4)次差值：Columns1through1200-0.00040.0002-0.02

30、76-0.00580.00060.0016-0.0003-0.0053形成的第(4)次Jacobi矩陣：Columns1through1100000000000004.19800001.191802.7273-1.188501.18852.72730-1.12174.24800-1.1039-4.9984000000000000-43.7449-1.191800-20.03024.1980003.4389000013.7881000001.4705-1.28060001.28061.4705004.271223.5927-45.9892-13.78813.438900-1.03372.8637

31、-1.0684-79.373711.143870.4222001.03374.2712-0.0002004.645516.3504-4.9899-13.04671.75021.470600-0.0012-16.35044.645516.7680-6.4292-1.47061.7502001.06842.8637-4.6945-82.22825.56080000000000-2.399800000001.503176.1844-5.5608Columns12through1300000-0.000400.00020-0.00020-0.00120-0.02760-0.0058-5.56080.0

32、00670.42220.00160.0473-0.0003-70.4222-0.0053Jacobi矩陣第(4)次回代運算Columns1through120.00210.00140.0056-0.00200.0259-0.00180.00160.00480.00020.0053各個節(jié)點電壓模0 1.10290.91780.6857 1.11141.2000功率方程第(5)次差值：Columns1through12-0.00000.00000.0000-0.0001-0.0032-0.0007-0.00000.0001-0.0000-0.0005形成的第(5)次Jacobi矩陣：Columns

33、1through11000000000000000000000-43.630023.62104.6137-16.3286004.1889-1.195800001.195802.7114-1.177601.17762.71140-0.85644.21620-1.0344-4.984300-20.05614.1889003.4337000013.6963000001.4034-1.25640001.25641.4034004.2621-45.8841-13.69633.433700-1.05172.8559-1.0801-79.237411.491170.3226001.05174.26211.0

34、8012.8559-5.0505-82.12235.86870000000000-2.3993000000001.841976.1690-5.8687Columns12through1316.3286-4.9688-12.93121.67111.4437004.613716.6803-6.4137-1.44371.67110000000-0.000000.000000.00000-0.00010-0.00320-0.0007-5.8687-0.000070.32260.00010.0580-0.0000-70.3226-0.0005Jacobi矩陣第(5)次回代運算Columns1throug

35、h12000.00030.00020.00080.0039-0.00030.00020.00070.00000.0008各個節(jié)點電壓模01.10260.91690.68201.11131.2000-0.0003功率方程第(6)次差值：1.0e-004*Columns1through120-0.00260.00070.0032-0.0182-0.7419-0.1647-0.00160.0184-0.0063-0.1036形成的第(6)次Jacobi矩陣：Columns1through1100000000000004.18760001.196402.7090-1.176001.17602.7090

36、0-0.81614.21140-1.0238-4.982200000000000000-43.6128-1.196400-20.05994.1876003.4330000013.6825000001.3932-1.25280001.25281.3932004.260723.6253-45.8683-13.68253.433000-1.05432.8548-1.0818-79.217011.543170.3078001.05434.26071.08182.8548-5.1036-82.10675.91480000000000-2.3993000000001.892476.1670-5.91484

37、.608916.3254-4.9658-12.91391.65911.439600-16.32544.608916.6670-6.4113-1.43961.659100Columns12through1300000-0.000000.000000.00000-0.00000-0.00010-0.0000-5.9148-0.000070.30780.00000.0596-0.0000-70.3078-0.0000Jacobi矩陣第（6）次回代運算1.0e-004*Columns1through120.2048-0.0728000.07640.04960.9503-0.06410.05640.17

38、530.00740.1923各個節(jié)點電壓模01.10260.91690.68191.11131.2000功率方程第（7）次差值：1.0e-007*Columns1through120-0.00150.00040.0019-0.0106-0.4343-0.0965-0.00100.0108-0.0037-0.06064.2 計算機法潮流計算結(jié)果迭代次數(shù)：6沒有達到精度要求的個數(shù)：9101010830各節(jié)點的電壓復(fù)數(shù)值E為（節(jié)點號從小到大排列）：1.20001.0947-0.1316i0.9068-0.1361i0.6128-0.2991i1.1073-0.0932i1.1996-0.0298i各

39、節(jié)點的電壓模值大小V為（節(jié)點號從小到大排列）：1.20001.10260.91690.68191.11131.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.3000i5.0000+2.5555i各條支路的首端功率Si為（順序同您輸入B1時一致）：S(1,2)=5.0135+1.8332iS(2,3)=3.0895+0.48941iS(

40、2,5)=-0.176-0.31075iS(3,4)=0.78738+0.39432iS(4,5)=-0.90639-0.43087iS(6,5)=5+2.5555i各條支路的末端功率Sj為(順序同您輸入B1時一致)：S(2,1)=-5.0135-1.1787iS(3,2)=-2.5874-0.79432iS(5,2)=0.17703-0.29553iS(4,3)=-0.69361-0.36913iS(5,4)=1.123+1.1889iS(5,6)=-5-2.1934i各條支路的功率損耗DS為(順序同您輸入B1時一致)：DS(1,2)=-8.8818e-016+0.65451iDS(2,3)

41、=0.50212-0.3049iDS(2,5)=0.0010207-0.60628iDS(3,4)=0.093766+0.025185iDS(4,5)=0.21659+0.75805iDS(6,5)=0+0.36211i總網(wǎng)損S=0.8135+0.88867i4.3 系統(tǒng)功率分布圖1a&42o.do,o.4321012345-3 5251so 5 2 LO'Q滓生需5 251.50 s 2.10-立o1234s6-5051s25o4-1.-2.sKYfc:驍XS35nlQ502126S4321Oo,o.o.o,0,6震蜜娜ew-BA.4J0rM4o0fc5ob。q本譯目晨e以g

42、01個人心得在這次課程設(shè)計中，我強烈感覺到自己在很多方面的不足，對別人的依賴性比較強。我想我會在以后的學(xué)習(xí)中不斷去發(fā)現(xiàn)自己在很多方面的不足，并一一改正，希望在以后的工作中不要犯同樣的錯誤。在這次課程設(shè)計中，我們盡量按照老師的要求做，但在具體的操作過程中，還是出現(xiàn)了很多的問題。搞完這個課程設(shè)計讓我感覺電力系統(tǒng)分析是一門很有用的課程。因為我對它的學(xué)到的知識比較少。在很多時候我很多東西都不了解。并且走了很多的彎路。而且我感覺自己的知識不夠連貫。好些時候都出現(xiàn)了卡殼的情況。這次課程設(shè)計后，我一定要重新對電力系統(tǒng)分析這門課程做進一步的了解。對在此過程中遺留下的問題做好好的研究。爭取早點對電力系統(tǒng)分析這門

43、課程有個全方位的了解。為在以后的畢業(yè)課程設(shè)計中多些方案。也為我子以后走上工作崗位，提升自己的專業(yè)技能，打下扎實的基礎(chǔ)。還有就是在幾天的課程設(shè)計中，使我養(yǎng)成了很好的學(xué)習(xí)習(xí)慣，和對學(xué)習(xí)知識的嚴謹?shù)膽B(tài)度，同時也養(yǎng)成了積極查閱相關(guān)資料的好習(xí)慣，好習(xí)慣的養(yǎng)成是來之不易的，我相信在以后的學(xué)習(xí)和工作中，我將繼續(xù)保持這些良好的習(xí)慣，并積極努力的學(xué)習(xí)。讓自己更上一層樓。同時在此也感謝一直指導(dǎo)我的老師，此次課程設(shè)計的完成與老師的指導(dǎo)師分不開的，終在我們的一起努力下，完成了這門課程設(shè)計。在此對指導(dǎo)老師致以衷心的感謝。還要感謝學(xué)院圖書館為我提供豐富的參考資料，也感謝班上同學(xué)給了我很多寶貴的意見和參考，使我獲益很多。附

44、錄：源程序>>%本程序的功能是用牛頓拉夫遜法進行潮流計算%B1矩陣：1、支路首端號；2、末端號；3、支路阻抗；4、線路對地電納（或變壓器導(dǎo)納）；%5、支路的變比；6、支路首端處于K側(cè)為1，1側(cè)為0；%7、線路/變壓器標識（0/1）變壓器參數(shù)當(dāng)支路首端處于K側(cè)標識為1時歸算至末端側(cè)，0歸算至首端側(cè)%B2矩陣：1、該節(jié)點發(fā)電機功率；2、該節(jié)點負荷功率；%3、PQ節(jié)點電壓初始值，或平衡節(jié)點及PV節(jié)點電壓的給定值%4、節(jié)點所接無功補償并聯(lián)電容（感）的電納%5、節(jié)點分類標號：1為平衡節(jié)點（應(yīng)為1號節(jié)點）；2為PQ節(jié)點；3為PV節(jié)clear;isb=1;%input('請輸入平衡母線節(jié)

45、點號：isb=');pr=0.00001;%input('請輸入誤差精度：pr=');%n=6;%input('請輸入節(jié)點數(shù)：n=');nl=6;%input('請輸入支路數(shù)：nl=');B1=120+0.03i01.0511;230.06+0.025i0+0.5i100;250.04+0.25i0+0.5i100;340.08+0.30i0+0.5i100;450.1+0.35i0100;650+0.015i01.0511B2=001.2001;02.1+1i1.0002;01.8+0.40i1.0002;01.6+0.8i1.0002

46、;03.7+1.3i1.0002;0-5+0i1.2003%input('請輸入各節(jié)點參數(shù)形成的矩陣：B2=');%X=10;20;30;40;50;60%n=4;%input('請輸入節(jié)點數(shù)：n=');nl=4;%input('請輸入支路數(shù)：nl=');%B1=124+16i0100;134+16i0100;232+8i0100;241.49+48.02i011/11001%input('請輸入由支路參數(shù)形成的矩陣：B1=');%B2=0011501;0011002;020+4i11002;010+6i1002%input(&#

47、39;請輸入各節(jié)點參數(shù)形成的矩陣：B2=');%Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl);%求導(dǎo)納矩陣%fori=1:n%ifX(i,2)=0;%p=X(i,1);%Y(p,p)=1/X(i,2);%end%endfor i=1:nl%從1到n1（總支路數(shù)）if B1(i,7)=1%if B1(i,6)=0p=B1(i,1);q=B1(i,2);elsep=B1(i,2);q=B1(i,1); endY(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5);Y(q,p)=Y(p,q);如果是變壓器支路%左節(jié)點（首端）處于1側(cè)%左節(jié)點