【正文】 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。 %打開(kāi)數(shù)據(jù)文件if pathname == 0 error(39。global line。 %打開(kāi)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ī)劃、基建部門(mén)提出改進(jìn)網(wǎng)架結(jié)構(gòu),加快基建進(jìn)度的建議。通過(guò)對(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)容，寫(xiě)出節(jié)點(diǎn)導(dǎo)納矩陣； （直角坐標(biāo)系形式）初值后，求解不平衡量；； ，重新修改初值，從2開(kāi)始重新循環(huán)計(jì)算； ，再計(jì)算各支路功率分布、功率損耗和平衡節(jié)點(diǎn)功率。除了文中特別加以標(biāo)注引用的內(nèi)容外，本論文不包含任何其他個(gè)人或集體已經(jīng)發(fā)表或撰寫(xiě)的成果作品。對(duì)本研究提供過(guò)幫助和做出過(guò)貢獻(xiàn)的個(gè)人或集體，均已在文中作了明確的說(shuō)明并表示了謝意。涉密論文按學(xué)校規(guī)定處理。計(jì)算精度要求各節(jié)點(diǎn)電壓修正量不大于。同時(shí)，為了提高供電的可靠性以及資源利用的綜合經(jīng)濟(jì)性，又把許多分散的各種形式的發(fā)電站，通過(guò)送電線路和變電所聯(lián)系起來(lái)。因此，潮流計(jì)算是電力系統(tǒng)中應(yīng)用最廣泛、最基本和最重要的一種電氣運(yùn)算。clear allformat 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 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。endfclose(myf)。 2 3 0。 此次課程設(shè)計(jì)的完成，還要感謝老師的幫助以及同學(xué)的協(xié)力合作，團(tuán)結(jié)才能出成績(jī)。對(duì)本文的研究做出重要貢獻(xiàn)的個(gè)人和集體均已在文中以明確方式標(biāo)明。首先非常感謝學(xué)校開(kāi)設(shè)這個(gè)課題，為本人日后從事計(jì)算機(jī)方面的工作提供了經(jīng)驗(yàn)，奠定了基礎(chǔ)。從他身上，我學(xué)到了許多能受益終生的東西。感謝老師四年來(lái)對(duì)我孜孜不倦的教誨，對(duì)我成長(zhǎng)的關(guān)心和愛(ài)護(hù)。在論文的撰寫(xiě)過(guò)程中老師們給予我很大的幫助，幫助解決了不少的難點(diǎn)，使得論文能夠及時(shí)完成，這里一并表示真誠(chéng)的感謝。四年的風(fēng)風(fēng)雨雨，我們一同走過(guò)，充滿著關(guān)愛(ài)，給我留下了值得珍藏的最美好的記憶。其次，我要感謝大學(xué)四年中所有的任課老師和輔導(dǎo)員在學(xué)習(xí)期間對(duì)我的嚴(yán)格要求，感謝他們對(duì)我學(xué)習(xí)上和生活上的幫助，使我了解了許多專(zhuān)業(yè)知識(shí)和為人的道理，能夠在今后的生活道路上有繼續(xù)奮斗的力量。本次畢業(yè)設(shè)計(jì)是對(duì)我大學(xué)四年學(xué)習(xí)下來(lái)最好的檢驗(yàn)。 盡我所知，除文中已經(jīng)特別注明引用的內(nèi)容和致謝的地方外，本論文不包含任何其他個(gè)人或集體已經(jīng)發(fā)表或撰寫(xiě)過(guò)的研究成果。 4 2 ]。 2 1。 %輸出線路計(jì)算結(jié)果到文件fprintf(myf, 39。節(jié)點(diǎn) 節(jié)點(diǎn)電壓 節(jié)點(diǎn)相角（角度） 節(jié)點(diǎn)注入功率\n39。global myf。 end end tem(I,:)=bus(nodenum(J,1),:)。 %保存到lPQ矩陣end%******************恢復(fù)節(jié)點(diǎn)原編號(hào)*************function ret global nl。 Iji=(UjUi)*Yt/K+Uj*Yt*(1K)/K^2。 %非接地節(jié)點(diǎn)，則表示Uj end if (K==0)amp。 lPQ(k,1)=I。 sum=sum+bus(I,2)*bus(J,2)*B。 for I=nPQ+1:nb %求PV結(jié)點(diǎn)的Q，平衡節(jié)點(diǎn)的P,Q if bus(I,6)==3 %對(duì)于平衡節(jié)點(diǎn)，求其P sum=0。出錯(cuò)指示 error(39。ang%d %10f \n39。 %由于求得的是delU/U，故還需作相應(yīng)處理 end x=xdx。,i)。 for I=1:nb+nPQ1 fprintf(myf, 39。%10f 39。 %得到功率偏差列向量 dx=Jac\del。\n39。 %收斂精度eps2=。 endend %***************牛頓拉夫遜法解非線性方程序組***************function NR global nb。 end if K0 % 變壓器線路:Zt和Ym為折算到i側(cè)的值,K在j側(cè) Y(I,I)=Y(I,I)+Yt+Ym。(J~=0) % 普通線路: K=0。 Y=zeros(nb,nb)。 %形成節(jié)點(diǎn)號(hào)對(duì)應(yīng)表bus(:,1)=newbus。 PV(nPV,:)=bus(I,:)。global nodenum。 %輸出文件名形成myf=fopen(outfile,39。39。)。 %計(jì)算結(jié)果輸出%*************************主程序結(jié)束*****************%******************打開(kāi)bus line的數(shù)據(jù)文件*****************function openfile global nb nl。 %輸出文件句柄%*******************子程序調(diào)用***************openfile。 %平衡節(jié)點(diǎn)，PQ節(jié)點(diǎn)，PV節(jié)點(diǎn)個(gè)數(shù)global nb。 —拉夫遜法潮流計(jì)算基本步驟牛頓—拉夫遜法潮流計(jì)算的基本步驟不外乎有如下幾步。 （1)在電網(wǎng)規(guī)劃階段,通過(guò)潮流計(jì)算,合理規(guī)劃電源容量及接入點(diǎn),合理規(guī)劃網(wǎng)架,選擇無(wú)功補(bǔ)償方案,滿足規(guī)劃水平的大、小方式下潮流交換控制、調(diào)峰、調(diào)相、調(diào)壓的要求。然后確定節(jié)點(diǎn)類(lèi)型，由于一般無(wú)發(fā)電設(shè)備的變電所、功率固定的發(fā)電廠為PQ節(jié)點(diǎn)；有可調(diào)無(wú)功設(shè)備的變電所、有勵(lì)磁儲(chǔ)備的發(fā)電廠為PV節(jié)點(diǎn)；主調(diào)頻電廠或出線多的發(fā)電廠為平衡節(jié)點(diǎn)。 設(shè)計(jì)要求1. 程序源代碼； ，輸出文件； 3. 程序說(shuō)明； 4. 選定算例的程序計(jì)算過(guò)程； 5. 選定算例的手算過(guò)程（至少迭代2次）（可選）。作者簽名：　　 　 　　 日　 期：　　 　 　　 學(xué)位論文原創(chuàng)性聲明本人鄭重聲明：所呈交的論文是本人在導(dǎo)師的指導(dǎo)下獨(dú)立進(jìn)行研究所取得的研究成果。作 者 簽 名：　　 　 　　日　 期：　　 　 　　 指導(dǎo)教師簽名：　　 　 　　 日　　期：　　 　 　　 使用授權(quán)說(shuō)明本人完全了解 大學(xué)關(guān)于收集、保存、使用畢業(yè)設(shè)計(jì)（論文）的規(guī)定，即：按照學(xué)校要求提交畢業(yè)設(shè)計(jì)（論文）的印刷本和電子版本；學(xué)校有權(quán)保存畢業(yè)設(shè)計(jì)（論文）的印刷本和電子版，并提供目錄檢索與閱覽服務(wù)；學(xué)校可以采用影印、縮印、數(shù)字化或其它復(fù)制手段保存論文；在不以贏利為目