电力系统潮流计算完整c语言程序(含网损计算的最终版)

for(j=0;j<Bus_Num;j++)
{
ia[i]=ia[i]+gY_G[n][j]*ge[j]-gY_B[n][j]*gf[j];
ib[i]=ib[i]+gY_G[n][j]*gf[j]+gY_B[n][j]*ge[j];
}
}
for(i=0,n=1;i<Bus_Num-1;i++,n++)
{
gDelta_PQ[2*i]=gDelta_P[i];
gDelta_PQ[2*i+1]=gDelta_Q[i];
}
if((fp=fopen("C:\\Documents and Settings\\Zorro\\桌面\\1\\data\\unbalance.txt","w"))==NULL)
if(gBus[n].Type==1)
gDelta_Q[i]=gDelta_Q[i]-gf[n]*(gY_G[n][j]*ge[j]-gY_B[n][j]*gf[j])+ge[n]*(gY_G[n][j]*gf[j]+gY_B[n][j]*ge[j]);
}
}
for(i=0;i<Bus_Num-1;i++)
{
gY_G[i][j]=0.0;
gY_B[i][j]=0.0;
}
for(l=0;l<Line_Num;l++)
{
i=gLine[l].No_I-1;
j=gLine[l].No_J-1;
r=gLine[l].R;
x=gLine[l].X;
d1=r*r+x*x;
{
printf("Can not open the file named 'unbalance.txt' \n");
return ;
}
fprintf(fp,"---Unbalance---\n");
for(i=0;i<=2*Bus_Num-3;i++)
{
fprintf(fp,"Unbalance[%d]=%10.5f\n",i+1,gDelta_PQ[i]);
#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;
// flow.cpp : Defines the entry point for the console application.
//
#include "stdafx.h"
#include <iostream.h>
#include <fstream.h>
#include "NEquation.h"
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);
if(gBus[n].Type==2)
{
gJaccobi[2*i+1][2*j]=0;
gJaccobi[2*i+1][2*j+1]=0;
}
else
gY_B[j][i]=gY_B[j][i]-b;
}
}
// output the Y matrix
if((fp=fopen("C:\\Documents and Settings\\Zorro\\桌面\\1\\data\\ymatrix.txt","w"))==NULL)
FILE *fp;
for(i=0;i<Bus_Num;i++)
{
if(gBus[i].Type==3)
{
a=gBus[i].Voltage;
gBus[i].Voltage=gBus[0].Voltage;
gBus[0].Voltage=a;
a=gBus[i].Phase;
gBus[i].Phase=gBus[0].Phase;
gBus[0].Phase=a;
a=gBus[i].GenP;
gBus[i].GenP=gBus[0].GenP;
gBus[0].GenP=a;
a=gBus[i].GenQ;
gBus[i].GenQ=gBus[Байду номын сангаас].GenQ;
}
{
if(gLine[i].No_J==bal_num)
gLine[i].No_J=1;
else if(gLine[i].No_J==1)
gLine[i].No_J=bal_num;
}
}
for(i=0;i<Line_Num-1;i++)
for(j=0;j<Line_Num-1;j++)
}
fclose(fp);
}
void GetYMatrix()
{
int i,j;
FILE *fp;
// calculate the Y matrix
int l,bal_num;
float r,x,d1,g,b;
for(i=0;i<Bus_Num-1;i++)
{
g=r/d1;
b=-x/d1;
if(gLine[l].k==0)
{
gY_G[i][i]=gY_G[i][i]+g;
gY_G[j][j]=gY_G[j][j]+g;
gY_B[i][i]=gY_B[i][i]+b+gLine[l].B;
gY_B[j][j]=gY_B[j][j]+b+gLine[l].B;
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()//Read the data
{
FILE *fp;
if(gBus[i].Type==3)
bal_num=i+1;
}
for(i=0;i<Line_Num;i++)
{
{
if(gLine[i].No_I==bal_num)
gLine[i].No_I=1;
else if(gLine[i].No_I==1)
gLine[i].No_I=bal_num;
gY_G[i][j]=gY_G[i][j]-g;
gY_G[j][i]=gY_G[j][i]-g;
gY_B[i][j]=gY_B[i][j]-b;
gY_B[j][i]=gY_B[j][i]-b;
}
else
{
gY_G[i][i]=gY_G[i][i]+g/gLine[l].k+g*(1-gLine[l].k)/gLine[l].k/gLine[l].k;
a=gBus[i].Type;
gBus[i].Type=gBus[0].Type;
gBus[0].Type=a;
}
}
for(i=0,n=1;i<Bus_Num-1;i++,n++)
{
gDelta_P[i]=gBus[n].GenP-gBus[n].LoadP;
if(gBus[n].Type==1)
{
printf("Can not open the file named 'ymatrix.txt' \n");
return ;
}
fprintf(fp,"---Y Matrix---\n");
for(i=0;i<=Bus_Num-1;i++)
{
for(j=0;j<=Bus_Num-1;j++)
}
fclose(fp);
}
void GetJaccobi()
{
FILE *fp;
int i,j,n;
float ia[Bus_Num-1],ib[Bus_Num-1];
for(i=0,n=1;i<Bus_Num-1;i++,n++)
{
ia[i]=0;
ib[i]=0;
{