【導(dǎo)讀】中提出的各種問題。對規(guī)劃中的電力系統(tǒng)，通過潮流計算可以檢驗所提出。的電力系統(tǒng)規(guī)劃方案能否滿足各種運行方式的要求；對運行中的電力系統(tǒng)，除它自身的重要作用之外，網(wǎng)絡(luò)原始數(shù)據(jù)輸入工作大量且易于出錯。WINDOWS操作系統(tǒng)的潮流計算軟件。而采用MATLAB界面直觀，運行穩(wěn)定，計算。的，各發(fā)電站孤立運行。隨著工農(nóng)業(yè)生產(chǎn)和城市的發(fā)展，電能的需要量迅速增加，距離輸送給電力用戶。同時，為了提高供電的可靠性以及資源利用的綜合經(jīng)濟性，現(xiàn)代電力系統(tǒng)提出了“靈活交流輸電和新型直流輸電”的概念。力和電力系統(tǒng)的穩(wěn)態(tài)水平。新型直流輸電技術(shù)是指應(yīng)用現(xiàn)電力電子技術(shù)的最新成。果，改善和簡化變流站的造價等。網(wǎng)在預(yù)想操作預(yù)想下的電網(wǎng)的潮流分布以及校驗運行的可靠性。潮流是確定電力網(wǎng)咯運行狀態(tài)的基本因。素，潮流問題是研究電力系統(tǒng)穩(wěn)態(tài)問題的基礎(chǔ)和前提。

　　

【正文】 2 5 3 6 注意上面兩式的區(qū)別。 2) 下面三條命令可以產(chǎn)生一個行矢量 1 a=linspace(x,y,n) 2 a=logspace(x,y,n) 3 a=[x:n:y] 第一條命令可以在線性空間產(chǎn)生一個值在 10x至 10y之間間隔點數(shù)為 n 的行矢量 (一組數(shù)據(jù) )。 第二條命令可以在對數(shù)空間產(chǎn)生一個值在 x至 y 之間等間隔的行矢量 (一組數(shù)據(jù) )。其行矢量的起始值是 x, 終值為 y, 點數(shù)為 n。 第三條命令產(chǎn)生 X 至 y 步長為 n 的行矢量。 但是 , 三個命令之間存在差別，下面的例子可以說明這一點。 例一 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 通過上面三個例子可以看出例一 , 例二中 n 代表選取的點數(shù)。而在例三 中 n 則表示步長 . 我們應(yīng)當(dāng) 注意它們的區(qū)別。 3) 矩陣的加 , 減 , 乘 , 除等 , 和其它語言書寫一樣。但要注意的是在運算符 前面加有 (.)則表示是元素對元素的操作 . 4)以下是 常用 的 運算命令 運算命令名 功能 Angle 求復(fù)數(shù)的角 Min 求最小值 Max 求最大值 Sum 求和 Roots 求多項式的根 Poly 由多項式的根求多項式的系數(shù) Polyval 求給定點多項式的值 Polyder 多項式求導(dǎo) 在進行潮流分布計算時， 實際 上是 由多個簡單系統(tǒng)構(gòu)成 的 復(fù)雜系統(tǒng)， 在求節(jié)點導(dǎo)納矩陣時的要用到反饋的指令，所以在 MATLAB 中有下面幾種命令可以解決兩個系統(tǒng)間的連接問題。 1） 系統(tǒng)的并聯(lián) parallel 命令可以實現(xiàn)兩個系統(tǒng)的并聯(lián)。示意圖如下： u1 y1 u + y u2 y2 并聯(lián)后的系統(tǒng)傳遞函數(shù)表示式為： 1 2 2 112 12( ) ( ) ( ) n d n dg s g s g s dd? ? ?? ? ? ?其中 n d1和 n d2分別為 g1（ s）、 g2（ s）的傳遞函數(shù)分子、分母系數(shù)行矢量。 命令格式： [ 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)，示意圖如下： u1 y u2 y2 串聯(lián)后系統(tǒng)的傳遞函數(shù)為 系統(tǒng) 1 系統(tǒng) 2 g1（ s） g2（ s） 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） 3） 系統(tǒng)的反饋 feedback 命令實現(xiàn)兩個系統(tǒng)的反饋連接，示意圖如下： u1 + y1 177。 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連接的符號，缺省時默認為負（即 sign = 1）。 4） 系統(tǒng)的閉環(huán) cloop 命令可以將系統(tǒng)的輸出反饋到系統(tǒng)的輸入構(gòu)成閉環(huán)系統(tǒng)，示意圖如下： u y 177。 正、負反饋后閉環(huán)系統(tǒng)為： 111 1 1()() ()g s ngs m g s n d?? ? 命令格式： [ n， d ] = cloop （ n1， d1， sign） [ ac， bc， cc， dc] = cloop（ a， b， c， d， sign） 通過以上對 MATLAB 基本指令的了解，我們就可以對所求的電力系統(tǒng)網(wǎng)絡(luò)的節(jié)點導(dǎo)納矩陣進行畫編程框架圖。 第 六 章 潮流計算流程圖 g1（ s） g2（ s） g1（ s） 在上一章節(jié)，我們對 MATLAB 的基本指令有了初步的了解，我們就可以畫出程序編程框圖，按照程序編程框圖，我們可以通過 MATLAB 軟件進行編程，實現(xiàn)我們所求的節(jié)點導(dǎo)納矩陣。 程序編程框圖如下： 潮流計算源程序圖 節(jié)點導(dǎo)納矩陣 程序如下 : n=input(39。請輸入節(jié)點數(shù)： n=39。)。 n1=input(39。請輸入支路數(shù)： n1=39。)。 B=input(39。請輸入由支路參數(shù)形成的矩陣： B=39。)。 X=input(39。請輸入由節(jié)點號及其對地阻抗形成的矩陣： X=39。)。 Y=zeros(n)。 for i=1:n if X(i,2)~=0。 p=X(i,1)。 Y(p,p)=1./X(i,2)。 end end for i=1:n1 if B(i,6)==0 p=B(i,1)。q=B(i,2)。 else p=B(i,2)。q=B(i,1)。 end Y(p,q)=Y(p,q)1./(B(i,3)*B(i,5))。 Y(q,p)=Y(p,q)。 Y(q,q)=Y(q,q)+1./(B(i,3)*B(i,5)^2)+B(i,4)./2。 Y(p,p)=Y(p,p)+1./B(i,3)+B(i,4)./2。 end disp(39。導(dǎo)納矩陣 Y=39。)。 disp(Y) 按照程序調(diào)試的如下 結(jié) 果： 請輸 入 節(jié)點數(shù) ： n=6 請輸 入支路 數(shù) ： n1=7 請輸 入由支路 參數(shù) 形成的矩 陣 ： B=[1 2 0 1。2 5 + 0 1 0。5 3 + 0 1 0。4 3 0 1。4 6 + 0 1 0。1 6 + 0 1 0。1 4 + 0 1 0] 請輸 入由 節(jié)點號 及其 對 地阻抗形成的矩 陣 ： X=[1 0。2 0。3 0。4 0。5 0。6 0] 導(dǎo)納 矩 陣 Y= 總結(jié) 此次課程設(shè)計我在 MATLAB 編程，潮流計算， WORD 文檔的編輯方面均有所提高，但也暴露出了一些問題；理論知識的儲備還是不足，對 MATLAB 的性能和特點還不能有一個全面的把握，對 WORD 軟件也不是很熟練，相信通過以后的學(xué)習(xí)能彌補這些不足，從而達到一個新的層次。 此次課程設(shè)計的 完 成， 還 要感謝老師的幫助以及同學(xué)的協(xié)力合作，團結(jié)才能出成績。 參考文獻 [1] 孟祥萍 . 《電力系統(tǒng)分析》 ． 高等教育出版社 .2021 [1] 李維波 . 《 MATLAB 在電氣工程中應(yīng)用》 ． 中國電力出版社 .2021 [1] 邱曉燕 .劉天琪 《電力系統(tǒng)分析的計算機算法》 ． 中國電力出版社 .2021 [1] 吳天明 .趙新力 .劉建存 《 MATLAB 電力系統(tǒng)設(shè)計與分析》 ． 國防工業(yè)出版社 .2021