华电潮流上机程序完整版197.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()//数据读入 { FILE *fp; int i; if((fp=fopen(E:\\新建文件夹\\flow\\data\\data.txt,r))==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); } fclose(fp); } void GetYMatrix() { int i,j; FILE *fp; // calculate the Y matrix int bus1,bus2; float r,x,d,g,b; 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++) { 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; //阻抗转化为导纳表示时的电纳 if(gLine[i].k==0) //不带变压器的普通支路 { gY_G[bus1][bus1]=gY_G[bus1][bus1]+g; //下标相同的为自导纳,不同的为互导纳 gY_G[bus2][bus2]=gY_G[bus2][bus2]+g; gY_G[bus1][bus2]=-g; //由于互导纳只涉及一个导纳,故直接等于支路导纳负值 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]=-b;