华电潮流上机计算实验报告.docVIP

《电力系统潮流上机》课程设计报告 院 系：电气与电子工程学院 班 级： 电气1205 学 号： 1121181223 学生姓名： 王城 指导教师： 孙英云 设计周数： 两周 成 绩： 日期：2015年7月7日 一、课程设计的目的与要求 培养学生的电力系统潮流计算机编程能力，掌握计算机潮流计算的相关知识 二、设计正文（详细内容见附录） 手算 计算机计算 思考题 三、课程设计总结或结论（详细内容见附录） 四、参考文献 《电力系统计算：电子数字计算机的应用》，西安交通大学等合编。北京：水利电力出版社； 《现代电力系统分析》，王锡凡主编，科学出版社； 《电力系统稳态分析》，陈珩，中国电力出版社，1995年，第三版； 附录（设计流程图、程序、表格、数据等） 机算潮流程序及结果 // dierti.cpp : 定义控制台应用程序的入口点。 // #include stdafx.h struct 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 #includestdlib.h #define NBUS 4 #define NLINE 4 /* Global variables */ int nL,nB,nVA,nSH; float X[NBUS]; int L; double def[2*NBUS]; double mn[50]; void Gauss(double a[50][50],double b[50], int n) /*定义高斯法 */ { int JS[50]; int i,j,k; float d,t,x[50]; FILE *fp; int L=1; for(i=0;i50;i++) JS[i]=0; for(k=0;kn;k++){ d=0.0; for(j=k;jn;j++) if(fabs(a[k][j])d){ /*在一行中找到一个最大值赋值d,并用JS[K]记住这个最大值所在的列号*/ d=fabs(a[k][j]); JS[k]=j; } if(fabs(d)0.000001) /*如果d的数值太小，做为被除数将带来很大的误差 */ L=0; else { if(JS[k]!=k) for(i=0;in;i++) { t=a[i][k]; a[i][k]=a[i][JS[k]]; /*进行列交换，让最大值始终在对角元上*/ a[i][JS[k]]=t; } } if(L==0)break; for(j=k+1;jn;j++)a[k][j]=a[k][j]/a[k][k]; /*对角元上的元素消为1*/ b[k]=b[k]/a[k][k]; for(i=k+1;in;i++){ for(j=k+1;jn;j++)a[i][j]=a[i][j]-a[i][k]*a[k][j]; /*使下三角阵的元素为0*/ b[i]=b[i]-a[i][k]*b[k]; } } if(fabs(a[n-1][n-1])0.00001){ /*用追赶法，解方程组，求未知数x*/ x[n-1]=b[n-1]; for(i=n-2;i=0;i--){ t=0.0; for(j=i+1;jn;j++)t=t+a[i][j]*x[