电力系统通用潮流计算c语言程序

《电力系统通用潮流计算c语言程序》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、#include#include#include#includeusingnamespacestd;//节点号类型负荷有功负荷无功母线数据（类型1＝PV节点，2＝PQ节点，3＝平衡节点）structBUS{intbusno;inttype;floatPd;floatQd;};//发电机数据节点号有功发电电压幅值structGenerator{intbusno;floatPg;floatVg;};//支路信息节点I节点JRXB/2kstructLine{intbusi;intbusj;floatR;floatX;floa

2、tB;floatk;};//deltaPdeltaQdeltaV^2//voidfun1(doubleYG[][50],doubleYB[][50],doublee[],doublef[],inttype[],intN,doubleW[],doubleP[],doubleQ[],doubleV[]){doubledP=0,dQ=0,dV=0;inti,j;for(i=0;i

3、]*V[i]-e[i]*e[i]-f[i]*f[i];dP=P[i]-e[i]*A-f[i]*B;W[2*i]=dP;dQ=Q[i]-f[i]*A+e[i]*B;if(type[i]==1)W[2*i+1]=dQ;elseW[2*i+1]=dV;}}//Jacobi矩阵//voidJacobi(doubleYG[][50],doubleYB[][50],doublee[50],doublef[50],inttype[50],intN,doubleJa[100][101]){inti,j;for(i=0;i

4、]==1){Ja[2*i][2*j]=-(YG[i][j]*e[i]+YB[i][j]*f[i]);Ja[2*i][2*j+1]=YB[i][j]*e[i]-YG[i][j]*f[i];Ja[2*i+1][2*j]=Ja[2*i][2*j+1];Ja[2*i+1][2*j+1]=-Ja[2*i][2*j];}else{Ja[2*i][2*j]=-YG[i][j]*e[i]+YB[i][j]*f[i];Ja[2*i][2*j+1]=YB[i][j]*e[i]-YG[i][j]*f[i];Ja[2*i+1][2*j+1]=Ja[2*i+1][2*j]=0;}}else{doublea[50]=

5、{0},b[50]={0};for(intk=0;k

6、[i]+YB[i][i]*f[i];}else{Ja[2*i+1][2*j]=-2*e[i];Ja[2*i+1][2*j+1]=-2*f[i];}}}}}}//高斯消元法解方程组函数//voidgauss(doublea[][101],intn){inti,j,k;doublec;for(k=0;k

7、k>=0;k--)for(j=k+1;j