电力系统稳态潮流机完整可运行程序.docVIP

// flow.cpp : Defines the entry point for the console application. // #include stdafx.h #include iostream.h #include fstream.h #include NEquation.h #include math.h #include config.h void test() { NEquation ob1; ob1.SetSize(2); ob1.Data(0,0)=1; ob1.Data(0,1)=2; ob1.Data(1,0)=2; ob1.Data(1,1)=1; ob1.Value(0)=5; ob1.Value(1)=4; ob1.Run(); printf(x1=%f\n,ob1.Value(0)); printf(x2=%f\n,ob1.Value(1)); } void GetData(int nBus, int nBranch)//Read the data { FILE *fp; int i; if((fp=fopen(data\\data5.txt,r))==NULL) { printf(Can not open the file named data.txt \n); return; } fscanf(fp,%d,%d,nBus,nBranch); for(i=0;inBus;i++) { fscanf(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); } for(i=0;inBranch;i++) { fscanf(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); } /*半容纳*/ fclose(fp); } void GetYMatrix(int nBus, int nBranch) { int i,j,bus1,bus2; float r,x,d,g,b; 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++) { if(gBranch[i].k==0) { 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(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;/*heng zong*/ 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 { gY_B[bus1][bus1]=gY_B[bus1][bus1]+b; // add your codes } } else {bus1=gBranch[i].No_I-1; bus2=gBranch[i].No_J-1; r=gBranch[i].R; x=gBranch[