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