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