《华电潮流上机课程设计报告(程序最终无误版)》由会员分享，可在线阅读，更多相关《华电潮流上机课程设计报告(程序最终无误版)（16页珍藏版）》请在金锄头文库上搜索。
1、.电力系统潮流上机课程设计报告院 系：电气与电子工程学院班 级：电气1108 学 号：1111550112 学生_ 龙日尚 指导刘宝柱设计周数：两周成 绩：日期：20XX1月10日.一、课程设计的目的与要求培养学生的电力系统潮流计算机编程能力,掌握计算机潮流计算的相关知识二、设计正文详细内容见附录1 手算2 计算机计算3 思考题三、课程设计总结或结论详细内容见附录四、参考文献1. 电力系统计算：电子数字计算机的应用,XX交通大学等合编。北京：水利电力出版社；2. 现代电力系统分析,王锡凡主编,科学出版社；3. 电力系统稳态分析,陈珩,中国电力出版社,1995年,第三版；附录设计流程图、程序、表
2、格、数据等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#includestring.h#includemath.h#inclu
3、destdlib.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 /*定义高斯法 */ int JS50; int i,j,k; float d,t,x50; FILE *fp; int L=1; fori=0;i JSi=0; fork=0;k d=0.0; forj=k;j iffabsd /*在一行中找到一个最大值赋值d,并用JSK记住这个最大值所在的列号*/ d=fabs; JSk
4、=j; iffabs /*如果d的数值太小,做为被除数将带来很大的误差 */ L=0; else iffori=0;i t=aik; aik=aiJSk; /*进行列交换,让最大值始终在对角元上*/ aiJSk=t; ifbreak; forj=k+1;jakj=akj/akk; /*对角元上的元素消为1*/ bk=bk/akk; fori=k+1;i forj=k+1;jaij=aij-aik*akj; /*使下三角阵的元素为0*/ bi=bi-aik*bk; iffabs0.00001 /*用追赶法,解方程组,求未知数x*/ xn-1=bn-1; for=0;i- t=0.0; forj=
5、i+1;jt=t+aij*xj; xi=; iffp=fopen=NULL /*将结果写到TXT文件中*/ printf;exit; fori=0;i fprintf; mni=xi; fprintf; fclose; if fclose;int _tmainFILE *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,
6、YBNBUS+1NBUS+1;double u502;i1=i2=i3=0;d1=d2=d3=d4=d5=d6=ps=qs=0.0;fori=0;iiffp=fopen=NULL printf;exit; fscanf;fori=0;isBi.Num=sBi.Type=0;sBi.Volt=1.0;sBi.Phase=sBi.GenP=sBi.GenQ=sBi.LoadP=sBi.LoadQ=0.0; fscanf;sBi.Num=i1;sBi.Volt=d1;sBi.Phase=d2;sBi.GenP=d3;sBi.GenQ=d4;sBi.LoadP=d5,sBi.LoadQ=d6;sBi.
7、Type=i2; ;fori=0;isLi.Num=sLi.NumI=sLi.NumJ=0;sLi.R=sLi.X=sLi.B=0.0;sLi.K=1.0;fscanf; sLi.Num=i1;sLi.NumI=i2;sLi.NumJ=i3;sLi.R=d1;sLi.X=d2;sLi.B=d3;sLi.K=d4; if fclose;/*Make Y Matrix*/fori=1;iforj=1;j YGij=0.0; YBij=0.0; ;forl=0; li=sLl.NumI; j=sLl.NumJ; r=sLl.R; x=sLl.X; d1=r*r+x*x; g=r/d1; b=-x/d
8、1; m=sLl.K; iffabs /普通支路 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; else /变压器支路YGii=YGii+g/m+g*/m; YGjj=YGjj+g/m+g*/m/m; YBii=YBii+b/m+b*/m; YBjj=YBjj+b/m+b*/m/m; YGij=YGij-g/m; YGji=YGji-g/m; YBij=YBij-b/m; YBji=YBji-b/m; /* C
9、heck the Y matrix */iffp=fopen=NULLprintf;exit;fprintf;fori=1;iforj=1;jiffabs0.000001fprintffp,Y=n,i,j,YGij,YBij;if fclose;/* 节点电压附初值 */fori=1;iifui0=0.0;ui1=1.0;else ifui1=sBi-1.Volt;ui0=0.0;else ifui1=sBi-1.Volt;ui0= sBi-1.Phase;for/* 迭代次数可以无限大 */* 节点电压附初值 */printf;fori=1;iprintf;printf;printf;/* 求偏移量 */double P_P10;double P_Q10;double P_UU
