【正文】 f(n)。endfor n=1:5 I(3,1)=I(3,1)+G(3,n)*e(n)B(3,n)*f(n)。endfor n=1:5 I(5,1)=I(5,1)+G(5,n)*e(n)B(5,n)*f(n)。M=[]。 %有功功率的不平衡量Q2=(2)*I(2,1)+e(2)*I(2,2)。Q4=(4)*I(4,1)+e(4)*I(4,2)。 N(m,m)=G(m,m)*e(m)+B(m,m)*f(m)+I(m,1)。 N(m,n)=G(m,n)*e(m)+B(m,n)*f(m)。M(2,2),L(2,2),M(2,3),L(2,3),M(2,4),L(2,4),M(2,5),L(2,5)。M(4,2),L(4,2),M(4,3),L(4,3),M(4,4),L(4,4),M(4,5),L(4,5)。雅克比矩陣J:39。C=[P2。P4。%解修正方程式disp(39。wucha=max(abs(A))。e(3)=e(3) +A(4,1)。e(5)=e(5) +A(8,1)。 disp(e)disp(39。u=e+f*i。 disp(u)。)。最后電壓誤差:39。enddisp(39。%計算各支路功率 endenddisp(39。（1）Clc，clear清屏，清除以前顯示結(jié)果；（2）disp 顯示待輸出的內(nèi)容，如disp(39。 矩陣的轉(zhuǎn)置。通過這次課程設(shè)計使我懂得了理論與實際相結(jié)合是很重要的，只有理論知識是遠遠不夠的，只有把所學(xué)的理論知識與實踐相結(jié)合起來，從理論中得出結(jié)論，才能真正為社會服務(wù)，從而提高自己的實際動手能力和獨立思考的能力。潮流計算靈活性和方便性的要求，對數(shù)字計算機的應(yīng)用也是一個很關(guān)鍵的問題。當(dāng)利用數(shù)字計算機進行潮流計算時，就失去了這種直觀性。在這個意義上，我們在編制潮流計算程序時，對使用的方便性和靈活性必須予以足夠的重視。同時讓我們意識到了團結(jié)的力量和重要性！希望在以后的學(xué)習(xí)中，同學(xué)之間都能互相幫助! 參考文獻及資料：《電力系統(tǒng)穩(wěn)態(tài)分析》，陳珩，中國電力出版社，2007，第三版《電力系統(tǒng)分析》，韓禎祥，浙江大學(xué)出版社，2005，第三版《電力系統(tǒng)分析》，蘇小林 閻曉霞 ，中國電力出版社《電力系統(tǒng)分析課程實際設(shè)計與綜合實驗》，祝書萍，中國電力出版社，2007，第一版—電力系統(tǒng)計算[M].北京:水利電力出版社，1978[M].北京：清華大學(xué)出版社，2006《火車發(fā)電廠、變電所二次接線技術(shù)規(guī)定（強電部分）》SDGJ─78《MALAB 科學(xué)運算完整解決法案》，精銳創(chuàng)組編，人民郵電大學(xué)出版社《MALAB 應(yīng)用集錦》，林雪松 周婧 林德新[M].北京：清華大學(xué)出版社，200628 / 29。由于以前沒學(xué)過Matlab，剛開始一直用的才C語言進行此編程，覺得還是C語言簡單些，所以花了幾天去上網(wǎng)查資料編程，在編程過程中遇到很多問題，如矩陣球逆矩陣，轉(zhuǎn)置一些大的循環(huán)，在老師的幫助下雖然有些進展，但是還是花費了大量時間沒搞出來，上網(wǎng)查的資料都是低階的，如矩陣球逆矩陣最多四階。電力系統(tǒng)潮流計算問題并不是單純的計算問題，把它當(dāng)作一個運行方式的調(diào)整問題可能更為確切。交流臺模擬了電力系統(tǒng)，因此在交流計算臺上計算潮流時，計算人員可以隨時監(jiān)視系統(tǒng)各部分運行狀態(tài)是否滿足要求，如發(fā)現(xiàn)某些部分運行不合理，則可以立即進行調(diào)整。在電力系統(tǒng)調(diào)度運行的多個領(lǐng)域都涉及到電網(wǎng)潮流計算。現(xiàn)在學(xué)到了不懂了要自己先找資料翻閱有關(guān)書籍。)，顯示：節(jié)點總數(shù)為：；（3）while，end；，for，end；循環(huán)語句，用于要循環(huán)的地方；（4）/ 正除；\ 反除；如DY=JJ\DW39。)。)。disp(wucha)。elsebreak。 disp(39。節(jié)點電壓的第C(k)次近似值:39。)。各點的電壓實部e(單位:V)為(節(jié)點號從小到大排列):39。e(4)=e(4) +A(6,1)。e(2)=e(2) +A(2,1)。)。P5。P3。 disp(J)。M(5,2),L(5,2),M(5,3),L(5,3),M(5,4),L(5,4),M(5,5),L(5,5)]。M(3,2),L(3,2),M(3,3),L(3,3),M(3,4),L(3,4),M(3,5),L(3,5)。 L(m,n)=H(m,n)。 L(m,m)=B(m,m)*e(m)+G(m,m)*f(m)I(m,2)。Q5=(5)*I(5,1)+e(5)*I(5,2)。Q3=(3)*I(3,1)+e(3)*I(3,2)。J=[]。 endH=[]。endfor n=1:5 I(4,1)=I(4,1)+G(4,n)*e(n)B(4,n)*f(n)。endfor n=1:5 I(2,1)=I(2,1)+G(2,n)*e(n)B(2,n)*f(n)。0,0。disp(B)k=0。)。0,0,]。B=[,0,0。,0。P5=。P3=。P1=0。PQ節(jié)點為：39。N=5disp(39。如：位陣的生成 eye(n) 、eye(m,n) ；零矩陣的生成 zeros(n) zeros(m,n) ；全1矩陣的生成 ones(n) ones(m,n) ；隨機元素矩陣函數(shù)rand(n,m) 、rand(n) ；對角矩陣diag(V) % V=[1 2 3 4]。第二節(jié) 矩陣的簡單基本運算矩陣是MATLAB數(shù)據(jù)存儲的基本單元，而矩陣的運算是MATLAB語言的核心，在MATLAB語言系統(tǒng)中幾乎一切運算均是以對矩陣的操作為基礎(chǔ)的。另外，MATLAB提供了一種特殊的工具：工具箱（TOOLBOXES）.這些工具箱主要包括：信號處理（SIGNAL PROCESSING）、控制系統(tǒng)（CONTROL SYSTEMS）、神經(jīng)網(wǎng)絡(luò)（NEURAL NETWORKS）、模糊邏輯(FUZZY LOGIC)、小波(WAVELETS)和模擬（SIMULATION）等等。MATLAB程序設(shè)計語言結(jié)構(gòu)完整，且具有優(yōu)良的移植性，它的基本數(shù)據(jù)元素是不需要定義的數(shù)組。假定系統(tǒng)中的第1,2,m號為P—Q節(jié)點,第m+1,m+2,n1為P—V節(jié)點,根據(jù)節(jié)點性質(zhì)的不同,得到如下迭代推算式: ⑴對于PQ節(jié)點 (324)⑵對于PV節(jié)點 (325)⑶對于平衡節(jié)點 平衡節(jié)點只設(shè)一個,電壓為已知,不參見迭代,其電壓為: (326)②. 修正方程式(235)和(236)兩組迭代式工包括2(n1)(235)和(236),并將其按泰勒級數(shù)展開,略去二次方程及以后各項,得到修正方程如下: (327) ③.雅可比矩陣各元素的算式式(328)中, 雅可比矩陣中的各元素可通過對式(324)和(325)時, 雅可比矩陣中非對角元素為 (328)當(dāng)時,雅可比矩陣中對角元素為: (329)第四章 Matlab概述第一節(jié) Matlab簡介目前電子計算機已廣泛應(yīng)用于電力系統(tǒng)的分析計算，潮流計算是其基本應(yīng)用軟件之一。當(dāng)采用直角坐標(biāo)時，潮流問題的待求量為各節(jié)點電壓的實部和虛部兩個分量由于平衡節(jié)點的電壓向量是給定的，因此待求共需要2(n1)個方程式。但若系統(tǒng)因無功緊張或其它原因?qū)е码妷嘿|(zhì)量很差或有重載線路而節(jié)點間角差很大時，仍用上述初始電壓就有可能出現(xiàn)問題。牛頓法所需的內(nèi)存量及每次迭代所需時間均較高斯法多。牛頓法當(dāng)初始估計值和方程的精確解足夠接近時，收斂速度非常快，具有平方收斂特性。接著就從出發(fā)，重復(fù)上述計算過程。其要點是把非線性方程式的求解過程變成反復(fù)地對相應(yīng)的線性方程式進行求解的過程。對現(xiàn)有電力系統(tǒng)的運行和擴建，對新的電力系統(tǒng)進行規(guī)劃設(shè)計以及對電力系統(tǒng)進行靜態(tài)和暫態(tài)穩(wěn)定分析都是以潮流計算為基礎(chǔ)。第三章 牛頓－拉夫遜法概述第一節(jié) 牛頓拉夫遜法基本原理電力系統(tǒng)潮流計算是電力系統(tǒng)分析中的一種最基本的計算，是對復(fù)雜電力系統(tǒng)正常和故障條件下穩(wěn)態(tài)運行狀態(tài)的計算。這一約束的主要意義就在于此。PU節(jié)點電壓幅值必須按上述條件給定?？梢赃x擇出線數(shù)多或者靠近電網(wǎng)中心的發(fā)電廠母線作平衡節(jié)點。對該節(jié)點，給定其電壓值，并在計算中取該節(jié)點電壓向量的方向作為參考軸，相當(dāng)于給定該點電壓向量的角度為零。這類節(jié)點在運行中往往要有一定可調(diào)節(jié)的無功電源。Q給定時，也作為PQ節(jié)點。然而在電力系統(tǒng)中，給出發(fā)電機或負荷連接母線上電壓或電流(都是向量)的情況是很少的，一般是給出發(fā)電機母線上發(fā)電機的有功功率(P)和母線電壓的幅值(U)，給出負荷母線上負荷消耗的有功功率(P)和無功功率(Q)。取 ，得到潮流方程的極坐標(biāo)形式： (215)b。它具有如下特點：（1）它是一組代數(shù)方程，因而表征的是電力系統(tǒng)的穩(wěn)定運行特性。二．非標(biāo)準(zhǔn)變比變壓器等值電路變壓器型等值電路更便于計算機反復(fù)計算,那就是變比。j=i)稱互導(dǎo)納，由此可得互導(dǎo)納數(shù)值上就等于在節(jié)點i施加單位電壓，其他節(jié)點全部接地時，經(jīng)節(jié)點j注入網(wǎng)絡(luò)的電流，因此可定義為：