【正文】 6=0。p=2*i1。q=2*j11。J(p,q)=X5。J(p,N)=DV。 m=p+1。J(m,q)=X1。J(m,N)=DP。q=q+1。J(p,q)=X6。J(m,q)=X2。elseif j1==iamp。j1~=isbX1=C(i)G(i,i)*e(i)B(i,i)*f(i)。X2=D(i)+B(i,i)*e(i)G(i,i)*f(i)。X5=2*e(i)。X6=2*f(i)。p=2*i1。q=2*j11。J(p,q)=X5。J(p,N)=DV。 m=p+1。J(m,q)=X1。J(m,N)=DP。q=q+1。J(p,q)=X6。J(m,q)=X2。endendendendend %求雅可比矩陣for k=3:N0k1=k+1。N1=N。for k2=k1:N1J(k,k2)=J(k,k2)./J(k,k)。endJ(k,k)=1。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。沒有達到精度要求的個數(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é)點的實際電壓標么值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 計算結果分析該算例的主接線圖來自《電力系統(tǒng)穩(wěn)態(tài)分析中》的一道例題，書中采用手算的方法，經(jīng)過3次迭代后得到如下的結果：========而仿真結果中====經(jīng)對比分析知Matlab仿真結果（見附錄B）與書中迭代的得出的結果基本相同。因此，可以證明改程序的正確性。河南城建學院本科畢業(yè)設計（論文） 結論 結 論電力系統(tǒng)潮流計算是研究電力系統(tǒng)穩(wěn)態(tài)運行的一項基本運算，它根據(jù)給定系統(tǒng)的網(wǎng)絡結構及運行條件來確定整個系統(tǒng)的運行狀態(tài):主要是各節(jié)點電壓(幅值和相角)，網(wǎng)絡中功率分布和功率損耗等狀態(tài)。它既是對電力系統(tǒng)規(guī)劃和運行方式的合理性、可靠性及經(jīng)濟性進行定量分析的依據(jù)又是電力系統(tǒng)穩(wěn)態(tài)和暫態(tài)穩(wěn)定計算的基礎，是電力系統(tǒng)一種非常重要和基本的計算。本設計是根據(jù)牛頓—拉夫遜法潮流計算原理，利用MATLAB的M語言編制使用計算機計算電力系統(tǒng)潮流的程序設計方案。該設計有計算方法可靠、收斂性好、具有通用性，占用較少的計算機內存，計算速度快等優(yōu)點。通過本次設計，我更好的理解和掌握電力系統(tǒng)潮流計算的原理與方法，也更加系統(tǒng)的理解程序設計語言M語言的編程技巧，培養(yǎng)了我綜合利用程序設計語言M語言進行科學計算，使所學知識轉化為分析和設計數(shù)學中實際問題的能力。提高建立程序文檔、歸納總結的能力。河南城建學院本科畢業(yè)設計（論文） 參考文獻 參考文獻[1]傅知蘭，電力系統(tǒng)電氣設備選擇與實用計算. 北京：中國電力出版社，2004[2]，北京：北京，中國水利水電出版社，1999[3]，中國電力出版社，2006[4]孟祥萍，電力系統(tǒng)分析，北京：高等教育出版社，2005[5] 徐政. ，機械工業(yè)出版社，2003[6].何仰贊. ，2002[3] 何仰贊，溫增銀. 電力系統(tǒng)分析. 武漢：華中科技大學出版社，2002[7] 宋文南，李樹鴻，. 天津：天津大學出版社，1990[8] 尹澤明，[M].北京，清華大學出版社，2002.[9] Enery System Theory：An York：McGrawHill Book Co.,1971（1st Edition），1982（2nd Edition）[10] :水利電力出版社，1985.[11] :水利電力出版社，1987.[12] 周全仁，:湖南科學技術出版社，1983.[13] 許主平，周少武，鄒軍安。電力系統(tǒng)計算機輔助設計。北京：中國電力出版社，2001。[14] [M].北京：中國電力出版社，2007. [15] 李維波. MATLAB 在電氣工程中的應用［M］. 北京: 中國電力出版社，2007.[16] et practical application to system ：Pitman Publishing ，1976[17] Power Control in Electric York：John Wileyamp。Sons，1982河南城建學院本科畢業(yè)設計（論文） 附錄A附錄A 程序流程圖開 始輸入原始數(shù)據(jù)形成導納矩陣給定電壓初值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（γ） 河南城建學院本科畢業(yè)設計（論文） 附錄B 附錄B Matlab仿真輸入的節(jié)點數(shù)，支路數(shù)和誤差精度形成的支路參數(shù)矩陣B1，節(jié)點參數(shù)矩陣B2 對地阻抗矩陣X及迭代結果 迭代結果河南城建學院本科畢業(yè)設計（論文） 致謝致謝本設計的完成是在我們的導師程漢蓬老師的細心指導下進行的。在每次設計遇到問題時老師不辭辛苦的講解才使得我的設計順利的進行。從設計的選題到資料的搜集直至最后設計的修改的整個過程中，花費了程老師很多的寶貴時間和精力，在此向導師表示衷心地感謝！導師嚴謹?shù)闹螌W態(tài)度，開拓進取的精神和高度的責任心都將使學生受益終生！還要感謝和我同一設計小組的幾位同學，是你們在我平時設計中和我一起探討問題，并指出我設計上的誤區(qū)，使我能及時的發(fā)現(xiàn)問題把設計順利的進行下去，沒有你們的幫助我不可能這樣順利地結稿，在此表示深深的謝意。