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

《电力系统稳态潮流上机完整可运行程序.doc》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、//flow.cpp:Definestheentrypointfortheconsoleapplication.//#include"stdafx.h"#include#include#include"NEquation.h"#include"math.h"#include"config.h"voidtest(){NEquationob1;ob1.SetSize(2);ob1.Data(0,0)=1;ob1.Data(0,1)=2;ob1.Data(1,0

2、)=2;ob1.Data(1,1)=1;ob1.Value(0)=5;ob1.Value(1)=4;ob1.Run();printf("x1=%f",ob1.Value(0));printf("x2=%f",ob1.Value(1));}voidGetData(int&nBus,int&nBranch)//Readthedata{FILE*fp;inti;if((fp=fopen("data\data5.txt","r"))==NULL){printf("Cannotopenthefilen

3、amed'data.txt'");return;}fscanf(fp,"%d,%d",&nBus,&nBranch);for(i=0;i

4、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);}voidGetYMatrix(int&nBus,int&nBranch){inti,j,bus1,bus2;floatr,x,d,g,b;FILE*fp;for(i=0;i

5、us;i++){for(j=0;j

6、[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;/*hengzong*/gY_B[bus2][bus2]=gY_B[bus2][bus2]+b+gBranch[i].B;gY_B[bus1][bus2]=gY_B[bus1][bus2]-b;gY_

7、B[bus2][bus1]=gY_B[bus2][bus1]-b;}else{gY_B[bus1][bus1]=gY_B[bus1][bus1]+b;//addyourcodes}}else{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;gY_G[bus1][bus1]=gY_G[bus1][bus1]+gBranch[i].k*gBranch[i].

8、k*g;gY_G[bus2][bus2]=gY_G[bus2][bus2]+g;gY_G[bus1][bus2]=gY_G[bus1][bus2]-g*gBranch[i].k;gY_G[bus2][bus1]=gY_G[bus2][bus1]-g*gBranch[i].k;gY_B[bus1][bus1]=gY_B[bus1][bus1]+gBranch[i].k*gBranch[i].k*b;gY_B[bus2][bus2]=gY_B[bus2][b