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

《华电潮流上机课程设计报告(程序最终无误版).doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、《电力系统潮流上机》课程设计报告院系：电气与电子工程学院班级：电气1108学号：1111550112学生：龙日尚指导教师：宝柱设计周数：两周成绩：日期：2014年1月10日...一、课程设计的目的与要求培养学生的电力系统潮流计算机编程能力，掌握计算机潮流计算的相关知识二、设计正文（详细容见附录）1．手算2．计算机计算3．思考题三、课程设计总结或结论（详细容见附录）四、参考文献1.《电力系统计算：电子数字计算机的应用》，交通大学等合编。北京：水利电力；2.《现代电力系统分析》，王锡凡主编，科学；3.《电力系统稳态分析》，珩，中国电力，

2、1995年，第三版；附录（设计流程图、程序、表格、数据等）4.机算潮流程序及结果...//dierti.cpp:定义控制台应用程序的入口点。//#include"stdafx.h"structLine//线路结构体{intNum,NumI,NumJ;//线路号左节点名右节点名floatR,X,B,K;//电阻电抗电纳变比（K等于1为普通支路，不等于1为变压器支路的变比）};structBus//节点结构体{intNum;floatVolt,Phase,GenP,GenQ,LoadP,LoadQ;intType;};#include"

3、stdio.h"#include"string.h"#include"math.h"#include"stdlib.h"#defineNBUS4#defineNLINE4/*Globalvariables*/intnL,nB,nVA,nSH;floatX[NBUS];intL;doubledef[2*NBUS];doublemn[50];voidGauss(doublea[50][50],doubleb[50],intn)/*定义高斯法*/{...intJS[50];inti,j,k;floatd,t,x[50];FILE*fp;i

4、ntL=1;for(i=0;i<50;i++)JS[i]=0;for(k=0;kd){/*在一行中找到一个最大值赋值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;i

5、进行列交换，让最大值始终在对角元上*/a[i][JS[k]]=t;}}if(L==0)break;for(j=k+1;j0.00001){/*用追赶

6、法，解方程组，求未知数x*/x[n-1]=b[n-1];for(i=n-2;i>=0;i--){t=0.0;for(j=i+1;j

7、ULL)fclose(fp);}int_tmain(intargc,_TCHAR*argv[]){FILE*fp;FILE*fpout;inti,j,k,l,h,n,v;inti1,i2,i3,kp,kq;floatd1,d2,d3,d4,d5,d6,r,x,g,b,tt,LL,e,ps,qs,shsh,m;structLinesL[NLINE];...structBussB[NBUS];floatYG[NBUS+1][NBUS+1],YB[NBUS+1][NBUS+1];doubleu[50][2];i1=i2=i3=0;d1=d

8、2=d3=d4=d5=d6=ps=qs=0.0;for(i=0;i