【正文】 【J】．繼電器，2004，32(21)：7l一75．wANG Hui，CHEN Kai，PENG zhe，et a1．Application of digital simulation眥hniques棚d severaJ simulation tools in power system[J】．Relay，2004，32(21)：71N)==0){return 0。}Bus::ScanfBusData(bus)。double **JacAug。//計算節(jié)點功率Bus::JisuanNodeScal(X,bus,N)。//修正節(jié)點電壓// Bus::PrintfNodeV(bus,N)。statemp=line[i].start。bb=Complex::subComplex(Complex::getconj(bus[endtemp1].V), Complex::getconj(bus[statemp1].V))。bus[i].V=Complex::Rec2Polar(bus[i].V)。endtemp=line[i].end。//輸入支路個數if((Lamp。ibus[i].=0。//計算節(jié)點導納矩陣Bus::PrintfNodeDnz(X,N)。LOOP:Bus::JisuanNodeI(X,bus,N)。gauss::gauss_slove(JacAug,x,NN)。=0。line[i].stoe=Complex::productComplex(bus[statemp1].V,dd)。i{bus[i]. = bus[i]. + bus[i].。i{int statemp,endtemp。double image。//求一個復數共軛static Complex getinverse(Complex c1)。} }。=+。=*+*。=(r1)。=0。cout”。static void gauss_output()。ipivot= fabs(a[k][k])。//行pivcol[k][1]=j。i{al=a[i][pivcol[k][1]]。//解方程for(i=(n2)。k){al=x[pivcol[k][1]]。00 0 節(jié)點4。00 節(jié)點8 節(jié)點9 0 節(jié)點10 0 節(jié)點11 節(jié)點12 0 節(jié)點13 節(jié)點14 0 節(jié)點15 0 節(jié)點16 0 節(jié)點17 0 節(jié)點18 0 節(jié)點19 0 節(jié)點20 0 節(jié)點21 0 節(jié)點22 0 節(jié)點23 0 節(jié)點24 0 。這個方法的原理比較簡單，要求的數字計算機的內存量也比較小，適應當時的電子數字計算機制作水平和電力系統(tǒng)理論水平，于是電力系統(tǒng)計算人員轉向以阻抗矩陣為主的逐次代入法（以下簡稱阻抗法）。為了克服阻抗法在內存和速度方面的缺點，后來發(fā)展了以阻抗矩陣為基礎的分塊阻抗法。牛頓法的特點是將非線性方程線性化。由于電力系統(tǒng)規(guī)模的不斷擴大，對計算速度的要求不斷提高，計算機的并行計算技術也將在潮流計算中得到廣泛的應用，成為重要的研究領域。在牛頓法的基礎上，根據電力系統(tǒng)的特點，抓住主要矛盾，對純數學的牛頓法進行了改造，得到了PQ分解法。但是，阻抗法的主要缺點就是占用計算機的內存很大，每次迭代的計算量很大。這種情況促使電力系統(tǒng)的研究人員不斷尋求新的更可靠的計算方法。00 0 節(jié)點7。節(jié)點3。} for(k=(n2)。j{a[i][j]=aik*a[k][j]/a[k][k]。a[k][j]=al。j{if(pivot{pivot=fabs(a[i][j])。pivcol=new int *[n]。=0。(8)。}Complex Complex::getinverse(Complex c1)//取倒數 { Complex Node。} Complex Complex::ComDivRea(Complex c1,double r1)//復數除數 { Complex Node。}Complex Complex::productComplex(Complex c1,Complex c2)//復數乘法 {Complex Node。}Complex Complex::CaddC(Complex c1,Complex c2)//復數加法 {Complex Node。//將復數復零static Complex Rec2Polar(Complex c1)。//求兩個復數商static Complex ComDivRea(Complex c1,double r2)。}include using namespace std。coutt(bus[i].)。}coutBus::JisuanNodeScal(X,bus,N)。cc=Complex::productComplex(bb , Complex::getconj(line[i].Y))。i{int statemp,endtemp。//whether converbence看迭代是否結束if(icon==1){coutBus::JisuanJacAug(JacAug,X,bus,N)。x=new double[NN]。for(i=0。//輸出支路參數Bus *bus=new Bus[N]。//i作為整個程序的循環(huán)變量int N=Bus::ScanfBusNo()。i{int statemp,endtemp。i{bus[i]. = bus[i]. + bus[i].。line[i].stoe=Complex::productComplex(bus[statemp1].V,dd)。=0。gauss::gauss_slove(JacAug,x,NN)。LOOP:Bus::JisuanNodeI(X,bus,N)。//計算節(jié)點導納矩陣Bus::PrintfNodeDnz(X,N)。ibus[i].=0。//輸入支路個數if((Lamp。End2實例驗證 【例題】設有一系統(tǒng)網絡結線見圖1，各支路阻抗和各節(jié)點功率均已以標幺值標示于圖1中，其中節(jié)點2連接的是發(fā)電廠，設節(jié)點1電壓保持U1＝，試計算其中的潮流分布，請輸入節(jié)點數：n=5 請輸入支路數：m=7 請輸入平衡母線的節(jié)點號：ph=l 請輸入支路信息：BI=[ l 2 0．02+0．06i O l 1；1 3 0．08+0．24i 0 1 2；2 3 0．06+0．18i 0 l 3： 2 4 0．06+0．18i O l 4： 2 5 0．04+0．12i 0 l 5： 3 4 0．01+0．03i 0 l 6： 4 5 0．08+0．24i O 1 7] 請輸入節(jié)點信息：B2=[ 0 0 1．06 3；0．2+0．20i 0 1 1；一O．45一O．15i 0 l l；一0．40．05i 0 l 1；一0．6—0．1i 0 1 l] 請輸入節(jié)點號及對地阻抗： A=[l 0；2 0；3 0；4 0；5 O ] 請輸入修正值：ip=0．000 0l參考文獻[1]陳珩．電力系統(tǒng)穩(wěn)定分析[M]．北京：中國電力出版社，2002：139—187．[2]鄭阿奇．MATLAB實用教程[M]．北京：電子工業(yè)出版社，2005：1243．[3] 束洪春，孫士云，等．云電送粵交商流混聯系統(tǒng)全過 程動態(tài)電壓研究[J】．中國電力，2008，4l(10)：l4． SHU Hong—ch吼，SUN Shiyun，et a1．Research on fun prc39。end J(m,q)=X3。X3= C(i)+G(i,j)*e(i)B(i,j)*f(i)。J(m,q)=X4。q=2*j1。X2=G(i,j)*e(i)+B(i,j)*f(i)。case 3 P=2*i1。[x,y]=size(DX)。endfunction [De,Df]=hxf(J,Df,ph,n,no)%該子函數是為求取De Df DX=JDf。Df(p)=P(i)P1。type resultm function [C,D,Df]=xxf(G,B,e,f,P,Q,n,B2,ph,V,no)%該子程序是用來求取Df for i=1:nIfi=phC(i)=0。V=abs(v)。whilemax(abs(De))ipamp。V(i)=B2(i3)。e=zeros（1，n）；Y(p,q)=Y(p,q)。牛頓拉夫遜法是求解非線性方程的迭代過程，其計算公式為DF=JDX，式中J為所求函數的雅可比矩陣；DX為需要求的修正值；DF為不平衡的列向量。U(i)=U(i)+oU(i)。M(i,i)=M(i,i)U(i)*U(j)*(G(i,j)*cos(a(i)a(j))+B(i,j)*sin(a(i)a(j)))。oQ(i)=oQ(i)+Q(i)。N(i,j)=0。Q=imag(S)。endY% Y為導納矩陣G=real(Y)。y(1,5)=1/(+)。為了便于用迭代法解方程組，需要將上述功率方程改寫成功率平衡方程，并對功率平衡方程求偏導，得出對應的雅可比矩陣，給未知節(jié)點賦電壓初值，將初值帶入功率平衡方程，得到功率不平衡量，這樣由功率不平衡量、雅可比矩陣、節(jié)點電壓不平衡量（未知的）構成了誤差方程，解誤差方程，得到節(jié)點電壓不平衡量，節(jié)點電壓加上節(jié)點電壓不平衡量構成節(jié)點電壓新的初值，將新的初值帶入原來的功率平衡方程，并重新形成雅可比矩陣，然后計算新的電壓不平衡量，這樣不斷迭代，不斷修正，一般迭代三到五次就能收斂。潮流計算的計算機算法包含高斯—賽德爾迭代法、牛頓拉夫遜法和P—Q分解法等，其中牛拉法計算原理較簡單、計算過程也不復雜，而且由于人們引入泰勒級數和非線性代數方程等在算法里從而進一步提高了算法的收斂性和計算速度。進行潮流計算前需要在文本方式下潮流計算信息窗口中定義作業(yè)，如圖2所示。具有電壓和功率的自動控制功能； 在解決電力系統(tǒng)潮流計算問題時，是以導納矩陣為基礎的，因此，只要我們能在迭代過程中盡可能保持方程式系數矩陣的稀疏性，就可以大大提高牛頓法潮流程序的效率。在用數字解算計算機解電力系統(tǒng)潮流問題的開始階段，普遍采取以節(jié)點導納矩陣為基礎的逐次代入法?？煽啃院徒洕?。同時，為了提高供電的可靠性以及資源利用的綜合經濟性，又把許多分散的各種形式的發(fā)電站，通過送電線路和變電所聯系起來。按電壓的不同表示方法，牛頓拉夫遜潮流計算分為直角坐標形式和極坐標形式兩種。第一篇：基于psasp的電力系統(tǒng)潮流計算基于psasp的電力系統(tǒng)潮流計算摘要：電力系統(tǒng)潮流計算是研究電力系統(tǒng)穩(wěn)態(tài)運行情況的一種計算，它根據給定的運行條件及系統(tǒng)接線情況確定整個電力系統(tǒng)各部分的運行狀態(tài)。實際電力系統(tǒng)的潮流計算主要采用牛頓拉夫遜法。隨著工農業(yè)生產和城市的發(fā)展，電能的需要量迅速增加，而熱能資源和水能資源豐富的地區(qū)又往往遠離用電比較集中的城市和工礦區(qū)，為了解決這個矛盾，就需要在動力資源豐富的地區(qū)建立大型發(fā)電站，然后將電能遠距離輸送給電力用戶。在電力系統(tǒng)規(guī)劃的設計和現有電力系統(tǒng)運行方式的研究中，都需要利用潮流計算來定量的分析比較供電方案或運行方式的合理性。實際電力系統(tǒng)的潮流技術那主要采用牛頓—拉夫遜法。這是數學中解決非線性方程式的典型方法，有較好的收斂性。提供PQ分解、牛頓、最佳乘了等五種計算方法，以保證良好的收斂性； 并且進行區(qū)域設置。參考文獻：[1]孟祥萍,[M].北京:高等教育出版社,2013,11.[2][M].北京:人民郵電出版社,2012,12.[3][J].華中科技大學,2013,11.[4][J].,6.[5]薛振宇,房大中,[J