【正文】 em(I,1)=I。 global nb。 end if K0 %變壓器線路 : Zt和 Ym為折算到 K側(cè)的值 ,K在 i側(cè) K=K。(J~=0) %普通線路非接地支路 。 %用 lPQ矩陣第一、二列保存線路兩端節(jié)點(diǎn)編號(hào) lPQ(k,2)=J。 end bus(I,5)=sum。 for J=1:nb ang=bus(I,3)bus(J,3)。超過最大迭代次數(shù)，不收斂停機(jī)！ 39。, I,angl(I,1))。 %修正 fprintf(myf,39。 for I=1:nb+nPQ1 fprintf(myf, 39。%10e 39。, Jac(I,J))。 %求得 dx fprintf(myf,39。)。 fprintf(myf, 39。 global nPQ。 Y(J,J)=Y(J,J)+Yt/K/K。 Y(I,I)=Y(I,I)+Yt+Ym。 %對(duì)導(dǎo)納矩陣賦初值 0 for k=1:nl I=line(k,1)。 for I=1:nl %根據(jù)節(jié)點(diǎn)號(hào)對(duì)應(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)。 else %求 PQ節(jié)點(diǎn)數(shù) nPQ = nPQ + 1。 nSW = 0。w39。 for I=1:length(dfile) if dfile(I)==39。 %打開數(shù)據(jù)文件 if pathname == 0 error(39。 global line。 %打開 bus line數(shù)據(jù)文件 change。 %節(jié)點(diǎn)個(gè)數(shù) global nl。 形成節(jié)點(diǎn)導(dǎo)納矩陣。 (2)在編制年運(yùn)行方式時(shí) ,在預(yù)計(jì)負(fù)荷增長(zhǎng)及新設(shè)備投運(yùn)基礎(chǔ)上 ,選擇典型方式進(jìn)行潮流計(jì)算 ,發(fā)現(xiàn)電網(wǎng)中薄弱環(huán)節(jié) ,供調(diào)度員日常調(diào)度控制參考 ,并對(duì)規(guī)劃、基建部門提出改進(jìn)網(wǎng)架結(jié)構(gòu) ,加快基建進(jìn)度的建議。通過對(duì)此電力系統(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)。 設(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)功率 。除了文中特別加以標(biāo)注引用的內(nèi)容外，本論文不包含任何其他個(gè)人或集體已經(jīng)發(fā)表或撰寫的成果作品。對(duì) 本研究提供過幫助和做出過貢獻(xiàn)的個(gè)人或集體，均已在文中作了明確的說明并表示了謝意。 涉密論文按學(xué)校規(guī)定處理。計(jì)算精度要求各節(jié)點(diǎn)電壓修正量不大于 510? 。同時(shí)，為了提 高供電的可靠性以及資源利用的綜合經(jīng)濟(jì)性，又把許多分散的各種形式的發(fā)電站，通過送電線路和變電所聯(lián)系起來。因此，潮流計(jì)算是電力系統(tǒng)中應(yīng)用最廣泛、最基本和最重要的一種電氣運(yùn)算。 clear all format long。 %節(jié)點(diǎn)編號(hào)對(duì)應(yīng)表 global lPQ。 %計(jì)算線路功率及損耗 ret。,39。 %求線路數(shù) nl [nb,mb]=size(bus)。.txt39。 global bus。 SW(nSW,:)=bus(I,:)。 %按 PQ PV SW順序排列的節(jié)點(diǎn)數(shù)據(jù)矩陣 newbus=[1:nb]39。 global Y。 %求 Ym K=line(k,7)。(J==0) % 對(duì)地支路 : K=0,J=0,R=X=0。 Y(I,J)=Y(I,J)+K*Yt。 max1=100。, real(Y(I,J)),imag(Y(I,J)))。 %從 bus矩陣中得到 PQ PV節(jié)點(diǎn)的相角、 PQ節(jié)點(diǎn)的電壓形成的列向量作為 x Jac=form_jac(bus,Y)。,i)。第 %d次迭代的功率偏差 dP和 dQ\n39。 end fprintf(myf,39。)。 for I=1:nb1 %輸出修正的 bus矩陣中 PQ、 PV結(jié)點(diǎn)的相角 bus(I,3)=angl(I,1)。 end 貴州大學(xué)電氣工程學(xué)院《電力系統(tǒng)分析》課程設(shè)計(jì) 14 if (max(abs(dx))eps1)amp。 global line。 for J=1:nb ang=bus(I,3)bus(J,3)。 for k=1:nl I=line(k,1)。 Ui=bus(I,2)*(cos(bus(I,3))+j*sin(bus(I,3)))。 Iji=0。 delS=Sij+Sji。 global lPQ。 global bus。)。\n 線路計(jì)算結(jié)果： \n39。 end fclose(myf)。 2 3 0。 此次課程設(shè)計(jì)的完成，還要感謝老師的幫助以及同學(xué)的協(xié)力合作，團(tuán)結(jié)才能貴州大學(xué)電氣工程學(xué)院《電力系統(tǒng)分析》課程設(shè)計(jì) 26 出成績(jī)。對(duì)本文的研究做出重要貢獻(xiàn)的個(gè)人和集體均已在文中以明確方式標(biāo)明。經(jīng)過這次畢業(yè)設(shè)計(jì)，我的能力有了很大的提高，比如操作能力、分析問題的能力、合作精神、嚴(yán)謹(jǐn)?shù)墓ぷ髯黠L(fēng)等方方面面都有很大的進(jìn)步。 另外，我還要感謝大學(xué) 四年和我一起走過的同學(xué)朋友對(duì)我的關(guān)心與支持，與他們一起學(xué)習(xí)、生活，讓我在大學(xué)期間生活的很充實(shí)，給我留下了很多難忘的回憶。 在我的十幾年求學(xué)歷程里，離不開父母的鼓勵(lì)和支持，是他們辛勤的勞作，無私的付出，為我創(chuàng)造良好的學(xué)習(xí)條件，我才能順利完成完成學(xué)業(yè)，感激他們一直以來對(duì)我的撫養(yǎng)與培育。 。 學(xué)友情深，情同兄妹。再次對(duì)周巍老師表示衷心的感謝。本次畢業(yè)設(shè)計(jì)大概持續(xù)了半年，現(xiàn)在終于到結(jié)尾了。 論文密級(jí)： □ 公開 □ 保密 （ ___年 __月至 __年 __月） (保密的學(xué)位論文在解密后應(yīng)遵守此協(xié)議 ) 作者簽名： _______ 導(dǎo)師簽名： _______ _______年 _____月 _____日 _______年 _____月 _____日 貴州大學(xué)電氣工程學(xué)院《電力系統(tǒng)分析》課程設(shè)計(jì) 28 獨(dú) 創(chuàng) 聲 明 本人鄭重聲明：所呈交的畢業(yè)設(shè)計(jì) (論文 )，是本人在指導(dǎo)老師的指導(dǎo)下，獨(dú)立進(jìn)行研究工作所取得的成果，成果不存在知識(shí)產(chǎn)權(quán)爭(zhēng)議。但從中我們學(xué)會(huì)了對(duì)這些實(shí)用性很強(qiáng)的軟件有深刻的了解和思考方法。 % b1 b2 ( R ) ( X ) (G) ( B ) ( K ) line = [ 1 2 0。 %2d %2d %10f+j%10f %10f+j%10f %10f+j%10f\n39。, 貴州大學(xué)電氣工程學(xué)院《電力系統(tǒng)分析》課程設(shè)計(jì) 17 bus(I,1),bus(I,2),bus(I,3)*180/pi,bus(I,4),bus(I,5))。 fprintf(myf, 39。 break end end end end ***********************將計(jì)算結(jié)果輸出 ************ function output global nb。 global Y。 end Sij=Ui*conj(Iij)。 end if (K==0)amp。 end Ym=line(k,5)+j*line(k,6)。 global bus。 end bus(I,4)=sum。 global nPQ。U%d %10f \n39。 angl=x(1:nb1,1)。 %輸出每次迭代的節(jié)點(diǎn)相角和電壓的偏差 fprintf(myf, 39。\n39。)。 fprintf(myf,39。 x=[angl。 %輸出節(jié)點(diǎn)導(dǎo)納矩陣 for I=1:nb for J=1:nb fprintf(myf, 39。 global Y。 end if K0 % 變壓器線路 :Zt和 Ym為折算到 K所在側(cè)的值 ,K在 i側(cè) Y(I,I)=Y(I,I)+Yt+Ym。 Y(J,I)=Y(I,J)。 if J~=0 %接地支路分母為零 Yt=1/Zt。 global bus。PV。 for I = 1:nb, % nb為總節(jié)點(diǎn)數(shù) type= bus(I,6)。 global nl。 end outfile=strcat(outfile,dfile(I))。 % strip off .m 貴州大學(xué)電氣工程學(xué)院《電力系統(tǒng)分析》課程設(shè)計(jì) 10 eval(dfile(1:lfile2))。 [dfile,pathname]=uigetfile(39。 %牛頓拉夫遜法解非線性方程序組 PQ。 %line:支路數(shù)據(jù) global Y。 計(jì)算線路功率分布， PV 節(jié)點(diǎn)無功功率和平衡節(jié)點(diǎn)注入功率。 總結(jié)為在 電力系統(tǒng)運(yùn)行方式 和規(guī)劃方案的研究中，都需要進(jìn)行潮流計(jì)算以比較運(yùn)行方式或規(guī)劃供電方案的可行性、可靠性和經(jīng)濟(jì)性。 第三章 電力系統(tǒng)潮流計(jì)算概述 電力系統(tǒng)概述 電力工業(yè)發(fā)展初期，電能是直接在用戶附近的發(fā)電站（或稱發(fā)電廠）中生產(chǎn)的，各發(fā)電站孤立運(yùn)行。 設(shè)計(jì)時(shí)間 2020年春季第 17周至第 18周 第二章 潮流計(jì)算設(shè)計(jì)題目及思路 潮流計(jì)算題目 題目：設(shè)計(jì)原始資料網(wǎng)絡(luò)接線如下圖所示： 其 中 已 知 ： 各 支 路 阻 抗 表 幺 值 參 數(shù) 分 別 為Z12=+,Z13=+,Z23=+,ZT1=,ZT2=,K=,L L2 電納為 。 作者簽名： 日期： 年 月 日 學(xué)位論 文版權(quán)使用授權(quán)書 本學(xué)位論文作者完全了解學(xué)校有關(guān)保留、使用學(xué)位論文的規(guī)定，同意學(xué)校保留并向國家有關(guān)部門或機(jī)構(gòu)送交論文的復(fù)印件和電子版，允許論文被查閱和借閱。貴州大學(xué)電氣工程學(xué)院 電力系統(tǒng)分析課程設(shè)計(jì)報(bào)告書 題目 : 高壓輸電網(wǎng)潮流的計(jì)算機(jī)算法程序設(shè)計(jì) 畢業(yè)設(shè)計(jì)（論文）原創(chuàng)性聲明和使