【文章內(nèi)容簡介】 的層次。第三篇：電力系統(tǒng)潮流計(jì)算程序電力系統(tǒng)潮流計(jì)算c語言程序，兩行，大家可以看看，仔細(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{coutComplex::PrintfComplex(bus[i].V)。coutt(bus[i].)。coutt(bus[i].)。cout}coutfor(i=0。i{int statemp,endtemp。statemp=line[i].start。endtemp=line[i].end。coutComplex::PrintfComplex(Complex::ComDivRea(line[i].stoe,))。Complex::PrintfComplex(Complex::ComDivRea(line[i].etos,))。cout} }return 0。}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{coutComplex::PrintfComplex(bus[i].V)。coutt(bus[i].)。coutt(bus[i].)。cout}coutfor(i=0。i{int statemp,endtemp。statemp=line[i].start。endtemp=line