《PQ分解法潮流计算程序》.doc

#include stdio.h #include math.h #define M 10 int pq;/*PQ节点数*/ int jd;/*系统节点总数*/ int br;/*系统支路数*/ int flag1[M],flag2[M];/*节点标志*/ int N;/*迭代次数*/ float r[M],x[M],g[M][M],b[M][M],p[M],q[M],e[M],f[M],u[M];/*系统参数*/ float w[2*M-2];/*功率误差*/ float jj[2*M-2][2*M];/*雅可比矩阵*/ float a[M],c[M];/*节点注入电流*/ float du[2*M-2];/*电压增量*/ float eps,eps1;/*精度*/ float g1[M],b1[M]; main() { FILE *fp; FILE *wp; /*读取系统参数*/ int i,j,t,dd=0; int m,n; float sum,l,k; if((fp=fopen(a:/liqiang/chaoliu/true/input.txt,r))==NULL) { printf(打不开文件!\n); return; } fscanf(fp,%d,%d,%d\n,jd,br,pq); printf(%d,%d,%d\n,jd,br,pq); fscanf(fp,%f\n,eps); printf(%f\n,eps); for(i=1;i=pq;i++) { fscanf(fp,%f,%f,%f,%f\n,e[i],f[i],p[i],q[i]); /*PQ节点参数*/ printf(%f,%f,%f,%f\n,e[i],f[i],p[i],q[i]); } for(i=pq+1;ijd;i++) { fscanf(fp,%f,%f,%f\n,e[i],f[i],p[i]); /*PV节点参数*/ printf(%f,%f,%f\n,e[i],f[i],p[i]); u[i]=e[i]*e[i]+f[i]*f[i]; } fscanf(fp,%f,%f\n,e[jd],f[jd]); /*平衡节点参数*/ printf(%f,%f\n,e[jd],f[jd]); for(i=1;i=br;i++) { fscanf(fp,%d,%d,%f,%f\n,flag1[i],flag2[i],r[i],x[i]); printf(%d,%d,%f,%f\n,flag1[i],flag2[i],r[i],x[i]); } fclose(fp);wp=fopen(a:/liqiang/chaoliu/true/input1.txt,w+); for(i=1;i=jd;i++) for(j=1;j=jd;j++) g[i][j]=b[i][j]=0; for(t=1;t=br;t++) { i=flag1[t]; j=flag2[t]; sum=r[t]*r[t]+x[t]*x[t]; l=r[t]/sum; g1[t]=l; k=-x[t]/sum; b1[t]=k; g[i][i]+=l;g[j][j]+=l; b[i][i]+=k;b[j][j]+=k; g[i][j]-=l;g[j][i]-=l; b[i][j]-=k;b[j][i]-=k; fprintf(wp,%d,%d,%3.5f,%3.5f\n,flag1[t],flag2[t],g[i][j],b[i][j]); }fclose(wp); for(i=1;i=jd;i++) { for(j=1;j=jd;j++) printf(%.5f+j%.5f ,g[i][j],b[i][j]); printf(\n); } ds(); compare(); printf(eps=%.5f\n,eps1); if(eps1eps) out();N=0; do {N=N+1;printf(\n^^^^^^^第%d次迭代^^^^^^^\n,N); for(i=1;i=pq;i++) for(j=1;jjd;j++) { if(i==j) {