【正文】 nPV = 0。 SW(nSW,:)=bus(I,:)。 PQ(nPQ,:)=bus(I,:)。 %按PQ PV SW順序排列的節(jié)點(diǎn)數(shù)據(jù)矩陣newbus=[1:nb]39。 break end end endend %******************建立節(jié)點(diǎn)導(dǎo)納矩陣*****************function Y = ybus(bus,line)global nl。global Y。 %讀入線路參數(shù) J=line(k,2)。 %求Ym K=line(k,7)。 Y(J,J)=Y(J,J)+Yt+Ym。(J==0) % 對(duì)地支路: K=0,J=0,R=X=0。 Y(I,J)=Y(I,J)Yt/K。 Y(I,J)=Y(I,J)+K*Yt。global bus。 max1=100。節(jié)點(diǎn)導(dǎo)納矩陣Y\n39。, real(Y(I,J)),imag(Y(I,J)))。end for i=1:max1 angl(:,1)=bus(1:nb1,3)。 %從bus矩陣中得到PQ PV節(jié)點(diǎn)的相角、PQ節(jié)點(diǎn)的電壓形成的列向量作為x Jac=form_jac(bus,Y)。第%d次迭代結(jié)果\n39。,i)。 %輸出每次迭代的雅可比矩陣 end fprintf(myf, 39。第%d次迭代的功率偏差dP和dQ\n39。, del(I,1))。 end fprintf(myf,39。%10e 39。)。第%d次迭代的節(jié)點(diǎn)相角delta和電壓U\n39。 for I=1:nb1 %輸出修正的bus矩陣中PQ、PV結(jié)點(diǎn)的相角 bus(I,3)=angl(I,1)。 end for I=1:nPQ %輸出修正的bus矩陣中PQ結(jié)點(diǎn)的電壓 bus(I,2)=u(I,1)。 end if (max(abs(dx))eps1)amp。)。global line。 A=real(Y(I,J))*cos(ang)+imag(Y(I,J))*sin(ang)。 for J=1:nb ang=bus(I,3)bus(J,3)。 %求取的PV、平衡節(jié)點(diǎn)的Q存回bus矩陣end%***************計(jì)算線路功率及損耗******************function flow global nl。 for k=1:nl I=line(k,1)。 if J~=0 Zt=line(k,3)+j*line(k,4)。 Ui=bus(I,2)*(cos(bus(I,3))+j*sin(bus(I,3)))。 Iij=Ui*(Yt+Ym)Uj*Yt。 Iji=0。 Iij=(UiUj)*Yt*K+Ui*(Ym+Yt*(1K))。delS=Sij+Sji。global bus。global lPQ。endbus=tem。global bus。\n牛頓－拉夫遜法潮流計(jì)算結(jié)果\n39。)。for I=1:nb, fprintf(myf, 39。\n 線路計(jì)算結(jié)果：\n39。)。endfclose(myf)。 4 2。 2 3 0。若線路功率轉(zhuǎn)化為標(biāo)幺值與matlab運(yùn)算的結(jié)果進(jìn)行比較。 此次課程設(shè)計(jì)的完成，還要感謝老師的幫助以及同學(xué)的協(xié)力合作，團(tuán)結(jié)才能出成績(jī)。本人完全意識(shí)到本聲明的法律結(jié)果由本人承擔(dān)。對(duì)本文的研究做出重要貢獻(xiàn)的個(gè)人和集體均已在文中以明確方式標(biāo)明。畢業(yè)設(shè)計(jì)（論文）使用授權(quán)聲明本人完全了解濱州學(xué)院關(guān)于收集、保存、使用畢業(yè)設(shè)計(jì)（論文）的規(guī)定。首先非常感謝學(xué)校開(kāi)設(shè)這個(gè)課題，為本人日后從事計(jì)算機(jī)方面的工作提供了經(jīng)驗(yàn)，奠定了基礎(chǔ)。這期間凝聚了很多人的心血，在此我表示由衷的感謝。從他身上，我學(xué)到了許多能受益終生的東西。最后，我要感謝我的父母對(duì)我的關(guān)系和理解，如果沒(méi)有他們?cè)谖业膶W(xué)習(xí)生涯中的無(wú)私奉獻(xiàn)和默默支持，我將無(wú)法順利完成今天的學(xué)業(yè)。感謝老師四年來(lái)對(duì)我孜孜不倦的教誨，對(duì)我成長(zhǎng)的關(guān)心和愛(ài)護(hù)。最后，我要特別感謝我的導(dǎo)師趙達(dá)睿老師、和研究生助教熊偉麗老師。在論文的撰寫(xiě)過(guò)程中老師們給予我很大的幫助，幫助解決了不少的難點(diǎn)，使得論文能夠及時(shí)完成，這里一并表示真誠(chéng)的感謝。老師們認(rèn)真負(fù)責(zé)的工作態(tài)度，嚴(yán)謹(jǐn)?shù)闹螌W(xué)精神和深厚的理論水平都使我收益匪淺。四年的風(fēng)風(fēng)雨雨，我們一同走過(guò)，充滿著關(guān)愛(ài)，給我留下了值得珍藏的最美好的記憶。從這里走出，對(duì)我的人生來(lái)說(shuō)，將是踏上一個(gè)新的征程，要把所學(xué)的知識(shí)應(yīng)用到實(shí)際工作中去。其次，我要感謝大學(xué)四年中所有的任課老師和輔導(dǎo)員在學(xué)習(xí)期間對(duì)我的嚴(yán)格要求，感謝他們對(duì)我學(xué)習(xí)上和生活上的幫助，使我了解了許多專業(yè)知識(shí)和為人的道理，能夠在今后的生活道路上有繼續(xù)奮斗的力量。首先，我要特別感謝我的知道郭謙功老師對(duì)我的悉心指導(dǎo)，在我的論文書(shū)寫(xiě)及設(shè)計(jì)過(guò)程中給了我大量的幫助和指導(dǎo)，為我理清了設(shè)計(jì)思路和操作方法，并對(duì)我所做的課題提出了有效的改進(jìn)方案。本次畢業(yè)設(shè)計(jì)是對(duì)我大學(xué)四年學(xué)習(xí)下來(lái)最好的檢驗(yàn)。（保密論文在解密后遵守此規(guī)定） 論文密級(jí)：□公開(kāi) □保密（___年__月至__年__月）(保密的學(xué)位論文在解密后應(yīng)遵守此協(xié)議)作者簽名：_______ 導(dǎo)師簽名：______________年_____月_____日 _______年_____月_____日獨(dú) 創(chuàng) 聲 明本人鄭重聲明：所呈交的畢業(yè)設(shè)計(jì)(論文)，是本人在指導(dǎo)老師的指導(dǎo)下，獨(dú)立進(jìn)行研究工作所取得的成果，成果不存在知識(shí)產(chǎn)權(quán)爭(zhēng)議。盡我所知，除文中已經(jīng)特別注明引用的內(nèi)容和致謝的地方外，本論文不包含任何其他個(gè)人或集體已經(jīng)發(fā)表或撰寫(xiě)過(guò)的研究成果。但從中我們學(xué)會(huì)了對(duì)這些實(shí)用性很強(qiáng)的軟件有深刻的了解和思考方法。 4 2 ]。 % b1 b2 ( R ) ( X ) (G) ( B ) ( K )line = [ 1 2 0。 2 1。 %2d %2d %10f+j%10f %10f+j%10f %10f+j%10f\n39。 %輸出線路計(jì)算結(jié)果到文件fprintf(myf, 39。, bus(I,1),bus(I,2),bus(I,3)*180/pi,bus(I,4),bus(I,5))。節(jié)點(diǎn) 節(jié)點(diǎn)電壓 節(jié)點(diǎn)相角（角度） 節(jié)點(diǎn)注入功率\n39。fprintf(myf, 39。global myf。 break end end endend ***********************將計(jì)算結(jié)果輸出************function output global nb。 end end tem(I,:)=bus(nodenum(J,1),:)。global Y。 %保存到lPQ矩陣end%******************恢復(fù)節(jié)點(diǎn)原編號(hào)*************function ret global nl。 end Sij=Ui*conj(Iij)。 Iji=(UjUi)*Yt/K+Uj*Yt*(1K)/K^2。 end if (K==0)amp。 %非接地節(jié)點(diǎn)，則表示Uj end if (K==0)amp。 end Ym=line(k,5)+j*line(k,6)。 lPQ(k,1)=I。global bus。 sum=sum+bus(I,2)*bus(J,2)*B。 end bus(I,4)=sum。 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。global nPQ。出錯(cuò)指示 error(39。U%d %10f \n39。ang%d %10f \n39。 angl=x(1:nb1,1)。 %由于求得的是delU/U，故還需作相應(yīng)處理 end x=xdx。 %輸出每次迭代的節(jié)點(diǎn)相角和電壓的偏差 fprintf(myf, 39。,i)。\n39。 for I=1:nb+nPQ1 fprintf(myf, 39。)。%10f 39