【正文】 【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。//計(jì)算節(jié)點(diǎn)功率Bus::JisuanNodeScal(X,bus,N)。//修正節(jié)點(diǎn)電壓// 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。//輸入支路個數(shù)if((Lamp。ibus[i].=0。//計(jì)算節(jié)點(diǎn)導(dǎo)納矩陣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。//求一個復(fù)數(shù)共軛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é)點(diǎn)4。00 節(jié)點(diǎn)8 節(jié)點(diǎn)9 0 節(jié)點(diǎn)10 0 節(jié)點(diǎn)11 節(jié)點(diǎn)12 0 節(jié)點(diǎn)13 節(jié)點(diǎn)14 0 節(jié)點(diǎn)15 0 節(jié)點(diǎn)16 0 節(jié)點(diǎn)17 0 節(jié)點(diǎn)18 0 節(jié)點(diǎn)19 0 節(jié)點(diǎn)20 0 節(jié)點(diǎn)21 0 節(jié)點(diǎn)22 0 節(jié)點(diǎn)23 0 節(jié)點(diǎn)24 0 。這個方法的原理比較簡單，要求的數(shù)字計(jì)算機(jī)的內(nèi)存量也比較小，適應(yīng)當(dāng)時(shí)的電子數(shù)字計(jì)算機(jī)制作水平和電力系統(tǒng)理論水平，于是電力系統(tǒng)計(jì)算人員轉(zhuǎn)向以阻抗矩陣為主的逐次代入法（以下簡稱阻抗法）。為了克服阻抗法在內(nèi)存和速度方面的缺點(diǎn)，后來發(fā)展了以阻抗矩陣為基礎(chǔ)的分塊阻抗法。牛頓法的特點(diǎn)是將非線性方程線性化。由于電力系統(tǒng)規(guī)模的不斷擴(kuò)大，對計(jì)算速度的要求不斷提高，計(jì)算機(jī)的并行計(jì)算技術(shù)也將在潮流計(jì)算中得到廣泛的應(yīng)用，成為重要的研究領(lǐng)域。在牛頓法的基礎(chǔ)上，根據(jù)電力系統(tǒng)的特點(diǎn)，抓住主要矛盾，對純數(shù)學(xué)的牛頓法進(jìn)行了改造，得到了PQ分解法。但是，阻抗法的主要缺點(diǎn)就是占用計(jì)算機(jī)的內(nèi)存很大，每次迭代的計(jì)算量很大。這種情況促使電力系統(tǒng)的研究人員不斷尋求新的更可靠的計(jì)算方法。00 0 節(jié)點(diǎn)7。節(jié)點(diǎn)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)//取倒數(shù) { Complex Node。} Complex Complex::ComDivRea(Complex c1,double r1)//復(fù)數(shù)除數(shù) { Complex Node。}Complex Complex::productComplex(Complex c1,Complex c2)//復(fù)數(shù)乘法 {Complex Node。}Complex Complex::CaddC(Complex c1,Complex c2)//復(fù)數(shù)加法 {Complex Node。//將復(fù)數(shù)復(fù)零static Complex Rec2Polar(Complex c1)。//求兩個復(fù)數(shù)商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看迭代是否結(jié)束if(icon==1){coutBus::JisuanJacAug(JacAug,X,bus,N)。x=new double[NN]。for(i=0。//輸出支路參數(shù)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)。//計(jì)算節(jié)點(diǎn)導(dǎo)納矩陣Bus::PrintfNodeDnz(X,N)。ibus[i].=0。//輸入支路個數(shù)if((Lamp。End2實(shí)例驗(yàn)證 【例題】設(shè)有一系統(tǒng)網(wǎng)絡(luò)結(jié)線見圖1，各支路阻抗和各節(jié)點(diǎn)功率均已以標(biāo)幺值標(biāo)示于圖1中，其中節(jié)點(diǎn)2連接的是發(fā)電廠，設(shè)節(jié)點(diǎn)1電壓保持U1＝，試計(jì)算其中的潮流分布，請輸入節(jié)點(diǎn)數(shù)：n=5 請輸入支路數(shù)：m=7 請輸入平衡母線的節(jié)點(diǎn)號：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é)點(diǎn)信息：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é)點(diǎn)號及對地阻抗： A=[l 0；2 0；3 0；4 0；5 O ] 請輸入修正值：ip=0．000 0l參考文獻(xiàn)[1]陳珩．電力系統(tǒng)穩(wěn)定分析[M]．北京：中國電力出版社，2002：139—187．[2]鄭阿奇．MATLAB實(shí)用教程[M]．北京：電子工業(yè)出版社，2005：1243．[3] 束洪春，孫士云，等．云電送粵交商流混聯(lián)系統(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)%該子函數(shù)是為求取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)。牛頓拉夫遜法是求解非線性方程的迭代過程，其計(jì)算公式為DF=JDX，式中J為所求函數(shù)的雅可比矩陣；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為導(dǎo)納矩陣G=real(Y)。y(1,5)=1/(+)。為了便于用迭代法解方程組，需要將上述功率方程改寫成功率平衡方程，并對功率平衡方程求偏導(dǎo)，得出對應(yīng)的雅可比矩陣，給未知節(jié)點(diǎn)賦電壓初值，將初值帶入功率平衡方程，得到功率不平衡量，這樣由功率不平衡量、雅可比矩陣、節(jié)點(diǎn)電壓不平衡量（未知的）構(gòu)成了誤差方程，解誤差方程，得到節(jié)點(diǎn)電壓不平衡量，節(jié)點(diǎn)電壓加上節(jié)點(diǎn)電壓不平衡量構(gòu)成節(jié)點(diǎn)電壓新的初值，將新的初值帶入原來的功率平衡方程，并重新形成雅可比矩陣，然后計(jì)算新的電壓不平衡量，這樣不斷迭代，不斷修正，一般迭代三到五次就能收斂。潮流計(jì)算的計(jì)算機(jī)算法包含高斯—賽德爾迭代法、牛頓拉夫遜法和P—Q分解法等，其中牛拉法計(jì)算原理較簡單、計(jì)算過程也不復(fù)雜，而且由于人們引入泰勒級數(shù)和非線性代數(shù)方程等在算法里從而進(jìn)一步提高了算法的收斂性和計(jì)算速度。進(jìn)行潮流計(jì)算前需要在文本方式下潮流計(jì)算信息窗口中定義作業(yè)，如圖2所示。具有電壓和功率的自動控制功能； 在解決電力系統(tǒng)潮流計(jì)算問題時(shí)，是以導(dǎo)納矩陣為基礎(chǔ)的，因此，只要我們能在迭代過程中盡可能保持方程式系數(shù)矩陣的稀疏性，就可以大大提高牛頓法潮流程序的效率。在用數(shù)字解算計(jì)算機(jī)解電力系統(tǒng)潮流問題的開始階段，普遍采取以節(jié)點(diǎn)導(dǎo)納矩陣為基礎(chǔ)的逐次代入法?？煽啃院徒?jīng)濟(jì)性。同時(shí)，為了提高供電的可靠性以及資源利用的綜合經(jīng)濟(jì)性，又把許多分散的各種形式的發(fā)電站，通過送電線路和變電所聯(lián)系起來。按電壓的不同表示方法，牛頓拉夫遜潮流計(jì)算分為直角坐標(biāo)形式和極坐標(biāo)形式兩種。第一篇：基于psasp的電力系統(tǒng)潮流計(jì)算基于psasp的電力系統(tǒng)潮流計(jì)算摘要：電力系統(tǒng)潮流計(jì)算是研究電力系統(tǒng)穩(wěn)態(tài)運(yùn)行情況的一種計(jì)算，它根據(jù)給定的運(yùn)行條件及系統(tǒng)接線情況確定整個電力系統(tǒng)各部分的運(yùn)行狀態(tài)。實(shí)際電力系統(tǒng)的潮流計(jì)算主要采用牛頓拉夫遜法。隨著工農(nóng)業(yè)生產(chǎn)和城市的發(fā)展，電能的需要量迅速增加，而熱能資源和水能資源豐富的地區(qū)又往往遠(yuǎn)離用電比較集中的城市和工礦區(qū)，為了解決這個矛盾，就需要在動力資源豐富的地區(qū)建立大型發(fā)電站，然后將電能遠(yuǎn)距離輸送給電力用戶。在電力系統(tǒng)規(guī)劃的設(shè)計(jì)和現(xiàn)有電力系統(tǒng)運(yùn)行方式的研究中，都需要利用潮流計(jì)算來定量的分析比較供電方案或運(yùn)行方式的合理性。實(shí)際電力系統(tǒng)的潮流技術(shù)那主要采用牛頓—拉夫遜法。這是數(shù)學(xué)中解決非線性方程式的典型方法，有較好的收斂性。提供PQ分解、牛頓、最佳乘了等五種計(jì)算方法，以保證良好的收斂性； 并且進(jìn)行區(qū)域設(shè)置。參考文獻(xiàn)：[1]孟祥萍,[M].北京:高等教育出版社,2013,11.[2][M].北京:人民郵電出版社,2012,12.[3][J].華中科技大學(xué),2013,11.[4][J].,6.[5]薛振宇,房大中,[J