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