【正文】 參考文獻(xiàn)1；中國(guó)電力出版社，2022 語(yǔ)言與應(yīng)用技術(shù)，北京；國(guó)防工業(yè)出版社 2022陳澤, MATLAB 在科學(xué)計(jì)算中的應(yīng)用，北京；電子工業(yè)出版社，2022：中國(guó)電力出版社，2022。最后感謝老師對(duì)我的悉心教導(dǎo)。什么是幸福？忙碌著，收獲著，簡(jiǎn)單著，快樂(lè)著，這就是幸福。早上 8 點(diǎn)多到教室，12 點(diǎn)下課了還渾然不知，很多次，晚上 10 點(diǎn)半時(shí)都是門(mén)衛(wèi)打著手電過(guò)來(lái)催我們走。軍令如山。俗話說(shuō)得好，不逼自己，你永遠(yuǎn)不知道自己有多強(qiáng)。其次，我們對(duì)穩(wěn)態(tài)潮流分析，掌握的更加深刻，不單單可以進(jìn)行簡(jiǎn)單的計(jì)算，還可以憑借計(jì)算機(jī)，對(duì)各種復(fù)雜網(wǎng)絡(luò)進(jìn)行潮流計(jì)算。一周內(nèi)，我們初步的掌握了對(duì) MATLAB 的應(yīng)用。首先，此次課程設(shè)計(jì)，是我們大學(xué)以來(lái)第一次嚴(yán)格進(jìn)行的課程設(shè)計(jì)，這是我們進(jìn)入社會(huì)前的一次小小的磨練，為我們今后參加工作打下一點(diǎn)基礎(chǔ)。各節(jié)點(diǎn)功率 sij 為：39。)S1=U(1)*sum(conj(I)) %電流for m=1:5 for n=1:5 S(m,n)=U(m)*(conj(U(m))conj(U(n)))*conj(Y(m,n))。enddisp(39。) u(m)=u(m)+ut(m1) endfor n=1:5 U(n)=u(n)*(cos(dt(n))+j*sin(dt(n)))。) utdisp( 39。disp( 39。 qq2=qq1(1:51,:)。各節(jié)點(diǎn)無(wú)功率不平衡量為：39。) dt(m)=dt(m)+dtt(m1)endfor m=2:5for n=1:5 qt(n)=u(m)*u(n)*(G(m,n)*sin(dt(m)dt(n))B(m,n)*cos(dt(m)dt(n)))。dtt=flipud(rot90(dtt))disp(39。dtt=(b1*dp2/u4)*u2。 %— ( )=U139。上下對(duì)稱(chēng)u4=u3(1:51,:)。 %取逆u3=flipud(rot90(u(2:5)))。) dp(m)=p(m)sum(pt) % = ??0iP?ip????njjiU10j????00sincosjjijijBG??endfor m=2:5u1=diag(u)。 % ????njjiU10j????00sincosjjijijBG??內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計(jì)說(shuō)明書(shū)20enddisp(39。dt(1)=dim。 %迭代次數(shù)，while wucha(1)pr。endk=0。請(qǐng)輸入 PQ 節(jié)點(diǎn)的無(wú)功功率： q=39。)。 %電壓實(shí)部初值endfor m=2:5p(m)=input(39。for m=2:5dt(m)=0。請(qǐng)輸入平衡節(jié)點(diǎn) 1 電壓初值的虛部，Im=39。)。ure=input(39。 b1=inv(B1)。 %real 實(shí)部，image 虛部B1=B(2:5, 2:5)。G=real(Y)。請(qǐng)輸入誤差精度：pr=39。)。1Y=input(39。平衡節(jié)點(diǎn)為：39。)。內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計(jì)說(shuō)明書(shū)18第三章 程序設(shè)計(jì) 設(shè)計(jì)流程3．2 程序設(shè)計(jì)針對(duì)兩機(jī)五節(jié)點(diǎn)網(wǎng)絡(luò)潮流計(jì)算基于 MATLAB 的程序如下：Clc %清除當(dāng)前屏幕上顯示的內(nèi)容，但并不清楚工作空間的數(shù)據(jù)Clear %清除工作空間內(nèi)的所有變量disp(39。 9265 和線路功率~1S~ijS 迭代收斂后，就可以計(jì)算平衡節(jié)點(diǎn)功率和線路功率。U?? = ； =50?內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計(jì)說(shuō)明書(shū)16表 1 迭代過(guò)程中各節(jié)點(diǎn)功率的不平衡量K )(2kP???kQ2)(3kP???kQ3)(4kP???kQ4)(5kP???kQ50 1 2 3 4 0,000024 5 6 表 2 迭代過(guò)程中各節(jié)點(diǎn)電壓的修正量k ?)(k???kU2?)(3k???kU3?)(4k???kU4?)(5k???kU50 1 2 3 4 5 6表 3 迭代過(guò)程中各節(jié)點(diǎn)電壓k )(2k???kU2)(3k???kU3)(4k???kU4)(5k???kU50 1 0,108571 2 0,089348 0,103974 3 4 內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計(jì)說(shuō)明書(shū)175 6 0,083906 = ； =每次迭代所得適于表 1~表3。)?B??0可得=—； =—； =—； =—?)0(2??)0(3??)0(4??)0(5?內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計(jì)說(shuō)明書(shū)15；所以=—； =—； =—； =—)1(2?)1(3?)1(4?)1(5?； iQ?按下式計(jì)算各節(jié)點(diǎn)無(wú)功功率不平衡量 = ??0ii????njjiU10j????11cossinijijjjBG??=—[—（—0）—（—0）]??02Q?—[+]—[（+）—（+）] —[（—+0,091810）—（0,036952+）]—[（—+）—（+）]=相似地可得 =—； =—； =—；??03Q???04Q???05Q? iU 由方程式 — ( )= U139。[ 0。他們就由導(dǎo)納矩陣的虛部部分中除第一行第一列外的各個(gè)元素所組成，即Y=[ 5+15j + 0 0 5+15j +5j +5j + + +5j 10+30j 0 0 +5j 10+30j +0 + 0 + ]Bˊ=B〞=[ 0 內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計(jì)說(shuō)明書(shū)14 0 ] 由此可見(jiàn)，網(wǎng)絡(luò)的節(jié)點(diǎn)導(dǎo)納矩陣雖為奇異矩陣，但它的虛數(shù)部分的子陣 Bˊ或 B〞則是非奇異矩陣，可以求逆，其逆陣為= =139。同。內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計(jì)說(shuō)明書(shū)12 PQ 分解法潮流計(jì)算流程圖 啟動(dòng)輸入原始數(shù)據(jù)形成節(jié)點(diǎn)導(dǎo)納矩陣形成 BˊB ＂并求其逆陣計(jì)算 計(jì)算電壓新值iP?iQiU新值是否滿(mǎn)足要求滿(mǎn) 足 停止置迭代次數(shù) k=0不滿(mǎn)足k+1內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計(jì)說(shuō)明書(shū)13 兩機(jī)五節(jié)點(diǎn)網(wǎng)絡(luò)潮流計(jì)算GG0 . 0 8 + j 0 . 2 4 0 . 0 1 + j 0 . 0 30.06+j0.180.06+j0.180 . 0 4 + j 0 . 1 20.08+j0.240.02+j0.060 . 4 5 + j 0 . 1 5 0 . 4 + j 0 . 0 50 . 6 + j 0 . 1 ( 0 . 2 + j 0 . 2 )123 45系統(tǒng)接線圖解： 1,.形成系數(shù)矩陣 Bˊ、B〞并求他們的逆陣由于節(jié)點(diǎn) 1 為平衡節(jié)點(diǎn)，其他的節(jié)點(diǎn)均為 PQ 節(jié)點(diǎn)，系數(shù)矩陣 ， 階數(shù)相39。??1i0i??0i（9） 運(yùn)用各節(jié)點(diǎn)電壓的新值自第三步開(kāi)始進(jìn)入下一次迭代。??0i???0Q（7） 解修正方程式( )=﹣B 〞 U（i=1,2 ，…，m ，i≠s ） 。???0i?（5） 求各節(jié)點(diǎn)電壓相位角的新值 = + （i=1,2，…，n，i ≠s） 。??0iU（3） 按式 = 計(jì)算有功功率的不平??0iP?ip????njji10jU????00sincojjijijBG??衡量 ，從而求出 （i=1，2，…，n≠s） 。內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計(jì)說(shuō)明書(shū)11第二章 理論計(jì)算 PQ 法潮流計(jì)算的基本步驟 運(yùn)用 PQ 分解法計(jì)算潮流分布時(shí)的基本步驟是：（1） 形成系數(shù)矩陣 Bˊ、B〞，并求其逆陣。 MATLAB 設(shè)計(jì)中，原始數(shù)據(jù)的填寫(xiě)格式是很關(guān)鍵的一個(gè)環(huán)節(jié)，它與程序使用的方便性和靈活性有著直接的關(guān)系。另外，MATLAB 提供了一種特殊的工具：工具箱（ TOOLBOXES）.這些工具箱主要包括：信號(hào)處理（SIGNAL PROCESSING） 、控制系統(tǒng)（CONTROL SYSTEMS） 、神經(jīng)網(wǎng)絡(luò)（NEURAL NETWORKS） 、模糊邏輯(FUZZY LOGIC)、小波(WAVELETS)和模擬（ SIMULATION）等等。MATLAB 與 C 語(yǔ)言和 FORTRAN 語(yǔ)言相比更容易掌握。 MATLAB 是一種交互式、面向?qū)ο蟮某绦蛟O(shè)計(jì)語(yǔ)言，廣泛的應(yīng)用于工業(yè)界和學(xué)術(shù)界，主要用于矩陣運(yùn)算，同時(shí)在數(shù)值分析、自動(dòng)控制模擬、MATLAB程序設(shè)計(jì)語(yǔ)言結(jié)構(gòu)完整，且具有優(yōu)良的移植性，它的基本數(shù)據(jù)元素是不需要定義的數(shù)組?，F(xiàn)有很多潮流計(jì)算方法。利 用 電 子 數(shù) 字 計(jì) 算 機(jī) 進(jìn) 行 電 力 系 統(tǒng) 潮 流 計(jì) 算 從 20 世 紀(jì) 50 年 代 中 期 就己 開(kāi) 始 ， 此 后 ， 潮 流 計(jì) 算 曾 采 用 了 各 種 不 同 的 方 法 ， 這 些 方 法 的 發(fā) 展 主 要 是圍 繞 著 對(duì) 潮 流 計(jì) 算 的 一 些 基 本 要 求 進(jìn) 行 的 ， 對(duì) 潮 流 計(jì) 算 的 要 求 可 以 歸 納 為 下面 幾 點(diǎn) ：　 　 (1)計(jì) 算 方 法 的 可 靠 性 或 收 斂 性　 　 (2)對(duì) 計(jì) 算 速 度 和 內(nèi) 存 量 的 要 求　 　 (3)計(jì) 算 的 方 便 性 和 靈 活 性一 般 情 況 下 ， 采 用 PQ 分 解 法 計(jì) 算 時(shí) 要 求 的 迭 代 次 數(shù) 較 采 用 牛 頓 拉 夫遜 法 時(shí) 多 ， 但 每 次 迭 代 所 需 的 時(shí) 間 則 較 采 用 牛 頓 拉 夫 遜 法 時(shí) 少 ， 以 致 總 的計(jì) 算 速 度 仍 是 PQ 分 解 法 快 。 原 則 上 講 ， 靜 態(tài) 安全 分 析 也 可 U 用 潮 流 計(jì) 算 來(lái) 代 替 。 電 力 系 統(tǒng) 穩(wěn) 態(tài) 分 析 潮 流 計(jì) 算 總 結(jié)潮 流 計(jì) 算 針 對(duì) 電 力 革 統(tǒng) 各 正 常 運(yùn) 行 方 式 ， 而 靜 態(tài) 安 全 分 析 則 要 研 究 各 種運(yùn) 行 方 式 下 個(gè) 別 系 統(tǒng) 元 件 退 出 運(yùn) 行 后 系 統(tǒng) 的 狀 況 。 因 此 ，潮 流 計(jì) 算 是 電 力 系 統(tǒng) 中 應(yīng) 用 最 廣 泛 、 最 基 本 和 最 重 要 的 一 種 電 氣 運(yùn) 算 。內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計(jì)說(shuō)明書(shū)9總 結(jié) 為 在 電 力 系 統(tǒng) 運(yùn) 行 方 式 和 規(guī) 劃 方 案 的 研 究 中 ， 都 需 要 進(jìn) 行 潮 流 計(jì)算 以 比 較 運(yùn) 行 方 式 或 規(guī) 劃 供 電 方 案 的 可 行 性 、 可 靠 性 和 經(jīng) 濟(jì) 性 。(3)正 常 檢 修 及 特 殊 運(yùn) 行 方 式 下 的 潮 流 計(jì) 算 ,用 于 日 運(yùn) 行 方 式 的 編 制 ,指導(dǎo) 發(fā) 電 廠 開(kāi) 機(jī) 方 式 ,有 功 、 無(wú) 功 調(diào) 整 方 案 及 負(fù) 荷 調(diào) 整 方 案 ,滿(mǎn) 足 線 路 、 變 壓器 熱 穩(wěn) 定 要 求 及 電 壓 質(zhì) 量 要 求 。具 體 表 現(xiàn) 在 以 下 方