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

《电力系统潮流上机》课程设计报告 院 系：电气与电子工程学院 班 级： 电气1205 学 号： 1121181223 学生姓名： 王城 指导教师： 孙英云 设计周数： 两周 成 绩： 日期：2015年7 月7 日 课程 课程设计报告 一、课程设计的目的与要求 培养学生的电力系统潮流计算机编程能力，掌握计算机潮流计算的相关知识 二、设计正文 （详细内容见附录） 1．手算 2．计算机计算 3．思考题 三、课程设计总结或结论 （详细内容见附录） 四、参考文献 1. 《电力系统计算：电子数字计算机的应用》，西安交通大学等合编。北京：水利电力出版社； 2. 《现代电力系统分析》，王锡凡主编，科学出版社； 3. 《电力系统稳态分析》，陈珩，中国电力出版社，1995年，第三版； 附录 （设计流程图、程序、表格、数据等） 1 课程 课程设计报告 4. 机算潮流程序及结果 // dierti.cpp : 定义控制台应用程序的入口点。 // #include stdafx.h struct Line //线路结构体 { int Num,NumI,NumJ; //线路号 左节点名 右节点名 floatR,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) /*定义高斯法 */ 2 课程 课程设计报告 { 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; }