【正文】 /Y(i,i)。 elseif(i=m+1) s_node(i)=P_PV(im)+Q_PV(im)*j。Q_PV=ones(1,nm)。%PQ節(jié)點(diǎn)無(wú)功功率不平衡量 elseif(i=m+1 amp。 a=m)%如果a號(hào)節(jié)點(diǎn)是平衡節(jié)點(diǎn)之后的PQ節(jié)點(diǎn),Ykb中的aa1。%PQ節(jié)點(diǎn)的Lia參數(shù) else if(a=s+1 amp。amp。%PQ節(jié)點(diǎn)Nia參數(shù) Ykb(2*(i1),2*(a1)1)=Ykb(2*(i1)1,2*(a1))。 a=m)%如果a號(hào)節(jié)點(diǎn)是平衡節(jié)點(diǎn)之后的PQ節(jié)點(diǎn) Ykb(2*(i1)1,2*(a1)1)=imag(Y(i,a))*real(U(i))real(Y(i,a))*imag(U(i))imag(I(i))。%PQ節(jié)點(diǎn)Nia參數(shù) Ykb(2*(i1),2*a1)=Ykb(2*(i1)1,2*a)。Ybk中的ii1。%PQ節(jié)點(diǎn)的Jia參數(shù)=Nia2real（I(i)）； Ykb(2*i1,2*(a1))=Ykb(2*i1,2*(a1)1)+2*imag(I(i))。%PV節(jié)點(diǎn)Lia參數(shù)=Hia end end elseif(a==i) if(a=s1)%如果a號(hào)節(jié)點(diǎn)是平衡節(jié)點(diǎn)之前的PQ節(jié)點(diǎn) Ykb(2*i1,2*a1)=imag(Y(i,a))*real(U(i))real(Y(i,a))*imag(U(i))imag(I(i))。 %PQ節(jié)點(diǎn)Nia參數(shù)Ykb(2*i,2*a1)=Ykb(2*i1,2*a)。%PQ節(jié)點(diǎn)無(wú)功功率不平衡量 else if(i=s+1 amp。%電壓不平衡量dumax=0。*U(i)39。A)/Y(i,i)。 ang=angle(Upv)。Ss=0。在進(jìn)行電力系統(tǒng)的其他分析，比如輸電線路的損耗，電力網(wǎng)絡(luò)的無(wú)功優(yōu)化，電力系統(tǒng)的故障分析，在這些分析中，潮流計(jì)算也可以得到廣泛應(yīng)用。（二）子程序的設(shè)計(jì)子程序主要包括：（1） 輸電線路數(shù)學(xué)模型計(jì)算子程序：line_para；（2） 變壓器數(shù)學(xué)模型計(jì)算子程序：transfer_para；（3） 支路導(dǎo)納矩陣計(jì)算子程序：Jmatrix；（4） 節(jié)點(diǎn)導(dǎo)納矩陣計(jì)算子程序：Ymatrix；（5） 牛頓——拉夫遜潮流計(jì)算子程序：DYnewton；（6） 高斯——賽德爾潮流計(jì)算子程序：DYgaosai；（7） 寫文件子程序：write；?！尽?，是一個(gè)電力網(wǎng)絡(luò)的拓?fù)浣Y(jié)構(gòu)。（2）框2中為了得到比較精確的節(jié)點(diǎn)電壓初值，要先調(diào)用高斯——賽德爾潮流計(jì)算子程序，以得到比較精確的節(jié)點(diǎn)電壓初值。如此不斷迭代，直至滿足對(duì)精確度的要求。 第一次迭代后的節(jié)點(diǎn)電壓節(jié)點(diǎn)編號(hào)123節(jié)點(diǎn)電壓++最大電壓變化量dumax為：dumax=。迭代過(guò)程中往往會(huì)出現(xiàn)越限，即按式（）求得的不能滿足的情況。再將、…代入第三式，又可解得。第三類節(jié)點(diǎn)稱平衡節(jié)點(diǎn)。其中對(duì)控制變量的約束條件是；對(duì)沒(méi)有電源的節(jié)點(diǎn)則為；這些限制條件取決于一系列技術(shù)經(jīng)濟(jì)因素，應(yīng)根據(jù)實(shí)際情況而定。其中、主要受、的控制；、主要受、的控制。潮流計(jì)算的額主要思想是迭代法。 三繞組變壓器的等值電路將三繞組變壓器化為兩個(gè)雙繞組變壓器后就可以按照雙繞組變壓器等值電路的求法進(jìn)行進(jìn)一步的化簡(jiǎn)。有的節(jié)點(diǎn)不與電流源相連，則其注入電流為零。牛頓—拉夫遜法，把非線性方程式的求解過(guò)程變成反復(fù)對(duì)相對(duì)應(yīng)的線性方程式的求解過(guò)程，通常稱為逐次線性化過(guò)程，這是牛頓—拉夫遜法的核心。目前，我國(guó)很多電力系統(tǒng)都采用了P—Q分解法潮流程序。阻抗法改善了系統(tǒng)潮流計(jì)算的收斂性問(wèn)題，解決了導(dǎo)納法無(wú)法求解的一些系統(tǒng)的潮流計(jì)算，在60年代獲得了廣泛的應(yīng)用，曾為我國(guó)的電力系統(tǒng)的設(shè)計(jì)、運(yùn)行和研究做出了很大的貢獻(xiàn)。電力系統(tǒng)的潮流計(jì)算的計(jì)算量非常巨大，通過(guò)人來(lái)計(jì)算是非常困難的。電力系統(tǒng)的潮流計(jì)算的電力系統(tǒng)穩(wěn)態(tài)分析、暫態(tài)分析和故障分析的基礎(chǔ)。因此，對(duì)潮流計(jì)算方法，首先要求它能可靠的收斂，并給出正確答案，由于電力系統(tǒng)結(jié)構(gòu)及參數(shù)的一些特點(diǎn)，并且隨著電力系統(tǒng)的不斷擴(kuò)大，潮流計(jì)算的方程式的階數(shù)越來(lái)越多。為了克服阻抗法在內(nèi)存和速度上的缺點(diǎn)，60年代中期發(fā)展了以阻抗矩陣為基礎(chǔ)的分塊阻抗法。后來(lái)又提出了根據(jù)直角坐標(biāo)形式的潮流方程是一個(gè)二次代數(shù)方程的特點(diǎn)，提出了采用直角坐標(biāo)的保留非線性快速潮流算法。牛頓——拉夫遜法具有很好的收斂性，計(jì)算速度快，計(jì)算結(jié)果準(zhǔn)確。在等值電路的求解過(guò)程中，輸電線路和變壓器的等值電路的求取是主要的工作。與節(jié)點(diǎn)2對(duì)應(yīng)的第二行非零非對(duì)角元數(shù)為2.（2）節(jié)點(diǎn)導(dǎo)納矩陣的對(duì)角元就等于與該節(jié)點(diǎn)所連接的導(dǎo)納的總和。于是， () 簡(jiǎn)單系統(tǒng)及其等值網(wǎng)絡(luò)（a）簡(jiǎn)單系統(tǒng)；（b）簡(jiǎn)單系統(tǒng)的等值網(wǎng)絡(luò)；（c）注入功率和注入電流 ()如令 ()并將式（）代入式（）展開，將有功功率和無(wú)功功率分列，可得 ()式（）就是這個(gè)簡(jiǎn)單系統(tǒng)的功率方程。換言之，擾動(dòng)向量、控制向量、狀態(tài)向量，都是階的列向量。這其實(shí)意味著讓這些電源調(diào)節(jié)它們發(fā)出的無(wú)功功率以保證與之連接的節(jié)點(diǎn)電壓為定值。進(jìn)行計(jì)算時(shí)，平衡節(jié)點(diǎn)是不可少的；PQ節(jié)點(diǎn)是大量的，PV節(jié)點(diǎn)是較少的，甚至可以沒(méi)有。先代入第一式，解得。 高斯——賽德爾程序框圖高斯——。迭代計(jì)算結(jié)束后，再根據(jù)式（）和式（）計(jì)算PV節(jié)點(diǎn)的無(wú)功功率和平衡節(jié)點(diǎn)的視在功率，并根據(jù)式組（）計(jì)算線路上的視在功率，根據(jù)式（）計(jì)算線路上的功率損耗。 ()式中，第一部分為給定的節(jié)點(diǎn)注入功率，第二部分為由節(jié)點(diǎn)電壓求得的節(jié)點(diǎn)注入功率，它們二者的差就是節(jié)點(diǎn)功率的不平衡量。 牛頓——拉夫遜潮流計(jì)算流程圖（三）算例分析【】。，+j0；PV節(jié)點(diǎn)5的有功功率為5，；PQ節(jié)點(diǎn)+、+、2+j1。)Y=zeros(n,n)。參考文獻(xiàn)參考文獻(xiàn)[1] 陳珩. 電力系統(tǒng)穩(wěn)態(tài)分析 [M]. 南京：南京工學(xué)院，1985051.[2] 王錫凡. 現(xiàn)代電力系統(tǒng)分析 [M]. 北京：科學(xué)出版社，20033.[3] 西安交通大學(xué). 電子數(shù)字計(jì)算機(jī)的應(yīng)用：電力系統(tǒng)計(jì)算 [M]. 西安：西安交通大學(xué)，1987101.[4] 南京工學(xué)院. 電力系統(tǒng) [M]. 南京：南京工學(xué)院，1980081.[5] 吳天明，彭彬. MATLAB電力系統(tǒng)設(shè)計(jì)與分析 [M]. 國(guó)防工業(yè)出版社，200401.[6] 林暉. VSCMTDC和交流混合的電力系統(tǒng)潮流計(jì)算 [D]. 湖南大學(xué)電氣與信息工程學(xué)院，.[7] 梁海平. 智能潮流分析系統(tǒng)的研究與開發(fā) [D]. 華北電力大學(xué)電力系統(tǒng)及其自動(dòng)化專業(yè)，.[8] 許建明. 經(jīng)典法和遺傳算法在電力系統(tǒng)有功優(yōu)化中的研究 [D]. 南昌大學(xué)信息工程學(xué)院，.[9] 趙芳 陳士方 張圣集. 保留非線性的電力系統(tǒng)概率潮流計(jì)算 [J]. 電力情報(bào)，2000年03期.[10] 成海濱 沈茂亞. 電力系統(tǒng)最優(yōu)潮流算法研究綜述 [J]. 電力電氣，2006年第25卷第11期.[11] 羅瑋. 電力市場(chǎng)環(huán)境下基于最優(yōu)潮流的輸電容量充裕度研究 [D]. 上海交通大學(xué)電力系統(tǒng)及其自動(dòng)化專業(yè)，[12] 郭金偉. 電力市場(chǎng)環(huán)境下基于最優(yōu)潮流的節(jié)點(diǎn)實(shí)時(shí)電價(jià)和購(gòu)電份額研究 [D]. 北京交通大學(xué)，2008.[13] 劉方. 電力系統(tǒng)動(dòng)態(tài)最優(yōu)潮流的模型與算法研究[D]. 重慶大學(xué)電氣工程學(xué)院，[14] 郗忠梅 李有安 趙法起 張博. 基于Matlab的電力系統(tǒng)潮流計(jì)算[J]. 山東農(nóng)業(yè)大學(xué)學(xué)報(bào)，2010年41卷第2期.[15] 樊宇璐 李世作 張志斌. 基于擬牛頓法的電力系統(tǒng)潮流計(jì)算[J]. 電氣開關(guān)，2009年第2期.[16] Yan Sun, NaiShan Flow Calculation Method for Islanded Power Network [J]. Department of Electrical Engineering Guangxi University ,Nanning, China,.[17] Andrey Pazderin, Sergey Flow Calculation by Combination ofNewtonRaphson Method and Newton’s Method in Optimization.[J].URAL STATE TECHNICAL UNIVERSITY – UPI,[18]D. A. Arzamastsev, P. I. Bartolomej, A. M. Holjan, Automatic control system and optimization of modes of electric power systems, Moscow:Vysshaya Shkola, 1983.[19]V. I. Tarasov, Nonlinear minimization methods for calculating steady states of electric power systems, – Novosibirsk: Nauka, 2001, p. 214.[20]A. Z. Gamm, Statistical methods of power system state estimation, Moscow: Nauka, 1976, p. 220.[21]V. M. Gornshtejn, B. P. Miroshnichenko, . Ponomarev, Methods of power system mode optimization, Moscow: Energiya, 1981, p. 336.附錄附錄附錄1 源程序1 高斯——賽德爾潮流計(jì)算源程序%高斯——賽爾德潮流計(jì)算%[s_node,u_node]=gaosai(Y,m,n,s,S_PQ,P_PV,U_PV,us,Qmax,Qmin,expl)%返回值：s_node: 節(jié)點(diǎn)功率% u_node： 節(jié)點(diǎn)電壓%輸入值： Y： 導(dǎo)納矩陣% m： PQ節(jié)點(diǎn)數(shù),包含平衡節(jié)點(diǎn)% n： 總結(jié)點(diǎn)數(shù)% s： 平衡節(jié)點(diǎn)編號(hào)，一般定為1號(hào)編號(hào)% S_PQ: PQ節(jié)點(diǎn)的視在功率[1:m],其中平衡節(jié)點(diǎn)的視在功率暫定為0；% P_PV：PV節(jié)點(diǎn)的有功功率[m+1:n]% U_PV: PV節(jié)點(diǎn)的電壓給定值,只有幅值，沒(méi)有相位角，是實(shí)數(shù)[m+1:n]% us: 平衡節(jié)點(diǎn)的電壓給定值，復(fù)數(shù)。 else %第i節(jié)點(diǎn)不是平衡節(jié)點(diǎn)，進(jìn)行以下處理 A=0。PV節(jié)點(diǎn)%gQ大于Qmax轉(zhuǎn)化為PQ節(jié)點(diǎn)\n39。 end end end du=abs(U(i)u)。 end end end s_node u_node=U2 牛頓——拉夫遜潮流計(jì)算程序%牛頓——拉夫遜潮流計(jì)算。%PV節(jié)點(diǎn)的無(wú)功功率；Q_PV=Q_PV.*。amp。i不變 Ykb(2*i1,2*(a1)1)=imag(Y(i,a))*real(U(i))real(Y(i,a))*imag(U(i))。amp。 i=m)%如果i號(hào)節(jié)點(diǎn)是平衡節(jié)點(diǎn)之后的PQ節(jié)點(diǎn)。%PQ節(jié)點(diǎn)的Jia參數(shù)=Nia Ykb(2*(i1),2*(a1))=Ykb(2*(i1)1,2*(a1)1)。%PQ節(jié)點(diǎn)Hia參數(shù) Ykb(2*(i1)1,2*(a1))=real(Y(i,a))*real(U(i))imag(Y(i,a))*imag(U(i))real(I(i))。%PQ節(jié)點(diǎn)Hia參數(shù) Ykb(2*(i1)1,2*a)=real(Y(i,a))*real(U(i))imag(Y(i,a))*imag(U(i))real(I(i))。 a=m)%如果a號(hào)節(jié)點(diǎn)是平衡節(jié)點(diǎn)之后的PQ節(jié)點(diǎn)。%PQ節(jié)點(diǎn)Nia參數(shù) Ykb(2*i,2*(a1)1)=Ykb(2*i1,2*(a1))2*real(I(i))。%PV節(jié)點(diǎn)Jia參數(shù)=NiaYkb(2*i1,2*(a1))=Ykb(2*i1,2*(a1)1)。 %PQ節(jié)點(diǎn)Hia參數(shù) Ykb(2*i1,2*a)=real(Y(i,a))*real(U(i))imag(Y(i,a))*imag(U(i))。%PQ節(jié)點(diǎn)有功功率不平衡量 dPQV2(2*i)=imag(S_PQ(i))imag(A)。dU=zeros(2*(n1),1)。)endfor