华电潮流上机程序(含网损计算)最终版2015111.docVIP

// flow.cpp: 主项目文件。 #include stdafx.h #include NEquation.h #include math.h #include stdio.h #include config.h using namespace System; 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)=4; ob1.Value(1)=6; ob1.Run(); printf(x1=%f\n,ob1.Value(0)); printf(x2=%f\n,ob1.Value(1)); printf(Test OK!\n); } void GetData() //读取数据 { int i; FILE *fp; fp=fopen(E:\\1121180617邵天赐\\data\\data.txt,r); if(fp==NULL) { printf(Can not open the file named data.txt \n); return; } for(i=0;i=Bus_Num-1;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;i=Line_Num-1;i++) { fscanf(fp,%d,%d,%d,%f,%f,%f,%f,gLine[i].No,gLine[i].No_I,gLine[i].No_J, gLine[i].R,gLine[i].X,gLine[i].B,gLine[i].k); } printf(Open the file OK! \n); fclose(fp); } void GetYMatrix()// 由Z阵计算Y阵 { int i,j,bus1,bus2; FILE *fp; float r,x,d,g,b,g1,b1,g2,b2,g3,b3; for(i=0;i=Bus_Num-1;i++) { for(j=0;j=Bus_Num-1;j++) { gY_G[i][j]=0; gY_B[i][j]=0; } } for(i=0;i=Line_Num-1;i++) { if(gLine[i].k==0) { bus1=gLine[i].No_I-1; bus2=gLine[i].No_J-1; r=gLine[i].R; x=gLine[i].X; d=r*r+x*x; g=r/d; b=-x/d; 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+gLine[i].B; gY_B[bus2][bus2]=gY_B[bus2][bus2]+b+gLine[i].B; gY_B[bus1][bus2]=gY_B[bus1][bus2]-b; gY_B[bus2][bus1]=gY_B[bus2][bus1]-b; } else { bus1=gLine[i].No_I-1; bus2=gLine[i].No_J-1; r=gLine[i].R; x=gLine[i].X; d=r*r+x*x; g=r/d; b=-x/d; g1=g/gLine[i].k; b1=b/gLine[i].k; g2=g*(gLine[i].k-1)/gLine[i].k; b2=b*(gLine[i].k-1)/gLine[i].k; g3=g*(1-gLine[i].k)/(gLine[i].k*gLine[i].k); b3=b*(1-gLine[i].k)/(gLine[i].k*gL