電力系統(tǒng)潮流計算,DOC

1、電力系統(tǒng)課程設(shè)計題目：電力系統(tǒng)潮流計算 院系名稱：電氣工程學(xué)院 專業(yè)班級：電氣F1206班 學(xué)生姓名：學(xué)號：指導(dǎo)教師：張孝遠(yuǎn)2.成績：目錄料指導(dǎo)老師簽名：.2日期：4卜算節(jié)點(diǎn)介紹4-變量的分類51概述5原始資2潮流2.1 TOC o 1-5 h z 3計算方法簡介6 HYPERLINK l bookmark41 o Current Document 3.1牛頓一拉夫遜法原理63.1.1牛頓一拉夫遜法概要6 HYPERLINK l bookmark44 o Current Document 3.1.2牛頓法的框圖及求解過程8 HYPERLINK l bookmark55 o Current Do

2、cument 3.2MATLAB 簡介94潮流分布計算10104.1系統(tǒng)的一次接線圖 HYPERLINK l bookmark61 o Current Document 4.2參數(shù)計算10 HYPERLINK l bookmark109 o Current Document 4.3豐大及枯大下地潮流分布情況144.3.1該地區(qū)變壓器的有功潮流分布數(shù)據(jù)154.3.2重、過載負(fù)荷元件統(tǒng)計表175設(shè)計心得17參考文獻(xiàn)18附錄:程序19原始資料一、系統(tǒng)接線圖見附件1。二、系統(tǒng)中包含發(fā)電廠、變電站、及其間的聯(lián)絡(luò)線路。500kV變電站以外的系統(tǒng)以一個等值發(fā)電機(jī)代替。各元件的參數(shù)見附件2。設(shè)計任務(wù)手動畫出該

3、系統(tǒng)的電氣一次接線圖，建立實(shí)際網(wǎng)絡(luò)和模擬網(wǎng)絡(luò)之間的聯(lián)系。1、2、根據(jù)已有資料，先手算出各元件的參數(shù)，后再用Matlab表格核算出各元件的參數(shù)。3、潮流計算1）對兩種不同運(yùn)行方式進(jìn)行潮流計算，注意110kV電網(wǎng)開環(huán)運(yùn)行。2）注意將電壓調(diào)整到合理的范圍110kV母線電壓控制在106kV117kV之間;220kV母線電壓控制在220kV242kV之間。附件一:500kV變電站附件二：1、變壓器：兩個220kV變電站均采用參數(shù)一致的三繞組變壓器，具體參數(shù) 如下。220kV變電站參數(shù)表高壓側(cè)繞組中壓側(cè)繞組低壓側(cè)繞組容量12012060(120)電壓220110(121)10.5110kV及以下的變電站的

4、變壓器省略，即可將負(fù)荷直接掛在110kV母線上。而110kV升壓變只計及以下參數(shù)。110kV變電站參數(shù)表序號變電 站名容量X1X01A162B283C1205D636E207G31.58F12.5+31.59水電站130MW、2*20MVA0.822500.5250010水電站272MW、3*31.5MVA0.344440.3333311水電站324MW、31.5MVA1.033330.3333312水電 站54*7.5MW、2*20MVA1.033330.3333313水電 站418MW、30MVA1.283331.0500014火電 廠50MW、31.5+40MVA0.482850.262

5、502、線路：具體參數(shù)如下。220kV線路參數(shù)表序號線路名稱導(dǎo)線牌號線路長度km1ML2x240102ML2x240103MH2x3005110kV線路參數(shù)表序號線路名稱導(dǎo)線牌號線路長度km1水電站1水電站2150302AB95803BC9514水電站4C150655BD95636CD240607水電站2C240758ED24014.79水電站3D15010.510水電站5T節(jié)點(diǎn)705.511T節(jié)點(diǎn)D70512T節(jié)點(diǎn)G240513HF1503.914HD185715水電站2L240/213016LG2401017DF150418火電廠D24013、發(fā)電機(jī)各發(fā)電機(jī)的參數(shù)如下:Xd Xq 裝機(jī) 容

6、量功率 因數(shù)水電 站10.1740.174300.85水電站20.1740.174720.85水電站30.1740.174240.85水電 站50.1740.174300.85水電 站40.1740.174180.85華鑫電廠0.1740.17450.85出力情況：水力發(fā)電機(jī)豐大出力70%,枯大出力20%?；鹆Πl(fā)電機(jī)豐大出力80%,枯大出力80%。4、負(fù)荷各110kV變電站豐大負(fù)荷按該站變電容量的50%估算，枯大負(fù)荷按該站變電 容量的60%估算。兩個220kV變電站的低壓側(cè)上各掛10MW的負(fù)荷，中壓側(cè)各掛20MW負(fù)荷。功率因素均為0.95。5、并聯(lián)電容器兩個220kV變電站的低壓側(cè)上均裝設(shè)并聯(lián)

7、補(bǔ)償。補(bǔ)償總量按該站變電容量 的20%裝設(shè)，分組原則以每組電容器的容量不超過10MVar且經(jīng)濟(jì)性較好為準(zhǔn)。1概述潮流計算是電力系統(tǒng)最基本最常用的計算。根據(jù)系統(tǒng)給定的運(yùn)行條件，網(wǎng) 絡(luò)接線及元件參數(shù)，通過潮流計算可以確定各母線的電壓，包括電壓的幅值和 相角，各元件流過的功率，整個系統(tǒng)的功率損耗等一系列系統(tǒng)中的潮流數(shù)據(jù)。近幾年，對潮流算法的研究仍然是如何改善傳統(tǒng)的潮流算法，即高斯-塞德 爾法、牛頓法和快速解耦法。牛頓法，由于其在求解非線性潮流方程時采用的 是逐次線性化的方法，為了進(jìn)一步提高算法的收斂性和計算速度，人們考慮采 用將泰勒級數(shù)的高階項或非線性項也考慮進(jìn)來，于是產(chǎn)生了二階潮流算法。后 來又提

8、出了根據(jù)直角坐標(biāo)形式的潮流方程是一個二次代數(shù)方程的特點(diǎn)，提出了 采用直角坐標(biāo)的保留非線性快速潮流算法。潮流計算在數(shù)學(xué)上是多元非線性方程組的求解問題，求解的方法有很多種， 牛頓一拉夫遜Newton-Raphson法是數(shù)學(xué)上解非線性方程組的有效方法，有較好 的收斂性。將N-R法用于潮流計算是以導(dǎo)納矩陣為基礎(chǔ)的，由于利用了導(dǎo)納矩 陣的對稱性，稀疏性及節(jié)點(diǎn)編號順序優(yōu)劃等技巧，使N-R法在收斂性，占用內(nèi) 存，計算速度等方面的優(yōu)點(diǎn)都超過了阻抗法總結(jié)為在電力系統(tǒng)運(yùn)行方式和規(guī)劃方案的研究中都需要進(jìn)行潮流計算以比 較運(yùn)行方式或規(guī)劃供電方案的可行性、可靠性和經(jīng)濟(jì)性。同時為了實(shí)時監(jiān)控電 力系統(tǒng)的運(yùn)行狀態(tài)也需要進(jìn)行大

9、量而快速的潮流計算。因此潮流計算是電力系 統(tǒng)中應(yīng)用最廣泛、最基本和最重要的一種電氣運(yùn)算。在系統(tǒng)規(guī)劃設(shè)計和安排系 統(tǒng)的運(yùn)行方式時采用離線潮流計算在電力系統(tǒng)運(yùn)行狀態(tài)的實(shí)時監(jiān)控中則采用在 線潮流計算。2潮流計算節(jié)點(diǎn)介紹常規(guī)的電力系統(tǒng)潮流計算中一般具有三種類型的節(jié)點(diǎn)：PQ、PV及平衡節(jié)點(diǎn)。 一個節(jié)點(diǎn)有四個變量，即注入有功功率、注入無功功率，電壓大小及相角。常 規(guī)的潮流計算一般給定其中的二個變量：PQ節(jié)點(diǎn)（注入有功功率及無功功率）， PV節(jié)點(diǎn)（注入有功功率及電壓的大小），平衡節(jié)點(diǎn)（電壓的大小及相角）。2.1變量的分類負(fù)荷消耗的有功、無功功率一一匕、qli、P22、ql2電源發(fā)出的有功、無功功率匕、qgi

10、、pg2、Qg2母線或節(jié)點(diǎn)的電壓大小和相位一一U 1、U2、8 1、5 2在這十二個變量中，負(fù)荷消耗的有功和無功功率無法控制，因它們?nèi)Q于 用戶，它們就稱為不可控變量或是擾動變量。電源發(fā)出的有功無功功率是可以 控制的自變量，因此它們就稱為控制變量。母線或節(jié)點(diǎn)電壓的大小和相位角一 是受控制變量控制的因變量。其中，U、U主要受Q、Q的控制，5、5主 1 2G1G 21 2要受PG1、PG2的控制。這四個變量就是簡單系統(tǒng)的狀態(tài)變量。為了保證系統(tǒng)的正常運(yùn)行必須滿足以下的約束條件:對控制變量對沒有電源的節(jié)點(diǎn)則為對狀態(tài)變量匕的約束條件則是對某些狀態(tài)變量5還有如下的約束條件i2.2節(jié)點(diǎn)的分類第一類稱PQ節(jié)點(diǎn)

11、。等值負(fù)荷功率P、Q和等值電源功率P、Q是給 Li LiGi Gi定的，從而注入功率P、Q是給定的，待求的則是節(jié)點(diǎn)電壓的大小氣和相位角 5,。屬于這類節(jié)點(diǎn)的有按給定有功、無功率發(fā)電的發(fā)電廠母線和沒有其他電源 的變電所母線。第二類稱pv節(jié)點(diǎn)。等值負(fù)荷和等值電源的有功功率pl匕.是給定的，從而注入有功功率P是給定的。等值負(fù)荷的無功功率q和節(jié)點(diǎn)電壓的大小U也 ini是給定的。待求的則是等值電源的無功功率Q，從而注入無功功率Qt和節(jié)點(diǎn)電 壓的相位角5。有一定無功功率儲備的發(fā)電廠和有一定無功功率電源的變電所 母線都可以作為PV節(jié)點(diǎn)；第三類平衡節(jié)點(diǎn)。潮流計算時一般只設(shè)一個平衡節(jié)點(diǎn)。等值負(fù)荷功率Pl、 Ql

12、?是給定的，節(jié)點(diǎn)電壓的大小和相位也是給定的。擔(dān)負(fù)調(diào)整系統(tǒng)頻率任務(wù)的發(fā) 電廠母線往往被選作為平衡節(jié)點(diǎn)。3計算方法簡介3.1牛頓一拉夫遜法原理3.1.1牛頓一拉夫遜法概要首先對一般的牛頓一拉夫遜法作一簡單的說明。已知一個變量X函數(shù)為：到此方程時，由適當(dāng)?shù)慕浦祒出發(fā)，根據(jù)：反復(fù)進(jìn)行計算，當(dāng)X滿足適當(dāng)?shù)氖諗織l件就是上面方程的根。這樣的方 法就是所謂的牛頓一拉夫遜法。這一方法還可以做下面的解釋，設(shè)第n次迭代得到的解語真值之差，即X 的誤差為時，則：把f （ X （n） +e ）在X（ n）附近對用泰勒級數(shù)展開上式省略去 2以后部分X（n）的誤差可以近似由上式計算出來。比較兩式，可以看出牛頓一拉夫遜法的

13、休整量和X（n）的誤差的一次項相等。用同樣的方法考慮，給出n個變量的n個方程：對其近似解X得修正量AX可以通過解下邊的方程來確定： TOC o 1-5 h z 11式中等號右邊的矩陣 f 都是對于X1，X；,., X：的值。這一矩陣稱為雅可比n（JACOBI）矩陣。按上述得到的修正向量AX, AX，,AX后，得到如下關(guān)系12n這比X。X，,X更接近真實(shí)值。這一步在收斂到希望的值以前重復(fù)進(jìn)行， 12n一般要反復(fù)計算滿足為預(yù)先規(guī)定的小正數(shù)，Xnn+1是第n次迭代Xn的近似值。3.1.2牛頓法的框圖及求解過程1、用牛頓法計算潮流時，有以下的步驟：（1）給這各節(jié)點(diǎn)電壓初始值e（o）, f（o）；（2）

14、將以上電壓初始值代入公式，求修正方程的常數(shù)項向量AP（0）, AQ （o） ,（AV 2 ）（o）；（3）將電壓初始值在帶入上述公式，求出修正方程中系數(shù)矩陣的各元素。（4）解修正方程式Ae（0）,Af （0）；(5)修正各節(jié)點(diǎn)電壓e=e(0) + Ae(0)，f =f (0) + Af (0)；（6）將e，f在帶入方程式，求出AP（i）,AQ（i）,（AV2）（1）；（7）檢驗(yàn)是否收斂，即maxAp（k）,AQ（k） Ksii（8）如果收斂，迭代到此結(jié)束，進(jìn)一步計算各線路潮流和平衡節(jié)點(diǎn)功率， 并打印輸出結(jié)果。如果不收斂，轉(zhuǎn)回（2）進(jìn)行下次迭代計算，直到收斂為止。2、程序框圖如下3.2MATLA

15、B 簡介MATLAB是用于算法開發(fā)、數(shù)據(jù)可視化、數(shù)據(jù)分析以及數(shù)值計算的高級技術(shù) 計算語言和交互式環(huán)境，主要包括MATLAB和Simulink兩大部分。是由美國 mathworks公司發(fā)布的主要面對科學(xué)計算、可視化以及交互式程序設(shè)計的高科技 計算環(huán)境。它將數(shù)值分析、矩陣計算、科學(xué)數(shù)據(jù)可視化以及非線性動態(tài)系統(tǒng)的 建模和仿真等諸多強(qiáng)大功能集成在一個易于使用的視窗環(huán)境中，為科學(xué)研究、 工程設(shè)計以及必須進(jìn)行有效數(shù)值計算的眾多科學(xué)領(lǐng)域提供了一種全面的解決方 案，并在很大程度上擺脫了傳統(tǒng)非交互式程序設(shè)計語言（如C、Fortran）的編 輯模式，代表了當(dāng)今國際科學(xué)計算軟件的先進(jìn)水平。MATLAB是一種交互式、

16、面向?qū)ο蟮某绦蛟O(shè)計語言廣泛應(yīng)用于工業(yè)界與學(xué) 術(shù)界主要用于矩陣運(yùn)算同時在數(shù)值分析、自動控制模擬、數(shù)字信號處理、 動態(tài)分析、繪圖等方面也具有強(qiáng)大的功能。MATLAB程序設(shè)計語言結(jié)構(gòu)完整 且具有優(yōu)良的移植性 它的基本數(shù)據(jù)元素 是不需要定義的數(shù)組。它可以高效率地解決工業(yè)計算問題 特別是關(guān)于矩陣和 矢量的計算。MATLAB與C語言和FORTRAN語言相比更容易被掌握。通過M語言 可以用類似數(shù)學(xué)公式的方式來編寫算法大大降低了程序所需的難度并節(jié)省了 時間，從而可把主要的精力集中在算法的構(gòu)思而不是編程上。目前電子計算機(jī)已廣泛應(yīng)用于電力系統(tǒng)的分析計算潮流計算是其基本應(yīng)用軟件之一?，F(xiàn)有很多潮流計算方法。對潮流計算

17、方法有五方面的要求（1）計算 速度快（2）內(nèi)存需要少（3）計算結(jié)果有良好的可靠性和可信（4）適應(yīng)性好 亦即能處理變壓器變比調(diào)整、系統(tǒng)元件的不同描述和與其它程序配合的能力強(qiáng)。4潮流分布計算4.1系統(tǒng)的一次接線圖圖4.1系統(tǒng)的一次連接圖4.2參數(shù)計算設(shè)定基準(zhǔn)值SB 100MVA , Ub二Uav.n,則各參數(shù)如下。（1）發(fā)電機(jī)的次暫態(tài)電抗：X=X*Sb/Sn, ZB=Ub2/Sn發(fā)電機(jī)參數(shù)單位（MW）電廠裝 機(jī) 容 量枯水 出力 比例豐水 出力 比例豐大有 功豐大無 功枯大有 功枯大無 功短路X*水電 站1300.20.721.00013.0156.0003.7180.493水電 站2720.20

18、.750.40031.23514.4008.9240.205417水電 站3240.20.716.80010.4124.8002.9750.61625水電 站5300.20.721.00013.0156.0003.7180.493華鑫 電廠500.80.840.00024.79040.00024.7900.2958（2）110KV升壓變壓器的參數(shù):電阻：R=P *U 2/1000S2；KNN電抗：X=Uk（%）*Un2/100Sn；電導(dǎo)：G=P0/1000Un2；電納：B=I0（%）*Sn/100Un2；式中U以KV為單位，S以MVA為單位，P、P以KW為單位。NN0 K110KV變壓器參數(shù)變

19、電站 名容量X1X0SbR*X*G*B*A16B28C120D63E20G31.5F12.5+31.5水電站130MW、2*20MVA0.82250.5251000.0373640.5812120.0001990.001445水電站272MW、3*31.5MVA0.344440.333331000.0211990.3690230.0002800.001992水電站324MW、31.5MVA1.033330.333331000.0211990.3690230.0002800.001992水電站54*7.5MW、2*20MVA1.033330.333331000.0373640.5812120.0

20、001990.001445火電廠50MW、31.5+40MVA0.482850.26251000.0211990.3690230.0002800.0019921000.0138380.2906060.0003790.002529（3）220KV三繞組變壓器的參數(shù):電阻：R=P *U 2/1000S2；KNN電抗：X=Uk(%)*Un2/100Sn；電導(dǎo)：G=P0/1000UN2；電納：B=I0(%)*Sn/100Un2220kV變電站參數(shù)表SBR*X*G*B*高壓側(cè)繞 組中壓側(cè)繞 組低壓側(cè)繞 組高壓 側(cè)繞 組 （1）1000.0029200.1166500.0016500.012000容量12

21、0120120中壓 側(cè)繞 組 （2）100-0.0021380.0716670.0004990.003630電壓22012110.5低壓側(cè)繞組（3）1000.003580-0.0045750.0000040.000027（4） 110KV線路參數(shù):r=p/S；二=0Sig； 一二5 Xl*=x1*l*SB/U2av.n。110KV線路參數(shù)標(biāo)么值序號線路名稱導(dǎo)線牌號線路長度 kmUB(KV)SB(MW)R*X*B*1水電 站1水電 站2150301151000.0476370.0907370.0115062AB95801151000.1996220.2522500.0297303BC951115

22、1000.0024950.0031530.0003724水電 站4C150651151000.1032140.1965970.0249295BD95631151000.1572020.1986470.0234126CD240601151000.0598870.1751230.0239647水電 站2C240751151000.0748580.2189040.0299558ED24014.71151000.0146720.0429050.0058719水電 站3D15010.51151000.0166730.0317580.00402710水電 站5T節(jié) 點(diǎn) 八、705.51151000.018

23、7150.0177580.00198611T節(jié) 點(diǎn) 八、D7051151000.0170130.0161440.00180512T節(jié) 點(diǎn) 八、G24051151000.0049910.0145940.00199713HF1503.91151000.0061930.0117960.00149614H0089980.0210660.00274015水電 站2L2401301151000.1297540.3794330.05192116LG240101151000.0099810.0291870.00399417D0063520.0120980.0

24、0153418火電 廠D24011151000.0009980.0029190.000399110KV變電站負(fù)荷參數(shù)變電站總?cè)?量豐變站比枯站變比豐大有功豐大無 功枯大有 功枯大無功A160.50.68.00002.62959.60003.1554B280.50.614.00004.601616.80005.5219C1200.50.660.000019.721072.000023.6653D630.50.631.500010.353537.800012.4243E200.50.610.00003.286812.00003.9442F440.50.622.00007.231126.40008.

25、6773G31.50.50.615.75005.176818.90006.2121各變電站補(bǔ)償電容參數(shù)、站別運(yùn)行 方式、ABCDEFGH低壓 側(cè)L低壓 側(cè)豐大34.609.805.901.92.0枯大34.40005.901.92.04.3豐大及枯大下地潮流分布情況電壓是衡量電力系統(tǒng)電能質(zhì)量的標(biāo)準(zhǔn)之一。電壓過高或過低，都將對人身 及其用電設(shè)備產(chǎn)生重大的影響。保證用戶的電壓接近額定值是電力系統(tǒng)調(diào)度的 基本任務(wù)之一。當(dāng)系統(tǒng)的電壓偏離允許值時，電力系統(tǒng)必須應(yīng)用電壓調(diào)節(jié)技術(shù) 調(diào)節(jié)系統(tǒng)電壓的大小，使其維持在允許值范圍內(nèi)。本文經(jīng)過手算形成了等值電 路圖，并編寫好了程序得出節(jié)點(diǎn)電壓標(biāo)幺值，使其滿足所要求的調(diào)

26、整范圍。我們首先對給定的程序輸入部分作了簡要的分析，程序開始需要我們確定 輸入節(jié)點(diǎn)數(shù)、支路數(shù)、平衡母線號、支路參數(shù)矩陣、節(jié)點(diǎn)參數(shù)矩陣。（1）為了保證整個系統(tǒng)潮流計算的完整性，我們把凡具有母線及發(fā)電機(jī)處 均選作節(jié)點(diǎn)，這樣，我們確定電廠一母線上的發(fā)電機(jī)作為平衡節(jié)點(diǎn)，節(jié)點(diǎn)號為 ，其它機(jī)組作為PV節(jié)點(diǎn)，節(jié)點(diǎn)號為，其余節(jié)點(diǎn)均為PQ節(jié)點(diǎn)，節(jié)點(diǎn)號 見等值電路圖。（2）確定完節(jié)點(diǎn)及編號后，各條支路也相應(yīng)確定了，我們對各支路參數(shù)進(jìn)行 了計算。根據(jù)所給實(shí)際電路圖和題中的已知條件，有以下公式計算各輸電線路的阻抗和對地支路電容的標(biāo)幺值和變壓器的阻抗標(biāo)幺值。4.3.1該地區(qū)變壓器的有功潮流分布數(shù)據(jù)該圖為豐大潮流模型圖：

27、豐大潮流模型圖該運(yùn)行方式的電壓合理，負(fù)荷分配也均勻，但是有些線路的負(fù)載率偏低。 比如水電廠2L站的負(fù)載率僅僅為21.1%，ML站的線路的負(fù)載率也只是 8.8%。該圖為枯大潮流模型圖：枯大潮流模型圖該運(yùn)行方式的電壓合理，負(fù)荷分配也均勻，有些線路的負(fù)載率偏低但是有 些線路的負(fù)載率則偏高。比如水電廠2L站的負(fù)載率僅僅為9.4%，ML站的線路的負(fù)載率也只是7.6%D H站的線路卻重載到84%。4.3.2重、過載負(fù)荷元件統(tǒng)計表類型負(fù)載率實(shí)際是在功率額定容量建議線路0.83870.486.2限D(zhuǎn)站負(fù)荷到4MW以下5設(shè)計心得通過這次課程設(shè)計，我發(fā)現(xiàn)自己有很多不足的地方，如基礎(chǔ)知識掌握不牢固, 很多知識點(diǎn)都忘

28、記了，計算速度慢及準(zhǔn)確性低，分析問題能力不夠全面等等。同 時，在設(shè)計的過程中遇到很多問題，如怎樣使用WORD的工具，計算公式輸入， 畫圖等。明白了有些東西看起來很簡單，但一旦做起來卻需要很多心思，要注 意到很多細(xì)節(jié)問題。要做到能好好理解課本的內(nèi)容，一定要認(rèn)認(rèn)真真做一次計 算。因此，完成課程設(shè)計使我對課本的內(nèi)容加深了理解?？傮w來說，這次的課 程設(shè)計不單在專業(yè)基礎(chǔ)方面反映了我的學(xué)習(xí)還要加倍努力，還在對一些軟件的應(yīng)用需要加強(qiáng)。計算在各種情況下的潮流分布，對于豐大和枯大情況下的潮流分布有了明 確的認(rèn)識。在本次課程設(shè)計過程當(dāng)中，鍛煉了自己實(shí)際操作分析能力，理論聯(lián) 系實(shí)際，對運(yùn)行中的電力系統(tǒng)，通過潮流計算

29、可以預(yù)知各種負(fù)荷變化和網(wǎng)絡(luò)結(jié) 構(gòu)的改變會不會危及系統(tǒng)的安全，系統(tǒng)中所有母線的電壓是否在允許的范圍以 內(nèi)，系統(tǒng)中各種元件（線路、變壓器等）是否會出現(xiàn)過負(fù)荷，以及可能出現(xiàn)過負(fù) 荷時應(yīng)事先采取哪些預(yù)防措施等。同時采用MATLAB進(jìn)行計算機(jī)的計算，在計算 時采用特殊算法使得潮流計算的過程更快，效率更高。而采用計算機(jī)的運(yùn)算應(yīng) 該是未來的一種趨勢，所以我會學(xué)習(xí)一定的編輯語言如C,C+等，以提高運(yùn)算準(zhǔn) 確性和快速性?？傮w而言，這次的課程設(shè)計對我們運(yùn)用所學(xué)知識，發(fā)現(xiàn)、提出、分析和解 決實(shí)際問題、鍛煉實(shí)踐能力的考察，使我們更清楚地知道不足之出，從而提高 我們。參考文獻(xiàn)于永源主編.電力系統(tǒng)分析湖南師范大學(xué)出版社M

30、.1992年7月陳珩編.電力系統(tǒng)穩(wěn)態(tài)分析水利電力出版社M.1995年1月第二版邱曉燕.劉天琪 電力系統(tǒng)分析的計算機(jī)算法 北京 中國電力出版200李光琦.電力系統(tǒng)暫態(tài)分析M.北京：水利電力出版社，2002.5?陸敏政主編.電力系統(tǒng)習(xí)題集水利電力出版社M.1990年附錄：程序%潮流計算MATLAB粗略程序？?？%creatanew_datat=0;s=0;r=0;w=0;number二input(Howmanynodearethere=)；%ConvertPqtoanewarrayforii=1:number?ifdata(ii,4)=1?t=t+1;?forjj=1:14?new_data1(t

31、,jj)=data(ii,jj);?end;?a(1,t)=ii;?s=s+1;?%recordthenumberofthePQnode?end;end;%Convertpvtoanewarrayforii=1:number?ifdata(ii,4)=2?t=t+1;?forjj=1:14?new_data1(t,jj)=data(ii,jj);?end;?a(1,t)=ii;?r=r+1;?%recordthenumberofthePVnode?end;end;%Convertset_vtoanewarrayforii=1:number?ifdata(ii,4)=3?t=t+1;?forjj

32、=1:14?new_data1(t,jj)=data(ii,jj);?end;?a(1,t)=ii;?w=w+1;?end;end;%creatanew_data2x,y=size(data2)?forii=1:x?forjj=1:2?formm=1:number?ifdata2(ii,jj)=a(1,mm)?new_data2(ii,jj)=mm;?end;?end;?end;end;forii=1:x?forjj=3:14?new_data2(ii,jj)=data2(ii,jj);?end;end;%creataYY=zeros(number,number);YY=zeros(numbe

33、r,number);yy=zeros(number,number);forii=1:x ?%forjj=1:14?iii=new_data2(ii,1);?jjj=new_data2(ii,2);?ifnew_data2(ii,5)=2?sub=new_data2(ii,6)./(new_data2(ii,7).*new_data2(ii,7)+new_data2(ii,6).*new_data2(ii,6)-new_data2(ii,7)./(new_data2(ii,7).*n ew_data2(ii,7)+new_data2(ii,6).*new_data2(ii,6)*i;?Y(iii

34、,jjj)=-sub./new_data2(ii,14);?YY(iii,jjj)=sub./new_data2(ii,14);?Y(jjj,iii)=-sub/new_data2(ii,14);?YY(jjj,iii)=sub./new_data2(ii,14);?yy(iii,jjj) = (1.-new_data2(ii,14)./(new_data2(ii,14).*new_data2(ii,14).*sub;?yy(jjj,iii) = (new_data2(ii,14)-1)./(new_data2(ii,14).*sub;?else?Y(iii,jjj)=-new_data2(i

35、i,6)./(new_data2(ii,7).*new_data2(+new_data2(ii,6).*new_data2(ii,6)+new_data2(ii,7)./(new_data2(.*new_data2(ii,7)+new_data2(ii,6).*new_data2(ii,6)*i;?YY(iii,jjj)=new_data2(ii,6)./(new_data2(ii,7).*new_data2+new_data2(ii,6).*new_data2(ii,6)-new_data2(ii,7)./(new_data2.*new_data2(ii,7)+new_data2(ii,6)

36、.*new_data2(ii,6)*i;?Y(jjj,iii)=-new_data2(ii,6)./(new_data2(ii,7).*new_data2(+new_data2(ii,6).*new_data2(ii,6)+new_data2(ii,7)./(new_data2(.*new_data2(ii,7)+new_data2(ii,6).*new_data2(ii,6)*i;?YY(jjj,iii)=new_data2(ii,6)./(new_data2(ii,7).*new_data2(+new_data2(ii,6).*new_data2(ii,6)-new_data2(ii,7)

37、./(new_data2(.*new_data2(ii,7)+new_data2(ii,6).*new_data2(ii,6)*i;?yy(iii,jjj)=new_data2(ii,8)./2.*i;?yy(jjj,iii)=new_data2(ii,8)./2.*i;?end;?%end;end;foriii=1:number?Y(iii,iii)=0;end;%forii=1:x?%?forjj=1:14?for?iii=1:number?forjj=1:number?%ifiii=jj?Y(iii,iii)=Y(iii,iii)+YY(iii,jj)+yy(iii,jj);?%end;

38、?end;end;%creatB,Gforii=1:number?forjj=1:number?G(ii,jj)=real(Y(ii,jj);?B(ii,jj)=imag(Y(ii,jj);?end;end;%creatInitial_PInitial_QInitial_Vforii=1:(s+r)?set_P(ii,1) = (new_data1(ii,9)-new_data1(ii,7)./100;end;forii=1:s;?set_Q(ii,1) = (new_data1(ii,10)-new_data1(ii,8)./100;end;forii=1:r?set_V(ii,1)=new

39、_data1(ii+s,12).*new_data1(ii+s,12);%trytomodifyforsikeofcorrectingend;Initial_p_q_v=set_P;set_Q;set_V;disp(Initial_p_q_v);%creatInitial_e,Initial_fforii=1:number-1?e(ii,1)=1;?f(ii,1)=0.0;%changeftotestusedtobe1.0end;?e(number,1)=new_data1(number,12);?f(number,1)=0;?%e(64,1)=0.88;%test118ieee?%f(64,

40、1)=0.39395826829394;?%f(14,1)=0;?%e(10,1)=1.045;?%e(11,1)=1.01;?%e(12,1)=1.07;?%e(13,1)=1.09;%StartNEWTOWNCALULATIONfortry_time=1:25%CreateverynodeconsumePQandUn=s;m=r;forii=1:(n+m)?sum1=0;?forjj=1:(n+m+1)?sum1=sum1+e(ii,1).*(G(ii,jj).*e(jj,1)-B(ii,jj).*f(jj,1)+f(ii,1).*(G(ii,jj).*f(jj,1)+B(ii,jj).*

41、e(jj,1);?end;?p(ii,1)=sum1;end;forii=1:n?sum2=0;?forjj=1:(n+m+1)?sum2=sum2+f(ii,1).*(G(ii,jj).*e(jj,1)-B(ii,jj).*f(jj,1)-e(ii,1).*(G(ii,jj).*f(jj,1)+B(ii,jj).*e(jj,1);?end;?q(ii,1)=sum2;end;disp(，q=，);disp(q);u=zeros(n+m),1);forii=(n+1):(n+m)?u(ii,1)=e(ii,1).*e(ii,1)+f(ii,1).*f(ii,1);end;forii=n+1:(

42、n+m)?extra_u(ii-n),1)=u(ii,1);end;disp(，extra_u=，);disp(extra_u);sum=p;q;extra_u;disp(sum)disp(s);disp(p);%creatJacobiandisp(n);disp(m);forii=1:(n+m)?forjj=1:(n+m)?if(ii=jj)?PF(ii,jj)=B(ii,jj).*e(ii,1)-G(ii,jj).*f(ii,1);?PE(ii,jj)=-G(ii,jj).*e(ii,1)-B(ii,jj).*f(ii,1);?else?ss=0;?qq=0;?fornum=1:(n+m+

43、1)?ss=ss+G(ii,num).*f(num,1)+B(ii,num).*e(num,1);?qq=qq+G(ii,num).*e(num,1)-B(ii,num).*f(num,1);?end;?PF(ii,jj)=-ss+B(ii,jj).*e(ii,1)-G(ii,jj).*f(ii,1);%TEST+1?PE(ii,jj)=-qq-G(ii,jj).*e(ii,1)-B(ii,jj).*f(ii,1);%TEST+1?end;?end;end;copy=3.14159;disp(=copy=)forii=1:n?forjj=1:m+n?if(ii=jj)?QE(ii,jj)=B(

44、ii,jj).*e(ii,1)-G(ii,jj).*f(ii,1);%TEST+1 ?QF(ii,jj)=G(ii,jj).*e(ii,1)+B(ii,jj).*f(ii,1);%TEST+1?else?ss=0;?qq=0;?fornum=1:(n+m+1)?ss=ss+G(ii,num).*f(num,1)+B(ii,num).*e(num,1);?qq=qq+G(ii,num).*e(num,1)-B(ii,num).*f(num,1);?end;?QF(ii,jj)=-qq+G(ii,jj).*e(ii,1)+B(ii,jj).*f(ii,1);%TEST+1?QE(ii,jj)=ss

45、+B(ii,jj).*e(ii,1)-G(ii,jj).*f(ii,1);%TEST+1 ?end;?end;end;%disp(QF);%disp(QF);%disp(QE);%disp(QE);UE=zeros(n+m),(n+m);UF=zeros(n+m),(n+m);forii=n+1:n+m?forjj=1:(n+m)?if(ii=jj)?UE(ii,jj)=0;?UF(ii,jj)=0;?else?ss=0;?qq=0;?fornum=1:(n+m+1)?ss=ss+G(ii,num).*f(num,1)+B(ii,num).*e(num,1);?qq=qq+G(ii,num).

46、*e(num,1)-B(ii,num).*f(num,1);?end;?UF(ii,jj)=-2.*f(ii,1);?UE(ii,jj)=-2.*e(ii,1);?end;?end;end;forii=(n+1):(n+m)?forjj=1:(n+m)?extra_UE(ii-n),jj)=UE(ii,jj);?extra_UF(ii-n),jj)=UF(ii,jj);?end;end;%disp(extra_UE);%disp(extra_UE);%disp(extra_Uf);%disp(extra_UF);Jacobian=PF,PE;QF,QE;extra_UF,extra_UE;%d

47、isp(Jacobian=)；%disp(Jacobian);%creatsubstractresultsubstract_result=Initial_p_q_v-sum;%disp(substract_result);%disp(substract_result);%calculatedelta_f_edelta_f_e=-inv(Jacobian)*substract_result;%disp(delta_f_e);forii=1:number-1;?f(ii,1)=f(ii,1)+delta_f_e(ii,1);?e(ii,1)=e(ii,1)+delta_f_e(ii+number-1,1);end;ifmax(substract_result)1e-4?break;end;end;%disp(substract_result);%disp(substract_result);%disp(e=)；%disp(e);%disp(f=)；%disp(f);forii=1:number?uuu(ii,1)=e(ii