【正文】 tabil時(shí)of hybrid AC／DC poWer tmnsmission System舶m Yu衄an proVince to 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。75．[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ǔ)上修改。int _tmain(int argc, _TCHAR* argv[]){int i。//輸入節(jié)點(diǎn)個(gè)數(shù)int L=Line::ScanflineNo()。amp。} //如果找不到兩個(gè)文件中的任意一個(gè)，退出Line *line=new Line[L]。//輸入支路參數(shù)Line::PrintfLineData(line,L)。//動(dòng)態(tài)分配結(jié)點(diǎn)結(jié)構(gòu)體for(int i=0。bus[i].=0。//輸入節(jié)點(diǎn)參數(shù)Bus::PrintfBusData(bus,N)。X=new Complex *[N]。iBus::JisuanNodeDnz(X,line,bus,L,N)。//輸出節(jié)點(diǎn)導(dǎo)納矩陣int NN=(N1)*2。JacAug=new double *[NN]。idouble *x。int count=1。//計(jì)算節(jié)點(diǎn)注入電流Bus::JisuanNodeScal(X,bus,N)。//計(jì)算節(jié)點(diǎn)功率Bus::JisuanNodeSdelta(bus,N)。//輸出節(jié)點(diǎn)功率差值int icon=wehcon1(bus,N)。//計(jì)算雅可比增廣矩陣 // Bus::PrintfJacAug(JacAug,N)。//解方程組求出電壓差值Bus::ReviseNodeV(bus,x,N)。count++。}else{for(i=0。Complex aa,bb,cc,dd,B。=line[i].B。endtemp=line[i].end。bb=Complex::subComplex(Complex::getconj(bus[statemp1].V), Complex::getconj(bus[endtemp1].V))。dd=Complex::CaddC(aa,cc)。aa=Complex::productComplex(Complex::getconj(bus[endtemp1].V), B)。cc=Complex::productComplex(bb , Complex::getconj(line[i].Y))。line[i].etos=Complex::productComplex(bus[endtemp1].V,dd)。//計(jì)算節(jié)點(diǎn)功率for(i=0。//發(fā)電機(jī)功率=注入功率+負(fù)荷功率bus[i].= bus[i].+ bus[i].。}coutfor(i=0。coutt(bus[i].)。cout}coutfor(i=0。statemp=line[i].start。coutComplex::PrintfComplex(Complex::ComDivRea(line[i].stoe,))。cout} }return 0。int _tmain(int argc, _TCHAR* argv[]){int i。//輸入節(jié)點(diǎn)個(gè)數(shù)int L=Line::ScanflineNo()。amp。} //如果找不到兩個(gè)文件中的任意一個(gè)，退出Line *line=new Line[L]。//輸入支路參數(shù)Line::PrintfLineData(line,L)。//動(dòng)態(tài)分配結(jié)點(diǎn)結(jié)構(gòu)體for(int i=0。bus[i].=0。//輸入節(jié)點(diǎn)參數(shù)Bus::PrintfBusData(bus,N)。X=new Complex *[N]。iBus::JisuanNodeDnz(X,line,bus,L,N)。//輸出節(jié)點(diǎn)導(dǎo)納矩陣int NN=(N1)*2。JacAug=new double *[NN]。idouble *x。int count=1。//計(jì)算節(jié)點(diǎn)注入電流Bus::JisuanNodeScal(X,bus,N)。//計(jì)算節(jié)點(diǎn)功率Bus::JisuanNodeSdelta(bus,N)。//輸出節(jié)點(diǎn)功率差值int icon=wehcon1(bus,N)。//計(jì)算雅可比增廣矩陣// Bus::PrintfJacAug(JacAug,N)。//解方程組求出電壓差值Bus::ReviseNodeV(bus,x,N)。count++。}else{for(i=0。Complex aa,bb,cc,dd,B。=line[i].B。endtemp=line[i].end。bb=Complex::subComplex(Complex::getconj(bus[statemp1].V), Complex::getconj(bus[endtemp1].V))。dd=Complex::CaddC(aa,cc)。aa=Complex::productComplex(Complex::getconj(bus[endtemp1].V), B)。cc=Complex::productComplex(bb , Complex::getconj(line[i].Y))。line[i].etos=Complex::productComplex(bus[endtemp1].V,dd)。//計(jì)算節(jié)點(diǎn)功率for(i=0。//發(fā)電機(jī)功率=注入功率+負(fù)荷功率bus[i].= bus[i].+ bus[i].。}coutfor(i=0。coutt(bus[i].)。cout}coutfor(i=0。statemp=line[i].start。coutComplex::PrintfComplex(Complex::ComDivRea(line[i].stoe,))。cout} }return 0。class Complex//定義復(fù)數(shù)類(lèi) { public: double real。intRecPolar。//求兩個(gè)復(fù)數(shù)和static Complex subComplex(Complex c1,Complex c2)。//求兩個(gè)復(fù)數(shù)積static Complex divideComplex(Complex c1,Complex c2)。//除數(shù)static Complex getconj(Complex c1)。//取倒數(shù)static double getComplexReal(Complex c1)。//求一個(gè)復(fù)數(shù)虛部static void PrintfComplex(Complex c1)。//顯示多個(gè)復(fù)數(shù)static void zeroComplex(Complex c1)。//取極坐標(biāo)Complex(){RecPolar=0。Complex Complex::Rec2Polar(Complex c1)//極坐標(biāo)表示 { Complex Node。=atan2(,)*180/。return Node。=+。return Node。=。return Node。=**。return Node。=(*+*)/(pow(,2)+pow(,2))。return Node。=(r1)。return Node。=。return Node。=1。Node=(Complex::divideComplex(Node,c1))。}double Complex::getComplexReal(Complex c1)//取實(shí)部 {return 。}voidComplex::PrintfComplex(Complex c1)//按直角坐標(biāo)輸出 { if(==0){ (6)。(ios::right)。(6)。(ios::left)。} else {coutComplex::zeroComplex(Complex c1)//清零 { =0。}class gauss { public: static void gauss_slove(double **a,double *x,int NN)。}。double eps,pivot,sum,aik,al。pivrow=new int[n]。for(i=0。pivrow[k]=k。pivcol[k][1]=k。i{for(j=k。pivrow[k]=i。//列}}}if(pivot{coutgetchar()。}if(pivrow[k]!=k)//行變換{for(j=k。a[pivrow[k]][j]=a[k][j]。}}if(pivcol[k][1]!=k)//列變換{for(i=0。a[i][pivcol[k][1]]=a[i][k]。}}if(k!=(n1))//將矩陣化為上三角形式{for(i=(k+1)。for(j=k。}}} } x[n1]=a[n1][n]/a[n1][n1]。i=0。for(j=(i+1)。 }x[i]=(a[i][n]sum)/a[i][i]。k=0。x[pivcol[k][1]]=x[pivcol[k][0]]。}coutcout}====節(jié)點(diǎn)功率計(jì)算值==== icon=1進(jìn)行第2次迭代 節(jié)點(diǎn)電壓修正量================= ====節(jié)點(diǎn)功率計(jì)算值==== icon=1進(jìn)行第3次迭代 節(jié)點(diǎn)電壓修正量================= ====節(jié)點(diǎn)功率計(jì)算值==== icon=0,迭代結(jié)束。節(jié)點(diǎn)2。0節(jié)點(diǎn)25 。0節(jié)點(diǎn)26 。節(jié) 點(diǎn)27 。0節(jié)點(diǎn)28 。0節(jié)點(diǎn)29 。。0 20節(jié) 點(diǎn)30