【正文】 ex::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é)點電壓// Bus::PrintfNodeV(bus,N)。gauss::gauss_slove(JacAug,x,NN)。//whether converbence看迭代是否結束if(icon==1){coutBus::JisuanJacAug(JacAug,X,bus,N)。//計算節(jié)點功率差值Bus::PrintfNodeScal(X,bus,N)。//計算節(jié)點功率Bus::JisuanNodeScal(X,bus,N)。LOOP:Bus::JisuanNodeI(X,bus,N)。x=new double[NN]。for(i=0。double **JacAug。//計算節(jié)點導納矩陣Bus::PrintfNodeDnz(X,N)。for(i=0。//輸出結點參數(shù)Complex **X。}Bus::ScanfBusData(bus)。ibus[i].=0。//輸出支路參數(shù)Bus *bus=new Bus[N]。//動態(tài)分配支路結構體Line::ScanfLineData(line)。N)==0){return 0。//輸入支路個數(shù)if((Lamp。//i作為整個程序的循環(huán)變量int N=Bus::ScanfBusNo()。謝謝include “” include include include include“” include“” include“” using namespace std。fout}foutfoutfout()。ki=L[k].busi1。ji=G[j].busno1。foutfoutfor(i=0。} }//ofstream fout。} if(sqrt(e[k]*e[k]+f[k]*f[k])Vl=sqrt(e[k]*e[k]+f[k]*f[k])。kVh){Vh=sqrt(e[k]*e[k]+f[k]*f[k])。double PG=0,PL=0,Prr=0,Vh=sqrt(e[0]*e[0]+f[0]*f[0]),Vl=sqrt(e[0]*e[0]+f[0]*f[0])。QZL[j][i]+=(e[j]*e[j]+f[j]*f[j])*(L[k].B)。pr[j][i]=I2*L[k].R。PZL[j][i]=(x1*x2+y1*y2)。y2=f[j]*YG[j][i]+e[j]*YB[j][i]。y1=f[j]*L[k].kf[i]。QZL[i][j]+=(e[i]*e[i]+f[i]*f[i])*(L[k].B)。pr[i][j]=I2*L[k].R。PZL[i][j]=(x1*x2+y1*y2)。y2=f[i]*YG[i][j]+e[i]*YB[i][j]。y1=f[i]/L[k].kf[j]。j=L[k].busj1。for(int k=0。}//支路//double PZL[100][101]={0},QZL[100][101]={0},pr[100][101]={0},qx[100][101]={0}。}P[i]=e[i]*a+f[i]*b。ja+=(YG[i][j]*e[j]YB[i][j]*f[j])。idouble a=0,b=0。C[i]+=D[i]。for(i=0。iJa[i][2*N2]=W[i]。}}//中間迭代過程//coutcoutif(itJacobi(YG,YB,e,f,type,N,Ja)。iif(itit=fabs(W[i])。double it=fabs(W[0])。f[i]=C[2*i+1]。tfor(i=0。}double W[100]={0},Ja[100][101]={0}。iC[2*i]=1。double C[100]={0},D[100]={0}。}}} //修正后// for(l=0。jif(i!=j){YG[i][i]=YG[i][i]+(YG[i][j]*K[i][j]*K[i][j])。K[j][i]=1。ifor(j=i。} } for(l=0。BB[i][j]=0。iYB[i][j]=0。//節(jié)點導納矩陣形成// double YB[50][50],YG[50][50],BB[50][50],K[50][50]。LD=NGS。} L[i].B=m5[i]。L[i].R=m3[i]。iL[i].busi=m1[i]。}L=new Line[ZLs]。finm2[i]m3[i]m4[i]m5[i]m6[i]。i++){finm1[i]。G[i].Vg=m3[i]。iG[i].busno=m1[i]。} G=new Generator[GS]。finm4[i]m3[i]。i++){finm1[i]。B[i].Qd=m4[i]。B[i].type=m2[i]。for(i=0。N++。if(m1[i]==0)break。支路數(shù)//for(i=0。float m3[50],m4[50],m5[50],m6[50]。if(!fin){coutgetchar()。Line*L。//節(jié)點數(shù)發(fā)電機數(shù)負荷數(shù) BUS*B。jvoid main(){ ifstream fin。k=0。j}}a[n1][n]/=a[n1][n1]。ic=a[i][k]。for(j=k。for(k=0。}}}}} }//高斯消元法解方程組函數(shù)// void gauss(double a[][101],int n){int i,j,k。}else {Ja[2*i+1][2*j]=2*e[i]。if(type[i]==1){Ja[2*i+1][2*j]=b[i]+YB[i][i]*e[i]YG[i][i]*f[i]。Ja[2*i][2*j]=a[i]YG[i][i]*e[i]YB[i][i]*f[i]。ka[i]+=(YG[i][k]*e[k]YB[i][k]*f[k])。}}else {double a[50]={0},b[50]={0}。Ja[2*i][2*j+1]=YB[i][j]*e[i]YG[i][j]*f[i]。Ja[2*i+1][2*j+1]=Ja[2*i][2*j]。Ja[2*i][2*j+1]=YB[i][j]*e[i]YG[i][j]*f[i]。ifor(j=0。} } //Jacobi矩陣// void Jacobi(double YG[][50],double YB[][50],double e[50],double f[50],int type[50],int N ,double Ja[100][101]){int i,j。if(type[i]==1)W[2*i+1]=dQ。W[2*i]=dP。}dV=V[i]*V[i]e[i]*e[i]f[i]*f[i]。jA+=YG[i][j]*e[j]YB[i][j]*f[j]。idouble A=0,B=0。int i,j。}。float B。float R。//支路信息 節(jié)點I 節(jié)點J R X B/2 k struct Line { int busi。float Vg。//發(fā)電機數(shù)據(jù) 節(jié)點號有功發(fā)電電壓幅值 struct Generator { int busno。float Qd。int type。第一篇：電力系統(tǒng)通用潮流計算C語言程序include include include include using namespace std。//節(jié)點號類型負荷有功負荷無功母線數(shù)據(jù)（類型 1＝PV節(jié)點，2＝PQ節(jié)點，3＝平衡節(jié)點）struct BUS { int busno。float Pd。}。float Pg。}。int busj。float X。float k。//deltaP deltaQ deltaV^2// void fun1(double YG[][50],double YB[][50],double e[],double f[],int type[],int N,double W[],double P[],double Q[],double V[]){ double dP=0,dQ=0,dV=0。for(i=0。for(j=0。B+=YG[i][j]*f[j]+YB[i][j]*e[j]。dP=P[i]e[i]*Af[i]*B。dQ=Q[i]f[i]*A+e[i]*B。else W[2*i+1]=dV。for(i=0。jif(i!=j){if(type[i]==1){Ja[2*i][2*j]=(YG[i][j]*e[i]+YB[i][j]*f[i])。Ja[2*i+1][2*j]=Ja[2*i][2*j+1]。}else {Ja[2*i][2*j]=YG[i][j]*e[i]+YB[i][j]*f[i]。Ja[2*i+1][2*j+1]=Ja[2*i+1][2*j]=0。for(int k=0。b[i]+=(YG[i][k]*f[k]+YB[i][k]*e[k])。Ja[2*i][2*j+1]=b[i]+YB[i][i]*e[i]YG[i][i]*f[i]。Ja[2*i+1][2*j+1]=a[i]+YG[i][i]*e[i]+YB[i][i]*f[i]。Ja[2*i+1][2*j+1]=2*f[i]。double c。kc=a[k][k]。jfor(i=k+1。for(j=k。for(k=n2。k)for(j=k+1。int N=0,GS=0,LD=0,ZLs=0。Generator*G。//從文本中讀入原始數(shù)據(jù)到數(shù)組中//(“C:”)。}int m1[50]={0},m2[50]={0}。int i,j,l。i++){finm1[i]。finm2[i]m3[i]m4[i]。} B=new BUS[N]。iB[i].busno=m1[i]。B[i].Pd=m3[i]。}for(i=0。if(m1[i]==0)break。GS++。for(i=0。G[i].Pg=m4[i]。}for(i=0。if(m1[i]==0)break。ZLs++。for(i=0。L[i].busj=m2[i]。L[i].X=m4[i]。L[i].k=m6[i]。()。for(i=0。YG[i][j]=0。K[i][j]=1。l} for(i=0。jK[i][j]=K[j][i]。}for(j=0。YB[i][i]=YB[i][i]+(YB[i][j]*K[i][j]*K[i][j]+BB[i][j])。l// 求A=e+f// double e[50]={0},f[50]={0}。for(i=0。} else C[2*i]=V[i]。//調用Jacobi函數(shù)和高斯函數(shù)//for(int t=1。ie[i]=C[2*i]。}fun1(YG,YB,e,f,type,N,W,P,Q,V)。for(i=1。j=i。for(i=0。}//高斯消元法解方程//gauss(Ja,2*N2)。iD[i]=Ja[i][2*(N1)]。}}//平衡節(jié)點//for(i=0。for(int j=0。b+=(YB[i][j]*e[j]+YG[i][j]*f[j])。Q[i]=f[i]*ae[i]*b。double x1=0,x2=0,y1=0,y2=0,I2=0。ki=L[k].busi1。x1=e[i]/L[k].ke[j]。x2=e[i]*YG[i][j]f[i]*YB[i][j]。QZL[i][j]=(x1*y2x2*y1)。I2=(PZL[i][j]*PZL[i][j]+QZL[i][j]*QZL[i][j])/(e[i]*e[i]+f[i]*f[i])。qx[i][j]=I2*L[k].X(e[i]*e[i]+f[i]*f[i]+e[j]*e[j]+f[j]*f[j])*L[k].B。x1=e[j]*L[k].ke[i]。x2=e[j]*YG[j][i]f[j]*YB[j][i]。QZL[j][i]=(x1*y2x2*y1)。I2=(PZL[j][i]*PZL[j][i]+QZL[j][i]*QZL[j][i])/(e[j]*e[j]+f[j]*f[j])。qx[j][i]=I2*L[k].X(e[i]*e[i]+f[i]*f[i]+e[j]*e[j]+f[j]*f[j])*L[k].B。}//全網(wǎng)數(shù)據(jù)// int high=1,low=1。for(k=0。high=k+1。low=k+1。(“C:”)。ifout}for(j=0。fout}foutfoutfor(k=0。j=L[k].busj1。}第二篇：電力系統(tǒng)潮流計算程序電力系統(tǒng)潮流計算c語言程序，兩行，大家可以看看，仔細研究，然后在這個基礎上修改