【正文】 。電壓迭代次數(shù)曲線39。電壓39。以下是每次迭代后各節(jié)點(diǎn)的電壓值（如圖所示）‘）；plot(Cs,Dy),xlabel(39。 end for i=1:ICT1 Cs(i)=i。 end DS(i)=Si(p,q)+Sj(q,p)。 else p=B1(i,2)。各條支路的功率損耗DS為（順序同您輸入B1時(shí)一樣）：’）；for i=1:n1 if B1(i,6)==0 p=B1(i,1)。disp(Sj(q,p))。q=B1(i,1)。q=B1(i,2)。enddisp(39。 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))))。 else p=B1(i,2)。各條支路的首端功率Si為（順序同您輸入B1時(shí)一樣）：‘）； for i=1:n1 if B1（i,6）==0 p=B1(i,1)。各節(jié)點(diǎn)的功率S為（節(jié)點(diǎn)號從小到大排列）：‘）； disp(S)。 end S(p)=E(p)*C(p)。for p=1:n C(p)=0。)。disp(39。)。disp(39。)。disp(39。 O(k)=atan(f(k)./e(k))*180./pi。disp(ICT2)。沒有達(dá)到精度要求的個(gè)數(shù)39。disp(ICT1)。迭代次數(shù)39。 endfor i=1:n Dy(ICT1,i)=dy(i)。for k=1:n dy(k)=sqrt(e(k)^2+f(k)^2)。end for k=3:N0 DET=abs(J(k,N))。 k1=k+1。 end endendfor k=3:2:N01 L=(k+1)./2。 end end for k3=k1:N0 for k2=k1:N1 J(k3,k2)=J(k3,k2)J(k3,k)*J(k,k2)。 for k3=3:k4 for k2=k1:N1 J(k3,k2)=J(k3,k2)J(k3,k)*J(k,k2)。endJ(k,k)=1。N1=N。J(m,q)=X2。q=q+1。 J(m,q)=X1。J(p,N)=DV。q=2*j11。 X6=2*f(i)。 X2=D(i)+B(i,i)*e(i)G(i,i)*f(i)。 elseif j1==iamp。J(p,q)=X6。J(m,N)=DP。m=p+1。J(p,q)=X5。 p=2*i1。 X5=0。j1~=i X1=G(i,j1)*e(i)B(i,j1)*f(i)。 DV=V(i)~2V2。J(m,q)=X2。q=q+1。 J(m,q)=X1。J(p,N)=DQ。q=2*j11。 X4=X1。 X2=B(i,j1)*e(i)G(i,j1)*f(i)。 for j1=1:n if j1~=isbamp。 if B2(i,6)~=3 DP=P(i)P1。P,Q39。 Q1=f(i)*C(i)D(i)*e(i)。 D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1)。 D(i)=0。 for i=1:n。while IT2~=0 IT2=0。N=NO+1。IT2=1。Q=imag(S)。 B(i,i)=B(i,i)+B2(i,5)。 V(i)=B2(i,4)。 for i=1:n e(i)=real(B2(i,3))。 end %求導(dǎo)納矩陣 G=real(Y)。 Y(p,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2。 end Y(p,q)=Y(p,q)1./(B1(i,3)*B1(i,5)。 else p=B1(i,2)。 end end for i=1:n1 if B1(i,6)==0 p=B1(i,1)。 p=X(i,1)。S1=zeros(n1)。V=seros(1,n)。e=zeros(1,n)。)。X=input(39。請輸入各節(jié)點(diǎn)參數(shù)形成的矩陣：B2=39。)。B1=input(39。請輸入誤差精度：pr=39。)。isb=input(39。請輸入支路數(shù)：n1=39。)。附錄Matlab仿真程序n=input(39。在此，謹(jǐn)向?qū)熖评蠋熤乱猿绺叩木匆夂椭孕牡母兄x！學(xué)部領(lǐng)導(dǎo)孫親錫主任、電氣工程及其自動化系主任劉沛、輔導(dǎo)員張甸老師等為我提供了良好的研究條件，謹(jǐn)向各位表示誠摯的敬意和謝忱。給我提供了大量數(shù)據(jù)資料和建議，告訴我應(yīng)該注意的細(xì)節(jié)問題，細(xì)心的給我指出錯(cuò)誤，修改論文。rQinR, YH. Song, Gh. GrigoraT, Optimal Operation and Planning of Power Systems, Casa de EditurR VENUS, IaTi, 2003致 謝值此本科論文完成之際，首先要感謝我的導(dǎo)師唐萃老師。參考文獻(xiàn)[1] 于永源，楊綺雯．電力系統(tǒng)分析(第二版)．北京:中國電力出版社,2004[2] 劉從愛．電力工程基礎(chǔ)．濟(jì)南：山東科學(xué)技術(shù)出版社，1997[3] 《電氣工程師手冊》第二版編輯委員會．電氣工程師手冊．北京：機(jī)械工業(yè)出版社，2000[4] 何仰贊，溫增銀．：華中科技大學(xué)出版社，2002[5] 吳際舜,侯志撿. 電力系統(tǒng)潮流計(jì)算的計(jì)算機(jī)方法. 上海：上海交通大學(xué)出版社, 2000 [6] 陳珩. 電力系統(tǒng)穩(wěn)態(tài)分析（第三版）. 北京：中國電力出版社, 2007[7] 孫麗華．電力工程基礎(chǔ)．北京：機(jī)械工業(yè)出版社：2006[8] 李維波. ：中國電力出版社，2007[9] 陳懇．直角坐標(biāo)牛頓一拉夫遜法潮流計(jì)算新解法[J]．電力系統(tǒng)及其自動化學(xué)報(bào)，1999，[10] 張錚．MATLAB程序設(shè)計(jì)與實(shí)例應(yīng)用[M]．北京：中國鐵道出版社，2O03[11] 張志涌,徐彥琴. ：北京航空航天大學(xué)出版社,2001[12] A. Abur, F. Magnago, Y. Lu, Educational Toolboxes for Power System Analysis, IEEE Computer Applications in Power,2000,13(4):3135[13] F. Milano, An Open Source Power System Analysis Toolbox, IEEE Trans. on Power Systems, 2005, 20(3) :1199 1206.[14] . Zimmerman, . MurrilloS225。不管采用什么算法，所有的潮流計(jì)算都是基于矩陣的迭代運(yùn)算。通過實(shí)例計(jì)算分析，取得了比較滿意的效果。對于不能求逆的矩陣我們通過在電力系統(tǒng)中至少有一條支路上有接地支路來實(shí)現(xiàn)其求逆。結(jié)束語在電力系統(tǒng)的潮流計(jì)算算法中牛頓一拉夫遜法是得到電力系統(tǒng)研究人員認(rèn)可的算法之一，在本文中我們采用牛頓一拉夫遜法，主要是同時(shí)考慮到計(jì)算的準(zhǔn)確和程序的運(yùn)行速度。（5）邏輯關(guān)系運(yùn)算 邏輯運(yùn)算是MATLAB中數(shù)組運(yùn)算所特有的一種運(yùn)算形式，也是幾乎所有的高級語言普遍適用的一種運(yùn)算。有了“對應(yīng)元素”的規(guī)定，數(shù)組的運(yùn)算實(shí)質(zhì)上就是針對數(shù)組內(nèi)部的每個(gè)元素進(jìn)行的。在數(shù)組運(yùn)算中有了“對應(yīng)關(guān)系”的規(guī)定，數(shù)組與常數(shù)之間的除法運(yùn)算沒有任何限制。而乘除法運(yùn)算有相當(dāng)大的區(qū)別，數(shù)組的乘除法是指兩同維數(shù)組對應(yīng)元素之間的乘除法，它們的運(yùn)算符為“.*”和“./”或“.\”。這種運(yùn)算不同于前面講的數(shù)學(xué)運(yùn)算，為有所區(qū)別，我們稱之為數(shù)組運(yùn)算。但需注意進(jìn)行數(shù)除時(shí)，常數(shù)通常只能做除數(shù)。在MATLAB6中兩者的區(qū)別不太大。在傳統(tǒng)的MATLAB算法中，右除是先計(jì)算矩陣的逆再相乘，而左除則不需要計(jì)算逆矩陣直接進(jìn)行除運(yùn)算。矩陣的基本數(shù)學(xué)運(yùn)算包括矩陣的四則運(yùn)算、與常數(shù)的運(yùn)算、逆運(yùn)算、行列式運(yùn)算、秩運(yùn)算、特征值運(yùn)算等基本函數(shù)運(yùn)算，這里進(jìn)行簡單介紹。原始數(shù)據(jù)輸入格式的設(shè)計(jì)，主要應(yīng)從使用的角度出發(fā)，原則是簡單明了，便于修改。不同領(lǐng)域、不同層次的用戶通過相應(yīng)工具的學(xué)習(xí)和應(yīng)用，可以方便地進(jìn)行計(jì)算、分析及設(shè)計(jì)工作。通過M語言，可以用類似數(shù)學(xué)公式的方式來編寫算法，大大降低了程序所需的難度并節(jié)省了時(shí)間，從而可把主要的精力集中在算法的構(gòu)思而不是編程上。它可以高效率地解決工業(yè)計(jì)算問題，特別是關(guān)于矩陣和矢量的計(jì)算。 MATLAB是一種交互式、面向?qū)ο蟮某绦蛟O(shè)計(jì)語言，廣泛應(yīng)用于工業(yè)界與學(xué)術(shù)界，主要用于矩陣運(yùn)算，同時(shí)在數(shù)值分析、自動控制模擬、數(shù)字信號處理、動態(tài)分析、繪圖等方面也具有強(qiáng)大的功能?，F(xiàn)有很多潮流計(jì)算方法。 (314)雅可比矩陣各元素的算式式(314)中, 雅可比矩陣中的各元素可通過對式(310)和(311), 雅可比矩陣中非對角元素為 (315)當(dāng)時(shí),雅可比矩陣中對角元素為: (316)由式(315)和(316)看出,雅可比矩陣的特點(diǎn)：矩陣中各元素是節(jié)點(diǎn)電壓的函數(shù),在迭代過程中,這些元素隨著節(jié)點(diǎn)電壓的變化而變化；導(dǎo)納矩陣中的某些非對角元素為零時(shí),則必有；雅可比矩陣不是對稱矩陣。以直角坐標(biāo)系形式表示：迭代推算式采用直角坐標(biāo)時(shí)