【正文】 。cc=Complex::productComplex(bb , Complex::getconj(line[i].Y))。aa=Complex::productComplex(Complex::getconj(bus[statemp1].V), B)。statemp=line[i].start。=0。i{int statemp,endtemp。goto LOOP。//修正節(jié)點(diǎn)電壓// Bus::PrintfNodeV(bus,N)。gauss::gauss_slove(JacAug,x,NN)。//whether converbence看迭代是否結(jié)束if(icon==1){coutBus::JisuanJacAug(JacAug,X,bus,N)。//計(jì)算節(jié)點(diǎn)功率差值 Bus::PrintfNodeScal(X,bus,N)。//計(jì)算節(jié)點(diǎn)功率Bus::JisuanNodeScal(X,bus,N)。LOOP:Bus::JisuanNodeI(X,bus,N)。x=new double[NN]。for(i=0。double **JacAug。//計(jì)算節(jié)點(diǎn)導(dǎo)納矩陣Bus::PrintfNodeDnz(X,N)。for(i=0。//輸出結(jié)點(diǎn)參數(shù)Complex **X。}Bus::ScanfBusData(bus)。ibus[i].=0。//輸出支路參數(shù)Bus *bus=new Bus[N]。//動(dòng)態(tài)分配支路結(jié)構(gòu)體Line::ScanfLineData(line)。N)==0){return 0。//輸入支路個(gè)數(shù)if((Lamp。//i作為整個(gè)程序的循環(huán)變量int N=Bus::ScanfBusNo()。}include “” include include include include“” include“” include“” using namespace std。Complex::PrintfComplex(Complex::ComDivRea(line[i].etos,))。endtemp=line[i].end。i{int statemp,endtemp。coutt(bus[i].)。i{coutComplex::PrintfComplex(bus[i].V)。bus[i].V=Complex::Rec2Polar(bus[i].V)。i{bus[i]. = bus[i]. + bus[i].。}coutBus::JisuanNodeScal(X,bus,N)。dd=Complex::CaddC(aa,cc)。bb=Complex::subComplex(Complex::getconj(bus[endtemp1].V), Complex::getconj(bus[statemp1].V))。line[i].stoe=Complex::productComplex(bus[statemp1].V,dd)。cc=Complex::productComplex(bb , Complex::getconj(line[i].Y))。aa=Complex::productComplex(Complex::getconj(bus[statemp1].V), B)。statemp=line[i].start。=0。i{int statemp,endtemp。goto LOOP。//修正節(jié)點(diǎn)電壓// Bus::PrintfNodeV(bus,N)。gauss::gauss_slove(JacAug,x,NN)。//whether converbence看迭代是否結(jié)束if(icon==1){coutBus::JisuanJacAug(JacAug,X,bus,N)。//計(jì)算節(jié)點(diǎn)功率差值Bus::PrintfNodeScal(X,bus,N)。//計(jì)算節(jié)點(diǎn)功率Bus::JisuanNodeScal(X,bus,N)。LOOP:Bus::JisuanNodeI(X,bus,N)。x=new double[NN]。for(i=0。double **JacAug。//計(jì)算節(jié)點(diǎn)導(dǎo)納矩陣Bus::PrintfNodeDnz(X,N)。for(i=0。//輸出結(jié)點(diǎn)參數(shù)Complex **X。}Bus::ScanfBusData(bus)。ibus[i].=0。//輸出支路參數(shù)Bus *bus=new Bus[N]。//動(dòng)態(tài)分配支路結(jié)構(gòu)體Line::ScanfLineData(line)。N)==0){return 0。//輸入支路個(gè)數(shù)if((Lamp。//i作為整個(gè)程序的循環(huán)變量int N=Bus::ScanfBusNo()。謝謝include “” include include include include“” include“” include“” using namespace std。jj=ii。data[ii][7]=。data[ii][5]=。data[ii][3]=data[ii1][3]+data[ii1][1]+data[ii][7]。data[ii][1]=。} } coutU1。k++。data[ii][4]=Tf[k].GetQt(data[ii1][3]+data[ii1][1]+data[ii][7],data[ii1][4]+data[ii1][0]+data[ii][6],data[ii1][2])。data[ii][2]=Tf[k].GetUt(data[ii1][3]+data[ii1][1]+data[ii][7],data[ii1][4]+data[ii1][0]+data[ii][6],data[ii1][2])。data[ii][6]=。data[ii][0]=Tf[k].GetP0()。Tf[k].I0=e。Tf[k].Uk=c。Tf[k].SN=a。Transformer Tf[100]。ii++。data[ii][5]=。//參數(shù)分別為前節(jié)點(diǎn)的有功、無(wú)功、電壓data[ii][3]=Tm[j].GetPl(data[ii1][3]+data[ii1][1]+data[ii][7],data[ii1][4]+data[ii1][0]+data[ii][6],data[ii1][2])。data[ii][7]=Tm[j].GetQb(data[ii1][2])。data[ii][1]=Tm[j].GetQf(data[ii1][2])。Tm[j].Bl=c。Tm[j].Rl=a。Transmission_line Tm[100]。ii++。data[ii][7]=。data[ii][5]=。data[ii][3]=N[i].Pn。data[ii][1]=。N[i].Qn=b。int i=0。if(x==1){coutcinab。x!=3){coutcoutcinx。x!=2amp。x!=1amp。if(x!=0amp。int x。double e。double c。coutwhile(1){double a。data[ii][7]=。data[ii][5]=。data[ii][3]=Pd。data[ii][1]=。ii=0。double data[300][8]。第一篇：電力系統(tǒng)潮流計(jì)算程序include include include “” include “” include “”void main(){ coutPdQd。coutUd。//分別為0前面有功，1前面無(wú)功，2壓降，3中間有功，4中間無(wú)功,5變比,6后面有功，7后面無(wú)功int ii。data[ii][0]=。data[ii][2]=Ud。data[ii][4]=Qd。data[ii][6]=。ii++。double b。double d。double f。coutcinx。amp。amp。amp。}if(x==0)break。Node N[100]。N[i].Pn=a。data[ii][0]=。data[ii][2]=。data[ii][4]=N[i].Qn。data[ii][6]=。i++。}if(x==2){coutcinabc。int j=0。Tm[j].Xl=b。data[ii][0]=。data[ii][6]=。data[ii][2]=Tm[j].GetUl(data[ii1][3]+data[ii1][1]+data[ii][7],data[ii1][4]+data[ii1][0]+data[ii][6],data[ii1][2])。data[ii][4]=Tm[j].GetQl(data[ii1][3]+data[ii1][1]+data[ii][7],data[ii1][4]+data[ii1][0]+data[ii][6],data[ii1][2])。j++。}if(x==3){coutcinabcdef。int k=0。Tf[k].Pk=b。Tf[k].P0=d。Tf[k].Kk=f。data[ii][1]=Tf[k].GetQ0()。data[ii][7]=。//參數(shù)要改data[ii][3]=Tf[k].GetPt(data[ii1][3]+data[ii1][1]+data[ii][7],data[ii1][4]+data[ii1][0]+data[ii][6],data[ii1][2])。data[ii][5]=Tf[k].Kk。ii++。data[ii][0]=。} data[ii][2]=U1。data[ii][4]=data[ii1][4]+data[ii1][0]+data[ii][6]。data[ii][6]=。int jj。if(fabs(data[ii][2]data[ii1][2])第二篇：電力系統(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::getc