【正文】 0 0 0 0] 請(qǐng)輸入結(jié)點(diǎn)電壓矩陣： U=[1+0i 1+0i +0i +0i] 請(qǐng)輸入各結(jié)點(diǎn)的功率： S=[ 0] 結(jié)點(diǎn)導(dǎo)納矩陣 : + + 0 + + 0 + + + + 0 0 + + 迭代中的雅克比矩陣 : 電 力系統(tǒng)基礎(chǔ)課程設(shè)計(jì)說(shuō)明書 28 迭代中的△ P： 迭代中的△ Q： 迭代中相角 : 迭代中電壓的模 : 平衡結(jié)點(diǎn)的功率 : 電 力系統(tǒng)基礎(chǔ)課程設(shè)計(jì)說(shuō)明書 29 + 全部線路功率分布 : 0 + 0 0 0 0 + + 0 第七章 總 結(jié) 此次課程設(shè)計(jì)首先讓我明白了要使電力系統(tǒng)運(yùn)行的穩(wěn)定，必須經(jīng)過(guò)精密的設(shè)計(jì)和計(jì)算。0 0 0 。平衡結(jié)點(diǎn)的功率 :39。迭代中相角 :39。迭代中的△ P： 39。 end 電 力系統(tǒng)基礎(chǔ)課程設(shè)計(jì)說(shuō)明書 26 U=b+i*c。 UU=UU1*UD2。 PQ=[PP39。 else N(m,R)=U1(m)*U1(R)*(G(m,R)*cos(O(m)O(R))+B(m,R)*sin(O(m)O(R)))。 (m)=Q(m)QT1(m)。 Q=imag(S)。)。 end 電 力系統(tǒng)基礎(chǔ)課程設(shè)計(jì)說(shuō)明書 21 end end for m=1:n for R=1:n if m==R Y(m,m)=C(m,m)+Z(m)。L=zeros(n2,n2)。 U=input(39。請(qǐng)輸入變比矩陣看 :K=39。)。 clear n=input(39。 y2 u2 連接后系統(tǒng)的傳遞函數(shù)表示為： 1 1 22 1 2 1 2()() 1 ( )g s n dgs g s d d n n??? ?? 命令格式： [ n， d ] = feedback（ n1， d1， n2， d2） 或： [ n， d ] = feedback（ n1， d1， n2， d2， sign） [ a， b， c， d ] = feedback（ a1， b1， c1， d1， a2， b2， c2， d2， sign） 其中 sign 是指示 y2到 u1連接的符號(hào)，缺省時(shí)默認(rèn)為負(fù)（即 sign = 1）。 但是 , 三個(gè)命令之間存在差別，下面的例子可以說(shuō)明這一點(diǎn)。 如 x=[1 2 3。它可以記錄下 工作的過(guò)程。 運(yùn)行下面兩種格式可以看出他們的區(qū)別 a=[1 2 3。用 ← 、 → 鍵來(lái)移動(dòng)光標(biāo)進(jìn)行修改。而乘除法運(yùn)算有相當(dāng)大的區(qū)別，數(shù)組的乘除法是指兩同維數(shù)組對(duì)應(yīng)元素之間的乘除法，它們的運(yùn)算符為“ .*”和“ ./”或“ .\?？梢灾苯诱{(diào)用 ,用戶也可以將自己編寫的實(shí)用程序?qū)氲?MATLAB 函數(shù)庫(kù)中方便自己以后調(diào)用，此外許 多的 MATLAB 愛(ài)好者都編寫了一些經(jīng)典的程序，用戶可以直接進(jìn)行下載就可以用。 迭代中的△ Q： 解修正方程式，可得各節(jié)點(diǎn)電壓幅值的修正量為 迭代中電壓的模 : 到這里為止，第一輪有功迭代和無(wú)功迭代便做完了。 ② 雅 可 比 矩陣 具 有結(jié) 構(gòu)對(duì) 稱 性， 數(shù) 據(jù)不 對(duì) 稱。 牛頓 拉夫遜法 潮流計(jì)算的公式。 修正和迭代 第一步，給定 PQ節(jié)點(diǎn)初值和各節(jié)點(diǎn)電壓相角初值。自從 60 年代中期，牛頓法中利用了最佳順序消去法以后，牛頓法在收斂性。這個(gè)方法的原理比較簡(jiǎn)單，要求的數(shù)字計(jì)算機(jī)內(nèi)存量比較差下，適應(yīng) 50 年代電子計(jì)算機(jī)制造水平和當(dāng)時(shí)電力 系統(tǒng)理論水平，但它的收斂性較差，當(dāng)系統(tǒng)規(guī)模變大時(shí)，迭代次數(shù)急劇上升，在計(jì)算中往往出現(xiàn)迭代不收斂的情況。潮流計(jì)算的目標(biāo)是求取電力系統(tǒng)在給定運(yùn)行狀態(tài)的計(jì)算，即節(jié)點(diǎn)電壓和功率分布，用以檢查系統(tǒng)各元件是否過(guò)負(fù)荷。各元件中流過(guò)的功率，系統(tǒng)的功率損耗等等。這種由發(fā)電機(jī)、升壓和降壓 變電所，送電線路以及用電設(shè)備有機(jī)連接起來(lái)的整體，即稱為電力系統(tǒng)。它的強(qiáng)大的矩陣處理功能給電力系統(tǒng)的分析、計(jì)算帶來(lái)很多方便，而且采用 MATLAB 界面直觀，運(yùn)行穩(wěn)定，計(jì)算準(zhǔn)確。除它自身的重要作用之外 ， 潮流計(jì)算還是網(wǎng)損計(jì)算、靜態(tài)安全 分析、暫態(tài)穩(wěn)定計(jì)算、小干擾靜態(tài)穩(wěn)定計(jì)算、短路計(jì)算、靜態(tài)和動(dòng)態(tài)等值計(jì)算的基礎(chǔ)。 實(shí)際電力系統(tǒng)的潮流計(jì)算主要采用牛頓 拉夫遜法。所以本次課程設(shè)計(jì)程序設(shè)計(jì)采用 MATLAB 計(jì)算。 現(xiàn)代電力系統(tǒng)提出了“靈活交流輸電和新型直流輸電”的概念。在電力系統(tǒng)規(guī)劃的設(shè)計(jì)和現(xiàn)有電力系統(tǒng)運(yùn)行方式的研究中，都需要利用潮流計(jì)算來(lái)定量的分析比較供電方案或運(yùn)行方式的合理性。各點(diǎn)電壓是否滿足要求，功率的分布和分配是否合理以及功率損耗等。這就迫使電力系統(tǒng)的計(jì)算人員轉(zhuǎn)向以阻抗矩陣為基礎(chǔ)的逐次代入法。內(nèi)存要求。 第二步，作第一次有功迭代，按公式計(jì)算節(jié)點(diǎn)有功功率不平衡量。把牛頓法用于潮流計(jì)算，采用極坐標(biāo)形式表示的如式（ 13）所示的形式。 如 非對(duì) 角 ij jiHH? ，ij ij i ij iH B e G f??， ji ij j ij jH B e G f??。 3 按公式計(jì)算平衡節(jié)點(diǎn)功率，得： P1+jQ1= + 經(jīng)過(guò)四輪迭代，節(jié)點(diǎn)不平衡功率也下降到 105以下，迭代到此結(jié)束。 MALAB 產(chǎn)品族可以用來(lái)進(jìn)行以下各種工作： ● 數(shù)值分析 ● 數(shù)值和符號(hào)計(jì)算 ● 工程與科學(xué)繪圖 ● 控制系統(tǒng)的設(shè)計(jì)與仿真 ● 數(shù)字圖像處理 技術(shù) ● 數(shù)字信號(hào)處理 技術(shù) ● 通訊系統(tǒng)設(shè)計(jì)與仿真 ● 財(cái)務(wù)與金融工程 電 力系統(tǒng)基礎(chǔ)課程設(shè)計(jì)說(shuō)明書 14 MATLAB 的應(yīng)用范圍非常廣，包括信號(hào)和圖像處理、通訊、控制系統(tǒng)設(shè)計(jì)、測(cè)試和測(cè)量、財(cái)務(wù)建模和分析以及計(jì)算生物學(xué)等眾多應(yīng)用領(lǐng)域。”前面講過(guò)常數(shù)與矩陣的除法運(yùn)算中常數(shù)只能做除數(shù)。 3） 所有 MATLAB 命令都用小寫字母。4 5 6。在每個(gè)工作過(guò)程之前使用該命令，工作結(jié)束后使用 diary off 則能將整個(gè)工作過(guò)程記錄下來(lái)。4 5 6] x = 1 2 3 4 5 6 而 x=[1 2 3。 例一 x=logspace(0,5,6) x = 1 10 100 1000 10000 100000 例二 x=linspace(0,10,11) x = 0 1 2 3 4 5 6 7 8 9 10 例三 x=[0:1:10] x = 0 1 2 3 4 5 6 7 8 9 10 通過(guò)上面三個(gè)例子可以看出例一 , 例二中 n代表選取的點(diǎn)數(shù)。 4） 系統(tǒng)的閉環(huán) cloop 命令可以將系統(tǒng)的輸出反饋到系統(tǒng)的輸入構(gòu)成閉環(huán)系統(tǒng)，示意圖如下： u y 177。請(qǐng)輸入結(jié)點(diǎn)數(shù)： n=39。 isb=input(39。)。請(qǐng)輸入結(jié)點(diǎn)電壓矩陣： U=39。QT1=zeros(1,n1+n2)。 else Y(m,R)=1/C(m,R)。 G=real(Y)。 while PRpr 電 力系統(tǒng)基礎(chǔ)課程設(shè)計(jì)說(shuō)明書 22 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)^