【正文】 e(k,2)。 global Y。 break end end end end %******************建立節(jié)點(diǎn)導(dǎo)納矩陣 ***************** function Y = ybus(bus,line) global nl。 %按 PQ PV SW順序排列的節(jié)點(diǎn)數(shù)據(jù)矩陣 newbus=[1:nb]39。 PQ(nPQ,:)=bus(I,:)。 SW(nSW,:)=bus(I,:)。 nPV = 0。 global bus。)。.txt39。.39。 %求線路數(shù) nl [nb,mb]=size(bus)。 you must select a valid data file39。,39。 global bus。 %計(jì)算線路功率及損耗 ret。 %節(jié)點(diǎn)重新編號(hào) ybus。 %節(jié)點(diǎn)編號(hào)對應(yīng)表 global lPQ。 %支路數(shù) global bus。 clear all format long。 設(shè)各節(jié)點(diǎn)電壓的初始值 和 求出功率偏差 根據(jù)公式算出雅可比矩陣。因此，潮流計(jì)算是電力系統(tǒng)中應(yīng)用最廣泛、最基本和最重要的一種電氣運(yùn)算。 貴州大學(xué)電氣工程學(xué)院《電力系統(tǒng)分析》課程設(shè)計(jì) 7 (3)正常檢修及 特殊運(yùn)行方式 下的潮流計(jì)算 ,用于日運(yùn)行方式的編制 ,指導(dǎo)發(fā)電廠開機(jī)方式 ,有功、無功調(diào)整方案及負(fù)荷調(diào)整方案 ,滿足線路、變壓器熱穩(wěn)定要求及電壓質(zhì)量要求。同時(shí)，為了提 高供電的可靠性以及資源利用的綜合經(jīng)濟(jì)性，又把許多分散的各種形式的發(fā)電站，通過送電線路和變電所聯(lián)系起來。綜合比較牛頓拉夫遜法（直角坐標(biāo)、極坐標(biāo)）、 PQ 分解法等多種求解方法的特點(diǎn)，最后確定采用牛頓拉夫遜法（極坐標(biāo)）。計(jì)算精度要求各節(jié)點(diǎn)電壓修正量不大于 510? 。 ；連調(diào)； 。 涉密論文按學(xué)校規(guī)定處理。對本文的研究做出重要貢獻(xiàn)的個(gè)人和集體，均已在文中以明確方式標(biāo)明。對 本研究提供過幫助和做出過貢獻(xiàn)的個(gè)人或集體，均已在文中作了明確的說明并表示了謝意。盡我所知，除文中特別加以標(biāo)注和致謝的地方外，不包含其他人或組織已經(jīng)發(fā)表或公布過的研究成果，也不包含我為獲得 及其它教育機(jī)構(gòu)的學(xué)位或?qū)W歷而使用過的材料。除了文中特別加以標(biāo)注引用的內(nèi)容外，本論文不包含任何其他個(gè)人或集體已經(jīng)發(fā)表或撰寫的成果作品。本人授權(quán) 大學(xué)可以將本學(xué)位論文的全部或部分內(nèi)容編入有關(guān)數(shù)據(jù)庫進(jìn)行檢索，可以采用影印、縮印或掃描等復(fù)制手段保存和匯編本學(xué)位論文。 設(shè)計(jì)題目 高壓輸電網(wǎng)潮流的計(jì)算機(jī)算法程序設(shè)計(jì)（ PQ分解法、牛頓 拉夫遜法） 或中壓配電網(wǎng)潮流的計(jì)算機(jī)算法程序設(shè)計(jì)（前推后代法、同倫延拓法等） 或電力系統(tǒng)短路故障的計(jì)算機(jī)算法程序設(shè)計(jì)（要求不限） 設(shè)計(jì)內(nèi)容 ，寫出節(jié)點(diǎn)導(dǎo)納矩陣； （直角坐標(biāo)系形式）初值后，求解不平衡量； ； ，重新修改初值，從 2開始重新循環(huán)計(jì)算； 貴州大學(xué)電氣工程學(xué)院《電力系統(tǒng)分析》課程設(shè)計(jì) 5 ，再計(jì)算各支路功率分布、功率損 耗和平衡節(jié)點(diǎn)功率 。運(yùn)用以極坐標(biāo)表示的高斯 — 塞德爾法、牛頓 — 拉夫遜法、 PQ 分解法計(jì)算該系統(tǒng)的潮流分布。通過對此電力系統(tǒng)分析，我們很容易得以節(jié)點(diǎn) 3 為 PQ 節(jié)點(diǎn)，節(jié)點(diǎn) 4 為 PV節(jié)點(diǎn)，節(jié) 點(diǎn) 5 為平衡節(jié)點(diǎn)。隨著工農(nóng)業(yè)生產(chǎn)和城市的發(fā)展，電能的需要量迅速增加，而熱能資源和水能資源豐富的地區(qū)又往往遠(yuǎn)離用電比較集中的城市和工礦區(qū)，為了解決這個(gè)矛盾，就需要在動(dòng)力資源豐富的地區(qū)建立大型發(fā)電站，然后將電能遠(yuǎn)距離輸送給電力用戶。 (2)在編制年運(yùn)行方式時(shí) ,在預(yù)計(jì)負(fù)荷增長及新設(shè)備投運(yùn)基礎(chǔ)上 ,選擇典型方式進(jìn)行潮流計(jì)算 ,發(fā)現(xiàn)電網(wǎng)中薄弱環(huán)節(jié) ,供調(diào)度員日常調(diào)度控制參考 ,并對規(guī)劃、基建部門提出改進(jìn)網(wǎng)架結(jié)構(gòu) ,加快基建進(jìn)度的建議。同時(shí)，為了 實(shí)時(shí)監(jiān)控 電力系統(tǒng)的運(yùn)行狀態(tài)，也需要進(jìn)行大量而快速的潮流計(jì)算。 形成節(jié)點(diǎn)導(dǎo)納矩陣。 (0)?(0)U (0)P? (0)Q?( 0 ) ( 0 ) ( 0 )1( 0 ) ( 0 ) ( 0 )1( c o s sin )( sin c o s )ni is i j ij ij ij ijjni is i j ij ij ij ijjP P U U G BQ Q U U G B??????? ? ? ?? ? ? ???(0)?? (0)U?( 0 )( 0 ) ( 0 )1 ( 0 )( 0 ) DP J UUQ ?????? ?? ? ???? ???? ??( 1 ) ( ) ( )( 1 ) ( ) ( )k k kk k kU U U? ? ???? ? ? ? ? ????? ? ? ? ? ??? ? ? ? ? ?貴州大學(xué)電氣工程學(xué)院《電力系統(tǒng)分析》課程設(shè)計(jì) 8 第四章程序流程圖及程序代碼 程序流程圖 根據(jù)牛頓 — 拉夫遜法潮流計(jì)算的基本步驟可以得到如下程序流程圖 啟 動(dòng)輸 入 原 始 數(shù) 據(jù)形 成 導(dǎo) 納 矩 陣 Y給 定 電 壓 初 值 e 0 ， f 0置 K ＝ 0利 用 式 ( 4 4 6 ) ( 4 4 7 ) 計(jì) 算d e l t P ( k ) , d e l t Q ( k ) , d e l t U ( k )I d e l t P ( k ) , d e l t Q ( k ) I E ? 利 用 式 ( 4 4 9 ) ( 4 5 0 ) 求 雅 可 比 矩 陣 各 元 素解 修 正 式 ( 4 4 8 ) 求 d e l t e ( k ) , d e l t f ( k )修 正 各 節(jié) 點(diǎn) 電 壓e ( k + 1 ) = e ( k ) , f ( k + 1 ) = f ( k )K = K + 1計(jì) 算 線 路 功 率 及平 衡 節(jié) 點(diǎn) 功 率輸 出是 貴州大學(xué)電氣工程學(xué)院《電力系統(tǒng)分析》課程設(shè)計(jì) 9 運(yùn)用 MATLAB 編程程序代碼 %*******************主程序開始 ******************** clc。 %節(jié)點(diǎn)個(gè)數(shù) global nl。 %Y:節(jié)點(diǎn)導(dǎo)納矩陣 global nodenum。 %打開 bus line數(shù)據(jù)文件 change。 %計(jì)算節(jié)點(diǎn)注入功率 flow。 global line。*.m39。 %打開數(shù)據(jù)文件 if pathname == 0 error(39。 %執(zhí)行數(shù)據(jù)文件中語句 end [nl,ml]=size(line)。 for I=1:length(dfile) if dfile(I)==39。 end outfile=strcat(outfile,39。w39。 global nPQ。 nSW = 0。 if type == 3, %求平衡節(jié)點(diǎn)數(shù) nSW = nSW + 1。 else %求 PQ節(jié)點(diǎn)數(shù) nPQ = nPQ + 1。SW]。 for I=1:nl %根據(jù)節(jié)點(diǎn)號(hào)對應(yīng)表重新排列線路矩陣節(jié)點(diǎn)編號(hào) for J=1:2 for k=1:nb if line(I,J)==nodenum(k,2) line(I,J)=nodenum(k,1)。 global line。 %對導(dǎo)納矩陣賦初值 0 for k=1:nl I=line(k,1)。 %非接地支路為阻抗的倒數(shù) end Ym=line(k,5)+j*line(k,6)。 Y(I,I)=Y(I,I)+Yt+Ym。 end if (K==0)amp。 Y(J,J)=Y(J,J)+Yt/K/K。 Y(J,J)=Y(J,J)+K*K*Yt。 global nPQ。 global myf。 fprintf(myf, 39。%10f+j*(%10f) 39。)。u]。 %求得 dx fprintf(myf,39。第 %d次迭代的雅比矩陣 J\n39。, Jac(I,J))。 end fprintf(myf,39。%10e 39。)。 for I=1:nb+nPQ1 fprintf(myf, 39。\n39。 %修正 fprintf(myf,39。 %將合一起的 x，拆成相角、電壓兩個(gè)列向量 u=x(nb:nb+nPQ1,1)。, I,angl(I,1))。, I,u(I,1))。超過最大迭代次數(shù)，不收斂停機(jī)！ 39。 global bus。 for J=1:nb ang=bus(I,3)bus(J,3)。 %求取的平衡節(jié)點(diǎn)的 P存回 bus矩陣 end sum=0。 end bus(I,5)=sum。 貴州大學(xué)電氣工程學(xué)院《電力系統(tǒng)分析》課程設(shè)計(jì) 15 global line。 %用 lPQ矩陣第一、二列保存線路兩端節(jié)點(diǎn)編號(hào) lPQ(k,2)=J。 K=line(k,7)。(J~=0) %普通線路非接地支路 。(J==0) %接地支路 Iij=Ui*Ym。 end if K0 %變壓器線路 : Zt和 Ym為折算到 K側(cè)的值 ,K在 i側(cè) K=K。 %求取節(jié)點(diǎn)注入功率 Sji=Uj*conj(Iji)。 global nb。 global nodenum。%bus矩陣中的節(jié)點(diǎn)號(hào)恢復(fù)原來的編號(hào) tem(I,1)=I。 global nl。 fprintf(myf, 39。 節(jié)點(diǎn)計(jì)算結(jié)果： \n39。)。 end fprintf(myf, 39。節(jié)點(diǎn) I 節(jié)點(diǎn) J 線路功率 S(I,J) 線路功率 S(J,I) 線路損耗 dS(I,J)\n39。, lPQ(I,1),lPQ(I,2),real(lPQ(I,3)),imag(lPQ(I,3)),real(lPQ(I,4)),imag(lPQ(I,4)),real(lPQ(I,5)),imag(lPQ(I,5)))。 3 1。 1 3 0。 貴州大學(xué)電氣工程學(xué)院《電力系統(tǒng)分析》課程設(shè)計(jì) 18 第五章 運(yùn)行結(jié)果分析 運(yùn)用 MATLAB 運(yùn)行結(jié)果 節(jié)點(diǎn)導(dǎo)納矩陣 Y +j*( ) +j*( ) +j*( ) +j*( ) +j*( ) +j*( ) +j*() +j*( ) +j*( ) +j*( ) +j*( ) +j*( ) +j*() +j*( ) +j*( ) +j*( )