华电 通用潮流上机程序.docVIP

#include stdafx.h #include iostream #include fstream #include NEquation.h #include math.h #include config.h using namespace std; void GetData(int nBus, int nBranch)//Read the data { FILE *fp; int i; if((fopen_s(fp,data\\data5.txt,r))!=0) { printf(Can not open the file named data.txt \n); return; } fscanf_s(fp,%d,%d,nBus,nBranch,100); for(i=0;inBus;i++) { fscanf_s(fp,%d,%f,%f,%f,%f,%f,%f,%d,gBus[i].No,gBus[i].Voltage,gBus[i].Phase, gBus[i].GenP,gBus[i].GenQ,gBus[i].LoadP,gBus[i].LoadQ,gBus[i].Type,100); } for(i=0;inBranch;i++) { fscanf_s(fp,%d,%d,%d,%f,%f,%f,%f,gBranch[i].No,gBranch[i].No_I,gBranch[i].No_J, gBranch[i].R,gBranch[i].X,gBranch[i].B,gBranch[i].k,100); } fclose(fp); } void GetYMatrix(int nBus, int nBranch) { int i,j,bus1,bus2; float r,x,d,g,b,k; FILE *fp; for(i=0;inBus;i++) { for(j=0;jnBus;j++) { gY_G[i][j]=0; gY_B[i][j]=0; } } for(i=0; inBranch; i++) { bus1=gBranch[i].No_I-1; bus2=gBranch[i].No_J-1; r=gBranch[i].R; x=gBranch[i].X; d=r*r+x*x; g=r/d; b=-x/d; if(gBranch[i].k==0) { if(bus1 != bus2) { gY_G[bus1][bus1]=gY_G[bus1][bus1]+g; gY_G[bus2][bus2]=gY_G[bus2][bus2]+g; gY_G[bus1][bus2]=gY_G[bus1][bus2]-g; gY_G[bus2][bus1]=gY_G[bus2][bus1]-g; gY_B[bus1][bus1]=gY_B[bus1][bus1]+b+gBranch[i].B; gY_B[bus2][bus2]=gY_B[bus2][bus2]+b+gBranch[i].B; gY_B[bus1][bus2]=gY_B[bus1][bus2]-b; gY_B[bus2][bus1]=gY_B[bus2][bus1]-b; } else { break; } } else { // add your codes } } // output the Y matrix if((fopen_s(fp,data\\ymatrix1.txt,w))!=0) { printf(Can not open the file named ymatrix.txt \n); return ; } fprintf(fp,---Y Matrix---\n); for(i=0;inBus;i++) { for(j=0;jnBus;j++) fprintf_s(fp,Y(%d,%d)=(%10.5f,%10.5f)\n,i+1,j+1,gY_G[i][j],gY_B[i][j]); } fclose(fp); } void SetInitial(int nBus) { int i; for(i=0;inBus;i++) { if(gBus[i].Type==2) { gf[i]=gBus[i].Voltage*sin(gBus[i].Phase); ge[i]