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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

27、><p>　　YN，d11 10MVA 兩個變壓器并聯(lián)</p><p>　　YN，Y，d11 10MVA 兩個三相變壓器并聯(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變壓器兩個：</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é)點編號。節(jié)點編號由小到大編寫。</p><p>　　2、對于含有變壓

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

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

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

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

38、lt;p>　　Times=1; %置迭代次數(shù)為初始值</p><p><b>　　%創(chuàng)建節(jié)點導納矩陣</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，用于存儲初始功率參數(shù)</p><p><b>　　h=0;</b></p>

43、;<p><b>　　j=0;</b></p><p>　　for i=1:n %對PQ節(jié)點的處理</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é)點的處理，注意這時不可再將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 用于存儲PV節(jié)點的初始電壓</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，用于存儲有功功率、無功功率和電壓幅值的不平衡量</p><p><b>　　h=0;</b></p><p>　　for i=1:n %對PQ節(jié)點的處理</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é)點的處理，注意這時不可再將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，用于存儲節(jié)點電流參數(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é)點的處理</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é)點的處理</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é)點電壓的不平衡量</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é)點電壓</b></p><p><b>　　j=0;</b></p><p>　　for i=1:n

74、%對PQ節(jié)點處理</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é)點的處理</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é)點電壓</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è)計，使我更加扎實的掌握了有關(guān)潮流計算方面的知識，在設(shè)計過程中雖然遇到了一些問題，但經(jīng)過一次又一次的思考，一遍又一遍的檢查終于找出了原因所在，也暴露出了前期我在這方面的知識欠缺和經(jīng)驗不足。實踐出真知，通過親自動手計算，使我們掌握的知識不再是紙上談兵。</p><p>

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

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

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

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

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

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

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

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