【文章內(nèi)容簡(jiǎn)介】
G啪gdong province【J】.Electric Power,2008,4l(10): l4.[4] 朱新立,湯涌,等.大電網(wǎng)安全分析的全過(guò)程動(dòng)態(tài)仿 真技術(shù)[J】.電網(wǎng)技術(shù),2008,32(22):23—28. SONG Xin—Ii,TANG Yof唱,et a1. Full dyn鋤ic simulation for the stabilhy a眥lysis of large power system【J】.Power System融IlrIolo影,2008,32(22): 23.28.[5]Roytelm鋤I,Shallidehpour S M.A prehsivc long teml dynaIIlic simulation for powcr system recoVery【J】. IEEE Transactions 0n Power Systems,1994,9(3). [6] 石雩梅,汪志宏,等.發(fā)電機(jī)勵(lì)磁系統(tǒng)數(shù)學(xué)模型及參 數(shù)對(duì)電網(wǎng)動(dòng)態(tài)穩(wěn)定性分析結(jié)果影響的研究[J】.繼電 器,2007,35(21):2227.SHI Xue.mei,WANG Zllihon舀et a1.Iksearch on the innuence of g鋤e翰to璐baScd ∞de詛iled excitation system models柚d parameterS t0 power鏟id dyn鋤ic stabil時(shí)【J】.Relay,2007,35(2 1):2227.[7] 方思立,朱方.快速勵(lì)磁系統(tǒng)對(duì)系統(tǒng)穩(wěn)定的影響[J】.中 國(guó)電機(jī)工程學(xué)報(bào),1986,6(1):20.28.FANG Si.1i,ZHU Fang.The effbct of f弧t.respon∞excitation system on the stability of power netwofk【J】. Proceedings ofthe CSEE,1986,6(1):2028.[8] 劉?。娏ο到y(tǒng)穩(wěn)定性及發(fā)電機(jī)勵(lì)磁控制[M】.北京: 中國(guó)電力出版社,2007.LIU Qu.Power system S詛bility鋤d generator excitation control【M】.BeUing:ChiIla Electric Powef Press,2007. [9] Dallachy J L,Anderson T.EXperience with rcplacing ro詛ting exciters wim static exciters【J】.1k InStitution of Electrical Engineers,1 996.[10] 陳利芳,陳天祿.淺談自并勵(lì)勵(lì)磁系統(tǒng)在大容量機(jī)組 中的應(yīng)用【J】.繼電器,2007,35(1):8l培4. CHEN Lif抽島CHEN Tian—lIL Application of 辯l仁exci組tion mode in large capacity髫memtor unit【J】. ReIay39。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 system【J】.Electric Powef,1994,27(12): 61.63.[12]梁小冰,黃方能.利用EMTDC進(jìn)行長(zhǎng)持續(xù)時(shí)間過(guò)程 的仿真研究【J】.電網(wǎng)技術(shù),2002,26(9):55.57. LIANG Xiaobing,HUANG Fan爭(zhēng)眥ng.How to cany out simulalion of long dul‘a(chǎn)tion processes by use of EMTDC【J】.Power System 11echnology,2002,26(9): 5557.[13]王卉,陳楷,彭哲,等.?dāng)?shù)字仿真技術(shù)在電力系統(tǒng)中 的應(yīng)用及常用的幾種數(shù)字仿真工具【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):7175.[14]IEEE Power Engmeering Socie哆.IEEE std 421.5.2005 IEEE玎ccOmmended practice for excitation system models for power system stabiI時(shí)studies【s】.第四篇:電力系統(tǒng)潮流計(jì)算程序電力系統(tǒng)潮流計(jì)算c語(yǔ)言程序,兩行,大家可以看看,仔細(xì)研究,然后在這個(gè)基礎(chǔ)上修改。謝謝include “” include include include include“” include“” include“” using namespace std。int _tmain(int argc, _TCHAR* argv[]){int i。//i作為整個(gè)程序的循環(huán)變量int N=Bus::ScanfBusNo()。//輸入節(jié)點(diǎn)個(gè)數(shù)int L=Line::ScanflineNo()。//輸入支路個(gè)數(shù)if((Lamp。amp。N)==0){return 0。} //如果找不到兩個(gè)文件中的任意一個(gè),退出Line *line=new Line[L]。//動(dòng)態(tài)分配支路結(jié)構(gòu)體Line::ScanfLineData(line)。//輸入支路參數(shù)Line::PrintfLineData(line,L)。//輸出支路參數(shù)Bus *bus=new Bus[N]。//動(dòng)態(tài)分配結(jié)點(diǎn)結(jié)構(gòu)體for(int i=0。ibus[i].=0。bus[i].=0。}Bus::ScanfBusData(bus)。//輸入節(jié)點(diǎn)參數(shù)Bus::PrintfBusData(bus,N)。//輸出結(jié)點(diǎn)參數(shù)Complex **X。X=new Complex *[N]。for(i=0。iBus::JisuanNodeDnz(X,line,bus,L,N)。//計(jì)算節(jié)點(diǎn)導(dǎo)納矩陣Bus::PrintfNodeDnz(X,N)。//輸出節(jié)點(diǎn)導(dǎo)納矩陣int NN=(N1)*2。double **JacAug。JacAug=new double *[NN]。for(i=0。idouble *x。x=new double[NN]。int count=1。LOOP:Bus::JisuanNodeI(X,bus,N)。//計(jì)算節(jié)點(diǎn)注入電流Bus::JisuanNodeScal(X,bus,N)。//計(jì)算節(jié)點(diǎn)功率Bus::JisuanNodeScal(X,bus,N)。//計(jì)算節(jié)點(diǎn)功率Bus::JisuanNodeSdelta(bus,N)。//計(jì)算節(jié)點(diǎn)功率差值Bus::PrintfNodeScal(X,bus,N)。//輸出節(jié)點(diǎn)功率差值int icon=wehcon1(bus,N)。//whether converbence看迭代是否結(jié)束if(icon==1){coutBus::JisuanJacAug(JacAug,X,bus,N)。//計(jì)算雅可比增廣矩陣 // Bus::PrintfJacAug(JacAug,N)。gauss::gauss_slove(JacAug,x,NN)。//解方程組求出電壓差值Bus::ReviseNodeV(bus,x,N)。//修正節(jié)點(diǎn)電壓// Bus::PrintfNodeV(bus,N)。count++。goto LOOP。}else{for(i=0。i{int statemp,endtemp。Complex aa,bb,cc,dd,B。=0。=line[i].B。statemp=line[i].start。endtemp=line[i].end。aa=Complex::productComplex(Complex::getconj(bus[statemp1].V), B)。bb=Complex::subComplex(Complex::getconj(bus[statemp1].V), Complex::getconj(bus[endtemp1].V))。cc=Complex::productComplex(bb , Complex::getconj(line[i].Y))。dd=Complex::CaddC(aa,cc)。line[i].stoe=Complex::productComplex(bus[statemp1].V,dd)。aa=Complex::productComplex(Complex::getconj(bus[endtemp1].V), B)。bb=Complex::subComplex(Complex::getconj(bus[endtemp1].V), Complex::getconj(bus[statemp1].V))。cc=Complex::productComplex(bb , Complex::getconj(line[i].Y))。dd=Complex::CaddC(aa,cc)。line[i].etos=Complex::productComplex(bus[endtemp1].V,dd)。}coutBus::JisuanNodeScal(X,bus,N)。//計(jì)算節(jié)點(diǎn)功率for(i=0。i{bus[i]. = bus[i]. + bus[i].。//發(fā)電機(jī)功率=注入功率+負(fù)荷功率bus[i].= bus[i].+ bus[i].。bus[i].V=Complex::Rec2Polar(bus[i].V)。}coutfor(i=0。i