【正文】 for R=1:n C(m,m)=C(m,m)+y(m,R)。 C=input(39。)。而在例三 中 n則表示步長 . 我們應當 注意它們的區(qū)別。 格式 diary(文件名 )（擴展名） . . diary off 一般來說擴展名可取， m這樣就可在 MATLAB＼ BIN＼目錄下存入該文件。大寫字 母和小寫字母分別表示不同的變量。附加的工具箱（單獨提供的專用 MATLAB 函數(shù)集）擴展了 MATLAB 環(huán)境，以解決這些應用領(lǐng)域內(nèi)特定類型的問題。 ③ 由式（ 17）可以看出，當導納矩陣中非對角元素 ijY 為零時。 第三步，做第一次無功迭代，按公式計算無功功率不平衡量，計算時電壓相角最新的修正值。阻抗法改善了系統(tǒng)潮流計算問題的收斂性，解決了導納無法求解的一些系統(tǒng)的潮流計算，在 60年代獲得了廣泛的應用，阻抗法德主要缺點是占用計算機內(nèi)存大，每次迭代的計算量大?？煽啃院徒?jīng)濟性 。 關(guān)鍵詞 ： 電力系統(tǒng)潮流計算 牛頓 — 拉夫遜法潮流計算 MATLAB 電 力系統(tǒng)基礎課程設計說明書 2 目 錄 一、 概述 設計目的與要求 ..................................................3 設計目的 .....................................................3 設計要求 .....................................................3 設計題目 .....................................................3 設計內(nèi)容 .....................................................3 二 電力系統(tǒng)潮流計算概述 ......................................4 電力系統(tǒng)簡介 ...................................................4 潮流計算簡介 ...................................................4 潮流計算的意義及其發(fā)展 ..................... ...................5 三 潮流計算設計題目 ...........................................6 潮流計算題目 ................................................. 6 對課題的分析及求解思路 ....................................... 7 四 潮流計算算法及手工計算 ....................................7 極坐標下 PQ法的算法 ..........................................7 節(jié)點電壓方程 ..................................................8 節(jié)點導納矩陣 ..................................................9 導納矩陣在潮流計算 ............................................10 潮流計算的手工計算 ............................................12 五 Matlab 概述 ..................................................13 Matlab 簡介 ................................................... 14 Matlab 的應用 ..................................................14 矩陣的運算 .................................................... 14 與常數(shù)的運算 .................................................. 14 基本數(shù)學運算 ...................................................14 邏輯關(guān)系運算 ...................................................14 Matlab 中的一些命令 ............................................. 15 六 潮流計算流程圖及源程序 ......................................18 潮流計算流程圖 ..................................................18 潮流計算源程序圖 ................................................19 運行計算結(jié)果 ....................................................27 七 總結(jié) ...................................................29 八 參考文獻 ...............................................29 電 力系統(tǒng)基礎課程設計說明書 3 第一章 系統(tǒng)概述 設計目的與要求 設計目的 1. 掌握電力系統(tǒng)潮流計算的基本原理； 2. 掌握并能熟練運用一門計算機語言（ MATLAB 語言或 C 語言或 C++語言）； 3. 采用計算機語言對潮流計算進行計算機編程。 潮流計算是電力系統(tǒng)分析最基本的計算。同時，為了提高供電的可靠性以及資源利用的綜合經(jīng)濟性，又把許多分散的各種形式的發(fā)電站，通過送電線路和變電所聯(lián)系起來。 電力系統(tǒng)潮流計算是電力系統(tǒng)分析中的一種最基本的計算，是對復雜電力系統(tǒng)正常和故障條件下穩(wěn)態(tài)運行狀態(tài)的計算。在解決電力系統(tǒng)潮流計算問題時，是以導納矩陣為基礎的，因此，只要我們能在迭代過程中盡可能保持方程式系數(shù)矩陣的稀疏性，就可以大大提高牛頓法潮流程序的效率。 潮流計算算法 本題采用了題目要求的牛頓－拉夫遜潮流計算的方法。Y11= Y12= + Y13= + Y14=0 Y21= + Y22= Y23= + Y24= 0 + Y31= + Y32= + Y33= Y34= + Y41=0 Y42= 0 + Y43= + Y44= 電 力系統(tǒng)基礎課程設計說明書 11 簡化雅可比矩陣 迭代中的雅克比矩陣 : 修正、迭代 給定 PQ 節(jié)點初值和各 節(jié)點電壓相角初值 V1= ， V2(0)=V3(0)=， V4= δ 2(0)=δ 3(0)=0， δ 4(0)=0 1 作第一次有功迭代，按公式計算節(jié)點有功功率不平衡量 迭代中的△ P： 電 力系統(tǒng)基礎課程設計說明書 12 2 做第一次無功迭代，按公式計算無功功率不平衡量，計算時電壓相角最新的修正值。 數(shù)組的加、減與矩陣的加、減運算完全相同。7 8 9] a = 1 2 3 4 5 6 7 8 9 7） 需要顯示命令的計算結(jié)果時，則語句后面不加 “ ； ” 號，否則要加 “ ； ” 號。 1) 在 MATLAB 中表示一個矢量要用方括號， 而列矢量的輸入只需在行矢量輸入格式基礎上加轉(zhuǎn)置符 (‘) 即可。 命令格式： [ 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ù)為 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） 系統(tǒng) 1 系統(tǒng) 2 g1（ s） g2（ s） 電 力系統(tǒng)基礎課程設計說明書 18 3） 系統(tǒng)的反饋 feedback 命令實現(xiàn)兩個系統(tǒng)的反饋連接，示意圖如下： u1 + y1 177。請輸入 PQ 結(jié)點數(shù)： n2=39。)。 end end end for m=1:n for R=1:n if m~=R Z(m)=Z(m)+1/C(m,R)。 P=real(S)。