【正文】 um(Pt)。F(5)=0。Q=imag(S)。 F=imag(U)。 end end end %導(dǎo)納矩陣 G=real(Y)。 %求互導(dǎo)納 for m=1:5 for n=1:5 if Z(m,n)==0 y(m,n)=0。+,0,+,+,+。 5） 分為優(yōu)、良、中、合格 、不合格五個等級。 。 當(dāng) j=i 時，對角元素為 111122()()()()22niij j ij j ii i ii i iijiniij j ij j ii i ii i iijiniij j ij j ii i ii i iijiniij j ij j ii i ii ii iijiiiiiiiPG e B f G e B f NePG f B e B e G f HfQG f B e B e G f LeQG e B f G e B f JfUeeUff?????? ?? ? ? ? ? ????? ? ? ? ? ????? ? ? ? ?????? ? ? ? ? ????????????????????????????????? （ 14） 當(dāng) ji? 時 ,矩陣非對角元素為： 22()0iiij i ij i ij ijijiiij i ij i ij ijjjiijjPQG e B f N JefPQB e G f H LfeUUef?? ? ? ?? ? ? ? ? ? ? ? ??? ??? ? ? ? ?? ? ? ? ? ??? ??? ? ? ?????? ?? （ 15） 三、極坐標(biāo)計算 對于潮流計算中待于求出功率的節(jié)點功率方程組，在某個近似解附近用泰勒級數(shù)展開略去二階及以上的高階得到已矩陣表示的修正方程： 式中節(jié)點 n 為節(jié)點數(shù)， m為 PV節(jié)點，雅閣比矩陣是（ 2n2m2）階非奇異矩陣，雅可比矩陣各元素表示如下： 牛頓拉夫遜極坐標(biāo)潮流計算的修正方程的迭代方程為： 設(shè)計資料及參數(shù)： 課程名稱 電力系統(tǒng)穩(wěn)態(tài)分析課程 設(shè)計 設(shè)計題目 兩機(jī)五節(jié)點網(wǎng)絡(luò)潮流計算 — 牛拉法 指導(dǎo)教師 劉景霞 時間 1周 一、教學(xué)要求 電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計 以設(shè)計和優(yōu)化電力系統(tǒng)的潮流分析為重點，提高學(xué)生綜合能力為目標(biāo)，盡可能結(jié)合實際工程進(jìn)行。把牛頓法用于潮流計算，采用直角坐標(biāo)形式。有了“對應(yīng)元素”的規(guī)定，數(shù)組的運算實質(zhì)上就是針對數(shù)組內(nèi)部的每個元素進(jìn)行的。而乘除法運算有相 當(dāng)大的區(qū)別，數(shù)組的乘除法是指 兩同維數(shù)組對應(yīng)元素之間的乘除法，它們的運算符為“ .*”和“ ./”或“ .\”。但需注意進(jìn)行數(shù)除時，常數(shù)通常只能做除數(shù)。不同領(lǐng)域、不 同層次的用戶通過相應(yīng)工具的學(xué)習(xí)和應(yīng)用，可以方便地進(jìn)行計算、分析及設(shè)計工作。它可以高效率地解決工業(yè)計算問題，特別是關(guān)于矩陣和矢量的計算。 P Q 分解法是為了改進(jìn)牛頓 拉夫遜法在內(nèi)存占用量及計算速度方面的不足， P Q 分解法根據(jù)電力系統(tǒng)實際運行狀態(tài)的物理特點，對極坐標(biāo)形式的牛頓 拉夫遜法修正方程式進(jìn)行了合理的簡化。這是因為雅可比 陣元素的數(shù)目約為 2(n 1) 2(n 1)個（直角坐標(biāo)），且其數(shù)值在迭代過程中不斷變化。但是高斯 賽德爾潮流計算法在牛頓法以及各種解耦法出現(xiàn)以后似乎成了一種邊緣性的方法。潮流計算是研究電力系統(tǒng)穩(wěn)態(tài)運行情況的一種計算，是根據(jù)給定的運行條件及系統(tǒng)接線情況確定整個電力系統(tǒng)各個部分的運行狀態(tài)，如各母線的電壓、各元件中流過的功率、系統(tǒng)的功率損耗等等。 交流計算臺。最初，電力系統(tǒng)潮流計算是通過人工手算的，后來為了適應(yīng)電力系統(tǒng)日益發(fā)展的需要，采用了 錯誤 !未指定書簽。經(jīng)過幾十年的時間，電力系統(tǒng)潮流計算已經(jīng)發(fā)展得十分成熟。 潮流計算方法分析比較 高斯 賽德爾潮流計算法原理簡單，編程實現(xiàn)容易，特別是對于配網(wǎng)潮流有其獨特優(yōu)勢。牛頓法的缺點是每次迭代的計算量和所需的內(nèi)存量較大。對于正常運行的系統(tǒng)，各節(jié)電電壓一般均在額定值附近，偏移不會太大，并且各節(jié)電的相角差也不大，所以對各節(jié)電可以采用統(tǒng)一的電壓初值。 MATLAB 程序設(shè)計語言結(jié)構(gòu)完整，且具有優(yōu)良的移植性，它的基本數(shù)據(jù)元素是不需要定義的數(shù)組。另外， MATLAB 提供了一種特殊的工具：工具箱（ TOOLBOXES ） . 這些工具箱主要包括：信號處理（ SIGNAL PROCESSING）、控制系統(tǒng)（ CONTROL SYSTEMS）、神經(jīng)網(wǎng)絡(luò)（ NEURAL NETWORKS ）、 模 糊 邏 輯 (FUZZY LOGIC) 、 小 波 (WAVELETS) 和 模 擬（ SIMULATION）等等。 矩陣的運算 常數(shù)與矩陣的運算即是同該矩陣的每一元素進(jìn)行運算。 數(shù)組的加、減與矩陣的加、減運算完全相同。另外，矩陣的數(shù)組運算中還有冪運算（運算符為 .^ ）、指數(shù)運算（ exp）、對數(shù)運算（ log）、和開方運算（ sqrt）等。 二、直角坐標(biāo)系計算 (1)牛頓 拉 夫遜法潮流計算的公式。假定系統(tǒng)中的第m+ m+?、 n1節(jié)點為 PV節(jié)點，則對其中每一 PV 節(jié)點可以列寫方程 ： 112 2 2 2 2 2( ) ( ) 0() nni is i is i ij j ij j i ij j ij jjjis i is i iP P P P e G e B f f G f B eU U U U e f???? ? ? ? ? ? ? ? ????? ? ? ? ? ? ???（ 13） i =（ m+ m+?、 n1） (4)形成雅可比矩陣。 。 四、進(jìn)度安排 根據(jù)給定的參數(shù)或工程具體要求，收集和查閱資料（半天） 學(xué)習(xí)軟件（ MATLAB 或 C 語言等） （一天半） 編程計算復(fù)雜系統(tǒng)潮流計算（三天） 編寫計算設(shè)計書（一天） 五、評分標(biāo)準(zhǔn) 課程設(shè)計成績評定依據(jù)包括以下幾點： 1） 工作態(tài)度（占 10%）； 2） 基本技能的掌握程度（占 20%）； 3） 程序編寫是否合理是否有運行結(jié)果 (40%)； 4） 課程設(shè)計說明書編寫水平 (占 30%)。 %節(jié)點阻抗矩陣 Z=[0,+,0,+,0。0,+,+,0,0]。 else Y(m,n)=sum(y(m,:))。 E=real(U)。0] P=real(S)。 while C E(5)=。 Qt(n)=(F(m)*(G(m,n)*E(n)B(m,n)*F(n))E(m)*(G(m,n)*F(n)+B(m,n)*E(n)))。 for m=1:4 for n=1:5 Bi(n)=G(m,n)*F(n)+B(m,n)*E(n)。 J(m,m)=2*B(m,m)*F(m)+sum(Ai)(G(m,m)*E(m)B(m,m)*F(m))。 JJ(2*m,2*m1)=J(m,m)。 J(m,n)=B(m,n)*F(m)G(m,n)*E(m)。 JJ(2*m,2*n1)=J(m,n)。 end dU=inv(JJ)*PQ39。 end for n=1:5 U(n)=E(n)+(F(n))*i。 end S(5)=U(5)*sum(conj(I)) Y = + 0