华电-通用潮流上机程序.doc

include "stdafx.h"
#include
#include
#include "NEquation.h"
#include "math.h"
#include "config.h"
using namespace std;

void GetData(int& nBus, int& nBranch)//Read the data
{
FILE *fp;
int i;

if((fopen_s(&fp,"data\\data5.txt","r"))!=0)
{
printf("Can not open the file named 'data.txt' \n");
return;
}
fscanf_s(fp,"%d,%d",&nBus,&nBranch,100);
for(i=0;ii=0;i

for(i=0;ibus1 != 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;
gY_B[bus2][bus2]=gY_B[bus2][bus2]+b+gBranch[i].B;
gY_B[bus1][bus2]=gY_B[bus1][bus2]-b;
gY_B[bus2][bus1]=gY_B[bus2][bus1]-b;
}
else
{
break;
}

}
else
{
// add your codes
}
}

// output the Y matrix
if((fopen_s(&fp,"data\\ymatrix1.txt","w"))!=0)
{
printf("Can not open the file named 'ymatrix.txt' \n");
return ;
}

fprintf(fp,"---Y Matrix---\n");
for(i=0;i
{
if(gBus[i].Type==2)
{
gf[i]=gBus[i].Voltage*sin(gBus[i].Phase);
ge[i]=gBus[i].Voltage*cos(gBus[i].Phase);
}
else
{
gf[i]=0;
ge[i]=1;
}
}
}

void GetUnbalance(int& nBus)

{ int i,j;
float sum1=0,sum2=0;
FILE *fp;

for(i=0;i
if (gBus[i].Type!=2)
{
sum1=0;
sum2=0;
for(j=0;jge[i]*(gY_G[i][j]*gf[j]+gY_B[i][j]*ge[j]);
}
gDelta_P[i]= gBus[i].GenP-gBus[i].LoadP-sum1;
gDelta_Q[i]= gBus[i].GenQ-gBus[i].LoadQ-sum2;

if (gBus[i].Type==1)
gDelta_Q[i]= gBus[i].Voltage* gBus[i].Voltage-(gf[i]*gf[i]+ge[i]*ge[i]);
}
}

for (i=0;ij/2) // i is or not=j.
{ s1=0;
s2=0;
for (k=0;kobi[i][j+1]=gY_G[i/2][j/2]*ge[j/2]+gY_B[i/2][j/2]*gf[j/2]+s2;

if (gBus[i/2].Type==0) //PQ pioint and PV Point difference.
{
gJaccobi[i+1][j]=-(gY_G[i/2][j/2]*ge[j/2]+gY_B[i/2][j/2]*gf[j/2])+s2;
gJaccobi[i+1][j+1]=gY_G[i/2][j/2]*gf[j/2]-gY_B[i/2][j/2]*ge[j/2]-s1;
}

else //pv bus.
{
gJaccobi[i+1][j]=2*gf[i/2];
gJaccobi[i+1][j+1]=2*ge[i/2];
}
}
else
{
gJaccobi[i][j]=gY_G[i/2][j/2]*gf[i/2]-gY_B[i/2][j/2]*ge[i/2];
gJaccobi[i][j+1]=gY_G[i/2][j/2]*ge[i/2]+gY_B[i/2][j/2]*gf[i/2];
if (gBus[i/2].Type==0)
{
gJaccobi[i+1][j]=-gJaccobi[i][j+1];
gJaccobi[i+1][j+1]=gJaccobi[i][j];
}

else
{
gJaccobi[i+1][j]=0;
gJaccobi[i+1][j+1]=0;
}

}

// output Jaccobimatrix.
if((fopen_s(&fp,"data\\Jaccobimatrix.txt","w"))!=0)
{
printf("Can not open the file named 'Jaccobimatrix.txt' \n");
return ;
}

fprintf(fp,"---JaccobiMatrix---\n");

for(i=0;i<2*nBus;i++)
{
fprintf_s(fp,"%d ,",i+1);

for(j=0;j<2*nBus;j++)
{
fprintf_s(fp,"[%10.5f]",gJaccobi[i][j]);
}
fprintf_s(fp,"\n");

}
fclose(fp);

}

void GetRevised(int& nBus)
{
int i,j;
NEquation ob1;
FILE *fp;

ob1.SetSize(2*nBus);

for(i=0;i<2*nBus;i++)
{
for(j=0;j<2*nBus;j++)
ob1.Data(i,j)=gJaccobi[i][j];

ob1.Value(i)=gDelta_PQ[i];
}

ob1.Run();

for(i=0;i<2*nBus;i++)
{
gDelta_fe[i]=ob1.Value(i);
}

for (i=0;ielta_e matrix.txt","w"))!=0)
{
printf("Can not open the file named 'Delta_f and Delta_e matrix.txt' \n");
return ;
}

fprintf(fp,"---Delta_f and Delta_e Matrix---\n");

for(i=0;i<2*nBus;i++)
{
fprintf(fp," Delta %d=%10.8f \n",i+1,gDelta_fe[i]);

}
fprintf_s(fp,"\n");
fclose(fp);
}

void GetNewValue(int& nBus)
{
int i;
FILE *fp;
for (i=0;i"Can not open the file named 'gf and ge matrix.txt' \n");
return ;
}

fprintf(fp,"---gf and ge Matrix---\n");

for(i=0;iint argc, char* argv[])
{
int i,Count_Num,j,n,nBus,nBranch;
float maxValue,s1,s2,k,su1,su2,deltalinep[12][12]={0,0},deltalineq[12][12]={0,0},deltap[12][12]={0,0},deltaq[12][12]={0,0};
FILE *fp;
float g,b,d,r,x;
GetData(nBus,nBranch);
GetYMatrix(nBus,nBranch);
SetInitial(nBus);

for(Count_Num=0;Count_Num<=20;Count_Num++)
{
GetUnbalance(nBus);
GetJaccobi(nBus);
GetRevised(nBus);
GetNewValue(nBus);

maxValue=fabs(gDelta_fe[0]);
for(i=1;i<=2*(nBus-1)-1;i++)
{
if(maxValuemaxValueut zhuru power.
for (i=0;i][j]*gf[j])-ge[i]*(gY_G[i][j]*gf[j]+gY_B[i][j]*ge[j]);
}
gBus[i].GenP=su1; gBus[i].GenQ=su2;
printf("zhurup[%d]=%10.7f zhuruq[%d]=%10.7f \n",i+1,gBus[i].GenP,i+1,gBus[i].GenQ);
}
}

for(n=0; nge[j])+gf[i]*(gf[i]-gf[j]))*g+(ge[i]*gf[j]-ge[j]*gf[i])*b;

deltaq[i][j]=-(ge[i]*ge[i]+gf[i]*gf[i])*gBranch[n].B+(ge[i]*gf[j]-ge[j]*gf[i])*g-(ge[i]*(ge[i]-ge[j])+gf[i]*(gf[i]-gf[j]))*b;

deltap[j][i]=(ge[j]*(ge[j]-ge[i])+gf[j]*(gf[j]-gf[i]))*g+(ge[j]*gf[i]-ge[i]gf[j])*b;

deltaq[j][i]=-(ge[j]*ge[j]+gf[j]*gf[j])*gBranch[n].B+(ge[j]*gf[i]-ge[i]*gf[j])*g-(ge[j]*(ge[j]-ge[i])+gf[j]*(gf[j]-gf[i]))*b;

}
else break;
}
else
{
k=gBranch[n].k;
deltap[i][j]= (ge[i]*ge[i]+gf[i]*gf[i])*k*(k-1)*g+(ge[i]*(ge[i]-ge[j])+gf[i]*(gf[i]-gf[j]))*g*k+(ge[i]*gf[j]-ge[j]*gf[i])*b*k;

deltaq[i][j]=-(ge[i]*ge[i]+gf[i]*gf[i])*k*(k-1)*b+(ge[i]*gf[j]-ge[j]*gf[i])*g*k-(ge[i]*(ge[i]-ge[j])+gf[i]*(gf[i]-gf[j]))*b*k;

deltap[j][i]=(ge[j]*ge[j]+gf[j]*gf[j])1-k)*g+(ge[j]*(ge[j]-ge[i])+gf[j]*(gf[j]-gf[i]))*g*k+(ge[j]*gf[i]-ge[i]*gf[j])*b*k;

deltaq[j][i]=-(ge[j]*ge[j]+gf[j]*gf[j])*(1-k)*b+(ge[j]*gf[i]-ge[i]*gf[j])*g*k-(ge[j]*(ge[j]-ge[i])+gf[j]*(gf[j]-gf[i]))*b*k;
}

}

for (i0;iineq[i][j]);
}

s1=0;s2=0;
for(i=0;i

if((fopen_s(&fp,"data\\bus voltage and zhuru power.txt","w"))!=0)
{
printf("Can not open the file named 'bus voltage and zhuru power.txt' \n");
return 0;
}

fprintf(fp,"---bus voltage and zhuru power---\n");

for(i=0;izhuruq[%d]=%10.7f \n",i+1,gBus[i].GenP-gBus[i].LoadP,i+1,gBus[i].GenQ-gBus[i].LoadQ);

}
fprintf_s(fp,"\n");
for(i=0;in");

for(i=0;i