【正文】 if k~=3。k4=k1。for k3=3:k4for k2=k1:N1J(k3,k2)= J(k3,k2)J(k3,k)*J(k,k2)。endJ(k3,k)=0。endif k==N0,break。endfor k3=k1:N0for k2=k1:N1J(k3,k2)=J(k3,k2)J(k3,k)*J(k,k2)。endJ(k3,k)=0。endelsefor k3=k1:N0for k2=k1:N1J(k3,k2)= J(k3,k2)J(k3,k)*J(k,k2)。endJ(k3,k)=0。endendend for k=3:2:N01L=(k+1)./2。e(L)=e(L)J(k,N)。k1=k+1。f(L)=f(L)J(k1,N)。endfor k=3:N0DET=abs(J(k,N))。if DET=prIT2=IT2+1。endendICT2(a)=IT2。ICT1=ICT1+1。for k=1:ndy(k)=sqrt(e(k)^2+f(k)^2)。endfor i=1:nDy(ICT1,i)=dy(i)。endend %用高斯消去法解“w=J*V”disp(39。迭代次數(shù)39。)。disp(ICT1)。disp(39。沒有達(dá)到精度要求的個數(shù)39。)。disp(ICT2)。for k=1:nV(k)=sqrt(e(k)^2+f(k)^2)。O(k)=atan(f(k)./e(k))*180/pi。end E=e+f*j。disp(39。各節(jié)點的實際電壓標(biāo)么值E為(節(jié)點號從小到大排列):39。)。disp(E)。disp(39。各節(jié)點的電壓大小V為(節(jié)點號從小到大排列):39。)。disp(V)。disp(39。各節(jié)點的電壓相角時shita為(節(jié)點號從小到大排列):39。)。disp(O)。for p=1:nC(p)=0。for q=1:nC(p)=C(p)+conj(Y(p,q))*conj(E(q))。endS(p)=E(p)*C(p)。end disp(39。各節(jié)點的功率S為(節(jié)點號從小到大排列):39。)。disp(S)。disp(39。各條支路的首端功率Si為(順序同輸入B1時一樣):39。)。for i=1:nlif B1(i,6)==0p=B1(i,1)。q=B1(i,2)。else p=B1(i,2)。q=B1(i,1)。endSi(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5))conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))))。disp(Si(p,q))。enddisp (39。各條支路的末端功率Sj為(順序同輸入B1時一樣):39。)。for i=1:nlif B1(i,6)==0p=B1(i,1)。q=B1(i,2)。else p=B1(i,2)。q=B1(i,1)。endSj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5))conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))))。disp(Sj(q,p))。enddisp(39。各條支路的功率損耗DS為(順序同輸入B1時一樣):39。 )。for i=1:nlif B1(i,6)==0p=B1(i,1)。q=B1(i,2)。else p=B1(i,2)。q=B1(i,1)。endDS(i)=Si(p,q)+Sj(q,p)。disp(DS(i))。end 計算結(jié)果分析該算例的主接線圖來自《電力系統(tǒng)穩(wěn)態(tài)分析中》的一道例題，書中采用手算的方法，經(jīng)過3次迭代后得到如下的結(jié)果：========而仿真結(jié)果中====經(jīng)對比分析知Matlab仿真結(jié)果（見附錄B）與書中迭代的得出的結(jié)果基本相同。因此，可以證明改程序的正確性。河南城建學(xué)院本科畢業(yè)設(shè)計（論文） 結(jié)論 結(jié) 論電力系統(tǒng)潮流計算是研究電力系統(tǒng)穩(wěn)態(tài)運行的一項基本運算，它根據(jù)給定系統(tǒng)的網(wǎng)絡(luò)結(jié)構(gòu)及運行條件來確定整個系統(tǒng)的運行狀態(tài):主要是各節(jié)點電壓(幅值和相角)，網(wǎng)絡(luò)中功率分布和功率損耗等狀態(tài)。它既是對電力系統(tǒng)規(guī)劃和運行方式的合理性、可靠性及經(jīng)濟(jì)性進(jìn)行定量分析的依據(jù)又是電力系統(tǒng)穩(wěn)態(tài)和暫態(tài)穩(wěn)定計算的基礎(chǔ)，是電力系統(tǒng)一種非常重要和基本的計算。本設(shè)計是根據(jù)牛頓—拉夫遜法潮流計算原理，利用MATLAB的M語言編制使用計算機計算電力系統(tǒng)潮流的程序設(shè)計方案。該設(shè)計有計算方法可靠、收斂性好、具有通用性，占用較少的計算機內(nèi)存，計算速度快等優(yōu)點。通過本次設(shè)計，我更好的理解和掌握電力系統(tǒng)潮流計算的原理與方法，也更加系統(tǒng)的理解程序設(shè)計語言M語言的編程技巧，培養(yǎng)了我綜合利用程序設(shè)計語言M語言進(jìn)行科學(xué)計算，使所學(xué)知識轉(zhuǎn)化為分析和設(shè)計數(shù)學(xué)中實際問題的能力。提高建立程序文檔、歸納總結(jié)的能力。河南城建學(xué)院本科畢業(yè)設(shè)計（論文） 參考文獻(xiàn) 參考文獻(xiàn)[1]傅知蘭，電力系統(tǒng)電氣設(shè)備選擇與實用計算. 北京：中國電力出版社，2004[2]，北京：北京，中國水利水電出版社，1999[3]，中國電力出版社，2006[4]孟祥萍，電力系統(tǒng)分析，北京：高等教育出版社，2005[5] 徐政. ，機械工業(yè)出版社，2003[6].何仰贊. ，2002[3] 何仰贊，溫增銀. 電力系統(tǒng)分析. 武漢：華中科技大學(xué)出版社，2002[7] 宋文南，李樹鴻，. 天津：天津大學(xué)出版社，1990[8] 尹澤明，[M].北京，清華大學(xué)出版社，2002.[9] Enery System Theory：An York：McGrawHill Book Co.,1971（1st Edition），1982（2nd Edition）[10] :水利電力出版社，1985.[11] :水利電力出版社，1987.[12] 周全仁，:湖南科學(xué)技術(shù)出版社，1983.[13] 許主平，周少武，鄒軍安。電力系統(tǒng)計算機輔助設(shè)計。北京：中國電力出版社，2001。[14] [M].北京：中國電力出版社，2007. [15] 李維波. MATLAB 在電氣工程中的應(yīng)用［M］. 北京: 中國電力出版社，2007.[16] et practical application to system ：Pitman Publishing ，1976[17] Power Control in Electric York：John Wileyamp。Sons，1982河南城建學(xué)院本科畢業(yè)設(shè)計（論文） 附錄A附錄A 程序流程圖開 始輸入原始數(shù)據(jù)形成導(dǎo)納矩陣給定電壓初值f( 0 ) ， e( 0)置 γ=0對于PQ節(jié)點，計算ΔP（γ）,ΔQ（γ），對于PU節(jié)點，計算ΔP（γ），ΔU2 （γ）按系統(tǒng)的潮流分布計算平衡節(jié)點功率 和線路功率是否 | ΔPj ,ΔQi|ε？ 是 計算雅克比矩陣個元素輸出以（γ+1）→γ解修正方程，求Δe（γ）, Δf（γ）用e（γ+1）= e（γ）+eΔ（γ），f（γ+1）=f（γ）+Δf（γ）修正節(jié)點電壓以e（γ+1）→ e（γ） f（γ+1）→f（γ） 河南城建學(xué)院本科畢業(yè)設(shè)計（論文） 附錄B 附錄B Matlab仿真輸入的節(jié)點數(shù)，支路數(shù)和誤差精度形成的支路參數(shù)矩陣B1，節(jié)點參數(shù)矩陣B2 對地阻抗矩陣X及迭代結(jié)果 迭代結(jié)果河南城建學(xué)院本科畢業(yè)設(shè)計（論文） 致謝致謝本設(shè)計的完成是在我們的導(dǎo)師程漢蓬老師的細(xì)心指導(dǎo)下進(jìn)行的。在每次設(shè)計遇到問題時老師不辭辛苦的講解才使得我的設(shè)計順利的進(jìn)行。從設(shè)計的選題到資料的搜集直至最后設(shè)計的修改的整個過程中，花費了程老師很多的寶貴時間和精力，在此向?qū)煴硎局孕牡馗兄x！導(dǎo)師嚴(yán)謹(jǐn)?shù)闹螌W(xué)態(tài)度，開拓進(jìn)取的精神和高度的責(zé)任心都將使學(xué)生受益終生！還要感謝和我同一設(shè)計小組的幾位同學(xué)，是你們在我平時設(shè)計中和我一起探討問題，并指出我設(shè)計上的誤區(qū)，使我能及時的發(fā)現(xiàn)問題把設(shè)計順利的進(jìn)行下去，沒有你們的幫助我不可能這樣順利地結(jié)稿，在此表示深深的謝意。