【正文】 utfile,39。 global nl。 global nodenum。 for I = 1:nb, % nb為總節(jié)點數(shù) type= bus(I,6)。 PV(nPV,:)=bus(I,:)。PV。 %形成節(jié)點號對應(yīng)表 貴州大學(xué)電氣工程學(xué)院《電力系統(tǒng)分析》課程設(shè)計 11 bus(:,1)=newbus。 global bus。 Y=zeros(nb,nb)。 if J~=0 %接地支路分母為零 Yt=1/Zt。(J~=0) % 普通線路 : K=0。 Y(J,I)=Y(I,J)。 end 貴州大學(xué)電氣工程學(xué)院《電力系統(tǒng)分析》課程設(shè)計 12 if K0 % 變壓器線路 :Zt和 Ym為折算到 i側(cè)的值 ,K在 j側(cè) Y(I,I)=Y(I,I)+Yt+Ym。 end if K0 % 變壓器線路 :Zt和 Ym為折算到 K所在側(cè)的值 ,K在 i側(cè) Y(I,I)=Y(I,I)+Yt+Ym。 end end %***************牛頓拉夫遜法解非線性方程序組 *************** function NR global nb。 global Y。 %收斂精度 eps2=。 %輸出節(jié)點導(dǎo)納矩陣 for I=1:nb for J=1:nb fprintf(myf, 39。\n39。 x=[angl。 %得到功率偏差列向量 貴州大學(xué)電氣工程學(xué)院《電力系統(tǒng)分析》課程設(shè)計 13 dx=Jac\del。 fprintf(myf,39。%10f 39。)。 for I=1:nb+nPQ1 fprintf(myf, 39。\n39。,i)。 %輸出每次迭代的節(jié)點相角和電壓的偏差 fprintf(myf, 39。 %由于求得的是 delU/U，故還需作相應(yīng)處理 end x=xdx。 angl=x(1:nb1,1)。ang%d %10f \n39。U%d %10f \n39。出錯指示 error(39。 global nPQ。 for I=nPQ+1:nb %求 PV結(jié)點的 Q，平衡節(jié)點的 P,Q if bus(I,6)==3 %對于平衡節(jié)點，求其 P sum=0。 end bus(I,4)=sum。 sum=sum+bus(I,2)*bus(J,2)*B。 global bus。 lPQ(k,1)=I。 end Ym=line(k,5)+j*line(k,6)。 %非接地節(jié)點，則表示 Uj end if (K==0)amp。 end if (K==0)amp。 Iji=(UjUi)*Yt/K+Uj*Yt*(1K)/K^2。 end Sij=Ui*conj(Iij)。 %保存到 lPQ矩陣 end %******************恢復(fù)節(jié)點原編號 ************* 貴州大學(xué)電氣工程學(xué)院《電力系統(tǒng)分析》課程設(shè)計 16 function ret global nl。 global Y。 end end tem(I,:)=bus(nodenum(J,1),:)。 break end end end end ***********************將計算結(jié)果輸出 ************ function output global nb。 global myf。 fprintf(myf, 39。節(jié)點 節(jié)點電壓 節(jié)點相角（角度） 節(jié)點注入功率 \n39。, 貴州大學(xué)電氣工程學(xué)院《電力系統(tǒng)分析》課程設(shè)計 17 bus(I,1),bus(I,2),bus(I,3)*180/pi,bus(I,4),bus(I,5))。 %輸出線路計算結(jié)果到文件 fprintf(myf, 39。 %2d %2d %10f+j%10f %10f+j%10f %10f+j%10f\n39。 2 1。 % b1 b2 ( R ) ( X ) (G) ( B ) ( K ) line = [ 1 2 0。 4 2 ]。但從中我們學(xué)會了對這些實用性很強的軟件有深刻的了解和思考方法。盡我所知，除文中已經(jīng)特別注明引用的內(nèi)容和致謝的地方外，本論文不包含任何其他個人或集體已經(jīng)發(fā)表或撰寫過的研究成果。 論文密級： □ 公開 □ 保密 （ ___年 __月至 __年 __月） (保密的學(xué)位論文在解密后應(yīng)遵守此協(xié)議 ) 作者簽名： _______ 導(dǎo)師簽名： _______ _______年 _____月 _____日 _______年 _____月 _____日 貴州大學(xué)電氣工程學(xué)院《電力系統(tǒng)分析》課程設(shè)計 28 獨 創(chuàng) 聲 明 本人鄭重聲明：所呈交的畢業(yè)設(shè)計 (論文 )，是本人在指導(dǎo)老師的指導(dǎo)下，獨立進(jìn)行研究工作所取得的成果，成果不存在知識產(chǎn)權(quán)爭議。 作者簽名 : 二〇 一 〇年 九 月 二十 日 畢業(yè)設(shè)計（論文）使用授權(quán)聲明 本人完全了解 濱州學(xué)院 關(guān)于收集、保存、使用畢業(yè)設(shè)計（論文）的規(guī)定。本次畢業(yè)設(shè)計大概持續(xù)了半年，現(xiàn)在終于到結(jié)尾了。沒有他們的幫助，我將無法順利完成這次設(shè)計。再次對周巍老師表示衷心的感謝。 四年的大學(xué)生活就快走入尾聲，我們的校園生活就要劃上句號，心中是無盡的難舍與眷戀。 學(xué)友情深，情同兄妹。是他們在我畢業(yè)的最后關(guān)頭給了我們巨大的幫助與鼓勵，給了我很多解決問題的思路，在此表示衷心的感激。 。他無論在理論上還是在 實踐中，都給與我很大的幫助，使我得到不少的提高這對于我以后的工作和學(xué)習(xí)都有一種巨大的幫助，感謝他耐心的輔導(dǎo)。 在我的十幾年求學(xué)歷程里，離不開父母的鼓勵和支持，是他們辛勤的勞作，無私的付出，為我創(chuàng)造良好的學(xué)習(xí)條件，我才能順利完成完成學(xué)業(yè)，感激他們一直以來對我的撫養(yǎng)與培育。 回首四年，取得了些許成績，生活中有快樂也有艱辛。 另外，我還要感謝大學(xué) 四年和我一起走過的同學(xué)朋友對我的關(guān)心與支持，與他們一起學(xué)習(xí)、生活，讓我在大學(xué)期間生活的很充實，給我留下了很多難忘的回憶。郭謙功老師淵博的知識、嚴(yán)謹(jǐn)?shù)淖黠L(fēng)和誨人不倦的態(tài)度給我留下了深刻的印象。經(jīng)過這次畢業(yè)設(shè)計，我的能力有了很大的提高，比如操作能力、分析問題的能力、合作精神、嚴(yán)謹(jǐn)?shù)墓ぷ髯黠L(fēng)等方方面面都有很大的進(jìn)步。 （保密論文在解密后遵守此規(guī)定） 作者簽名 : 二〇 一 〇年 九 月 二十 日 貴州大學(xué)電氣工程學(xué)院《電力系統(tǒng)分析》課程設(shè)計 29 致 謝 時間 飛逝，大學(xué)的學(xué)習(xí)生活很快就要過去，在這四年的學(xué)習(xí)生活中，收獲了很多，而這些成績的取得是和一直關(guān)心幫助我的人分不開的。對本文的研究做出重要貢獻(xiàn)的個人和集體均已在文中以明確方式標(biāo)明。本人完全意識到本聲明的法律結(jié)果由本人承擔(dān)。 此次課程設(shè)計的完成，還要感謝老師的幫助以及同學(xué)的協(xié)力合作，團(tuán)結(jié)才能貴州大學(xué)電氣工程學(xué)院《電力系統(tǒng)分析》課程設(shè)計 26 出成績。若線路功率轉(zhuǎn)化為標(biāo)幺值與 matlab運算的結(jié)果進(jìn)行比較。 2 3 0。 4 2。 end fclose(myf)。)。\n 線路計算結(jié)果： \n39。 for I=1:nb, fprintf(myf, 39。)。\n牛頓－拉夫遜法潮流計算結(jié)果 \n39。 global bus。 end bus=tem。 global lPQ。 global bus。 delS=Sij+Sji。 Iij=(UiUj)*Yt*K+Ui*(Ym+Yt*(1K))。 Iji=0。 Iij=Ui*(Yt+Ym)Uj*Yt。 Ui=bus(I,2)*(cos(bus(I,3))+j*sin(bus(I,3)))。 if J~=0 Zt=line(k,3)+j*line(k,4)。 for k=1:nl I=line(k,1)。 %求取的 PV、平衡節(jié)點的 Q存回 bus矩陣 end %***************計算線路功率及損耗 ****************** function flow global nl。 for J=1:nb ang=bus(I,3)bus(J,3)。 A=real(Y(I,J))*cos(ang)+imag(Y(I,J))*sin(ang)。 global line。)。 end 貴州大學(xué)電氣工程學(xué)院《電力系統(tǒng)分析》課程設(shè)計 14 if (max(abs(dx))eps1)amp。 end for I=1:nPQ %輸出修正的 bus矩陣中 PQ結(jié)點的電壓 bus(I,2)=u(I,1)。 for I=1:nb1 %輸出修正的 bus矩陣中 PQ、 PV結(jié)點的相角 bus(I,3)=angl(I,1)。第 %d次迭代的節(jié)點相角 delta和電壓 U\n39。)。%10e 39。 end fprintf(myf,39。, del(I,1))。第 %d次迭代的功率偏差 dP和 dQ\n39。 %輸出每次迭代的雅可比矩陣 end fprintf(myf, 39。,i)。第 %d次迭代結(jié)果 \n39。 %從 bus矩陣中得到 PQ PV節(jié)點的相角、 PQ節(jié)點的電壓形成的列向量作為 x Jac=form_jac(bus,Y)。 end for i=1:max1 angl(:,1)=bus(1:nb1,3)。, real(Y(I,J)),imag(Y(I,J)))。節(jié)點導(dǎo)納矩陣 Y\n39。 max1=100。 global bus。 Y(I,J)=Y(I,J)+K*Yt。 Y(I,J)=Y(I,J)Yt/K。(J==0) % 對地支路 : K=0,J=0,R=X=0。 Y(J,J)=Y(J,J)+Yt+Ym。 %求 Ym K=line(k,7)。 %讀入線路參數(shù) J=lin