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