潮流计算实训心得体会.docxVIP

潮流计算实训心得体会 　　篇一：华电潮流上机计算实验报告 　　《电力系统潮流上机》课程设计报告  　　院 系：电气与电子工程学院  　　班 级： 电气1205  　　学 号：  　　学生姓名： 王城  　　指导教师：孙英云  　　设计周数：  　　成 绩：  　　日期：XX年7月7 　　日 　　一、课程设计的目的与要求  　　培养学生的电力系统潮流计算机编程能力，掌握计算机潮流计算的相关知识  　　二、设计正文（详细内容见附录）  　　1． 手算  　　2． 计算机计算  　　3． 思考题  　　三、课程设计总结或结论（详细内容见附录）  　　四、参考文献  　　1. 《电力系统计算：电子数字计算机的应用》，西安交通大学等合编。北京：水利电力出版社；  　　2. 《现代电力系统分析》，王锡凡主编，科学出版社；  　　3. 《电力系统稳态分析》，陈珩，中国电力出版社，1995年，第三版；  　　附录（设计流程图、程序、表格、数据等） 　　4. 机算潮流程序及结果 　　// : 定义控制台应用程序的入口点。  　　//  　　#include  　　struct Line //线路结构体  　　{  　　int Num,NumI,NumJ;//线路号 左节点名右节点名  　　float R,X,B,K; //电阻 电抗 电纳  　　};  　　struct Bus //节点结构体  　　{  　　int Num ;  　　float Volt,Phase,GenP,GenQ,LoadP,LoadQ;  　　int Type;  　　};  　　#include  　　#include  　　#include  　　#include  　　#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) /*变比（K等于1为普通支路，定义高斯法 */  　　不等于1为变压器支路的变比）  　　{  　　int JS[50];  　　int i,j,k;  　　float d,t,x[50];  　　FILE *fp;  　　int L=1;  　　for(i=0;i　　for(k=0;k　　d=;  　　for(j=k;j　　if(fabs(a[k][j])d){/*在一行中找到一个最大值赋值d,并用JS[K]记住这个最大值所在的列号*/  　　d=fabs(a[k][j]);  　　JS[k]=j;  　　}  　　if(fabs(d)　　L=0;  　　else {  　　if(JS[k]!=k)  　　for(i=0;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;j　　b[k]=b[k]/a[k][k];  　　for(i=k+1;i　　for(j=k+1;j　　b[i]=b[i]-a[i][k]*b[k];  　　}  　　}  　　if(fabs(a[n-1][n-1])){/*用追赶法，解方程组，求未知数x*/ x[n-1]=b[n-1];  　　for(i=n-2;i=0;i--){  　　t=;  　　for(j=i+1;j　　x[i]=(b[i]-t);  　　}  　　}  　　if((fp=fopen(,w))==NULL) /*将结果写到TXT文件中*/  　　{printf(err);exit(0);}  　　for(i=0;i　　fprintf(fp,%lf,x[i]);  　　mn[i]=x[i];  　　fprintf(fp,\n);}  　　fclose(fp);  　　if(fp!=NULL) fclose(fp);  　　}  　　int _tmain(int argc, _TCHAR* argv[])  　　{  　　FILE *fp;  　　FILE *fpout;  　　int i,j,k,l,h,n,v;  　　int i1,i