電力系統(tǒng)潮流計(jì)算課程設(shè)計(jì)

1、<p>　　信息工程學(xué)院課程設(shè)計(jì)任務(wù)書</p><p>　學(xué)生姓名學(xué) 號成 績</p><p>　設(shè)計(jì)題目電力系統(tǒng)潮流計(jì)算</p><p>　設(shè)計(jì)內(nèi)容1. 對給定的網(wǎng)絡(luò)查找潮流計(jì)算所需的各元件等值參數(shù)，畫出等值電路圖。2.輸入各支路數(shù)據(jù)，各節(jié)點(diǎn)數(shù)據(jù)利用給定的程序進(jìn)行在變電所在某一負(fù)荷情況下的潮流計(jì)算，并對計(jì)算結(jié)果進(jìn)行分析。3.跟隨變電所負(fù)荷按一定比例發(fā)

2、生變化，進(jìn)行潮流計(jì)算分析。1）4個(gè)變電所的負(fù)荷同時(shí)以2%的比例增大；2）4個(gè)變電所的負(fù)荷同時(shí)以2%的比例下降1和4號變電所的負(fù)荷同時(shí)以2%的比例下降，而2和3號變電所的負(fù)荷同時(shí)以2%的比例上升；4.在不同的負(fù)荷情況下，分析潮流計(jì)算的結(jié)果，如果各母線電壓不滿足要求，進(jìn)行電壓的調(diào)整。（變電所低壓母線電壓10KV要求調(diào)整范圍在9.5-10.5之間；電壓35KV要求調(diào)整范圍在35-36之間）5.輪流斷開環(huán)網(wǎng)一回線，分析潮流的分布。6.利用仿真軟

3、件，進(jìn)行繪制系統(tǒng)圖進(jìn)行上述各種情況潮流的分析，并進(jìn)行結(jié)果的比較。7.最終形成課程設(shè)計(jì)成品說明書。</p><p>　設(shè)計(jì)要求1. 在讀懂程序的基礎(chǔ)上畫出潮流計(jì)算基本流程圖。2. 通過輸入數(shù)據(jù)，進(jìn)行潮流計(jì)算輸出結(jié)果</p><p>　時(shí)間安排第3周：對系統(tǒng)模型分析和資料收集，找到計(jì)算的思路，總結(jié)計(jì)算方法；理論計(jì)算第4周：對整個(gè)模型進(jìn)行詳細(xì)計(jì)算，MATLAB仿真并完成word文檔。</p

4、><p>　參考資料[1]何仰贊, 溫增銀《電力系統(tǒng)分析》（第三版）［M]. 華中科技大學(xué)，2002[1]康華光，陳大欽. 電子技術(shù)基礎(chǔ)—模擬部分（第五版）［M]. 北京：高等教育出版社，2005[3]鐘鍵偉 《電路基礎(chǔ)》中國電力出版社[4]湯蘊(yùn)繆 羅應(yīng)立《電機(jī)學(xué)》機(jī)械工業(yè)出版社</p><p><b>　　目 錄</b></p><p>　　

5、1 任務(wù)提出與方案論證3</p><p>　　1.1 潮流計(jì)算的意義和用途3</p><p>　　1.2 設(shè)計(jì)計(jì)算思路4</p><p><b>　　2 總體設(shè)計(jì)5</b></p><p>　　2.1潮流計(jì)算設(shè)計(jì)題目5</p><p><b>　　2.2課題分析6</b

6、></p><p><b>　　3 詳細(xì)設(shè)計(jì)8</b></p><p>　　3.1初步數(shù)據(jù)計(jì)算8</p><p><b>　　3.2仿真12</b></p><p><b>　　4 總結(jié)21</b></p><p><b>　　參

7、考文獻(xiàn)22</b></p><p>　　1 任務(wù)提出與方案論證</p><p>　　1.1 潮流計(jì)算的意義和用途</p><p>　　1.1.1什么是潮流計(jì)算</p><p>　　潮流計(jì)算是在給定電力系統(tǒng)網(wǎng)絡(luò)結(jié)構(gòu)、參數(shù)和決定系統(tǒng)運(yùn)行狀態(tài)的邊界條件的情況下確定系統(tǒng)穩(wěn)態(tài)運(yùn)行狀態(tài)的一種基本方法,是電力系統(tǒng)規(guī)劃和運(yùn)營中不可缺少的一個(gè)重要

8、組成部分?？梢哉f,它是電力系統(tǒng)分析中最基本、 最重要的計(jì)算,是系統(tǒng)安全、 經(jīng)濟(jì)分析和實(shí)時(shí)控制與調(diào)度的基礎(chǔ)。常規(guī)潮流計(jì)算的任務(wù)是根據(jù)給定的運(yùn)行條件和網(wǎng)路結(jié)構(gòu)確定整個(gè)系統(tǒng)的運(yùn)行狀態(tài)，如各母線上的電壓（幅值及相角）、網(wǎng)絡(luò)中的功率分布以及功率損耗等。潮流計(jì)算的結(jié)果是電力系統(tǒng)穩(wěn)定計(jì)算和故障分析的基礎(chǔ)。在電力系統(tǒng)運(yùn)行方式和規(guī)劃方案的研究中，都需要進(jìn)行潮流計(jì)算以比較運(yùn)行方式或規(guī)劃供電方案的可行性、可靠性和經(jīng)濟(jì)性。同時(shí)，為了實(shí)時(shí)監(jiān)控電力系統(tǒng)的運(yùn)行狀態(tài)，

9、也需要進(jìn)行大量而快速的潮流計(jì)算。因此，潮流計(jì)算是電力系統(tǒng)中應(yīng)用最廣泛、最基本和最重要的一種電氣運(yùn)算。</p><p>　　1.1.2潮流計(jì)算的意義</p><p>　　電力系統(tǒng)已經(jīng)與我們的生活息息相關(guān)，不可分割。進(jìn)行電力系統(tǒng)潮流計(jì)算是保證電力系統(tǒng)正常運(yùn)行的必要計(jì)算。具體來講電力系統(tǒng)潮流計(jì)算具有以下意義：</p><p>　?。?）在電網(wǎng)規(guī)劃階段,通過潮流計(jì)算,合理規(guī)

10、劃電源容量及接入點(diǎn),合理規(guī)劃網(wǎng)架,選擇無功補(bǔ)償方案,滿足規(guī)劃水平的大、小方式下潮流交換控制、調(diào)峰、調(diào)相、調(diào)壓的要求。 　　</p><p>　?。?）在編制年運(yùn)行方式時(shí),在預(yù)計(jì)負(fù)荷增長及新設(shè)備投運(yùn)基礎(chǔ)上,選擇典型方式進(jìn)行潮流計(jì)算,發(fā)現(xiàn)電網(wǎng)中薄弱環(huán)節(jié),供調(diào)度員日常調(diào)度控制參考,并對規(guī)劃、基建部門提出改進(jìn)網(wǎng)架結(jié)構(gòu),加快基建進(jìn)度的建議。 　</p><p>　　（3）正常檢修及特殊運(yùn)行方式下的潮

11、流計(jì)算,用于日運(yùn)行方式的編制,指導(dǎo)發(fā)電廠開機(jī)方式,有功、無功調(diào)整方案及負(fù)荷調(diào)整方案,滿足線路、變壓器熱穩(wěn)定要求及電壓質(zhì)量要求。 　</p><p>　?。?）預(yù)想事故、設(shè)備退出運(yùn)行對靜態(tài)安全的影響分析及作出預(yù)想的運(yùn)行方式調(diào)整方案?？偨Y(jié)為在電力系統(tǒng)運(yùn)行方式和規(guī)劃方案的研究中，都需要進(jìn)行潮流計(jì)算以比較運(yùn)行方式或規(guī)劃供電方案的可行性、可靠性和經(jīng)濟(jì)性。同時(shí)，為了實(shí)時(shí)監(jiān)控電力系統(tǒng)的運(yùn)行狀態(tài)，也需要進(jìn)行大量而快速的潮流計(jì)算。

12、</p><p>　　1.1.3潮流計(jì)算的用途</p><p>　　潮流計(jì)算的是電力系統(tǒng)中非常重要的分析計(jì)算，用以研究系統(tǒng)規(guī)劃和運(yùn)行中提出的各種問題。對規(guī)劃中電力系統(tǒng)，通過潮流計(jì)算可以檢驗(yàn)所提出的電力系統(tǒng)規(guī)劃方案能否滿足各種運(yùn)行方式的要求；對運(yùn)行中的電力系統(tǒng)，同過潮流計(jì)算可以預(yù)知各種負(fù)荷變化和網(wǎng)絡(luò)結(jié)構(gòu)的改變會(huì)不會(huì)危機(jī)系統(tǒng)的安全，系統(tǒng)中所有母線的電壓是否在允許范圍內(nèi)，，系統(tǒng)中各種元件（線路、

13、變壓器等）是否會(huì)出現(xiàn)過負(fù)荷，以及可能出現(xiàn)過負(fù)荷時(shí)應(yīng)事先采取哪些預(yù)防措施等。潮流計(jì)算是電力系統(tǒng)分析中最基本的計(jì)算。除它自身的重要作用之外，在《電力系統(tǒng)分析綜合程序》中，潮流計(jì)算、網(wǎng)損計(jì)算、靜態(tài)安全分析、暫態(tài)穩(wěn)定計(jì)算、小干擾靜態(tài)穩(wěn)定計(jì)算、短路計(jì)算、靜態(tài)和動(dòng)態(tài)等值計(jì)算的基礎(chǔ)。</p><p>　　(1)在電網(wǎng)規(guī)劃階段,通過潮流計(jì)算,合理規(guī)劃電源容量及接入點(diǎn),合理規(guī)劃網(wǎng)架,選擇無功補(bǔ)償方案,滿足規(guī)劃水平的大、小方式下潮流

14、交換控制、調(diào)峰、調(diào)相、調(diào)壓的要求?！　?2)在編制年運(yùn)行方式時(shí),在預(yù)計(jì)負(fù)荷增長及新設(shè)備投運(yùn)基礎(chǔ)上,選擇典型方式進(jìn)行潮流計(jì)算,發(fā)現(xiàn)電網(wǎng)中薄弱環(huán)節(jié),供調(diào)度員日常調(diào)度控制參考,并對規(guī)劃、基建部門提出改進(jìn)網(wǎng)架結(jié)構(gòu),加快基建進(jìn)度的建議?！　?3)正常檢修及特殊運(yùn)行方式下的潮流計(jì)算,用于日運(yùn)行方式的編制,指導(dǎo)發(fā)電廠開機(jī)方式,有功、無功調(diào)整方案及負(fù)荷調(diào)整方案,滿足線路、變壓器熱穩(wěn)定要求及電壓質(zhì)量要求?！　?4)預(yù)想事故、設(shè)備退出運(yùn)行對靜態(tài)安全的

15、影響分析及作出預(yù)想的運(yùn)行方式調(diào)整方案?！　】偨Y(jié)為在電力系統(tǒng)運(yùn)行方式和規(guī)劃方案的研究中，都需要進(jìn)行潮流計(jì)算以比較運(yùn)行方式或規(guī)劃供電方案的可行性、可靠性和經(jīng)濟(jì)性。同時(shí)，為了實(shí)時(shí)監(jiān)控電力系統(tǒng)的運(yùn)行狀態(tài)，也需要進(jìn)行大量而快速的潮流計(jì)算。因此，潮流計(jì)算是電力系統(tǒng)中應(yīng)用最廣泛、最基本和最重要的一種電氣運(yùn)算。在系統(tǒng)規(guī)劃設(shè)計(jì)和安排系統(tǒng)的運(yùn)行方式時(shí)，采用離線潮流計(jì)算；在電力系統(tǒng)運(yùn)行狀態(tài)的實(shí)時(shí)監(jiān)控中，則采用在線潮流計(jì)算。</p><p

16、>　　1.2 設(shè)計(jì)計(jì)算思路</p><p>　　運(yùn)用MATLAB軟件進(jìn)行潮流計(jì)算，對給定的題目進(jìn)行分析計(jì)算，再應(yīng)用DDRTS軟件，構(gòu)建系統(tǒng)圖進(jìn)行仿真，最終得到合理的系統(tǒng)潮流！</p><p>　　首先，畫出系統(tǒng)的等效電路圖，在計(jì)算出各元件參數(shù)的基礎(chǔ)上，應(yīng)用牛頓—拉夫遜法以及MATLAB軟件進(jìn)行計(jì)算給定系統(tǒng)圖進(jìn)行了四種不同負(fù)荷下的潮流計(jì)算，經(jīng)過調(diào)節(jié)均得到符合電壓限制及功率限制的潮流分布

17、。</p><p>　　其次，輪流斷開環(huán)網(wǎng)的三條支路，在新的系統(tǒng)結(jié)構(gòu)下進(jìn)行次潮流計(jì)算，結(jié)果亦均滿足潮流分布要求，牛頓—拉夫遜法具有良好的收斂性，上訴計(jì)算過程經(jīng)過四到五次迭代后均能收斂。</p><p>　　最后應(yīng)用Simulink軟件，構(gòu)建系統(tǒng)圖，對給定負(fù)荷重新進(jìn)行分析，潮流計(jì)算后的結(jié)果也能滿足相應(yīng)的參數(shù)要求。</p><p><b>　　2設(shè)計(jì)思路<

18、;/b></p><p>　　2.1 潮流計(jì)算設(shè)計(jì)題目</p><p>　　1、系統(tǒng)圖：兩個(gè)發(fā)電廠分別通過變壓器和輸電線路與四個(gè)變電所相連。如圖2-1所示。</p><p><b>　　圖2-1系統(tǒng)圖</b></p><p><b>　　2、發(fā)電廠資料：</b></p><

19、p>　　母線1和2為發(fā)電廠高壓母線，發(fā)電廠一總裝機(jī)容量為（300MW），母線3為機(jī)壓母線，機(jī)壓母線上裝機(jī)容量為（100MW），最大負(fù)荷和最小負(fù)荷分別為40MW和20MW；發(fā)電廠二總裝機(jī)容量為（200MW）。</p><p><b>　　3、變電所資料：</b></p><p>　　變電所1、2、3、4低壓母線的電壓等級分別為：10KV 35KV 35KV

20、 10KV</p><p>　　變電所的負(fù)荷分別為：</p><p>　　40MW 60MW 70MW 50MW </p><p>　　每個(gè)變電所的功率因數(shù)均為cosφ=0.85；</p><p>　　變電所1和變電所2分別配有兩臺(tái)容量為75MVA的變壓器，短路損耗414KW，短路電壓（%）=16.7；變電所3和變電所4分別

21、配有兩臺(tái)容量為63MVA的變壓器，短路損耗為245KW，短路電壓（%）=10.5；</p><p><b>　　4、輸電線路資料：</b></p><p>　　發(fā)電廠和變電所之間的輸電線路的電壓等級及長度標(biāo)于圖中，單位長度的電阻為，單位長度的電抗為，單位長度的電納為。</p><p><b>　　2.2數(shù)據(jù)計(jì)算</b>&l

22、t;/p><p><b>　　課題分析</b></p><p><b>　　1潮流計(jì)算等值電路</b></p><p><b>　　2建立電力系統(tǒng)模型</b></p><p>　　在Simulink中按照電力系統(tǒng)原型選擇元件進(jìn)行建模。</p><p>　　在

23、電力系統(tǒng)模型的建立工程中主要涉及到的是：元器件的選擇及其參數(shù)的設(shè)置；發(fā)電機(jī)選型；變壓器選擇；線路的選擇；負(fù)荷模型的選擇；母線選擇。</p><p><b>　　3模型的調(diào)試與運(yùn)行</b></p><p>　　建立系統(tǒng)模型，并設(shè)置好參數(shù)以后，就可以在Simulink環(huán)境下進(jìn)行仿真運(yùn)行。</p><p>　　4計(jì)算數(shù)據(jù)：首先計(jì)算各點(diǎn)的運(yùn)算負(fù)荷，其次

24、根據(jù)運(yùn)算負(fù)荷以及各個(gè)線路及變壓器的功率損耗計(jì)算各節(jié)點(diǎn)功率，然后用各節(jié)點(diǎn)的功率來計(jì)算電壓降落，從而計(jì)算各個(gè)節(jié)點(diǎn)的電壓。</p><p>　　1）變壓器的阻抗、，勵(lì)磁損耗。</p><p>　　2）線路的阻抗、，無功損耗。</p><p>　　3）變壓器及線路的功率損耗、</p><p>　　4)變壓器及線路的電壓降落、</p>&

25、lt;p>　　50各節(jié)點(diǎn)的功率和電壓</p><p>　　牛頓-拉夫遜潮流計(jì)算程序框圖如圖3-2所示</p><p><b>　　3 詳細(xì)設(shè)計(jì)</b></p><p><b>　　3.1初步數(shù)據(jù)計(jì)算</b></p><p>　　3.3.1變壓器的阻抗、，勵(lì)磁損耗</p><

26、p>　　YN，d11 16MVA 一個(gè)變壓器</p><p>　　YN，d11 20MVA 兩個(gè)變壓器并聯(lián)</p><p>　　(3) YN，d11 63MVA 四個(gè)變壓器并聯(lián)</p><p>　　YN，d11 16MVA 兩個(gè)變壓器并聯(lián)</p><p>　　YN，d11 63MVA 一個(gè)變壓器</p

27、><p>　　YN，d11 10MVA 兩個(gè)變壓器并聯(lián)</p><p>　　YN，Y，d11 10MVA 兩個(gè)三相變壓器并聯(lián)</p><p>　　3.3.2 線路阻抗及無功損耗</p><p><b>　　32km線路：</b></p><p><b>　　40km線路：<

28、/b></p><p><b>　　80km線路：</b></p><p><b>　　30km線路：</b></p><p><b>　　70km線路：</b></p><p>　　3.3.3變壓器及線路功率損耗</p><p>　　Ynd11

29、16Mw 變壓器</p><p>　　SE1=92.0458+j78.1295--17.7.91—j33.7316=74.3367+j44.3979 </p><p>　　YN，d11 63MVA變壓器：</p><p>　　SE2=50+j30.99--0.1055--j5.7677=49.8945+j25.2223</p><p>

30、　　SE3= SE1+ SE2=74.3367+j44.3979+49.8945+j25.2223=124.2312+j69.6202</p><p>　　YN，d11 10MVA變壓器兩個(gè)：</p><p>　　SE= 15+6.890625+j14.4703=21.890625+j14.4703</p><p><b>　　30km線路</b&g

31、t;</p><p>　　SG= 21.890625+j14.4703+1.7703+j3.372=23.2209+j17.8423</p><p><b>　　三相變壓器：</b></p><p><b>　　中壓側(cè)：</b></p><p>　　SH1= 23.6609+j17.8423+0.0

32、2246—j0.009=23.6834+j17.833</p><p><b>　　低壓側(cè)：</b></p><p>　　SH2= 50--35+j37.5--0.5138—j8.7679=14.4862+j28.7321</p><p><b>　　高壓側(cè)：</b></p><p>　　3.3.4

33、變壓器及線路電壓降落</p><p><b>　　3.2仿真</b></p><p>　　3.2.1MATLAB仿真</p><p>　　相關(guān)的原始數(shù)據(jù)輸入格式如下：</p><p>　　1、B1是支路參數(shù)矩陣，第一列和第二列是節(jié)點(diǎn)編號。節(jié)點(diǎn)編號由小到大編寫。</p><p>　　2、對于含有變壓

34、器的支路，第一列為低壓側(cè)節(jié)點(diǎn)編號，第二列為高壓側(cè)節(jié)點(diǎn)編號，將變壓器的串聯(lián)阻抗置于低壓側(cè)處理，第三列為支路的串列阻抗參數(shù)，第四列為支路的對地導(dǎo)納參數(shù)，第五烈為含變壓器支路的變壓器的變比，第六列為變壓器是否是否含有變壓器的參數(shù)，其中“1”為含有變壓器，“0”為不含有變壓器。</p><p>　　3、B2為節(jié)點(diǎn)參數(shù)矩陣，其中第一列為節(jié)點(diǎn)注入發(fā)電功率參數(shù)；第二列為節(jié)點(diǎn)負(fù)荷功率參數(shù)；第三列為節(jié)點(diǎn)電壓參數(shù)；第六列為節(jié)點(diǎn)類型參

35、數(shù)，其中“1”為平衡節(jié)點(diǎn)，“2”為PQ節(jié)點(diǎn)，“3”為PV節(jié)點(diǎn)參數(shù)。</p><p>　　4、X為節(jié)點(diǎn)號和對地參數(shù)矩陣。其中第一列為節(jié)點(diǎn)編號，第二列為節(jié)點(diǎn)對地參數(shù)。</p><p>　　3.2.2潮流計(jì)算源程序</p><p><b>　　close all</b></p><p><b>　　clear all

36、</b></p><p>　　n=input('請輸入節(jié)點(diǎn)數(shù):n=');</p><p>　　n1=input('請輸入支路數(shù):n1=');</p><p>　　isb=input('請輸入平衡節(jié)點(diǎn)號:isb=');</p><p>　　pr=input('請輸入誤差精度:p

37、r=');</p><p>　　B1=input('請輸入支路參數(shù):B1=');</p><p>　　B2=input('請輸入節(jié)點(diǎn)參數(shù):B2=');</p><p>　　X=input('節(jié)點(diǎn)號和對地參數(shù):X=');</p><p>　　Y=zeros(n);</p>&

38、lt;p>　　Times=1; %置迭代次數(shù)為初始值</p><p><b>　　%創(chuàng)建節(jié)點(diǎn)導(dǎo)納矩陣</b></p><p>　　for i=1:n1</p><p>　　if B1(i,6)==0 %不含變壓器的支路</p><p>　　p=B1(i,1);</p&

39、gt;<p>　　q=B1(i,2);</p><p>　　Y(p,q)=Y(p,q)-1/B1(i,3);</p><p>　　Y(q,p)=Y(p,q);</p><p>　　Y(p,p)=Y(p,p)+1/B1(i,3)+0.5*B1(i,4);</p><p>　　Y(q,q)=Y(q,q)+1/B1(i,3)+0.5*

40、B1(i,4);</p><p>　　else %含有變壓器的支路</p><p>　　p=B1(i,1);</p><p>　　q=B1(i,2);</p><p>　　Y(p,q)=Y(p,q)-1/(B1(i,3)*B1(i,5));</p><p>　　Y(q,p)=Y(p,q);

41、</p><p>　　Y(p,p)=Y(p,p)+1/B1(i,3);</p><p>　　Y(q,q)=Y(q,q)+1/(B1(i,5)^2*B1(i,3));</p><p><b>　　end</b></p><p><b>　　end</b></p><p><

42、;b>　　Y</b></p><p>　　OrgS=zeros(2*n-2,1);</p><p>　　DetaS=zeros(2*n-2,1); %將OrgS、DetaS初始化</p><p>　　%創(chuàng)建OrgS，用于存儲(chǔ)初始功率參數(shù)</p><p><b>　　h=0;</b></p>

43、;<p><b>　　j=0;</b></p><p>　　for i=1:n %對PQ節(jié)點(diǎn)的處理</p><p>　　if i~=isb&B2(i,6)==2 </p><p><b>　　h=h+1;</b></p><p><b>　　

44、for j=1:n</b></p><p>　　OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));</p>

45、<p>　　OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));</p><p><b>　　end</b&g

46、t;</p><p><b>　　end</b></p><p><b>　　end</b></p><p>　　for i=1:n %對PV節(jié)點(diǎn)的處理，注意這時(shí)不可再將h初始化為0</p><p>　　if i~=isb&B2(i,6)==3</p>&l

47、t;p><b>　　h=h+1;</b></p><p><b>　　for j=1:n</b></p><p>　　OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*

48、(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));</p><p>　　OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(

49、Y(i,j))*real(B2(j,3)));</p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　OrgS</b></p><

50、p>　　%創(chuàng)建PVU 用于存儲(chǔ)PV節(jié)點(diǎn)的初始電壓</p><p>　　PVU=zeros(n-h-1,1);</p><p><b>　　t=0;</b></p><p><b>　　for i=1:n</b></p><p>　　if B2(i,6)==3</p><p&

51、gt;<b>　　t=t+1;</b></p><p>　　PVU(t,1)=B2(i,3);</p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　PVU</b></p><

52、;p>　　%創(chuàng)建DetaS，用于存儲(chǔ)有功功率、無功功率和電壓幅值的不平衡量</p><p><b>　　h=0;</b></p><p>　　for i=1:n %對PQ節(jié)點(diǎn)的處理</p><p>　　if i~=isb&B2(i,6)==2</p><p><b>　　h=h

53、+1;</b></p><p>　　DetaS(2*h-1,1)=real(B2(i,2))-OrgS(2*h-1,1);</p><p>　　DetaS(2*h,1)=imag(B2(i,2))-OrgS(2*h,1);</p><p><b>　　end</b></p><p><b>　　en

54、d</b></p><p><b>　　t=0;</b></p><p>　　for i=1:n %對PV節(jié)點(diǎn)的處理，注意這時(shí)不可再將h初始化為0</p><p>　　if i~=isb&B2(i,6)==3</p><p><b>　　h=h+1;</b>&

55、lt;/p><p><b>　　t=t+1;</b></p><p>　　DetaS(2*h-1,1)=real(B2(i,2))-OrgS(2*h-1,1);</p><p>　　DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2-real(B2(i,3))^2-imag(B2(i,3))^2;</

56、p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　DetaS</b></p><p>　　%創(chuàng)建I，用于存儲(chǔ)節(jié)點(diǎn)電流參數(shù)</p><p>　　i=zeros(n-1,1);</p>

57、<p><b>　　h=0;</b></p><p><b>　　for i=1:n</b></p><p><b>　　if i~=isb</b></p><p><b>　　h=h+1;</b></p><p>　　I(h,1)=(OrgS(

58、2*h-1,1)-OrgS(2*h,1)*sqrt(-1))/conj(B2(i,3));</p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　I</b></p><p>　　%創(chuàng)建Jacbi(雅可比矩陣)<

59、/p><p>　　Jacbi=zeros(2*n-2);</p><p><b>　　h=0;</b></p><p><b>　　k=0;</b></p><p>　　for i=1:n %對PQ節(jié)點(diǎn)的處理</p><p>　　if B2(i,6)==2</p

60、><p><b>　　h=h+1;</b></p><p><b>　　for j=1:n</b></p><p><b>　　if j~=isb</b></p><p><b>　　k=k+1;</b></p><p>　　if i=

61、=j %對角元素的處理</p><p>　　Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1));</p><p>　　Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+r

62、eal(I(h,1));</p><p>　　Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k)+2*real(I(h,1));</p><p>　　Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-2*imag(I(h,1));</p><p>　　else %非對角元素的處理</p><p&g

63、t;　　Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3));</p><p>　　Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));</p><p>　　Jacbi(2*h,2*k-1)=-Jacbi(2

64、*h-1,2*k);</p><p>　　Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1);</p><p><b>　　end</b></p><p>　　if k==(n-1) %將用于內(nèi)循環(huán)的指針置于初始值，以確保雅可比矩陣換行</p><p><b>　　k=0;</b>&

65、lt;/p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　end</b

66、></p><p><b>　　k=0;</b></p><p>　　for i=1:n %對PV節(jié)點(diǎn)的處理</p><p>　　if B2(i,6)==3</p><p><b>　　h=h+1;</b></p><p><b>　　for j=

67、1:n</b></p><p><b>　　if j~=isb</b></p><p><b>　　k=k+1;</b></p><p>　　if i==j %對角元素的處理</p><p>　　Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i

68、,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1));</p><p>　　Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1));</p><p>　　Jacbi(2*h,2*k-1)=2*imag(B2(i,3));</p>

69、<p>　　Jacbi(2*h,2*k)=2*real(B2(i,3));</p><p>　　else %非對角元素的處理</p><p>　　Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3));</p><p>　　Jacbi(2*h-1,2*

70、k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));</p><p>　　Jacbi(2*h,2*k-1)=0;</p><p>　　Jacbi(2*h,2*k)=0;</p><p><b>　　end</b></p><p>　　if k==(n-1)

71、 %將用于內(nèi)循環(huán)的指針置于初始值，以確保雅可比矩陣換行</p><p><b>　　k=0;</b></p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　end</b></p>

72、<p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　Jacbi</b></p><p>　　%求解修正方程，獲取節(jié)點(diǎn)電壓的不平衡量</p><p>　　DetaU=zeros(2*n-2,1);</p&

73、gt;<p>　　DetaU=inv(Jacbi)*DetaS;</p><p><b>　　DetaU</b></p><p><b>　　%修正節(jié)點(diǎn)電壓</b></p><p><b>　　j=0;</b></p><p>　　for i=1:n

74、%對PQ節(jié)點(diǎn)處理</p><p>　　if B2(i,6)==2</p><p><b>　　j=j+1;</b></p><p>　　B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);</p><p><b>　　end</b></p&

75、gt;<p><b>　　end</b></p><p>　　for i=1:n %對PV節(jié)點(diǎn)的處理</p><p>　　if B2(i,6)==3</p><p><b>　　j=j+1;</b></p><p>　　B2(i,3)=B2(i,3)+DetaU(2*j,1

76、)+DetaU(2*j-1,1)*sqrt(-1);</p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　B2</b></p><p>　　%開始循環(huán)********************************

77、**************************************</p><p>　　while abs(max(DetaU))>pr</p><p>　　OrgS=zeros(2*n-2,1); %!!!初始功率參數(shù)在迭代過程中是不累加的，所以在這里必須將其初始化為零矩陣</p><p><b>　　h=0;</b&

78、gt;</p><p><b>　　j=0;</b></p><p><b>　　for i=1:n</b></p><p>　　if i~=isb&B2(i,6)==2</p><p><b>　　h=h+1;</b></p><p><

79、b>　　for j=1:n</b></p><p>　　OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));<

80、;/p><p>　　OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));</p><p><b>　　end&

81、lt;/b></p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　for i=1:n</b></p><p>　　if i~=isb&B2(i,6)==3</p><p>&l

82、t;b>　　h=h+1;</b></p><p><b>　　for j=1:n</b></p><p>　　OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i

83、,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));</p><p>　　OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*r

84、eal(B2(j,3)));</p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　OrgS</b></p><p><

85、b>　　%創(chuàng)建DetaS</b></p><p><b>　　h=0;</b></p><p><b>　　for i=1:n</b></p><p>　　if i~=isb&B2(i,6)==2</p><p><b>　　h=h+1;</b><

86、;/p><p>　　DetaS(2*h-1,1)=real(B2(i,2))-OrgS(2*h-1,1);</p><p>　　DetaS(2*h,1)=imag(B2(i,2))-OrgS(2*h,1);</p><p><b>　　end</b></p><p><b>　　end</b></

87、p><p><b>　　t=0;</b></p><p><b>　　for i=1:n</b></p><p>　　if i~=isb&B2(i,6)==3</p><p><b>　　h=h+1;</b></p><p><b>　　t

88、=t+1;</b></p><p>　　DetaS(2*h-1,1)=real(B2(i,2))-OrgS(2*h-1,1);</p><p>　　DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2-real(B2(i,3))^2-imag(B2(i,3))^2;</p><p><b>　　end&l

89、t;/b></p><p><b>　　end</b></p><p><b>　　DetaS</b></p><p><b>　　%創(chuàng)建I</b></p><p>　　i=zeros(n-1,1);</p><p><b>　　h=0;

90、</b></p><p><b>　　for i=1:n</b></p><p><b>　　if i~=isb</b></p><p><b>　　h=h+1;</b></p><p>　　I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqr

91、t(-1))/conj(B2(i,3));</p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　I</b></p><p><b>　　%創(chuàng)建Jacbi</b></p><

92、p>　　Jacbi=zeros(2*n-2);</p><p><b>　　h=0;</b></p><p><b>　　k=0;</b></p><p><b>　　for i=1:n</b></p><p>　　if B2(i,6)==2</p><

93、;p><b>　　h=h+1;</b></p><p><b>　　for j=1:n</b></p><p><b>　　if j~=isb</b></p><p><b>　　k=k+1;</b></p><p><b>　　if i=

94、=j</b></p><p>　　Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1));</p><p>　　Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real

95、(I(h,1));</p><p>　　Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k)+2*real(I(h,1));</p><p>　　Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-2*imag(I(h,1));</p><p><b>　　else</b></p><p>

96、;　　Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3));</p><p>　　Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));</p><p>　　Jacbi(2*h,2*k-1)=-Jacbi(2*

97、h-1,2*k);</p><p>　　Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1);</p><p><b>　　end</b></p><p>　　if k==(n-1)</p><p><b>　　k=0;</b></p><p><b>

98、;　　end</b></p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　end</b></p><p><

99、;b>　　k=0;</b></p><p><b>　　for i=1:n</b></p><p>　　if B2(i,6)==3</p><p><b>　　h=h+1;</b></p><p><b>　　for j=1:n</b></p>

100、<p><b>　　if j~=isb</b></p><p><b>　　k=k+1;</b></p><p><b>　　if i==j</b></p><p>　　Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*i

101、mag(B2(i,3))+imag(I(h,1));</p><p>　　Jacbi(2*h-1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1));</p><p>　　Jacbi(2*h,2*k-1)=2*imag(B2(i,3));</p><p>　　Jacbi(2*

102、h,2*k)=2*real(B2(i,3));</p><p><b>　　else</b></p><p>　　Jacbi(2*h-1,2*k-1)=-imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3));</p><p>　　Jacbi(2*h-1,2*k)=real(Y(i,j))*r

103、eal(B2(i,3))+imag(Y(i,j))*imag(B2(i,3));</p><p>　　Jacbi(2*h,2*k-1)=0;</p><p>　　Jacbi(2*h,2*k)=0;</p><p><b>　　end</b></p><p>　　if k==(n-1)</p><p&g

104、t;<b>　　k=0;</b></p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　end</b></p>&

105、lt;p><b>　　end</b></p><p><b>　　Jacbi</b></p><p>　　DetaU=zeros(2*n-2,1);</p><p>　　DetaU=inv(Jacbi)*DetaS;</p><p><b>　　DetaU</b><

106、/p><p><b>　　%修正節(jié)點(diǎn)電壓</b></p><p><b>　　j=0;</b></p><p><b>　　for i=1:n</b></p><p>　　if B2(i,6)==2</p><p><b>　　j=j+1;<

107、/b></p><p>　　B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);</p><p><b>　　end</b></p><p><b>　　end</b></p><p><b>　　for i=1:n</b&

108、gt;</p><p>　　if B2(i,6)==3</p><p><b>　　j=j+1;</b></p><p>　　B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);</p><p><b>　　end</b></p>

109、<p><b>　　end</b></p><p><b>　　B2</b></p><p>　　Times=Times+1; %迭代次數(shù)加1</p><p><b>　　end</b></p><p><b>　　Times</b>&

110、lt;/p><p><b>　　4 總結(jié)</b></p><p>　　通過此次課程設(shè)計(jì)，使我更加扎實(shí)的掌握了有關(guān)潮流計(jì)算方面的知識，在設(shè)計(jì)過程中雖然遇到了一些問題，但經(jīng)過一次又一次的思考，一遍又一遍的檢查終于找出了原因所在，也暴露出了前期我在這方面的知識欠缺和經(jīng)驗(yàn)不足。實(shí)踐出真知，通過親自動(dòng)手計(jì)算，使我們掌握的知識不再是紙上談兵。</p><p>

111、　　過而能改，善莫大焉。在課程設(shè)計(jì)過程中，我們不斷發(fā)現(xiàn)錯(cuò)誤，不斷改正，不斷領(lǐng)悟，不斷獲取。最終的檢測調(diào)試環(huán)節(jié)，本身就是在踐行“過而能改，善莫大焉”的知行觀。這次課程設(shè)計(jì)終于順利完成了，在設(shè)計(jì)中遇到了很多問題，最后在老師的指導(dǎo)下，終于游逆而解。在今后社會(huì)的發(fā)展和學(xué)習(xí)實(shí)踐過程中，一定要不懈努力，不能遇到問題就想到要退縮，一定要不厭其煩的發(fā)現(xiàn)問題所在，然后一一進(jìn)行解決，只有這樣，才能成功的做成想做的事，才能在今后的道路上劈荊斬棘，而不是知難而

112、退，那樣永遠(yuǎn)不可能收獲成功，收獲喜悅，也永遠(yuǎn)不可能得到社會(huì)及他人對你的認(rèn)可！</p><p>　　課程設(shè)計(jì)誠然是一門專業(yè)課，給我很多專業(yè)知識以及專業(yè)技能上的提升，同時(shí)又是一門講道課，一門辯思課，給了我許多道，給了我很多思，給了我莫大的空間。同時(shí)，設(shè)計(jì)讓我感觸很深。使我對抽象的理論有了具體的認(rèn)識。通過這次課程設(shè)計(jì)，我掌握了潮流計(jì)算的原理。</p><p>　　我認(rèn)為，在這學(xué)期的實(shí)驗(yàn)中，不僅培

113、養(yǎng)了獨(dú)立思考、動(dòng)手操作的能力，在各種其它能力上也都有了提高。更重要的是，在實(shí)驗(yàn)課上，我們學(xué)會(huì)了很多學(xué)習(xí)的方法。而這是日后最實(shí)用的，真的是受益匪淺。要面對社會(huì)的挑戰(zhàn)，只有不斷的學(xué)習(xí)、實(shí)踐，再學(xué)習(xí)、再實(shí)踐。這對于我們的將來也有很大的幫助。以后，不管有多苦，我想我們都能變苦為樂，找尋有趣的事情，發(fā)現(xiàn)其中珍貴的事情。就像中國提倡的艱苦奮斗一樣，我們都可以在實(shí)驗(yàn)結(jié)束之后變的更加成熟，會(huì)面對需要面對的事情。</p><p>

114、　　回顧起此課程設(shè)計(jì)，至今我仍感慨頗多，從理論到實(shí)踐，在這段日子里，可以說得是苦多于甜，但是可以學(xué)到很多很多的東西，同時(shí)不僅可以鞏固了以前所學(xué)過的知識，而且學(xué)到了很多在書本上所沒有學(xué)到過的知識。通過這次課程設(shè)計(jì)使我懂得了理論與實(shí)際相結(jié)合是很重要的，只有理論知識是遠(yuǎn)遠(yuǎn)不夠的，只有把所學(xué)的理論知識與實(shí)踐相結(jié)合起來，從理論中得出結(jié)論，才能真正為社會(huì)服務(wù)，從而提高自己的實(shí)際動(dòng)手能力和獨(dú)立思考的能力。在設(shè)計(jì)的過程中遇到問題，可以說得是困難重重，但

115、可喜的是最終都得到了解決。 </p><p>　　通過本次課程設(shè)計(jì)讓我有復(fù)習(xí)了一次潮流計(jì)算的相關(guān)知識，跟家清晰了什么事潮流計(jì)算以及潮流計(jì)算的在電力系統(tǒng)的重要性。電力系統(tǒng)的穩(wěn)定運(yùn)行狀況即是正常運(yùn)行狀況，是指電力系統(tǒng)在穩(wěn)定運(yùn)行條件下電壓、功率的分布，也稱為潮流分布。電力系統(tǒng)分析的潮流計(jì)算是電力系統(tǒng)分析的一個(gè)重要的部分。通過對電力系統(tǒng)潮流分布的分析和計(jì)算，可進(jìn)一步對系統(tǒng)運(yùn)行的安全性，經(jīng)濟(jì)性進(jìn)行分析、評估，提出改進(jìn)措施。

116、同時(shí)潮流分布也是電力系統(tǒng)規(guī)劃設(shè)計(jì)的一項(xiàng)基礎(chǔ)工作。</p><p>　　整個(gè)計(jì)算過程的模型建立并不是十分復(fù)雜，但計(jì)算過程十分繁瑣、計(jì)算量相當(dāng)?shù)拇?，而且由于枝?jié)太多很容易算錯(cuò)。不過在計(jì)算潮流計(jì)算的過程中卻對以往學(xué)過的電力系統(tǒng)分析的相關(guān)知識進(jìn)行了一次較為深入的復(fù)習(xí)。而且整個(gè)計(jì)算對計(jì)算量的要求很大，鍛煉了我們的計(jì)算能力。而且對細(xì)節(jié)的把握也得到了鍛煉，做題的精細(xì)程度得到了提高。此次設(shè)計(jì)也讓我明白了思路即出路，有什么不懂不明

117、白的地方要及時(shí)請教或上網(wǎng)查詢，只要認(rèn)真鉆研，動(dòng)腦思考，動(dòng)手實(shí)踐，就沒有弄不懂的知識，收獲頗豐。</p><p><b>　　參考文獻(xiàn)</b></p><p>　　[1]康華光，陳大欽. 電子技術(shù)基礎(chǔ)—模擬部分（第五版）［M]. 北京：高等教育出版社，2005</p><p>　　[2]皮文兵. 一種寬輸入范圍的Gillbert模擬乘法器設(shè)計(jì)［J

118、]. 電子設(shè)計(jì)應(yīng)用.2007.13（1）：88-90</p><p>　　[3] 焦春生.新型綠色能效D類音頻放大器設(shè)計(jì)應(yīng)用[EB/OL] </p><p>　　[4] 張筑生. 微分半動(dòng)力系統(tǒng)的不變集 [D]. 北京: 北京大學(xué)數(shù)學(xué)系數(shù)學(xué)研究所, 1983</p><p>　　[5]鐘鍵偉 《電路基礎(chǔ)》中國電力出版社</p><p>　　[