《华电潮流上机课程设计报告(程序最终无误版)》由会员分享，可在线阅读，更多相关《华电潮流上机课程设计报告(程序最终无误版)（19页珍藏版）》请在金锄头文库上搜索。
1、 电力系统潮流上机课程设计报告 院 系：电气与电子工程学院 班 级： 电气1108 学 号： 1111550112 学生姓名： 龙日尚 指导教师： 刘宝柱 设计周数： 两周 成 绩： 日期：2014年1月10日 课程 课程设计报告一、课程设计的目的与要求培养学生的电力系统潮流计算机编程能力，掌握计算机潮流计算的相关知识二、设计正文（详细内容见附录）1 手算2 计算机计算3 思考题三、课程设计总结或结论（详细内容见附录）四、参考文献1. 电力系统计算：电子数字计算机的应用，西安交通大学等合编。北京：水利电力出版社；2. 现代电力系统分析，王锡凡主编，科学出版社；3. 电力系统稳态分析，陈珩，中国
2、电力出版社，1995年，第三版；附录（设计流程图、程序、表格、数据等）4. 机算潮流程序及结果/ dierti.cpp : 定义控制台应用程序的入口点。/#include stdafx.hstruct Line /线路结构体 int Num,NumI,NumJ; /线路号 左节点名 右节点名 float R,X,B,K; /电阻 电抗 电纳 变比（K等于1为普通支路， 不等于1为变压器支路的变比） ;struct Bus /节点结构体 int Num ; float Volt,Phase,GenP,GenQ,LoadP,LoadQ; int Type; ;#includestdio.h#inc
3、ludestring.h#includemath.h#includestdlib.h#define NBUS 4#define NLINE 4/* Global variables */int nL,nB,nVA,nSH;float XNBUS;int L;double def2*NBUS;double mn50;void Gauss(double a5050,double b50, int n) /*定义高斯法 */ int JS50; int i,j,k; float d,t,x50; FILE *fp; int L=1; for(i=0;i50;i+) JSi=0; for(k=0;kn
4、;k+) d=0.0; for(j=k;jd) /*在一行中找到一个最大值赋值d,并用JSK记住这个最大值所在的列号*/ d=fabs(akj); JSk=j; if(fabs(d)0.000001) /*如果d的数值太小，做为被除数将带来很大的误差 */ L=0; else if(JSk!=k)for(i=0;in;i+) t=aik; aik=aiJSk; /*进行列交换，让最大值始终在对角元上*/ aiJSk=t; if(L=0)break; for(j=k+1;jn;j+)akj=akj/akk; /*对角元上的元素消为1*/ bk=bk/akk; for(i=k+1;in;i+) f
5、or(j=k+1;j0.00001) /*用追赶法，解方程组，求未知数x*/ xn-1=bn-1; for(i=n-2;i=0;i-) t=0.0; for(j=i+1;jn;j+)t=t+aij*xj; xi=(bi-t); if(fp=fopen(gauss.txt,w)=NULL) /*将结果写到TXT文件中*/ printf(err);exit(0); for(i=0;in;i+) fprintf(fp,%lf,xi); mni=xi; fprintf(fp,n); fclose(fp); if(fp!=NULL) fclose(fp);int _tmain(int argc, _TC
6、HAR* argv)FILE *fp;FILE *fpout;int i,j,k,l,h,n,v;int i1,i2,i3,kp,kq;float d1,d2,d3,d4,d5,d6,r,x,g,b,tt,LL,e,ps,qs,shsh,m;struct Line sLNLINE;struct Bus sBNBUS;float YGNBUS+1NBUS+1,YBNBUS+1NBUS+1;double u502;i1=i2=i3=0;d1=d2=d3=d4=d5=d6=ps=qs=0.0;for(i=0;iNBUS;i+)if(fp=fopen(in.txt,r)=NULL) printf(Ca
7、n not open the file named in.txt n);exit(0); fscanf(fp,%d,%d,%d,&nB,&nL,&nSH);for(i=0;inB;i+)sBi.Num=sBi.Type=0;sBi.Volt=1.0;sBi.Phase=sBi.GenP=sBi.GenQ=sBi.LoadP=sBi.LoadQ=0.0; fscanf(fp,%d,%f,%f,%f,%f,%f,%f,%d,&i1,&d1,&d2,&d3,&d4,&d5,&d6,&i2);sBi.Num=i1;sBi.Volt=d1;sBi.Phase=d2;sBi.GenP=d3;sBi.Gen
8、Q=d4;sBi.LoadP=d5,sBi.LoadQ=d6;sBi.Type=i2; ;for(i=0;inL;i+)sLi.Num=sLi.NumI=sLi.NumJ=0;sLi.R=sLi.X=sLi.B=0.0;sLi.K=1.0;fscanf(fp,%2d %3d %3d %f %f %f %f,&i1,&i2,&i3,&d1,&d2,&d3,&d4); sLi.Num=i1;sLi.NumI=i2;sLi.NumJ=i3;sLi.R=d1;sLi.X=d2;sLi.B=d3;sLi.K=d4; if(fp!=NULL) fclose(fp);/*Make Y Matrix*/for
9、(i=1;inB+1;i+)for(j=1;jnB+1;j+) YGij=0.0; YBij=0.0; ;for(l=0; lnL; l+)i=sLl.NumI; j=sLl.NumJ; r=sLl.R; x=sLl.X; d1=r*r+x*x; g=r/d1; b=-x/d1; m=sLl.K; if(fabs(sLl.K-1.0)0.000001) /普通支路 YGii=YGii+g; YGjj=YGjj+g; YBii=YBii+b+sLl.B; YBjj=YBjj+b+sLl.B; YGij=YGij-g; YGji=YGji-g; YBij=YBij-b; YBji=YBji-b;
10、else /变压器支路YGii=YGii+g/m+g*(m-1)/m; YGjj=YGjj+g/m+g*(1-m)/m/m; YBii=YBii+b/m+b*(m-1)/m; YBjj=YBjj+b/m+b*(1-m)/m/m; YGij=YGij-g/m; YGji=YGji-g/m; YBij=YBij-b/m; YBji=YBji-b/m; /* Check the Y matrix */if(fp=fopen(GGBB.txt,w)=NULL)printf(Can not open the file named GGBB.txt n);exit(0);fprintf(fp,-Y Matrix-n);for(i=1;inB+1;i+)for(j=1;j0.000001)fprintf(fp,Y(%3d,%-3d)=(%10.5f,%10.5f)n,i,j,YGij,YBij);if(fp!=NULL) fclose(fp);/* 节点电压附初值 */for(i=1;inB+1;i+)if(sBi-1.Type=0)ui0=0.0;ui1=1.0;else if(sBi-1.Type=1)ui1=sBi-1.Volt;ui0=0.0;else if(sBi-1.Type=2)ui1=sBi-1.Volt;ui0= sBi-1.Phase;for(v
