华电潮流上机课程设计报告程序最终无误版.docVIP

《电力系统潮流上机》课程设计报告 院 系：电气与电子工程学院 班 级： 电气1108 学 号： 1111550112 学生姓名： 龙日尚 指导教师： 刘宝柱 设计周数： 两周 成 绩： 日期：2014年1月10日一、课程设计附录（）// #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[j]; x[i]=(b[i]-t); } } if((fp=fopen(gauss.txt,w))==NULL) /*将结果写到TXT文件中*/ {printf(err);exit(0);} for(i=0;in;i++){ fprintf(fp,%lf,x[i]); mn[i]=x[i]; fprintf(fp,\n);} fclose(fp); if(fp!=NULL) fclose(fp); } in