【正文】 2*G(m,m)PT1(m)。39。 for m=1:n1+n2 O(m)=O(m)+OO(m)。 for R=1:n PH1(R)=U(isb)*conj(Y(isb,R))*conj(U(R))。)。)。)。0 0 0 0。在進(jìn)行課題設(shè)計(jì)的過程中，加深了我對(duì)潮流計(jì)算的認(rèn)識(shí)，尤其是對(duì)牛頓拉夫遜潮流計(jì)算的求解思路有 了比較透徹的理解。 0 0 0。 程序結(jié)果 請(qǐng)輸入結(jié)點(diǎn)數(shù)： n=4 請(qǐng)輸入 PV結(jié)點(diǎn)數(shù)： n1=1 請(qǐng)輸入 PQ結(jié)點(diǎn)數(shù)： n2=2 請(qǐng)輸入平衡結(jié)點(diǎn)： isb=4 請(qǐng)輸入精確度： pr= 請(qǐng)輸入變比矩陣看 :K=[0 0 0 0。 disp(39。 disp(39。 end end end disp(39。 c(m)=U1(m)*sin(O(m))。 end UD2=diag(UD)。 disp(JJ)。 end end end for m=1:n1+n2 for R=1:n2 if m==R N(m,m)=U1(m)^2*G(m,m)PT1(m)。 end QT1(m)=sum(QT)。 P=real(S)。迭代中的雅克比矩陣 :39。 end end end for m=1:n for R=1:n if m~=R Z(m)=Z(m)+1/C(m,R)。N=zeros(n1+n2,n2)。)。 K=input(39。請(qǐng)輸入 PQ 結(jié)點(diǎn)數(shù)： n2=39。)。 命令格式： [ n， d ] = paralltl（ n1， d1， n2， d2） [a， b， c， d] = paralltl（ a1， b1， c1， d1， a2， b2， c2， d2） 2） 系統(tǒng)的串聯(lián) series 命令實(shí)現(xiàn)兩個(gè)系統(tǒng)的串聯(lián)，示意圖如下： u1 y u2 y2 串聯(lián)后系統(tǒng)的傳遞函數(shù)為 121212( ) ( ) ( ) nng s g s g s dd?? ? ? ? 命令格式： [ n， d ] = series（ n1， d1， n2， d2） [ a， b， c， d ] = series（ a1， b1， c1， d1， a2， b2， c2， d2） 系統(tǒng) 1 系統(tǒng) 2 g1（ s） g2（ s） 電 力系統(tǒng)基礎(chǔ)課程設(shè)計(jì)說明書 18 3） 系統(tǒng)的反饋 feedback 命令實(shí)現(xiàn)兩個(gè)系統(tǒng)的反饋連接，示意圖如下： u1 + y1 177。 第三條命令產(chǎn)生 X至 y步長(zhǎng)為 n 的行矢量。 1) 在 MATLAB 中表示一個(gè)矢量要用方括號(hào)， 而列矢量的輸入只需在行矢量輸入格式基礎(chǔ)上加轉(zhuǎn)置符 (‘) 即可。 但這個(gè)命令不是存儲(chǔ)數(shù)據(jù)，而是存儲(chǔ)輸入與屏幕上輸出的內(nèi)容。7 8 9] a = 1 2 3 4 5 6 7 8 9 7） 需要顯示命令的計(jì)算結(jié)果時(shí)，則語(yǔ)句后面不加 “ ； ” 號(hào)，否則要加 “ ； ” 號(hào)。 2） 可用 ↑ 、 ↓ 鍵來(lái)重現(xiàn)已輸入的數(shù)據(jù)或命令。 數(shù)組的加、減與矩陣的加、減運(yùn)算完全相同。在新的版本中也加入了對(duì) C， FORTRAN， C++ ， JAVA 的支持。Y11= Y12= + Y13= + Y14=0 Y21= + Y22= Y23= + Y24= 0 + Y31= + Y32= + Y33= Y34= + Y41=0 Y42= 0 + Y43= + Y44= 電 力系統(tǒng)基礎(chǔ)課程設(shè)計(jì)說明書 11 簡(jiǎn)化雅可比矩陣 迭代中的雅克比矩陣 : 修正、迭代 給定 PQ 節(jié)點(diǎn)初值和各 節(jié)點(diǎn)電壓相角初值 V1= ， V2(0)=V3(0)=， V4= δ 2(0)=δ 3(0)=0， δ 4(0)=0 1 作第一次有功迭代，按公式計(jì)算節(jié)點(diǎn)有功功率不平衡量 迭代中的△ P： 電 力系統(tǒng)基礎(chǔ)課程設(shè)計(jì)說明書 12 2 做第一次無(wú)功迭代，按公式計(jì)算無(wú)功功率不平衡量，計(jì)算時(shí)電壓相角最新的修正值。 ① 雅可比矩陣中的諸元素都是節(jié)點(diǎn)電壓的函數(shù)，因此在迭代過程中，它們將隨著節(jié)點(diǎn)電壓的變化而不斷的變 化。 潮流計(jì)算算法 本題采用了題目要求的牛頓－拉夫遜潮流計(jì)算的方法。 第四章 潮流計(jì)算算法及手工計(jì)算 極坐標(biāo)下 PQ 法的算法 節(jié)點(diǎn)導(dǎo)納矩陣 Y 根據(jù)題目提供的各節(jié)點(diǎn)的參數(shù)，求得節(jié)點(diǎn)導(dǎo)納矩陣 Yii = ijyjiy ??0 yYikik ?? 簡(jiǎn)化雅可比矩陣 B/和 B// 通過上一步的導(dǎo)納矩陣，形成有功迭代和無(wú)功迭代的簡(jiǎn)化雅可比矩陣 B/和 B// 對(duì)雅可比矩陣進(jìn)行三角分解，形成因子表，為后面進(jìn)行 修正方程計(jì)算作好準(zhǔn)備。在解決電力系統(tǒng)潮流計(jì)算問題時(shí)，是以導(dǎo)納矩陣為基礎(chǔ)的，因此，只要我們能在迭代過程中盡可能保持方程式系數(shù)矩陣的稀疏性，就可以大大提高牛頓法潮流程序的效率。 在用數(shù)字解 算計(jì)算機(jī)解電力系統(tǒng)潮流問題的開始階段，普遍采取以節(jié)點(diǎn)導(dǎo)納矩陣為基礎(chǔ)的逐次代入法。 電力系統(tǒng)潮流計(jì)算是電力系統(tǒng)分析中的一種最基本的計(jì)算，是對(duì)復(fù)雜電力系統(tǒng)正常和故障條件下穩(wěn)態(tài)運(yùn)行狀態(tài)的計(jì)算。 電力系統(tǒng)潮流計(jì)算是研究電力系統(tǒng)穩(wěn)態(tài)運(yùn)行情況的一種計(jì)算，它根據(jù)給定的運(yùn)行條件及系統(tǒng)接線情況確定整個(gè)電力系統(tǒng)各部分的運(yùn)行狀態(tài)：各母線的電壓。同時(shí)，為了提高供電的可靠性以及資源利用的綜合經(jīng)濟(jì)性，又把許多分散的各種形式的發(fā)電站，通過送電線路和變電所聯(lián)系起來(lái)。目前 MATLAB 已成為國(guó)際控制界最流行、使 用最廣泛的語(yǔ)言了。 潮流計(jì)算是電力系統(tǒng)分析最基本的計(jì)算。按電壓的不同表示方法，牛頓 拉夫遜潮流計(jì)算分為直角坐標(biāo)形式和極坐標(biāo)形式兩種。 關(guān)鍵詞 ： 電力系統(tǒng)潮流計(jì)算 牛頓 — 拉夫遜法潮流計(jì)算 MATLAB 電 力系統(tǒng)基礎(chǔ)課程設(shè)計(jì)說明書 2 目 錄 一、 概述 設(shè)計(jì)目的與要求 ..................................................3 設(shè)計(jì)目的 .....................................................3 設(shè)計(jì)要求 .....................................................3 設(shè)計(jì)題目 .....................................................3 設(shè)計(jì)內(nèi)容 .....................................................3 二 電力系統(tǒng)潮流計(jì)算概述 ......................................4 電力系統(tǒng)簡(jiǎn)介 ...................................................4 潮流計(jì)算簡(jiǎn)介 ...................................................4 潮流計(jì)算的意義及其發(fā)展 ..................... ...................5 三 潮流計(jì)算設(shè)計(jì)題目 ...........................................6 潮流計(jì)算題目 ................................................. 6 對(duì)課題的分析及求解思路 ....................................... 7 四 潮流計(jì)算算法及手工計(jì)算 ....................................7 極坐標(biāo)下 PQ法的算法 ..........................................7 節(jié)點(diǎn)電壓方程 ..................................................8 節(jié)點(diǎn)導(dǎo)納矩陣 ..................................................9 導(dǎo)納矩陣在潮流計(jì)算 ............................................10 潮流計(jì)算的手工計(jì)算 ............................................12 五 Matlab 概述 ..................................................13 Matlab 簡(jiǎn)介 ................................................... 14 Matlab 的應(yīng)用 ..................................................14 矩陣的運(yùn)算 ...........