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