華北電力大學(xué)潮流上機課程設(shè)計報告

1、 課程設(shè)計報告( 20142015年度第一學(xué)期)名 稱：電力系統(tǒng)潮流上機 院 系：電氣與電子工程學(xué)院 班 級： 電氣1215 學(xué) 號： 1121960726 學(xué)生姓名： 指導(dǎo)教師： 王莉麗 設(shè)計周數(shù)： 兩周 成 績： 日期： 2015年1月4日 課程設(shè)計報告一、課程設(shè)計的目的與要求培養(yǎng)學(xué)生的電力系統(tǒng)潮流計算機編程能力，掌握計算機潮流計算的相關(guān)知識二、設(shè)計正文（詳細(xì)內(nèi)容見附錄）1 手算2 計算機計算3 思考題三、課程設(shè)計總結(jié)或結(jié)論四、參考文獻(xiàn)1. 電力系統(tǒng)分析基礎(chǔ)，李庚銀，機械工業(yè)出版社，2011年，第一版；2. 電力系統(tǒng)穩(wěn)態(tài)分析，陳珩，中國電力出版社，2007年，第三版；附錄（包括：1.手算過

2、程及結(jié)果；2.計算機計算流程圖、表格、數(shù)據(jù)；3.思考題答案）附錄1. 手算過程及結(jié)果2. 計算機計算流程圖3.完整的潮流程序頭文件（定義變量）：#define Bus_Num 9#define Line_Num 9#define Precision 1e-5struct Busint No ;float Voltage,Phase,GenP,GenQ,LoadP,LoadQ;int Type;gBusBus_Num;struct Lineint No,No_I,No_J;float R,X,B,k;gLineLine_Num;float gY_GBus_NumBus_Num,gY_BBus_N

3、umBus_Num;float gDelta_PBus_Num-1,gDelta_QBus_Num-1,gDelta_PQ2*(Bus_Num-1);float gJaccobi2*(Bus_Num-1)2*(Bus_Num-1);float gDelta_fBus_Num-1,gDelta_eBus_Num-1,gDelta_fe2*(Bus_Num-1);float gfBus_Num,geBus_Num;主程序：/ flow.cpp: 主項目文件#include "stdafx.h"#include "NEquation.h"#include &q

4、uot;math.h"#include "stdio.h"#include "config.h"using namespace System;void test()NEquation ob1;ob1.SetSize(2);ob1.Data(0,0)=1;ob1.Data(0,1)=2;ob1.Data(1,0)=2;ob1.Data(1,1)=1;ob1.Value(0)=4;ob1.Value(1)=6;ob1.Run();printf("x1=%fn",ob1.Value(0);printf("x2=%fn&q

5、uot;,ob1.Value(1);void GetData()/Read the data FILE *fp;int i;fp=fopen("E:1121960726flowdatadata.txt","r");if(fp=NULL) printf("Can not open the file named 'data.txt' n");return;for(i=0;i<=Bus_Num-1;i+)fscanf(fp,"%d,%f,%f,%f,%f,%f,%f,%d",&gBusi.N

6、o,&gBusi.Voltage,&gBusi.Phase,&gBusi.GenP,&gBusi.GenQ,&gBusi.LoadP,&gBusi.LoadQ,&gBusi.Type);for(i=0;i<=Line_Num-1;i+)fscanf(fp,"%d,%d,%d,%f,%f,%f,%f",&gLinei.No,&gLinei.No_I,&gLinei.No_J,&gLinei.R,&gLinei.X,&gLinei.B,&gLinei.k);fc

7、lose(fp);void GetYMatrix() int i,j,bus1,bus2; float r,x,d,g,b,g1,b1,g2,b2,g3,b3; FILE *fp; for(i=0;i<=Bus_Num-1;i+) for(j=0;j<=Bus_Num-1;j+) gY_Gij=0; gY_Bij=0; for(i=0; i<=Line_Num-1; i+) if(gLinei.k=0) bus1=gLinei.No_I-1; bus2=gLinei.No_J-1; r=gLinei.R; x=gLinei.X; d=r*r+x*x; g=r/d; b=-x/

8、d; gY_Gbus1bus1=gY_Gbus1bus1+g; gY_Gbus2bus2=gY_Gbus2bus2+g;gY_Gbus1bus2=gY_Gbus1bus2-g; gY_Gbus2bus1=gY_Gbus2bus1-g; gY_Bbus1bus1=gY_Bbus1bus1+b+gLinei.B; gY_Bbus2bus2=gY_Bbus2bus2+b+gLinei.B; gY_Bbus1bus2=gY_Bbus1bus2-b; gY_Bbus2bus1=gY_Bbus2bus1-b; else bus1=gLinei.No_I-1; bus2=gLinei.No_J-1; x=g

9、Linei.X; d=r*r+x*x; g=r/d; b=-x/d; g1=g/gLinei.k; b1=b/gLinei.k; g2=g*(1-gLinei.k)/(gLinei.k*gLinei.k); b2=b*(1-gLinei.k)/(gLinei.k*gLinei.k); g3=g*(gLinei.k-1)/gLinei.k; b3=b*(gLinei.k-1)/gLinei.k; gY_Gbus1bus1=gY_Gbus1bus1+g1+g2; gY_Gbus2bus2=gY_Gbus2bus2+g1+g3; gY_Gbus1bus2=gY_Gbus1bus2-g1; gY_Gb

10、us2bus1=gY_Gbus2bus1-g1; gY_Bbus1bus1=gY_Bbus1bus1+b1+b2; gY_Bbus2bus2=gY_Bbus2bus2+b1+b3; gY_Bbus1bus2=gY_Bbus1bus2-b1; gY_Bbus2bus1=gY_Bbus2bus1-b1; / output the Y matrix fp=fopen("E:1121960726flowdataymatrix.txt","w");if(fp=NULL)printf("Can not open the file named 'ym

11、atrix.txt' n");return ;fprintf(fp,"-Y Matrix-n");for(i=0;i<=Bus_Num-1;i+)for(j=0;j<=Bus_Num-1;j+)fprintf(fp,"Y(%d,%d)=(%10.5f,%10.5f)n",i+1,j+1,gY_Gij,gY_Bij);fclose(fp);void SetInitial()int i;for(i=0;i<=Bus_Num-1;i+)if(gBusi.Type=3)gfi=gBusi.Voltage*sin(gBusi.P

12、hase);gei=gBusi.Voltage*cos(gBusi.Phase);elsegfi=0;gei=1;void GetUnbalance() int i,j; FILE *fp; for(i=0;i<=Bus_Num-2;i+) gDelta_Pi=gBusi+1.GenP-gBusi+1.LoadP; if(gBusi+1.Type=2) /PV節(jié)¨²點Ì? gDelta_Qi=gBusi+1.Voltage*gBusi+1.Voltage-(gei+1*gei+1+gfi+1*gfi+1); else gDelta_Qi=gBusi+1.Ge

13、nQ-gBusi+1.LoadQ; for(j=0;j<=Bus_Num-1;j+) gDelta_Pi=gDelta_Pi-gei+1*(gY_Gi+1j*gej-gY_Bi+1j*gfj)-gfi+1*(gY_Gi+1j*gfj+gY_Bi+1j*gej); if(gBusi+1.Type=1) /PQ節(jié)¨²點Ì? gDelta_Qi=gDelta_Qi-gfi+1*(gY_Gi+1j*gej-gY_Bi+1j*gfj)+gei+1*(gY_Gi+1j*gfj+gY_Bi+1j*gej); for(i=0;i<=Bus_Num-2;i+) /合?并

14、¡é gDelta_PQ2*i=gDelta_Pi; gDelta_PQ2*i+1=gDelta_Qi; fp=fopen("E:1121960726flowdataunbalance.txt","w"); if(fp=NULL) printf("無法打開文件：”'unbalance.txt' n"); return ; fprintf(fp,"-Unbalance-n"); for(i=0;i<=2*Bus_Num-3;i+) fprintf(fp,"Unbal

15、ance%d=%10.5fn",i+1,gDelta_PQi); fclose(fp);void GetJaccobi() int i,j; float gaBus_Num-1,gbBus_Num-1; FILE *fp; for(i=0;i<=Bus_Num-2;i+) /計算注入電流 gai=0; gbi=0; for(j=0;j<=Bus_Num-1;j+) gai=gai+gY_Gi+1j*gej-gY_Bi+1j*gfj; gbi=gbi+gY_Gi+1j*gfj+gY_Bi+1j*gej; for(i=0;i<=Bus_Num-2;i+) for(j=0

16、;j<=Bus_Num-2;j+) if(i!=j) gJaccobi2*i2*j=-gY_Bi+1j+1*gei+1+gY_Gi+1j+1*gfi+1; gJaccobi2*i2*j+1=gY_Gi+1j+1*gei+1+gY_Bi+1j+1*gfi+1;if(gBusi+1.Type=2) /PV節(jié)¨²點Ì? gJaccobi2*i+12*j=0; gJaccobi2*i+12*j+1=0; else /PQ gJaccobi2*i+12*j=-gJaccobi2*i2*j+1; gJaccobi2*i+12*j+1=gJaccobi2*i2*j; el

17、se gJaccobi2*i2*j=-gY_Bi+1j+1*gei+1+gY_Gi+1j+1*gfi+1+gbi;gJaccobi2*i2*j+1=gY_Gi+1j+1*gei+1+gY_Bi+1j+1*gfi+1+gai;if(gBusi+1.Type=2) /PV節(jié)¨²點Ì? gJaccobi2*i+12*j=2*gfi+1; gJaccobi2*i+12*j+1=2*gei+1; else /PQ節(jié)點 gJaccobi2*i+12*j=-gY_Gi+1j+1*gei+1-gY_Bi+1j+1*gfi+1+gai; gJaccobi2*i+12*j+1=-gY

18、_Bi+1j+1*gei+1+gY_Gi+1j+1*gfi+1-gbi; fp=fopen("E:1121180825flowdatajaccobi.txt","w"); if(fp=NULL) printf("無法打開文件：” 'jaccobi.txt' n"); return ; fprintf(fp,"-Jaccobi Matrix-n"); for(i=0;i<=2*Bus_Num-3;i+) for(j=0;j<=2*Bus_Num-3;j+) fprintf(fp,"

19、;jaccobi(%d,%d)=%10.5fn",i+1,j+1,gJaccobiij); fclose(fp);void GetRevised() int i,j; FILE *fp; NEquation ob1; /解矩陣方程 ob1.SetSize(2*(Bus_Num-1); for(i=0;i<=2*Bus_Num-3;i+) for(j=0;j<=2*Bus_Num-3;j+) ob1.Data(i,j)=gJaccobiij; for(i=0;i<=2*Bus_Num-3;i+) ob1.Value(i)=gDelta_PQi; ob1.Run();

20、for(i=0;i<=Bus_Num-2;i+) gDelta_fi=ob1.Value(2*i); gDelta_ei=ob1.Value(2*i+1); gDelta_fe2*i=gDelta_fi; gDelta_fe2*i+1=gDelta_ei; fp=fopen("E:1121960726flowdatarevised.txt","w"); if(fp=NULL) printf("無法打開文件：” 'revised.txt' n"); return ; fprintf(fp,"-Revise

21、d-n"); for(i=0;i<=2*Bus_Num-3;i+) fprintf(fp,"revised%d=%10.5fn",i+1,gDelta_fei); fclose(fp);void GetNewValue() int i; FILE *fp; for(i=0;i<=Bus_Num-2;i+) gfi+1=gfi+1+gDelta_fi; gei+1=gei+1+gDelta_ei; fp=fopen("E:1121960726flowdatanewvalue.txt","w"); if(fp=NULL) printf("無法打開文件：”'newvalue.txt' n"); return ; fprintf(fp,"-New Value-n"); for(i=0;i<=Bus_Num-2;i+) fprintf(fp,"f(%d)=%10.5f,e(%d)=%10.5fn",