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