电力系统潮流计算C语言程序.doc

#include #include #define N 10/#define n 5 void main()double t; /中间变量int n; /统计节点个数int i,j,k; /循环变量int r; /PV节点个数 /int r=0; int m; /PQ节点个数double RNN;/节点间的电阻double XNN;/节点间的电抗double g1NN=0; /节点间的电导double b1NN=0; /节点间的电纳double GNN=0;double BNN=0;double RRNN=0;double XXNN=0;double KKNN=0;double gg1NN=0; /节点间对地的电导double bb1NN=0; /节点间对地的电纳double PN=0; /节点有功初值，平衡节点除外double QN=0; /节点无功初值，平衡节点除外double B1N-1N-1; /系数矩阵1double B2N-1N-1; /系数矩阵2double A1N-1N-1; /系数矩阵1的逆矩阵double A2N-1N-1; /系数矩阵2的逆矩阵double aN2*N; /求逆矩阵时用的矩阵double w=0,s=0,l,m1; /求逆矩阵时用的变量double UN; /节点电压幅值double radN; /节点电压幅角double cpN; /有功不平衡量double cqN; /无功不平衡量double ec; /容许误差cout 请将平衡节点设为零，接下来设PQ节点，最后是PV节点。; cout n ;cout 请输入PV节点个数(rn): r;m=n-r-1;cout PQ节点个数(n-r-1)为： m= m endl;cout 请输入节点间电阻：n; for( i=0 ; in ; i+ ) for( j=0 ; j Rj; cout n请输入节点间电抗：n; for( i=0 ; in ; i+ ) for( j=0 ; j Xj; cout n请输入节点对地容抗：n; for( i=0 ; in ; i+ ) for( j=0 ; j XXj; cout n请输入节点间变比：n; for( i=0 ; in ; i+ ) for( j=0 ; jKKj; for( i=0 ; in ; i+ ) for( j=0 ; jn ; j+ ) if(KKj!=1) t=KKj*Rj;Rj=Rj=t; t=KKj*Xj;Xj=Xj=t; RRj=KKj*KKj*Rj/(1-KKj); XXj=KKj*KKj*Xj/(1-KKj); RRj=KKj*Rj/(KKj-1); XXj=KKj*Xj/(KKj-1); for( i=0 ; in ; i+ ) for( j=0 ; jn ; j+ ) if( (Rj!=0)|(Xj!=0) ) g1j=Rj/(Rj*Rj+Xj*Xj); b1j=-Xj/(Rj*Rj+Xj*Xj); gg1j=RRj/(RRj*RRj+XXj*XXj); bb1j=-XXj/(RRj*RRj+XXj*XXj); else g1j=0; b1j=0; for( i=0 ; in ; i+ ) for( j=0 ; jn ; j+ ) if( (RRj!=0)|(XXj!=0) ) gg1j=RRj/(RRj*RRj+XXj*XXj); bb1j=-XXj/(RRj*RRj+XXj*XXj); else gg1j=0; bb1j=0; for( i=0 ; in ; i+ ) for( j=0 ; jn ; j+ ) if(j!=i) Gj=-g1j; else for( k=0 ; kn ; k+ ) Gj+=(g1k+gg1k); for( i=0 ; in ; i+ ) for( j=0 ; jn ; j+ ) if(j!=i) Bj=-b1j; else for( k=0 ; kn ; k+ ) Bj+=b1k+bb1k; for( i=0 ; in-1 ; i+ ) for( j=0 ; jn-1 ; j+ ) B1j=Bi+1j+1; for( i=0 ; in-r-1 ; i+ ) for( j=0 ; jn-r-1 ; j+ ) B2j=Bi+1j+1; int n1=n-1;for (i=0;in1;i+) for (j=0;jn1;j+) aj=B1j; for (i=0;in1;i+) for (j=n1;j2*n1;j+) if(i+n1)=j) aj=1; else aj=0; for (i=0;in1;i+) w=a; for(j=i;j2*n1;j+) aj=aj/w; for(k=i+1;kn1;k+) s=ak; for(j=i;j=0;i-) l=a; for(j=i;j=0;k-) m1=ak; for(j=i;j2*n1;j+) akj=akj-m1*aj; for(i=0;in1;i+) for(j=n1;j2*n1;j+) A1j-n1=aj ; int n2=n-r-1;for (i=0;in2;i+) for (j=0;jn2;j+) aj=B2j; for (i=0;in2;i+) for (j=n1;j2*n2;j+) if(i+n2)=j) aj=1; else aj=0; for (i=0;in2;i+) w=a; for(j=i;j2*n2;j+) aj=aj/w; for(k=i+1;kn2;k+) s=ak; for(j=i;j=0;i-) l=a; for(j=i;j=0;k-) m1=ak; for(j=i;j2*n2;j+) akj=akj-m1*aj; for(i=0;in2;i+) for(j=n1;j2*n2;j+) A2j-n2=aj ; cout U0; /1.05cout 请依次输入其余 n-1 个节点电压幅值初值： endl;for( i=1 ; i U; cout 请输入 n 个节点电压相位角初值: endl;for( i=0 ; i rad;cout 请输入 n-1 个节点注入有功初值：(平衡节点除外）n;for( i=1 ; i P; cout 请输入 m 个节点注入无功初值：(平衡节点除外）n;for( i=1 ; i Q; cout ec; / 0.00001 for( k=0 ; k=100 ; k+ )double cp1=0,cq1=0; /中间变量int kp=1,kq=1; for( i=1 ; in ; i+ ) double s1=0; for( j=0 ; jn ; j+ ) s1+=Uj*(Gj*cos(rad-radj)+Bj*sin(rad-radj); cp=P-U*s1; if(cp1abs(cp)cp1=cp; if(cp1ec)kp=0;if(kq=0)break; else double a1NN; /系数矩阵的逆 double bN; /常数 double xN; double c_radN; for( i=0 ; in-1 ; i+ ) for( j=0 ; jn-1 ; j+ ) a1j=A1j; for( i=0 ; in-1 ; i+ ) b=cpi+1/Ui+1; for( i=0 ; in-1 ; i+ ) double sum=0; for( j=0 ; jn-1 ; j+ ) sum+=a1j*bj; x=sum; for( i=0; in-1 ; i+ ) c_rad=-x/Ui+1; for( i=0 ; in-1 ; i+ ) radi+1=radi+1+c_rad; kq=1; for( i=1 ; i=m ; i+ ) double s1=0; for( j=0 ; jn ; j+ ) s1+=Uj*(Gj*sin(rad-radj)-Bj*cos(rad-radj); cq=Q-U*s1; if(cq1abs(cq)cq1=cq; if(cq1ec) kq=0; if(kp=0)break; else double a1NN; /系数矩阵 double bN; /常数 double xN; double CUN; for( i=0 ; im ; i+ ) for( j=0 ; jm ; j+ ) a1j=A2j; for( i=0 ; im ; i+ ) b=cqi+1/Ui+1; for( i=0 ; in-1 ; i+ ) double