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

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

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;ob1.Data(1,1)=1;ob1.Value(0)=5;ob1.Val

2、ue(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("Cannotopenthefilenamed'data.txt'");return;}fscanf(fp,"%d,%d",&nBus,&nBranch);for(i=0;i

3、anf(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

4、.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、h[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;/*hengzong*/gY_B[bus2][bus2]=gY_B[bus2][bus2]+b+gBranch[i

6、].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;//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]

7、.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