【正文】 lPQ(k,1)=I。 %非接地節(jié)點，則表示 Uj end if (K==0)amp。 Iji=(UjUi)*Yt/K+Uj*Yt*(1K)/K^2。 %保存到 lPQ矩陣 end %******************恢復(fù)節(jié)點原編號 ************* 貴州大學電氣工程學院《電力系統(tǒng)分析》課程設(shè)計 16 function ret global nl。 end end tem(I,:)=bus(nodenum(J,1),:)。 global myf。節(jié)點 節(jié)點電壓 節(jié)點相角（角度） 節(jié)點注入功率 \n39。 %輸出線路計算結(jié)果到文件 fprintf(myf, 39。 2 1。 4 2 ]。盡我所知，除文中已經(jīng)特別注明引用的內(nèi)容和致謝的地方外，本論文不包含任何其他個人或集體已經(jīng)發(fā)表或撰寫過的研究成果。 作者簽名 : 二〇 一 〇年 九 月 二十 日 畢業(yè)設(shè)計（論文）使用授權(quán)聲明 本人完全了解 濱州學院 關(guān)于收集、保存、使用畢業(yè)設(shè)計（論文）的規(guī)定。沒有他們的幫助，我將無法順利完成這次設(shè)計。 四年的大學生活就快走入尾聲，我們的校園生活就要劃上句號，心中是無盡的難舍與眷戀。是他們在我畢業(yè)的最后關(guān)頭給了我們巨大的幫助與鼓勵，給了我很多解決問題的思路，在此表示衷心的感激。他無論在理論上還是在 實踐中，都給與我很大的幫助，使我得到不少的提高這對于我以后的工作和學習都有一種巨大的幫助，感謝他耐心的輔導(dǎo)。 回首四年，取得了些許成績，生活中有快樂也有艱辛。郭謙功老師淵博的知識、嚴謹?shù)淖黠L和誨人不倦的態(tài)度給我留下了深刻的印象。 （保密論文在解密后遵守此規(guī)定） 作者簽名 : 二〇 一 〇年 九 月 二十 日 貴州大學電氣工程學院《電力系統(tǒng)分析》課程設(shè)計 29 致 謝 時間 飛逝，大學的學習生活很快就要過去，在這四年的學習生活中，收獲了很多，而這些成績的取得是和一直關(guān)心幫助我的人分不開的。本人完全意識到本聲明的法律結(jié)果由本人承擔。若線路功率轉(zhuǎn)化為標幺值與 matlab運算的結(jié)果進行比較。 4 2。)。 for I=1:nb, fprintf(myf, 39。\n牛頓－拉夫遜法潮流計算結(jié)果 \n39。 end bus=tem。 global bus。 Iij=(UiUj)*Yt*K+Ui*(Ym+Yt*(1K))。 Iij=Ui*(Yt+Ym)Uj*Yt。 if J~=0 Zt=line(k,3)+j*line(k,4)。 %求取的 PV、平衡節(jié)點的 Q存回 bus矩陣 end %***************計算線路功率及損耗 ****************** function flow global nl。 A=real(Y(I,J))*cos(ang)+imag(Y(I,J))*sin(ang)。)。 end for I=1:nPQ %輸出修正的 bus矩陣中 PQ結(jié)點的電壓 bus(I,2)=u(I,1)。第 %d次迭代的節(jié)點相角 delta和電壓 U\n39。%10e 39。, del(I,1))。 %輸出每次迭代的雅可比矩陣 end fprintf(myf, 39。第 %d次迭代結(jié)果 \n39。 end for i=1:max1 angl(:,1)=bus(1:nb1,3)。節(jié)點導(dǎo)納矩陣 Y\n39。 global bus。 Y(I,J)=Y(I,J)Yt/K。 Y(J,J)=Y(J,J)+Yt+Ym。 %讀入線路參數(shù) J=line(k,2)。 break end end end end %******************建立節(jié)點導(dǎo)納矩陣 ***************** function Y = ybus(bus,line) global nl。 PQ(nPQ,:)=bus(I,:)。 nPV = 0。)。.39。 you must select a valid data file39。 global bus。 %節(jié)點重新編號 ybus。 %支路數(shù) global bus。 設(shè)各節(jié)點電壓的初始值 和 求出功率偏差 根據(jù)公式算出雅可比矩陣。 貴州大學電氣工程學院《電力系統(tǒng)分析》課程設(shè)計 7 (3)正常檢修及 特殊運行方式 下的潮流計算 ,用于日運行方式的編制 ,指導(dǎo)發(fā)電廠開機方式 ,有功、無功調(diào)整方案及負荷調(diào)整方案 ,滿足線路、變壓器熱穩(wěn)定要求及電壓質(zhì)量要求。綜合比較牛頓拉夫遜法（直角坐標、極坐標）、 PQ 分解法等多種求解方法的特點，最后確定采用牛頓拉夫遜法（極坐標）。 ；連調(diào)； 。對本文的研究做出重要貢獻的個人和集體，均已在文中以明確方式標明。盡我所知，除文中特別加以標注和致謝的地方外，不包含其他人或組織已經(jīng)發(fā)表或公布過的研究成果，也不包含我為獲得 及其它教育機構(gòu)的學位或?qū)W歷而使用過的材料。本人授權(quán) 大學可以將本學位論文的全部或部分內(nèi)容編入有關(guān)數(shù)據(jù)庫進行檢索，可以采用影印、縮印或掃描等復(fù)制手段保存和匯編本學位論文。運用以極坐標表示的高斯 — 塞德爾法、牛頓 — 拉夫遜法、 PQ 分解法計算該系統(tǒng)的潮流分布。隨著工農(nóng)業(yè)生產(chǎn)和城市的發(fā)展，電能的需要量迅速增加，而熱能資源和水能資源豐富的地區(qū)又往往遠離用電比較集中的城市和工礦區(qū)，為了解決這個矛盾，就需要在動力資源豐富的地區(qū)建立大型發(fā)電站，然后將電能遠距離輸送給電力用戶。同時，為了 實時監(jiān)控 電力系統(tǒng)的運行狀態(tài)，也需要進行大量而快速的潮流計算。 (0)?(0)U (0)P? (0)Q?( 0 ) ( 0 ) ( 0 )1( 0 ) ( 0 ) ( 0 )1( c o s sin )( sin c o s )ni is i j ij ij ij ijjni is i j ij ij ij ijjP P U U G BQ Q U U G B??????? ? ? ?? ? ? ???(0)?? (0)U?( 0 )( 0 ) ( 0 )1 ( 0 )( 0 ) DP J UUQ ?????? ?? ? ???? ???? ??( 1 ) ( ) ( )( 1 ) ( ) ( )k k kk k kU U U? ? ???? ? ? ? ? ????? ? ? ? ? ??? ? ? ? ? ?貴州大學電氣工程學院《電力系統(tǒng)分析》課程設(shè)計 8 第四章程序流程圖及程序代碼 程序流程圖 根據(jù)牛頓 — 拉夫遜法潮流計算的基本步驟可以得到如下程序流程圖 啟 動輸 入 原 始 數(shù) 據(jù)形 成 導(dǎo) 納 矩 陣 Y給 定 電 壓 初 值 e 0 ， f 0置 K ＝ 0利 用 式 ( 4 4 6 ) ( 4 4 7 ) 計 算d e l t P ( k ) , d e l t Q ( k ) , d e l t U ( k )I d e l t P ( k ) , d e l t Q ( k ) I E ? 利 用 式 ( 4 4 9 ) ( 4 5 0 ) 求 雅 可 比 矩 陣 各 元 素解 修 正 式 ( 4 4 8 ) 求 d e l t e ( k ) , d e l t f ( k )修 正 各 節(jié) 點 電 壓e ( k + 1 ) = e ( k ) , f ( k + 1 ) = f ( k )K = K + 1計 算 線 路 功 率 及平 衡 節(jié) 點 功 率輸 出是 貴州大學電氣工程學院《電力系統(tǒng)分析》課程設(shè)計 9 運用 MATLAB 編程程序代碼 %*******************主程序開始 ******************** clc。 %Y:節(jié)點導(dǎo)納矩陣 global nodenum。 %計算節(jié)點注入功率 flow。*.m39。 %執(zhí)行數(shù)據(jù)文件中語句 end [nl,ml]=size(line)。 end outfile=strcat(outfile,39。 global nPQ。 if type == 3, %求平衡節(jié)點數(shù) nSW = nSW + 1。SW]。 global line。 %非接地支路為阻抗的倒數(shù) end Ym=line(k,5)+j*line(k,6)。 end if (K==0)amp。 Y(J,J)=Y(J,J)+K*K*Yt。 global myf。%10f+j*(%10f) 39。u]。第 %d次迭代的雅比矩陣 J\n39。 end fprintf(myf,39。)。\n39。 %將合一起的 x，拆成相角、電壓兩個列向量 u=x(nb:nb+nPQ1,1)。, I,u(I,1))。 global bus。 %求取的平衡節(jié)點的 P存回 bus矩陣 end sum=0。 貴州大學電氣工程學院《電力系統(tǒng)分析》課程設(shè)計 15 global line。 K=line(k,7)。(J==0) %接地支路 Iij=Ui*Ym。 %求取節(jié)點注入功率 Sji=Uj*conj(Iji)。 global nodenum。 global nl。 節(jié)點計算結(jié)果： \n39。 end fprintf(myf, 39。, lPQ(I,1),lPQ(I,2),real(lPQ(I,3)),imag(lPQ(I,3)),real(lPQ(I,4)),imag(lPQ(I,4)),real(lPQ(I,5)),imag(lPQ(I,5)))。 1 3 0。 第六章 總結(jié) 此次課程設(shè)計我在 MATLAB 編程，潮流計算， WORD 文檔的編輯方面均有所提高，但也暴露出了一 些問題；理論知識的儲備還是不足，對 MATLAB 的性能和特點還不能有一個全面的把握，對 WORD 軟件也不是很熟練，相信通過以后的學習能彌補這些不足，從而達到一個新的層次。盡我所知，除文中已經(jīng)注明引用的內(nèi)容外，本設(shè)計（論文）不含任何其他個人或集體已經(jīng)發(fā)表或撰寫過的作品成果。本次畢業(yè)設(shè)計是對我大學四年學習下來最好的檢驗。 其次，我要感謝大學四年中所有的任課老師和輔導(dǎo)員在學習期間對我的嚴格要求，感謝他們對我學習上和生活上的幫助，使我了解了許多專業(yè)知識和為人的道理，能夠在今后的生活道路上有繼續(xù)奮斗的力量。四年的風風雨雨，我們一同走過，充滿著關(guān)愛，給我留下了值得珍藏的最美好的記