【文章內(nèi)容簡介】 輯運(yùn)算是 MATLAB 中數(shù)組運(yùn)算所特有的一種運(yùn)算形式，也是幾乎所有的高級(jí)語言普遍適用的一種運(yùn)算。 Matlab中的一些命令 1） 一般 MATLAB命令格式為 信息工程學(xué)院課程設(shè)計(jì)（論文） 11 [輸出參數(shù) 1，輸出參數(shù) 2， ??]= （命令名）（輸入?yún)?shù) 1，輸入?yún)?shù) 2， ?? ） 輸出參數(shù)用方括號(hào)，輸入?yún)?shù)用圓括號(hào)如果輸出參數(shù)只有一個(gè)可不使用括號(hào)。 2） 可用 ↑ 、 ↓ 鍵來重現(xiàn)已輸入的數(shù)據(jù)或命令。用 ← 、 → 鍵來移動(dòng)光標(biāo)進(jìn)行修改。 3） 所有 MATLAB命令都用小寫字母。大寫字母和小寫字 母分別表示不同的變量。 4） 常用數(shù)有特定的名字，如 pi（ =） 、 Inf（ =∞ ）、 NaN則表示不定型求得的結(jié)果（如 0/0）。 5） 矩陣的輸入要一行一行的進(jìn)行，每行各元素用空格或（，）分開，每行用（；）分開 。 6） MATLAB書寫格式為 A=[1 2 3 ； 4 5 6 ； 7 8 9] 在 MATLAB中運(yùn)行如下程序可得到 A矩陣 a=[1 2 3。4 5 6。7 8 9] a = 1 2 3 4 5 6 7 8 9 7） 需要顯示命令的計(jì)算結(jié)果時(shí)，則語句后面不加 “ ； ” 號(hào)，否則要加 “ ； ” 號(hào)。 運(yùn)行下面兩種格式可以看出他們的區(qū)別 a=[1 2 3。4 5 6。7 8 9] a=[1 2 3。4 5 6。7 8 9]。 a = （不顯示計(jì)算結(jié)果） 1 2 3 4 5 6 8） 當(dāng)輸入語句過長需要換行時(shí)，應(yīng)加上 “ ? ” 后再回車，則可連續(xù)輸入。 9） diary命令 使用該命令可以在窗口中以 ASCII碼形式記錄所有的輸入和輸出 。 但這個(gè)命令不是存儲(chǔ)數(shù)據(jù)，而是存儲(chǔ)輸入與屏幕上輸出的內(nèi)容。它可以記錄下工作的過程 。在每個(gè)工作過程之前使用該命令，工作結(jié)束后使用 diary off 則能將整個(gè)工作過程記錄下來。 格式 diary(文件名 )（擴(kuò)展名） . . diary off 一般來說擴(kuò)展名可取， m這樣就可在 MATLAB＼ BIN＼目錄下存入該文件。 10） save 命令 該命令存儲(chǔ)定義的變量或演算結(jié)果，也可以用來存儲(chǔ)指定的變量。 命令格式為 save 文件名 . 擴(kuò)展名 11） what 命 令 信息工程學(xué)院課程設(shè)計(jì)（論文） 12 該命令可以在當(dāng)目錄下顯示 MATLAB文件和 MAT數(shù)據(jù)文件 12） dir 命令 顯示當(dāng)前目錄下的所有文件 . 13） clear命令 14） [d1,d2,d3,..]=size(a) 求矩 陣 的大小， 對(duì) m*n二 維 矩 陣 ，第一 個(gè)為 行 數(shù) m,第二個(gè)為 列 數(shù) n。 如果輸入 calear a b c ，則表示清除工作空間中指定變量 a,b,c；如果僅僅輸入 calear命令，則清除整個(gè)工作空間。 與此同時(shí)， MATLAB 具有強(qiáng)大的矩陣運(yùn)算功能 , 但由于我們 在求節(jié)點(diǎn)導(dǎo) 納矩陣時(shí)用 的不多 , 因此這里我們只作簡單介紹 。 1) 在 MATLAB中表示一個(gè)矢量要用方括號(hào)， 而列矢量的輸入只需在行矢量輸入格式基礎(chǔ)上加轉(zhuǎn)置符 (‘) 即可。 如 x=[1 2 3。4 5 6] x = 1 2 3 4 5 6 而 x=[1 2 3。4 5 6]39。 (加轉(zhuǎn)置符 ) x = 1 4 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] 第一條命令可以在線性空間產(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步長為 n的行矢量。 但是 , 三個(gè)命令之間存在差別，下面的例子可以說明這一點(diǎn)。 例一 x=logspace(0,5,6) x = 信息工程學(xué)院課程設(shè)計(jì)（論文） 13 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 通過上面三個(gè)例子可以看出例一 , 例二中 n代表選取的點(diǎn)數(shù)。而在例三 中 n則表示步長 . 我們應(yīng)當(dāng) 注意它們的區(qū)別。 3) 矩陣的加 , 減 , 乘 , 除等 , 和其它語言書寫一樣。但要注意的是在運(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è)簡單系統(tǒng)構(gòu)成 的 復(fù)雜系統(tǒng)， 在求節(jié)點(diǎn)導(dǎo)納矩陣時(shí)的要用到反饋的指令，所以在 MATLAB中有下面幾種命令可以解決兩個(gè)系統(tǒng)間的連 接問題。 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ù)行矢量。 系統(tǒng) 1 系統(tǒng) 2 信息工程學(xué)院課程設(shè)計(jì)（論文） 14 命令格式： [ 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） 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?? ? g1（ s） g2（ s） g1（ s） g2（ s） g1（ s） 信息工程學(xué)院課程設(shè)計(jì)（論文） 15 命令格式： [ n， d ] = cloop （ n1， d1， sign） [ ac， bc， cc， dc] = cloop（ a， b， c， d， sign） 通過以上對(duì) MATLAB基本指令的了解，我們就可以對(duì)所求的電力系統(tǒng)網(wǎng)絡(luò)的節(jié)點(diǎn)導(dǎo)納矩陣進(jìn)行畫編程框架圖。 潮流計(jì)算流程圖 在上一節(jié)，我們對(duì) MATLAB的基本指令有了初步的了解，我們就可以畫出程序編程框圖，按照程序編程框圖，我們可以通過 MATLAB軟件進(jìn)行編程，實(shí)現(xiàn)我們所求的節(jié)點(diǎn)導(dǎo)納矩陣。 程序編程框圖如下 ： 信息工程學(xué)院課程設(shè)計(jì)（論文） 16 潮流計(jì)算源程序圖及結(jié)果 節(jié)點(diǎn)導(dǎo)納矩陣程序如下 : n=input(39。請(qǐng)輸入節(jié)點(diǎn)數(shù)： n=39。)。 信息工程學(xué)院課程設(shè)計(jì)（論文） 17 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。