【正文】 fprintf(myf, 39。\n牛頓－拉夫遜法潮流計算結(jié)果\n39。endbus=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é)點導納矩陣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 endend %******************建立節(jié)點導納矩陣*****************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ù)公式算出雅可比矩陣。 (3)正常檢修及特殊運行方式下的潮流計算,用于日運行方式的編制,指導發(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ù)庫進行檢索，可以采用影印、縮印或掃描等復制手段保存和匯編本學位論文。運用以極坐標表示的高斯—塞德爾法、牛頓—拉夫遜法、PQ分解法計算該系統(tǒng)的潮流分布。隨著工農(nóng)業(yè)生產(chǎn)和城市的發(fā)展，電能的需要量迅速增加，而熱能資源和水能資源豐富的地區(qū)又往往遠離用電比較集中的城市和工礦區(qū)，為了解決這個矛盾，就需要在動力資源豐富的地區(qū)建立大型發(fā)電站，然后將電能遠距離輸送給電力用戶。同時，為了實時監(jiān)控電力系統(tǒng)的運行狀態(tài)，也需要進行大量而快速的潮流計算。第四章程序流程圖及程序代碼根據(jù)牛頓—拉夫遜法潮流計算的基本步驟可以得到如下程序流程圖%*******************主程序開始********************clc。 %Y:節(jié)點導納矩陣global nodenum。 %計算節(jié)點注入功率 flow。*.m39。 %執(zhí)行數(shù)據(jù)文件中語句end [nl,ml]=size(line)。 endoutfile=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。global line。 K=line(k,7)。(J==0) %接地支路 Iij=Ui*Ym。 %求取節(jié)點注入功率 Sji=Uj*conj(Iji)。global nodenum。global nl。 節(jié)點計算結(jié)果：\n39。endfprintf(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。第6章 總結(jié) 此次課程設(shè)計我在MATLAB編程，潮流計算，WORD文檔的編輯方面均有所提高，但也暴露出了一些問題；理論知識的儲備還是不足，對MATLAB的性能和特點還不能有一個全面的把握，對WORD軟件也不是很熟練，相信通過以后的學習能彌補這些不足，從而達到一個新的層次。盡我所知，除文中已經(jīng)注明引用的內(nèi)容外，本設(shè)計（論文）不含任何其他個人或集體已經(jīng)發(fā)表或撰寫過的作品成果。作者簽名: 二〇一〇年九月二十日致 謝時間飛逝，大學的學習生活很快就要過去，在這四年的學習生活中，收獲了很多，而這些成績的取得是和一直關(guān)心幫助我的人分不開的。郭謙功老師淵博的知識、嚴謹?shù)淖黠L和誨人不倦的態(tài)度給我留下了深刻的印象?；厥姿哪?，取得了些許成績，生活中有快樂也有艱辛。他無論在理論上還是在實踐中，都給與我很大的幫助，使我得到不少的提高這對于我以后的工作和學習都有一種巨大的幫助，感謝他耐心的輔導。在我的十幾年求學歷程里，離不開父母的鼓勵和支持，是他們辛勤的勞作，無私的付出，為我創(chuàng)造良好的學習條件，我才能順利完成完成學業(yè)，感激他們一直以來對我的撫養(yǎng)與培育。另外，我還要感謝大學四年和我一起走過的同學朋友對我的關(guān)心與支持，與他們一起學習、生活，讓我在大學期間生活的很充實，給我留下了很多難忘的回憶。經(jīng)過這次畢業(yè)設(shè)計，我的能力有了很大的提高，比如操作能力、分析問題的能力、合作精神、嚴謹?shù)墓ぷ髯黠L等方方面面都有很大的進步。作者簽名: 二〇一〇年九月二十日對本文的研究做出重要貢獻的個人和集體，均已在文中以明確方式注明并表示感謝。 第5章 運行結(jié)果分析節(jié)點導納矩陣Y +j*( ) +j*( ) +j*( ) +j*( ) +j*( ) +j*( ) +j*() +j*( ) +j*( ) +j*( ) +j*( ) +j*( ) +j*() +j*( ) +j*( ) +j*( ) +j*( ) +j*( ) +j*() +j*( ) +j*( ) +j*( ) +j*( ) +j*( ) +j*() 第1次迭代結(jié)果第1次迭代的雅比矩陣J 第1次迭代的功率偏差dP和dQ+000 +000 +000 +000 +000 +000 第1次迭代的節(jié)點相角和電壓的偏差dx 第1