【正文】 第6次迭代的節(jié)點(diǎn)相角delta和電壓Uang1 ang2 ang3 ang4 U1 U2 U3 牛頓－拉夫遜法潮流計算結(jié)果 節(jié)點(diǎn)計算結(jié)果：節(jié)點(diǎn) 節(jié)點(diǎn)電壓 節(jié)點(diǎn)相角（角度） 節(jié)點(diǎn)注入功率 1 +j 2 +j 3 +j 4 +j 5 +j 線路計算結(jié)果：節(jié)點(diǎn)I 節(jié)點(diǎn)J 線路功率S(I,J) 線路功率S(J,I) 線路損耗dS(I,J) 1 2 +j +j +j 1 3 +j +j +j 2 3 +j +j +j 5 3 +j +j +j 4 2 +j +j +j 此圖中的容量都是以100MW為基準(zhǔn)值。 1 3 0。 3 1。, 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)))。節(jié)點(diǎn)I 節(jié)點(diǎn)J 線路功率S(I,J) 線路功率S(J,I) 線路損耗dS(I,J)\n39。endfprintf(myf, 39。)。 節(jié)點(diǎn)計算結(jié)果：\n39。 fprintf(myf, 39。global nl。%bus矩陣中的節(jié)點(diǎn)號恢復(fù)原來的編號 tem(I,1)=I。global nodenum。global nb。 %求取節(jié)點(diǎn)注入功率 Sji=Uj*conj(Iji)。 end if K0 %變壓器線路: Zt和Ym為折算到K側(cè)的值,K在i側(cè) K=K。(J==0) %接地支路 Iij=Ui*Ym。(J~=0) %普通線路非接地支路。 K=line(k,7)。 %用lPQ矩陣第一、二列保存線路兩端節(jié)點(diǎn)編號 lPQ(k,2)=J。global line。 end bus(I,5)=sum。 %求取的平衡節(jié)點(diǎn)的P存回bus矩陣 end sum=0。 for J=1:nb ang=bus(I,3)bus(J,3)。global bus。超過最大迭代次數(shù)，不收斂停機(jī)！39。, I,u(I,1))。, I,angl(I,1))。 %將合一起的x，拆成相角、電壓兩個列向量 u=x(nb:nb+nPQ1,1)。 %修正 fprintf(myf,39。\n39。 for I=1:nb+nPQ1 fprintf(myf, 39。)。%10e 39。 end fprintf(myf,39。, Jac(I,J))。第%d次迭代的雅比矩陣J\n39。 %求得dx fprintf(myf,39。u]。)。%10f+j*(%10f) 39。 fprintf(myf, 39。global myf。global nPQ。 Y(J,J)=Y(J,J)+K*K*Yt。 Y(J,J)=Y(J,J)+Yt/K/K。 end if (K==0)amp。 Y(I,I)=Y(I,I)+Yt+Ym。 %非接地支路為阻抗的倒數(shù) end Ym=line(k,5)+j*line(k,6)。 %對導(dǎo)納矩陣賦初值0 for k=1:nl I=line(k,1)。global line。 for I=1:nl %根據(jù)節(jié)點(diǎn)號對應(yīng)表重新排列線路矩陣節(jié)點(diǎn)編號 for J=1:2 for k=1:nb if line(I,J)==nodenum(k,2) line(I,J)=nodenum(k,1)。SW]。else %求PQ節(jié)點(diǎn)數(shù) nPQ = nPQ + 1。 if type == 3, %求平衡節(jié)點(diǎn)數(shù) nSW = nSW + 1。 nSW = 0。global nPQ。w39。 endoutfile=strcat(outfile,39。for I=1:length(dfile) if dfile(I)==39。 %執(zhí)行數(shù)據(jù)文件中語句end [nl,ml]=size(line)。 %打開數(shù)據(jù)文件if pathname == 0 error(39。*.m39。global line。 %計算節(jié)點(diǎn)注入功率 flow。 %打開bus line數(shù)據(jù)文件 change。 %Y:節(jié)點(diǎn)導(dǎo)納矩陣global nodenum。 %節(jié)點(diǎn)個數(shù)global nl。第四章程序流程圖及程序代碼根據(jù)牛頓—拉夫遜法潮流計算的基本步驟可以得到如下程序流程圖%*******************主程序開始********************clc。 形成節(jié)點(diǎn)導(dǎo)納矩陣。同時，為了實(shí)時監(jiān)控電力系統(tǒng)的運(yùn)行狀態(tài)，也需要進(jìn)行大量而快速的潮流計算。 (2)在編制年運(yùn)行方式時,在預(yù)計負(fù)荷增長及新設(shè)備投運(yùn)基礎(chǔ)上,選擇典型方式進(jìn)行潮流計算,發(fā)現(xiàn)電網(wǎng)中薄弱環(huán)節(jié),供調(diào)度員日常調(diào)度控制參考,并對規(guī)劃、基建部門提出改進(jìn)網(wǎng)架結(jié)構(gòu),加快基建進(jìn)度的建議。隨著工農(nóng)業(yè)生產(chǎn)和城市的發(fā)展，電能的需要量迅速增加，而熱能資源和水能資源豐富的地區(qū)又往往遠(yuǎn)離用電比較集中的城市和工礦區(qū)，為了解決這個矛盾，就需要在動力資源豐富的地區(qū)建立大型發(fā)電站，然后將電能遠(yuǎn)距離輸送給電力用戶。通過對此電力系統(tǒng)分析，我們很容易得以節(jié)點(diǎn)3為PQ節(jié)點(diǎn)，節(jié)點(diǎn)4為PV節(jié)點(diǎn)，節(jié)點(diǎn)5為平衡節(jié)點(diǎn)。運(yùn)用以極坐標(biāo)表示的高斯—塞德爾法、牛頓—拉夫遜法、PQ分解法計算該系統(tǒng)的潮流分布。 設(shè)計題目高壓輸電網(wǎng)潮流的計算機(jī)算法程序設(shè)計（PQ分解法、牛頓拉夫遜法）或中壓配電網(wǎng)潮流的計算機(jī)算法程序設(shè)計（前推后代法、同倫延拓法等）或電力系統(tǒng)短路故障的計算機(jī)算法程序設(shè)計（要求不限） 設(shè)計內(nèi)容，寫出節(jié)點(diǎn)導(dǎo)納矩陣； （直角坐標(biāo)系形式）初值后，求解不平衡量；； ，重新修改初值，從2開始重新循環(huán)計算； ，再計算各支路功率分布、功率損耗和平衡節(jié)點(diǎn)功率。本人授權(quán)　　 　 　大學(xué)可以將本學(xué)位論文的全部或部分內(nèi)容編入有關(guān)數(shù)據(jù)庫進(jìn)行檢索，可以采用影印、縮印或掃描等復(fù)制手段保存和匯編本學(xué)位論文。除了文中特別加以標(biāo)注引用的內(nèi)容外，本論文不包含任何其他個人或集體已經(jīng)發(fā)表或撰寫的成果作品。盡我所知，除文中特別加以標(biāo)注和致謝的地方外，不包含其他人或組織已經(jīng)發(fā)表或公布過的研究成果，也不包含我為獲得 及其它教育機(jī)構(gòu)的學(xué)位或?qū)W歷而使用過的材料。對本研究提供過幫助和做出過貢獻(xiàn)的個人或集體，均已在文中作了明確的說明并表示了謝意。對本文的研究做出重要貢獻(xiàn)的個人和集體，均已在文中以明確方式標(biāo)明。涉密論文按學(xué)校規(guī)定處理。 ；連調(diào)； 。計算精度要求各節(jié)點(diǎn)電壓修正量不大于。綜合比較牛頓拉夫遜法（直角坐標(biāo)、極坐標(biāo)）、PQ分解法等多種求解方法的特點(diǎn)，最后確定采用牛頓拉夫遜法（極坐標(biāo)）。同時，為了提高供電的可靠性以及資源利用的綜合經(jīng)濟(jì)性，又把許多分散的各種形式的發(fā)電站，通過送電線路和變電所聯(lián)系起來。 (3)正常檢修及特殊運(yùn)行方式下的潮流計算,用于日運(yùn)行方式的編制,指導(dǎo)發(fā)電廠開機(jī)方式,有功、無功調(diào)整方案及負(fù)荷調(diào)整方案,滿足線路、變壓器熱穩(wěn)定要求及電壓質(zhì)量要求。因此，潮流計算是電力系統(tǒng)中應(yīng)用最廣泛、最基本和最重要的一種電氣運(yùn)算。 設(shè)各節(jié)點(diǎn)電壓的初始值 和 求出功率偏差根據(jù)公式算出雅可比矩陣。clear allformat long。 %支路數(shù) global bus。 %節(jié)點(diǎn)編號對應(yīng)表global lPQ。 %節(jié)點(diǎn)重新編號 ybus。 %計算線路功率及損耗 ret。global bus。,39。 you must select a valid data file39。 %求線路數(shù)nl[nb,mb]=size(bus)。.39。.txt39。)。global bus。