【正文】 em models for power system stabiI時(shí)studies【s】．第四篇：電力系統(tǒng)潮流計(jì)算程序電力系統(tǒng)潮流計(jì)算c語言程序，兩行，大家可以看看，仔細(xì)研究，然后在這個(gè)基礎(chǔ)上修改。//輸入節(jié)點(diǎn)參數(shù)Bus::PrintfBusData(bus,N)。//計(jì)算節(jié)點(diǎn)功率Bus::JisuanNodeSdelta(bus,N)。endtemp=line[i].end。}coutfor(i=0。amp。//輸出節(jié)點(diǎn)導(dǎo)納矩陣int NN=(N1)*2。//解方程組求出電壓差值Bus::ReviseNodeV(bus,x,N)。aa=Complex::productComplex(Complex::getconj(bus[endtemp1].V), B)。statemp=line[i].start。//取倒數(shù)static double getComplexReal(Complex c1)。return Node。return Node。(6)。pivrow[k]=k。a[i][pivcol[k][1]]=a[i][k]。x[pivcol[k][1]]=x[pivcol[k][0]]。0 20節(jié) 點(diǎn)30 。這個(gè)方法把一個(gè)大系統(tǒng)分割為幾個(gè)小的地區(qū)系統(tǒng)，在計(jì)算機(jī)內(nèi)只需存儲(chǔ)各個(gè)地區(qū)系統(tǒng)的阻抗矩陣及它們之間的聯(lián)絡(luò)線的阻抗，這樣不僅大幅度的節(jié)省了內(nèi)存容量，同時(shí)也提高了節(jié)省速度。但是，到目前為止這些新的模型和算法還不能取代牛頓法和PQ分解法的地位。阻抗法改善了電力系統(tǒng)潮流計(jì)算問題的收斂性，解決了導(dǎo)納法無法解決的一些系統(tǒng)的潮流計(jì)算，在當(dāng)時(shí)獲得了廣泛的應(yīng)用，曾為我國(guó)電力系統(tǒng)設(shè)計(jì)、運(yùn)行和研究作出了很大的貢獻(xiàn)。0節(jié)點(diǎn)28 。 }x[i]=(a[i][n]sum)/a[i][i]。a[pivrow[k]][j]=a[k][j]。pivrow=new int[n]。}voidComplex::PrintfComplex(Complex c1)//按直角坐標(biāo)輸出 { if(==0){ (6)。return Node。return Node。//求兩個(gè)復(fù)數(shù)積static Complex divideComplex(Complex c1,Complex c2)。coutt(bus[i].)。bb=Complex::subComplex(Complex::getconj(bus[statemp1].V), Complex::getconj(bus[endtemp1].V))。//輸出節(jié)點(diǎn)功率差值int icon=wehcon1(bus,N)。X=new Complex *[N]。int _tmain(int argc, _TCHAR* argv[]){int i。//計(jì)算節(jié)點(diǎn)功率for(i=0。Complex aa,bb,cc,dd,B。int count=1。//動(dòng)態(tài)分配結(jié)點(diǎn)結(jié)構(gòu)體for(int i=0。End J(:,no)=[]。q=q+1。X2=C(i)+G(i,j)*e(i)+B(i,j)*f(i)。switch B2(i4)X4=0。If ph=n Df(noJ=[]。printcut(result1,S,b,c,result2)。t=0。Y=zeros（n）；Y(p,q)=Y(p,q)1./(B1(i3)*B1(i5)。endfor i=1:4a(i)=a(i)+oa(i)。endoP(i)=oP(i)+P(i)。P=real(S)。y(1,4)=1/(+)。參考文獻(xiàn)：[1]孟祥萍,[M].北京:高等教育出版社,2013,11.[2][M].北京:人民郵電出版社,2012,12.[3][J].華中科技大學(xué),2013,11.[4][J].,6.[5]薛振宇,房大中,[J].天津大學(xué)學(xué)報(bào),2011.[6]王錫凡,方萬良,[M].北京:科學(xué)出版社,2013.[7]周曉峰,[J]..[8]王曉茹,[J]..[9]楊清,[J]..[10][M].北京:呼國(guó)電力科學(xué)研究院,2002.第二篇：電力系統(tǒng)潮流計(jì)算南 京 理 工 大 學(xué)《電力系統(tǒng)穩(wěn)態(tài)分析》課程報(bào)告姓名XX學(xué) 號(hào)： 5*** 自動(dòng)化學(xué)院 電氣工程基于牛頓拉夫遜法的潮流計(jì)算例題編程報(bào)學(xué)院(系): 專業(yè): 題目: 任課教師 碩士導(dǎo)師 告楊偉 XX2015年6月10號(hào)基于牛頓拉夫遜法的潮流計(jì)算例題編程報(bào)告摘要：電力系統(tǒng)潮流計(jì)算的目的在于：確定電力系統(tǒng)的運(yùn)行方式、檢查系統(tǒng)中各元件是否過壓或者過載、為電力系統(tǒng)繼電保護(hù)的整定提供依據(jù)、為電力系統(tǒng)的穩(wěn)定計(jì)算提供初值、為電力系統(tǒng)規(guī)劃和經(jīng)濟(jì)運(yùn)行提供分析的基礎(chǔ)。提供PQ分解、牛頓、最佳乘了等五種計(jì)算方法，以保證良好的收斂性； 實(shí)際電力系統(tǒng)的潮流技術(shù)那主要采用牛頓—拉夫遜法。隨著工農(nóng)業(yè)生產(chǎn)和城市的發(fā)展，電能的需要量迅速增加，而熱能資源和水能資源豐富的地區(qū)又往往遠(yuǎn)離用電比較集中的城市和工礦區(qū)，為了解決這個(gè)矛盾，就需要在動(dòng)力資源豐富的地區(qū)建立大型發(fā)電站，然后將電能遠(yuǎn)距離輸送給電力用戶。第一篇：基于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)接線情況確定整個(gè)電力系統(tǒng)各部分的運(yùn)行狀態(tài)。同時(shí)，為了提高供電的可靠性以及資源利用的綜合經(jīng)濟(jì)性，又把許多分散的各種形式的發(fā)電站，通過送電線路和變電所聯(lián)系起來。在用數(shù)字解算計(jì)算機(jī)解電力系統(tǒng)潮流問題的開始階段，普遍采取以節(jié)點(diǎn)導(dǎo)納矩陣為基礎(chǔ)的逐次代入法。具有電壓和功率的自動(dòng)控制功能； 潮流計(jì)算的計(jì)算機(jī)算法包含高斯—賽德爾迭代法、牛頓拉夫遜法和P—Q分解法等，其中牛拉法計(jì)算原理較簡(jiǎn)單、計(jì)算過程也不復(fù)雜，而且由于人們引入泰勒級(jí)數(shù)和非線性代數(shù)方程等在算法里從而進(jìn)一步提高了算法的收斂性和計(jì)算速度。y(1,5)=1/(+)。Q=imag(S)。oQ(i)=oQ(i)+Q(i)。U(i)=U(i)+oU(i)。e=zeros（1，n）；Y(p,q)=Y(p,q)。whilemax(abs(De))ipamp。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。endfunction [De,Df]=hxf(J,Df,ph,n,no)%該子函數(shù)是為求取De Df DX=JDf。case 3 P=2*i1。q=2*j1。X3= C(i)+G(i,j)*e(i)B(i,j)*f(i)。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ì)算其中的潮流分布，請(qǐng)輸入節(jié)點(diǎn)數(shù)：n=5 請(qǐng)輸入支路數(shù)：m=7 請(qǐng)輸入平衡母線的節(jié)點(diǎn)號(hào)：ph=l 請(qǐng)輸入支路信息：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] 請(qǐng)輸入節(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] 請(qǐng)輸入節(jié)點(diǎn)號(hào)及對(duì)地阻抗： A=[l 0；2 0；3 0；4 0；5 O ] 請(qǐng)輸入修正值：ip=0．000 0l參考文獻(xiàn)[1]陳珩．電力系統(tǒng)穩(wěn)定分析[M]．北京：中國(guó)電力出版社，2002：139—187．[2]鄭阿奇．MATLAB實(shí)用教程[M]．北京：電子工業(yè)出版社，2005：1243．[3] 束洪春，孫士云，等．云電送粵交商流混聯(lián)系統(tǒng)全過 程動(dòng)態(tài)電壓研究[J】．中國(guó)電力，2008，4l(10)：l4． SHU Hong—ch吼，SUN Shiyun，et a1．Research on fun prc39。ibus[i].=0。LOOP:Bus::JisuanNodeI(X,bus,N)。=0。i{bus[i]. = bus[i]. + bus[i].。//i作為整個(gè)程序的循環(huán)變量int N=Bus::ScanfBusNo()。for(i=0。//whether converbence看迭代是否結(jié)束if(icon==1){coutBus::JisuanJacAug(JacAug,X,bus,N)。cc=Complex::productComplex(bb , Complex::getconj(line[i].Y))。coutt(bus[i].)。//求兩個(gè)復(fù)數(shù)商static Complex ComDivRea(Complex c1,double r2)。}Complex Complex::CaddC(Complex c1,Complex c2)//復(fù)數(shù)加法 {Complex Node。} Complex Complex::ComDivRea(Complex c1,double r1)//復(fù)數(shù)除數(shù) { Complex Node。(8)。pivcol=new int *[n]。a[k][j]=al。} for(k=(n2)。00 0 節(jié)點(diǎn)7。但是，阻抗法的主要缺點(diǎn)就是占用計(jì)算機(jī)的內(nèi)存很大，每次迭代的計(jì)算量很大。由于電力系統(tǒng)規(guī)模的不斷擴(kuò)大，對(duì)計(jì)算速度的要求不斷提高，計(jì)算機(jī)的并行計(jì)算技術(shù)也將在潮流計(jì)算中得到廣泛的應(yīng)用，成為重要的研究領(lǐng)域。為了克服阻抗法在內(nèi)存和速度方面的缺點(diǎn)，后來發(fā)展了以阻抗矩陣為基礎(chǔ)的分塊阻抗法。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 。k){al=x[pivcol[k][1]]。i{al=a[i][pivcol[k][1]]。ipivot= fabs(a[k][k])。cout”。=(r1)。=+。//求一個(gè)復(fù)數(shù)共軛static Complex getinverse(Complex c1)。i{int statemp,endtemp。line[i].stoe=Complex::productComplex(bus[statemp1].V,dd)。gauss::gauss_slove(JacAug,x,NN)。//計(jì)算節(jié)點(diǎn)導(dǎo)納矩陣Bus::PrintfNodeDnz(X,N)。//輸入支路個(gè)數(shù)if((Lamp。bus[i].V=Complex::Rec2Polar(bus[i].V)。statemp=line[i].start。//計(jì)算節(jié)點(diǎn)功率Bus::JisuanNodeScal(X,bus,N)。}Bus::ScanfBusData(bus)。2007，35(1)：8184．[11] 方思立，劉增煌，孟慶和．大型汽輪發(fā)電機(jī)自并勵(lì)勵(lì) 磁系統(tǒng)的應(yīng)用條件【J]．中國(guó)電力，1994，27(12)：61．63． FANG Si．Ii，LIU Zenghu鋤g，MENG Qin爭(zhēng)hc．m application conditions of large turbine generator selfexcitation syst