华电 通用潮流上机程序.docx
《华电 通用潮流上机程序.docx》由会员分享,可在线阅读,更多相关《华电 通用潮流上机程序.docx(18页珍藏版)》请在冰豆网上搜索。
华电通用潮流上机程序
#include"stdafx.h"
#include
#include
#include"NEquation.h"
#include"math.h"
#include"config.h"
usingnamespacestd;
voidGetData(int&nBus,int&nBranch)//Readthedata
{
FILE*fp;
inti;
if((fopen_s(&fp,"data\\data5.txt","r"))!
=0)
{
printf("Cannotopenthefilenamed'data.txt'\n");
return;
}
fscanf_s(fp,"%d,%d",&nBus,&nBranch,100);
for(i=0;i{
fscanf_s(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,100);
}
for(i=0;i{
fscanf_s(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,100);
}
fclose(fp);
}
voidGetYMatrix(int&nBus,int&nBranch)
{
inti,j,bus1,bus2;
floatr,x,d,g,b,k;
FILE*fp;
for(i=0;i{
for(j=0;j{
gY_G[i][j]=0;
gY_B[i][j]=0;
}
}
for(i=0;i{
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;
if(gBranch[i].k==0)
{
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;
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
{
//addyourcodes
}
}
//outputtheYmatrix
if((fopen_s(&fp,"data\\ymatrix1.txt","w"))!
=0)
{
printf("Cannotopenthefilenamed'ymatrix.txt'\n");
return;
}
fprintf(fp,"---YMatrix---\n");
for(i=0;i{
for(j=0;jfprintf_s(fp,"Y(%d,%d)=(%10.5f,%10.5f)\n",i+1,j+1,gY_G[i][j],gY_B[i][j]);
}
fclose(fp);
}
voidSetInitial(int&nBus)
{
inti;
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;
}
}
}
voidGetUnbalance(int&nBus)
{inti,j;
floatsum1=0,sum2=0;
FILE*fp;
for(i=0;i{
if(gBus[i].Type!
=2)
{
sum1=0;
sum2=0;
for(j=0;j{
sum1=sum1+ge[i]*(gY_G[i][j]*ge[j]-gY_B[i][j]*gf[j])+gf[i]*(gY_G[i][j]*gf[j]+gY_B[i][j]*ge[j]);//chuqubalancepointaddyourcodes
sum2=sum2+gf[i]*(gY_G[i][j]*ge[j]-gY_B[i][j]*gf[j])-ge[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;i{
j=2*i;
gDelta_PQ[j]=gDelta_P[i];
gDelta_PQ[j+1]=gDelta_Q[i];
}
}
voidGetJaccobi(int&nBus)
{
inti,j,k;
floats1=0,s2=0;
FILE*fp;
for(i=0;i<2*nBus;i=i+2)//addyourcodes
for(j=0;j<2*nBus;j=j+2)
if(i/2==j/2)//iisornot=j.
{s1=0;
s2=0;
for(k=0;k{
s1=s1+gY_G[i/2][k]*gf[k]+gY_B[i/2][k]*ge[k];
s2=s2+gY_G[i/2][k]*ge[k]-gY_B[i/2][k]*gf[k];
}
gJaccobi[i][j]=gY_G[i/2][j/2]*gf[j/2]-gY_B[i/2][j/2]*ge[j/2]+s1;
gJaccobi[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)//PQpiointandPVPointdifference.
{
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//pvbus.
{
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;
}
}
//outputJaccobimatrix.
if((fopen_s(&fp,"data\\Jaccobimatrix.txt","w"))!
=0)
{
printf("Cannotopenthefilenamed'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);
}
voidGetRevised(int&nBus)
{
inti,j;
NEquationob1;
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;i{
gDelta_f[i]=gDelta_fe[2*i];
gDelta_e[i]=gDelta_fe[2*i+1];
}
//outputdelta_fanddelta_e.
if((fopen_s(&fp,"data\\Delta_fandDelta_ematrix.txt","w"))!
=0)
{
printf("Cannotopenthefilenamed'Delta_fandDelta_ematrix.txt'\n");
return;
}
fprintf(fp,"---Delta_fandDelta_eMatrix---\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);
}
voidGetNewValue(int&nBus)
{
inti;
FILE*fp;
for(i=0;iif(gBus[i].Type!
=2)
{
gf[i]=gf[i]+gDelta_f[i];
ge[i]=ge[i]+gDelta_e[i];}
if((fopen_s(&fp,"data\\gfandgematrix.txt","w"))!
=0)
{
printf("Cannotopenthefilenamed'gfandgematrix.txt'\n");
return;
}
fprintf(fp,"---gfandgeMatrix---\n");
for(i=0;i{
fprintf(fp,"gf[%d]=%10.6f,ge[%d]=%10.6f\n",i+1,gf[i],i+1,ge[i]);
}
fprintf_s(fp,"\n");
fclose(fp);
}
intmain(intargc,char*argv[])
{
inti,Count_Num,j,n,nBus,nBranch;
floatmaxValue,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;
floatg,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(maxValue{
maxValue=fabs(gDelta_fe[i]);
}
}
if(maxValue{
break;
}
}
printf("countnumber%d\n",Count_Num);
for(i=0;i{
printf("bus(%d).voltage=%10.6f,angle=%.6f,e=%.6f,f=%.6f\n",i+1,sqrt(ge[i]*ge[i]+gf[i]*gf[i]),atan(gf[i]/ge[i])/3.1415926*180,ge[i],gf[i]);
}
//outputzhurupower.
for(i=0;i{
if(gBus[i].Type!
=0)
{
su1=0;
su2=0;
for(j=0;j{
su1=su1+ge[i]*(gY_G[i][j]*ge[j]-gY_B[i][j]*gf[j])+gf[i]*(gY_G[i][j]*gf[j]+gY_B[i][j]*ge[j]);
su2=su2+gf[i]*(gY_G[i][j]*ge[j]-gY_B[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.7fzhuruq[%d]=%10.7f\n",i+1,gBus[i].GenP,i+1,gBus[i].GenQ);
}
}
for(n=0;n{
i=gBranch[n].No_I-1;
j=gBranch[n].No_J-1;
r=gBranch[n].R;
x=gBranch[n].X;
d=r*r+x*x;
g=r/d;
b=-x/d;
if(gBranch[n].k==0)
{
if(i!
=j)
{deltap[i][j]=(ge[i]*(ge[i]-ge[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;
}
elsebreak;
}
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(i=0;ifor(j=i+1;j{
deltalinep[i][j]=deltap[i][j]+deltap[j][i];
deltalineq[i][j]=deltaq[i][j]+deltaq[j][i];
printf("deltalinep[%d][%d]=%10.7f,deltalineq[%d][%d]=%10.7f\n",i+1,j+1,deltalinep[i][j],i+1,j+1,deltalineq[i][j]);
}
s1=0;s2=0;
for(i=0;ifor(j=0;j{
s1=s1+deltalinep[i][j];
s2=s2+deltalineq[i][j];
}
printf("totallineP=%10.7f,totallineQ=%10.7f\n",s1,s2);
if((fopen_s(&fp,"data\\busvoltageandzhurupower.txt","w"))!
=0)
{
printf("Cannotopenthefilenamed'busvoltageandzhurupower.txt'\n");
return0;
}
fprintf(fp,"---busvoltageandzhurupower---\n");
for(i=0;i{
fprintf_s(fp,"zhurup[%d]=%10.7f,zhuruq[%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;i{
fprintf_s(fp,"bus(%d).voltage=%10.5f,,angle=%10.5f\n",i+1,sqrt(ge[i]*ge[i]+gf[i]*gf[i]),atan(gf[i]/ge[i]));
}
fprintf_s(fp,"\n");
for(i=0;ifor(j=i+1;jfprintf_s(fp,"deltalinep[%d][%d]=%10.5f,deltalineq[%d][%d]=%10.5f\n",i+1,j+1,deltalinep[i][j],i+1,j+1,deltalineq[i][j]);
fprintf_s(fp,"\n");
fclose(fp);
return(0);
}
4,4
1,1.00,0.00,0.00,0.00,0.30,0.18,0
2,1.00,0.00,0.00,0.00,0.55,0.13,0
3,1.10,0.00,0.50,0.00,0.00,0.00,1
4,1.05,0.00,0.00,0.00,0.00,0.00,2
1,1,2,0.10,0.41,0.01528,0.0
2,1,3,0.00,0.30,0.00000,1.1
3,1,4,0.12,