【正文】 0MW。π，表 的 功角是角度的形式表示的。將 matlab 程序計(jì)算的結(jié)果整理如下： 1）各節(jié)點(diǎn)的電壓和注入功率如表 。 8） 基于 MATLAB 平臺(tái)，得出計(jì)算結(jié)果，采用誤差精度為 ε= 105。 7）當(dāng) 和 都滿足要求時(shí)，計(jì)算結(jié)束。計(jì)算不平衡功率 ΔQ 判斷 ΔQ 是否小于可接受誤差。 5） 初始化各節(jié)點(diǎn)電壓（幅值和相角）的存儲(chǔ)矩陣和各節(jié)點(diǎn)注入功率（實(shí)部和虛部）的存儲(chǔ)矩陣，以存放迭代過程中的中間值和輸出最終值； 6） 構(gòu)造各節(jié)點(diǎn)的潮流方程，進(jìn)入潮流計(jì)算程序 ， 計(jì)算不平衡功率 ΔP， 判斷 ΔP 是否小于可接受誤差。 基于 Matlab 的 PQ 分解法程序 編程思想如下： 1） 根據(jù)支路信息矩陣 B1 和節(jié)點(diǎn)信息矩陣 B2，先形成待求解電網(wǎng)的節(jié)點(diǎn)導(dǎo)納矩陣： 2） 先對(duì)空矩陣加入對(duì)角線元素，即先加入每個(gè)節(jié)點(diǎn)的對(duì)地導(dǎo)納 y0； 3） 根據(jù)支路信息矩陣 B1，對(duì)導(dǎo)納矩陣再加入互導(dǎo)納 Ypq； 4） 將互導(dǎo)納 Ypq 加入到導(dǎo)納矩陣的對(duì)角線元素，形成自導(dǎo)納。待求解網(wǎng)架中， 1 號(hào)節(jié)點(diǎn)為平衡節(jié)點(diǎn)，2~6 號(hào)節(jié)點(diǎn)為 PQ 節(jié)點(diǎn)， 7 號(hào)節(jié)點(diǎn)（即母線 G）為 PV 節(jié)點(diǎn)。本題形成的 B1 矩陣見附錄。首先利用線路的標(biāo)幺值參數(shù)整合成矩陣，使其包含線路的所有信息，再形成節(jié)點(diǎn)導(dǎo)納矩陣。本文編寫了基于 Matlab 平臺(tái)的 PQ 分解法計(jì)算程序，對(duì)電網(wǎng)進(jìn)行了潮流計(jì)算，得出了系統(tǒng)穩(wěn)定時(shí)各節(jié)點(diǎn)的電壓，以及線路損耗等，分析了系統(tǒng)穩(wěn)定運(yùn)行時(shí)各線路的輸送功率。 表 編號(hào) 1 2 3 4 5 6 G 節(jié)點(diǎn)類型 平衡節(jié)點(diǎn) PQ 節(jié)點(diǎn) PQ 節(jié)點(diǎn) PQ 節(jié)點(diǎn) PQ 節(jié)點(diǎn) PQ 節(jié)點(diǎn) PV 節(jié)點(diǎn) 有功功率 P/MW － 180 350 260 250 150 837 無功功率 Q/MVar － 100 180 120 120 70 － 電壓 V/kV 230 － － － － － 242 采用 Matlab 計(jì)算電網(wǎng)潮流 通過對(duì)待求解電網(wǎng)的等效線路參數(shù)的計(jì)算與標(biāo)幺化， 形成節(jié)點(diǎn)導(dǎo)納矩陣。節(jié)點(diǎn) 6 均為 PQ 節(jié)點(diǎn)，并且?guī)в幸欢ǖ呢?fù)荷。各線路參數(shù)的標(biāo)幺化計(jì)算如下： 線路 L12 長(zhǎng)度 15km，則： Z12= 15 （ +） =（ +） Ω B12=2 2π 50 15 106= 104S 標(biāo)幺值： Z12*=（ +） 100/2202= 104+ B12*= 104 2202/100= 線路 L13 長(zhǎng)度 18km，則： Z13= 18 （ +） =（ +） Ω B13=2 2π 50 18 106= 104S 標(biāo)幺值： Z13*=（ +） 100/2202= 103+ B12*= 104 2202/100= 線路 L24 長(zhǎng)度 12km，則： Z24= 12 （ +） =（ +） Ω B24=2 2π 50 12 106= 105S 標(biāo)幺值： Z24*=（ +） 100/2202= 104+ B24*= 105 2202/100= 線路 L45 長(zhǎng)度 10km，則： Z45= 10 （ +） =（ +） Ω B45=2 2π 50 10 106= 105S 標(biāo)幺值： Z45*=（ +） 100/2202= 104+ B45*= 105 2202/100= 線路 L46 長(zhǎng)度 12km，則： Z46=12 （ +） =（ +） Ω B46=2π 50 12 106= 105S 標(biāo)幺值： Z46*=（ +） 100/2202= 103+ B46*= 105 2202/100= 線路 L3G長(zhǎng)度 16km，則： Z3G= 16 （ +） =（ +） Ω B3G=2 2π 50 16 106= 104S 標(biāo)幺值： Z3G*=（ +） 100/2202= 104+ B3G*= 104 2202/100= 線路 L5G長(zhǎng)度 25km，則： Z5G= 25 （ +） =（ +） Ω B5G=2 2π 50 25 106= 104S 標(biāo)幺值： Z5G*=（ +） 100/2202= 103+ B5G*= 104 2202/100= 13 線路 L6G長(zhǎng)度 12km，則： Z6G= 12 （ +） =（ +） Ω B6G=2 2π 50 12 106= 105S 標(biāo)幺值： Z6G*=（ +） 100/2202= 104+ B6G*= 105 2202/100= 綜上，電網(wǎng)的線路參數(shù)的有名值和標(biāo)幺值如下表所示： 表 線路名 長(zhǎng)度 /km 有名值 標(biāo)幺值 阻抗 Z/Ω 導(dǎo)納 B/S 阻抗 Z*/Ω 導(dǎo)納 B*/S L12 2 15 + 104 104+ L13 2 18 + 104 103+ L24 2 12 + 105 104+ L45 2 10 + 105 104+ L46 12 + 105 103+ L3G 2 16 + 104 104+ L5G 2 25 + 104 103+ L6G 2 12 + 105 104+ 線路零序參數(shù)的標(biāo)幺值計(jì)算與正序相似，如表 。 參考資料以及說明 《電力系統(tǒng)分析》（上、下冊(cè)）華中科技大學(xué)出版 《發(fā)電廠電氣部分》高等學(xué)校教材 《電網(wǎng)調(diào)度運(yùn)用技術(shù)》東北大學(xué)出版社 PowerWorld 15 使用手冊(cè) 基于 MATLAB/Simulink 的系統(tǒng)仿真技術(shù)與應(yīng)用 本畢業(yè)設(shè)計(jì)（論文）任務(wù)書于 20xx 年 12 月 23 日發(fā)出，應(yīng)于 20xx 年 1 月 3 日前完成，然后提交課程考試委員會(huì)進(jìn)行答辯。 =1e8。 例如，求解如下微分方程： 12212sin( )xx?? 初值為 x1=1,x2=0 從 0s積分到 3s，步長(zhǎng)： 則，首先定義函數(shù) myfunc，計(jì)算微分方程右端項(xiàng)的值： function dx=myfunc(x) dx=[x(2) 2*sin(x(1) ] 。 說明： odefun：列向量 1*n，通過函數(shù)計(jì)算柄輸出的微分方程的右端項(xiàng)； t1,tf: 分別制定積分的時(shí)間起點(diǎn)和終點(diǎn)； x0：列向量 1*n，狀態(tài)變量初值 options：：微分優(yōu)化參數(shù)，是一個(gè)結(jié)構(gòu)體，使用 odeset可以設(shè)置其具體參數(shù)，詳細(xì)內(nèi)容查看幫助。這里可采用 ode45編程（大家也可選擇和對(duì)比其它函數(shù)，不同編號(hào)采用的數(shù)值積分算法不同）。 附注： ODE 函數(shù)說明 Matlab提供了一階常微分方程組求解的系列函數(shù)： ode**。 編制課程設(shè)計(jì)報(bào)告。 發(fā)電機(jī)采用三階模型，用 PowerWorld 作為分析工具，對(duì)問題 4 的兩種故障方案進(jìn)行穩(wěn)定計(jì)算，給出搖擺曲線，并計(jì)算故障的極限切除時(shí)間。 試手工計(jì)算序網(wǎng)以及用網(wǎng)絡(luò)變換法求解轉(zhuǎn)移阻抗； 針對(duì)問題 4，用 Matlab 編制穩(wěn)定計(jì)算程序（三臺(tái)機(jī)可 并聯(lián)等值成一臺(tái)機(jī)），發(fā)電機(jī)采用二階經(jīng)典模型（注：用 ode45 函數(shù)既可求解），用軟件計(jì)算出搖擺曲線，要求輸出發(fā)電機(jī)功角，角速度。 時(shí)三相重合閘。 3. 考慮該電廠開機(jī)三臺(tái)，均為額定容量輸出的方式，用編制的程序計(jì)算潮流。負(fù)荷采用恒阻抗模型，負(fù)序阻抗標(biāo)幺制取。發(fā)電廠廠用電均按出力的 7%考慮。不計(jì)內(nèi)阻和空載損耗。電廠 220kV 側(cè)采用單母分段接線，正常運(yùn)行時(shí)分段開關(guān)閉合。C 長(zhǎng)期運(yùn)行允許的最大電流： 1190A。F/km； 零序參數(shù)： r0 = , x0 = , C0 = 181。變電站參數(shù)和負(fù)荷如下表： 編號(hào) 類型 220kV 最大負(fù)荷， MVA 1 500kV 站 平衡節(jié)點(diǎn) 2 220kV 站 180+j110 3 220kV 站 350+j180 4 220kV 站 260+j120 5 220kV 站 250+j120 6 220kV 站 150+j70 各線路長(zhǎng)度如圖所示。圖表整潔，布局合理，文字注釋必須使用工程字書寫，不準(zhǔn)用徒手畫 3）畢業(yè)論文須用 A4 單面打印，論文 50 頁以上的雙面打印 4）圖表應(yīng)繪制于無格子的頁面上 5）軟件工程類課題應(yīng)有程序清單，并提供電子文檔 1）設(shè)計(jì)（論文） 2）附件：按照任務(wù)書、開題報(bào)告、外文譯文、譯文原文（復(fù)印件）次序裝訂 3）其它 5 目 錄 課 程 設(shè)計(jì)題 目 C： 電 力系 統(tǒng)穩(wěn) 定分析和 計(jì) 算 ............................................................................... 6 算 .......................................................................................................................... 11 線 路 參數(shù) 的 計(jì) 算和 標(biāo) 幺化 ................................................................................................ 11 節(jié) 點(diǎn)信息 統(tǒng)計(jì) .................................................................................................................. 14 潮流 計(jì) 算 .......................................................................................................................... 14 采用 Matlab 計(jì) 算 電網(wǎng) 潮流 .............................................................................................. 14 網(wǎng)絡(luò) 信息 處 理 ................................................................................................................................................... 14 基于 Matlab 的 PQ 分解法程序 ................................................................................................................. 15 PQ 分解法潮流 計(jì) 算的 結(jié) 果分析 ............................................................................................................. 15 采用 Power World 計(jì) 算潮流 ............................................................................................ 16 Power world 仿 真軟 件 簡(jiǎn) 介 ........................................................................................................................... 16 基于 Power World 軟 件 構(gòu) 建潮流模型 ..................................................................................................... 17 采用 power world 仿 真軟 件的潮流 計(jì) 算 結(jié) 果 .................................................................................