電力系統潮流計算的c語言實現(編輯修改稿)

　

【文章內容簡介】 n%5s%5s%10s%12s%12s%12sn ”, “I”, “J”, “PIJ”，“QIJ”, “PJI”, “QJI”)。for(i =1。i {j = idk[i]。dkk[i] = v[j] * v[j] / dkk[i]。} for(l = 1。l {ii = iabs(izl[l])。// izl[]: 支路左節點號。jj = iabs(jzl[l])。// jzl[]: 支路右節點號。i=nnew[ii]。j=nnew[jj]。//轉換為新節點號。ii = newsort[ii]。jj = newsort[jj]。//轉化為相應舊號 r = zr[l]。x = zx[l]。yk = zyk[l]。vi = v[i]。//v[]: 電壓幅值。ci = va[i]。//va[]: 電壓相角。vj = v[j]。cj = va[j]。//支路左、右節點電壓值由極坐標轉換為直角坐標 ei = vi * cos(ci)。fi = vi * sin(ci)。// ei: 支路左節點電壓實部，fi: 支路左節點電壓虛部。ej = vj * cos(cj)。fj = vj * sin(cj)。// ej: 支路右節點電壓實部，fj: 支路右節點電壓虛部。if((izl[l] 0)||(jzl[l] 0))// 變壓器支路。{ if(izl[l] 0){ei = ei / yk。fi = fi / yk。}// yk=zyk[l]else{ej = ej / yk。fj = fj / yk。}yk = 。}de = eifj。zf = r * r + x * x。fii =(de * r + df * x)/ zf。fir =(df * rfir * ei。pji =fii * ejvi * vi * yk。qji = qji1。i ++){ip = i。k1 = iabs(izl[i])。k3 = iabs(jzl[i])。for(j = i + 1。j {k2 = iabs(izl[j])。k4 = iabs(jzl[j])。if(k2 k1 ||(k2 == k1 amp。amp。 k4 k3)){ip = j。k1 = k2。k3 = k4。}}if(i!= ip){iswap(amp。izl[i], amp。izl[ip])。iswap(amp。jzl[i], amp。jzl[ip])。swap(amp。zr[i], amp。zr[ip])。swap(amp。zx[i], amp。zx[ip])。swap(amp。zyk[i], amp。zyk[ip])。}} for(i = 1。i k = nob[i]。nob[i] = newsort[k]。printf(“nob[%d] = %dn”, i, nob[i])。} for(i = 1。i for(j = i+1。j {if(nob[i] nob[j]){k = nob[i]。nob[i] = nob[j]。nob[j] = k。}} } printf(“nob[%d] = %dn”, i, nob[i])。}void yy1(){//**** 本函數形成節點導納陣(不包括接地支路)****// int j1。double r, x, yk, zf, gij, bij。int i, j, i1, l, ll。for(i = 1。i {gii[i] = 。bii[i] = 。}// 導納陣對角元(與節點一一對應)先清零。for(i = 1。i yg[i] = 。yb[i] = 。} //導納陣非零非對角元(與支路一一對應)先清零。ll = 1。for(l = 1。l i1 = izl[l]。j1 = jzl[l]。i = iabs(i1)。j = iabs(j1)。if(i == j)// 支路左節點號。// 支路右節點號。// 變壓器支路有一節點號為負值。continue。//排除左、右節點號相等的情況。r = zr[l]。x = zx[l]。yk = zyk[l]。// zr[],zx[],zyk[]:支路三參數。zf = r * r + x * x。gij = r / zf。//bij =x / zf。bij =1/x。yg[ll] = yg[ll]bij。gii[i] = gii[i] + gij。bii[i] = bii[i] + bij。gii[j] = gii[j] + gij。bii[j] = bii[j] + bij。if((i!= iabs(izl[l + 1]))||(j!= iabs(jzl[l + 1])))ll++。} // 打印導納矩陣。對角元實部為gii，虛部為bii，// 非零非對角元實部為yb[]，虛部為yb[]，列足碼為iy[]。fprintf(fp2, “*******GII(1),BII(1)********n”)。printf2(gii,bii,n)。}void y3(){//**** 本函數形成節點導納陣，追加接地支路 ****//int j1。double r, x, yk, zf, gij, bij。int i, j, i1, l, ll, kk = 0。for(i = 1。i {j = idk[i]。bii[j] = bii[j]gij / yk。yb[ll] = yb[ll]bij / yk。}if((i!= iabs(izl[l + 1]))||(j!= iabs(jzl[l + 1]))){ll++。kk = 0。}elsekk = 1。} // 打印導納矩陣。對角元實部為gii，虛部為bii，// 非零非對角元實部為yb[]，虛部為yb[]，列足碼為iy[]。fprintf(fp2, “*******GII,BII********”)。printf2(gii,bii,n)。fprintf(fp2, “n*******YYYYY********”)。for(i = 1。i {if(i % 2 == 1)fprintf(fp2, “n”)。fprintf(fp2, “%%%8d”, yg[i], yb[i], iy[i])。} }第二篇：電力系統通用潮流計算C語言程序include include include include using namespace std。//節點號類型負荷有功負荷無功母線數據（類型 1＝PV節點，2＝PQ節點，3＝平衡節點）struct BUS { int busno。int type。float Pd。float Qd。}。//發電機數據 節點號有功發電電壓幅值 struct Generator { int busno。float Pg。float Vg。}。//支路信息 節點I 節點J R X B/2 k struct Line { int busi。int busj。float R。float X。float B。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。int i,j。for(i=0。idouble A=0,B=0。for(j=0。jA+=YG[i][j]*e[j]YB[i][j]*f[j]。B+=YG[i][j]*f[j]+YB[i][j]*e[j]。}dV=V[i]*V[i]e[i]*e[i]f[i]*f[i]。dP=P[i]e[i]*Af[i]*B。W[2*i]=dP。dQ=Q[i]f[i]*A+e[i]*B。if(type[i]==1)W[2*i+1]=dQ。else W[2*i+1]=dV。} } //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。for(i=0。ifor(j=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][2*j+1]=YB[i][j]*e[i]YG[i][j]*f[i]。Ja[2*i+1][2*j]=Ja[2*i][2*j+1]。Ja[2*i+1][2*j+1]=Ja[2*i][2*j]。}else {Ja[2*i][2*j]=YG[i][j]*e[i]+YB[i][j]*f[i]。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+1][2*j]=0。}}else {double a[50]={0},b[50]={0}。for(int k=0。ka[i]+=(YG[i][k]*e[k]YB[i][k]*f[k])。b[i]+=(YG[i][k]*f[k]+YB[i][k]*e[k])。Ja[2*i][2*j]=a[i]YG[i][i]*e[i]YB[i][i]*f[i]。Ja[2*i][2*j+1]=b[i]+YB[i][i]*e[i]YG[i][i]*f[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+1][2*j+1]=a[i]+YG[i][i]*e[i]+YB[i][i]*f[i]。}else {Ja[2*i+1][2*j]=2*e[i]。Ja[2*i+1][2*j+1]=2*f[i]。}}}}} }//高斯消元法解方程組函數// void gauss(double a[][101],int n){int i,j,k。double c。for(k=0。kc=a[k][k]。for(j=k。jfor(i=k+1。ic=a[i][k]。for(j=k。j}}a[n1][n]/=a[n1][n1]。for(k=n2。k=0。k)for(j=k+1。jvoid main(){ ifstream fin。int N=0,GS=0,LD=0,ZLs=0。//節點數發電機數負荷數 BUS*B。Generator*G。Line*L。//從文本中讀入原始數據到數組中//(“C:”)。if(!fin){coutgetchar()。}int m1[50]={0},m2[50]={0}。float m3[50],m4[50],m5[50],m6[50]。int i,j,l。支路數//for(i=0。i++){finm1[i]。if(m1[i]==0)break。finm2[i]m3[i]m4[i]。N++。} B=new BUS[N]。for(i=0。iB[i].busno=m1[i]。B[i].type=m2[i]。B[i].Pd=m3[i]。B[i].Qd=m4[i]。}for(i=0。i++){finm1[i]。if(m1[i]==0)break。finm4[i]m3[i]。GS++。} G=new Generator[GS]。for(i=0。iG[i].busno=m1[i]。G[i].Pg=m4[i]。G[i].Vg=m3[i]。}for(i=0。i++){finm1[i]。if(m1[i]==0)break。finm2[i]m3[i]m4[i]m5[i]m6[i]。ZLs++。}L=new Line[ZLs]。for(i=0。iL[i].busi=m1[i]。L[i].busj=m2[i]。L[i].R=m3[i]。L[i].X=m4[i]。} L[i].B=m5[i]。L[i].k=m6[i]。LD=NGS。()。//節點導納矩陣形成// double YB[50][50],YG[50][50],BB[50][50],K[50][50]。for(i=0。iYB[i][j]=0。YG[i][j]=0。BB[i][j]=0。K[i][j]=1。} } for(l=0。l} for(i=0。ifor(j=i。jK[i][j]=K[j][i]。K[j][i]=1。}for(j=0。jif(i!=j){YG[i][i]=YG[i][i]+(YG[i][j]*K[i][j]*K[i][j])。YB[i][i]=YB[i][i]+(YB[i][j]*K[i][j]*K[i][j]+BB[i][j])。}}} //修正后// for(l=0。l// 求A=e+f// double e[50]={0},f[50]={0}。double C[100]={0},D[100]={0}。for(i=0。iC[2*i]=1。} else C[2*i]=V[i]。}double W[100]={0},Ja[100][101]={0}。//調用Jacobi函數和高斯函數//for(int t=1。t