潮流上机课程设计报告[1].华电.docVIP

《电力系统潮流上机》课程设计报告 院 系：电气与电子工程学院 班 级： 电气1108班 学 号： 1111180806 学生姓名： 陈子君 指导教师： 刘宝柱 设计周数： 两周 成 绩： 日期：2014年1月10日 课程课程设计报告 一、课程设计的目的与要求 培养学生的电力系统潮流计算机编程能力，掌握计算机潮流计算的相关知识 二、设计正文（详细内容见附录） 1.手算： 要求应用牛顿-拉夫逊法或P-Q分解法手算求解，要求精度为0.001MW。 节点1为平衡节点，电压，节点2为PQ节点，负荷功率，节点3是PV节点，，两条支路分别为，，对地支路 2.计算机计算：编写潮流计算程序，要求如下： 2.1据给定的潮流计算任务书整理潮流计算的基础数据：节点的分类，线路模型，等值变压器模型，电压等级的归算，标幺值的计算； 2.2基础数据的计算机存储：节点数据，支路数据（包括变压器）； 2.3用牛顿-拉夫逊法计算； 2.4根据所选潮流计算方法画流程图，划分出功能模块，有数据输入模块，导纳阵形成模块，解线性方程组模块，计算不平衡功率模块，形成雅可比矩阵模块，解修正方程模块，计算线路潮流，网损，PV节点无功功率和平衡节点功率，数据输出模块； 2.5据上述模块编制程序并上机调试程序，得出潮流计算结果； 2.6源程序及其程序中的符号说明集、程序流图 简单系统如下图所示，支路数据如下： ，，， ，， 节点数据如下： ，， ，， // 2014.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 awe[100]; int _tmain(int argc, _TCHAR* argv[]) { FILE *fp; FILE *fpout; int i,j,l,h,n,k; int i1=0,i2=0,i3=0,kp,kq; double a,b,k1; double H=0,N=0,J=0,L=0,R=0,S=0; float d1=0,d2=0,d3=0,d4=0,d5=0,d6=0,r,x,g,b1,b2,g1,g2,tt,LL,e,ps=0,qs=0,shsh,m; struct Line sL[NLINE]; struct Bus sB[NBUS]; float YG[NBUS+1][NBUS+1],YB[NBUS+1][NBUS+1]; double u[50][2],p[50],da,db,pp[50],ppl[50],Q[10][10],gJaccobi[50][50],sP[10][10],sQ[10][10],WP=0.0,WQ=0.0; for(i=0;iNBUS;i++) {if((fp=fopen(H:\\2014\\2014\\in.txt,r))==NULL) { printf(Can not open the file named in.txt \n); exit(0); } fscanf(fp,%d,%d,%d,nB,nL,nSH);} for(i=0;inB;i++){ sB[i].Num=sB[i].Type=0;sB[i].Volt=1.0; sB[i].Phase=sB[i].GenP=sB[i].GenQ=sB[i].LoadP=sB[i].LoadQ=0.0; fscanf(fp,%d,%f,%f,%f,%f,%f,%f,%d,i1,d1,d2,d3,d4,d5,d6,i2); sB[i].Num=i