【正文】 參考文獻(xiàn) [1] 孟祥萍 . 《電力系統(tǒng)分析》 ． 高等教育出版社 .2021 [1] 李維波 . 《 MATLAB 在電氣工程中應(yīng)用》 ． 中國電力出版社 .2021 [1] 邱曉燕 .劉天琪 《電力系統(tǒng)分析的計(jì)算機(jī)算法》 ． 中國電力出版社 .2021 [1] 吳天明 .趙新力 .劉建存 《 MATLAB 電力系統(tǒng)設(shè)計(jì)與分析》 ． 國防工業(yè)出版社 .2021 。6 0] 導(dǎo)納 矩 陣 Y= 總結(jié) 此次課程設(shè)計(jì)我在 MATLAB 編程，潮流計(jì)算， WORD 文檔的編輯方面均有所提高，但也暴露出了一些問題；理論知識(shí)的儲(chǔ)備還是不足，對(duì) MATLAB 的性能和特點(diǎn)還不能有一個(gè)全面的把握，對(duì) WORD 軟件也不是很熟練，相信通過以后的學(xué)習(xí)能彌補(bǔ)這些不足，從而達(dá)到一個(gè)新的層次。4 0。2 0。1 6 + 0 1 0。4 3 0 1。2 5 + 0 1 0。)。 end disp(39。 Y(q,q)=Y(q,q)+1./(B(i,3)*B(i,5)^2)+B(i,4)./2。 end Y(p,q)=Y(p,q)1./(B(i,3)*B(i,5))。 else p=B(i,2)。 end end for i=1:n1 if B(i,6)==0 p=B(i,1)。 p=X(i,1)。 Y=zeros(n)。請(qǐng)輸入由節(jié)點(diǎn)號(hào)及其對(duì)地阻抗形成的矩陣： X=39。)。 B=input(39。請(qǐng)輸入支路數(shù)： n1=39。)。 程序編程框圖如下： 潮流計(jì)算源程序圖 節(jié)點(diǎn)導(dǎo)納矩陣 程序如下 : n=input(39。 正、負(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） 通過以上對(duì) MATLAB 基本指令的了解，我們就可以對(duì)所求的電力系統(tǒng)網(wǎng)絡(luò)的節(jié)點(diǎn)導(dǎo)納矩陣進(jìn)行畫編程框架圖。 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）。示意圖如下： 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ù)行矢量。但要注意的是在運(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)間的連接問題。而在例三 中 n 則表示步長 . 我們應(yīng)當(dāng) 注意它們的區(qū)別。 但是 , 三個(gè)命令之間存在差別，下面的例子可以說明這一點(diǎn)。其行矢量的起始值是 x, 終值為 y, 點(diǎn)數(shù)為 n。 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ù) )。4 5 6]39。 如 x=[1 2 3。 與此同時(shí)， MATLAB 具有強(qiáng)大的矩陣運(yùn)算功能 , 但由于我們 在求節(jié)點(diǎn)導(dǎo)納矩陣時(shí)用 的不多 , 因此這里我們只作簡單介紹 。 命令格式為 save 文件名 . 擴(kuò)展名 11） what 命令 該命令可以在當(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。 格式 diary(文件名 )（擴(kuò)展名） . . diary off 一般來說擴(kuò)展名可取， m這樣就可在 MATLAB＼ BIN＼目錄下存入該文件。它可以記錄下工作的過程。 9） diary 命令 使用該命令可以在窗口中以 ASCII 碼形式記錄所有的輸入和輸出 。7 8 9]。7 8 9] a=[1 2 3。 運(yùn)行下面兩種格式可以看出他們的區(qū)別 a=[1 2 3。4 5 6。 5） 矩陣的輸入要一行一行的進(jìn)行，每行各元素用空格或（，）分開，每行用（；）分開 。大寫字母和小寫字母分別表示不同的變量。用 ← 、 → 鍵來移動(dòng)光標(biāo)進(jìn)行修改。 . Matlab中的一些命令 1） 一般 MATLAB 命令格式為 [輸出參數(shù) 1，輸出參數(shù) 2， ??]= （命令名）（輸入?yún)?shù) 1，輸入?yún)?shù) 2， ?? ） 輸出參數(shù)用方括號(hào)，輸入?yún)?shù)用圓括號(hào)如果輸出參數(shù)只有一個(gè)可不使用括號(hào)。矩陣的冪運(yùn)算與數(shù)組的冪運(yùn)算有很大的區(qū)別。在數(shù)組運(yùn)算中有了“對(duì)應(yīng)關(guān)系”的規(guī)定，數(shù)組與常數(shù)之間的除法運(yùn)算沒有任何限制。而乘除法運(yùn)算有相當(dāng)大的區(qū)別，數(shù)組的乘除法是指兩同維數(shù)組對(duì)應(yīng)元素之間的乘除法，它們的運(yùn)算符為“ .*”和“ ./”或“ .\。這種運(yùn)算不同于前面講的數(shù)學(xué)運(yùn)算，為有所區(qū)別，我們稱之為數(shù)組運(yùn)算。但需注意進(jìn)行數(shù)除時(shí)，常數(shù)通常只能做除量。附加 的工具箱（單獨(dú)提供的專用 MATLAB 函數(shù)集）擴(kuò)展 了 MATLAB 環(huán)境，以解決這些應(yīng)用領(lǐng)域內(nèi)特定類型的問題?？梢灾苯诱{(diào)用 ,用戶也可以將自己編寫的實(shí)用程序?qū)氲?MATLAB 函數(shù)庫中方便自己以后調(diào)用，此外許多的 MATLAB 愛好者都編寫了一些經(jīng)典的程序，用戶可以直接進(jìn)行下載就可以用。 Matlab 的應(yīng)用 MATLAB 的基本數(shù)據(jù)單位是矩陣，它的指令表達(dá)式與數(shù)學(xué)、工程中常用的形式十分相似，故用 MATLAB 來解算問題要比用 C， FORTRAN 等語言完成相同的事情簡捷得多，并且 mathwork 也吸收了像 Maple 等軟件的優(yōu)點(diǎn) ,使MATLAB 成為一個(gè)強(qiáng)大的數(shù)學(xué)軟件。它在數(shù)學(xué)類科技應(yīng)用軟件中在 數(shù)值計(jì)算 方面首屈一指。它將數(shù)值分析、矩陣計(jì)算、科學(xué)數(shù)據(jù)可視化以及非線性動(dòng)態(tài)系統(tǒng)的建模和仿真等諸多強(qiáng)大功能集成在一個(gè)易于使用的視窗環(huán)境中，為科學(xué)研 究、工程設(shè)計(jì)以及必須進(jìn)行有效數(shù)值計(jì)算的眾多科學(xué)領(lǐng)域提供了一種全面的解決方案，并在很大程度上擺脫了傳統(tǒng)非交互式程序設(shè)計(jì)語言（如 C、 Fortran）的編輯模式，代表了當(dāng)今國際科學(xué)計(jì)算軟件的先進(jìn)水平。本節(jié)只是對(duì)導(dǎo)納矩陣在潮流分布的計(jì)算機(jī)算法一些簡單的描述，我們將在下一章對(duì)其進(jìn)行詳細(xì)講解。 運(yùn) 用 計(jì) 算機(jī) 進(jìn) 行潮流 計(jì) 算，一般要完成以下幾 個(gè) 步 驟 ： 建立 數(shù)學(xué) 模型、確立 計(jì) 算方法、制定 計(jì) 算機(jī)流程并 編 制程序、上機(jī) 計(jì) 算及對(duì)計(jì) 算 結(jié) 果 進(jìn) 行分析。目前 這 種方法多半與 牛頓型 算法配合使用以 彌補(bǔ)后 者的不足。高斯 塞德 爾 迭代法由于其 簡單 而在早期的潮流 計(jì) 算程序中得以采用。 在計(jì)算潮流分布時(shí)，我們必須先導(dǎo)出該網(wǎng)絡(luò)的導(dǎo)納矩陣，而進(jìn)行潮流計(jì)算時(shí)解非