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