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

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

2、机算潮流程序及结果./ dierti.cpp :/#include "stdafx.h"定义控制台应用程序的入口点。struct Line/线路结构体int Num,NumI,NumJ;/线路号左节点名右节点名floatR,X,B,K;/电阻电抗电纳变比（ K 等于1 为普通支路，不等于1 为变压器支路的变比）;struct Bus/节点结构体int Num ;float Volt,Phase,GenP,GenQ,LoadP,LoadQ;int Type;#include"stdio.h"#include"string.h"#incl

3、ude"math.h"#include"stdlib.h"#define NBUS 4#define NLINE 4/* Global variables */int nL,nB,nVA,nSH;float XNBUS;int L;double def2*NBUS;double mn50;void Gauss(double a5050,double b50, int n)./*定义高斯法*/int JS50;int i,j,k;float d,t,x50;FILE *fp;int L=1;for(i=0;i<50;i+) JSi=0;for(k=0

4、;k<n;k+)d=0.0;for(j=k;j<n;j+)if(fabs(akj)>d)/*在一行中找到一个最大值赋值d, 并用JSK 记住这个最大值所在的列号*/d=fabs(akj);JSk=j;if(fabs(d)<0.000001)/*如果d 的数值太小，做为被除数将带来很大的误差*/L=0;else if(JSk!=k)for(i=0;i<n;i+)t=aik;aik=aiJSk;/*进行列交换，让最大值始终在对角元上*/aiJSk=t;if(L=0)break;for(j=k+1;j<n;j+)akj=akj/akk;/*对角元上的元素消为1*/

5、bk=bk/akk;for(i=k+1;i<n;i+)for(j=k+1;j<n;j+)aij=aij-aik*akj;/*使下三角阵的元素为0*/.bi=bi-aik*bk;.if(fabs(an-1n-1)>0.00001)/*用追赶法，解方程组，求未知数x*/xn-1=bn-1;for(i=n-2;i>=0;i-)t=0.0;for(j=i+1;j<n;j+)t=t+aij*xj;xi=(bi-t);if(fp=fopen("gauss.txt","w")=NULL)/*将结果写到TXT文件中 */printf(&qu

6、ot;err");exit(0);for(i=0;i<n;i+)fprintf(fp,"%lf",xi);mni=xi;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,i2,i3,kp,kq;float d1,d2,d3,d4,d5,d6,r,x,g,b,tt,LL,e,ps,qs,shsh,m;struct Line sLNL

7、INE;.struct Bus sBNBUS;float YGNBUS+1NBUS+1,YBNBUS+1NBUS+1;double u502;i1=i2=i3=0;d1=d2=d3=d4=d5=d6=ps=qs=0.0;for(i=0;i<NBUS;i+)if(fp=fopen("in.txt","r")=NULL) printf("Can not open the file named 'in.txt' n");exit(0);fscanf(fp,"%d,%d,%d",&nB,&am

8、p;nL,&nSH);for(i=0;i<nB;i+)sBi.Num=sBi.Type=0;sBi.Volt=1.0;sBi.Phase=sBi.GenP=sBi.GenQ=sBi.LoadP=sBi.LoadQ=0.0;fscanf(fp,"%d,%f,%f,%f,%f,%f,%f,%d",&i1,&d1,&d2,&d3,&d4,&d5,&d6,&i2);sBi.Num=i1;sBi.Volt=d1;sBi.Phase=d2;sBi.GenP=d3;sBi.GenQ=d4;sBi.LoadP=d

9、5,sBi.LoadQ=d6;sBi.Type=i2;for(i=0;i<nL;i+)sLi.Num=sLi.NumI=sLi.NumJ=0;sLi.R=sLi.X=sLi.B=0.0;sLi.K=1.0;fscanf(fp,"%2d %3d %3d %f %f %f %f",&i1,&i2,&i3,&d1,&d2,&d3,&d4);sLi.Num=i1;sLi.NumI=i2;sLi.NumJ=i3;sLi.R=d1;sLi.X=d2;sLi.B=d3;sLi.K=d4;if(fp!=NULL) fclose(

10、fp);./*Make Y Matrix*/for(i=1;i<nB+1;i+)for(j=1;j<nB+1;j+)YGij=0.0;YBij=0.0;for(l=0; l<nL; l+)i=sLl.NumI;j=sLl.NumJ;r=sLl.R;x=sLl.X;d1=r*r+x*x;g=r/d1;b=-x/d1;m=sLl.K;if(fabs(sLl.K-1.0)<0.000001)/普通支路YGii=YGii+g;YGjj=YGjj+g;YBii=YBii+b+sLl.B;YBjj=YBjj+b+sLl.B;YGij=YGij-g;YGji=YGji-g;YBij=

11、YBij-b;YBji=YBji-b;else / YGii=YGii+g/m+g*(m-1)/m; . . .变压器支路YGjj=YGjj+g/m+g*(1-m)/m/m;YBii=YBii+b/m+b*(m-1)/m;YBjj=YBjj+b/m+b*(1-m)/m/m;YGij=YGij-g/m;YGji=YGji-g/m;YBij=YBij-b/m;YBji=YBji-b/m;/* Check the Y matrix */if(fp=fopen("GGBB.txt","w")=NULL)printf("Can not open the

12、file named 'GGBB.txt' n");exit(0);fprintf(fp,"-Y Matrix-n");for(i=1;i<nB+1;i+)for(j=1;j<nB+1;j+)if(fabs(YBij-0.0)>0.000001)fprintf(fp,"Y(%3d,%-3d)=(%10.5f,%10.5f)n",i,j,YGij,YBij);if(fp!=NULL) fclose(fp);/*节点电压附初值*/for(i=1;i<nB+1;i+)if(sBi-1.Type=0)ui0=0.

13、0;ui1=1.0;else if(sBi-1.Type=1)ui1=sBi-1.Volt;ui0=0.0;else if(sBi-1.Type=2)ui1=sBi-1.Volt;ui0= sBi-1.Phase; . .for(v=1;v+)/*迭代次数可以无限大*/*节点电压附初值*/printf("迭代第 %d次赋予的电压初值为e+jf:n",v);for(i=1;i<nB+1;i+)printf("%lf,%lfn",ui1,ui0);printf("n");printf("n");/*求偏移量 */

14、double P_P10;double P_Q10;double P_UU10;for(i=1;i<nB+1;i+)if(sBi-1.Type=2)P_Pi=0.0;P_Qi=0.0;P_UUi=1.05;if(sBi-1.Type=0)double tempP=0.0;double tempQ=0.0;for(j=1;j<nB+1;j+)tempP+=YGij*uj1-YBij*uj0;tempQ+=YGij*uj0+YBij*uj1;.P_Pi=(sBi-1.GenP-sBi-1.LoadP)-tempP*ui1-tempQ*ui0;P_Qi=(sBi-1.GenQ-sBi-1

15、.LoadQ)-tempP*ui0+tempQ*ui1;P_UUi=0.0;if(sBi-1.Type=1)double tempP=0.0;double tempQ=0.0;for(j=1;j<nB+1;j+)tempP+=YGij*uj1-YBij*uj0;tempQ+=YGij*uj0+YBij*uj1;P_Pi=(sBi-1.GenP-sBi-1.LoadP)-tempP*ui1-tempQ*ui0;P_UUi=sBi-1.Volt*sBi-1.Volt-ui1*ui1-ui0*ui0;P_Qi=0.0;/*偏移量阵 */double P_PQ6;int a=0;for(i=1;

16、i<3;i+)P_PQa=P_Pi;a=a+2;a=1;for(i=1;i<3;i+)P_PQa=P_Qi;a=a+2;.P_PQ4=P_P3;P_PQ5=P_UU3;printf("迭代第 %d次的偏移量为：n",v);for(i=0;i<6;i+)printf("%f",P_PQi);printf("n");printf("n");printf("n");/*雅可比矩阵*/double H66,N66,J66,L66,R66,S66,aa6,bb6;for(i=1;i<

17、;5;i+) if(fabs(sBi-1.Type-2.0)<0.000001) continue;elsefor(j=1;j<5;j+)if(i!=j)Hij=-YBij*ui1+YGij*ui0;Nij=YGij*ui1+YBij*ui0;Jij=-Nij;Lij=Hij;Rij=0;Sij=0;elseaai=bbi=0.0;for(n=1;n<5;n+).aai+=YGin*un1-YBin*un0;bbi+=YGin*un0+YBin*un1;Hii=-YBii*ui1+YGii*ui0+bbi;Nii=YGii*ui1+YBii*ui0+aai;Jii=-YGii

18、*ui1-YBii*ui0+aai;Lii=YGii*ui0-YBii*ui1-bbi;Rii=2*ui0;Sii=2*ui1;double ss5050;for(i=0;i<6;i+)for(j=0;j<6;j+)ssij=0.0;for(i=1;i<3;i+)for(j=1;j<4;j+)ss2*i-22*j-2=Hij;ss2*i-22*j-1=Nij;ss2*i-12*j-2=Jij;ss2*i-12*j-1=Lij;i=3;for(j=1;j<4;j+)for(j=1;j<4;j+).ss2*i-22*j-2=Hij;ss2*i-22*j-1=Ni

19、j;ss2*i-12*j-2=Rij;ss2*i-12*j-1=Sij;printf("迭代第 %d次的雅可比矩阵为：n",v);for(i=0;i<6;i+)for(j=0;j<6;j+)printf("%10f",ssij);printf("n");printf("n");printf("n");Gauss(ss,P_PQ,6);for(i=1;i<nB;i+)ui0=ui0+mn2*(i-1);ui1=ui1+mn2*i-1;double max;max=fabs(P_P

20、Q0);for(i=0;i<=5;i+)if (max<fabs(P_PQi)max=fabs(P_PQi);if(fabs(max)<0.0001)printf("满足精度要求，迭代终止，迭代次数为%dn",v);printf("n");printf("n");break;./*叠代循环的括号*/printf("最终求得的节点电压值为e+jf:n");for(i=1;i<nB+1;i+)printf("%lf,%lfn",ui1,ui0);printf("n&

21、quot;);printf("n");double uu5,Phase5;for(i=1;i<nB+1;i+)uui=sqrt(ui1*ui1+ui0*ui0);Phasei=atan(ui0/ui1);for(i=1;i<nB+1;i+)printf("%lf,%lfn",uui,Phasei);* 计算线路功率和平衡节点PV 节点功率 */double P5,Q5;double tempP=0.0;double tempQ=0.0;for(i=1;i<nB+1;i+)for(j=1;j<nB+1;j+)tempP+=YGij*

22、uj1-YBij*uj0;tempQ+=YGij*uj0+YBij*uj1;Pi=tempP*ui1+tempQ*ui0;Qi=tempP*ui0-tempQ*ui1;tempP=0.0;tempQ=0.0;.for(i=1;i<nB+1;i+)printf("节点 %d注入功率为 %lf,%lfn",i,Pi,Qi);/*支路功率 */double V42;for(i=1;i<5;i+)for(j=0;j<3;j+)Vij=uij;double sP55,sQ55;double dsq,dsp,dp,sumgen;for(i=1;i<NBUS+1;

23、i+)for(j=1;j<NBUS+1;j+)sPij=0.0;sQij=0.0;for(l=0; l<nL; l+).i=sLl.NumI;j=sLl.NumJ;r=sLl.R;x=sLl.X;d1=r*r+x*x;g=r/d1;b=-x/d1;.if(fabs(sLl.K-1.0)<0.000001)/*Normal lines or transformers*/sPij=Vi1*Vi1*g-Vi1*Vj1*(g*cos(Vi0-Vj0)+b*sin(Vi0-Vj0);sQij=-(Vi1*Vi1*sLl.B+Vi1*Vi1*b+Vi1*Vj1*(g*sin(Vi0-Vj

24、0)-b*cos(Vi0-Vj0);sPji=Vj1*Vj1*g-Vi1*Vj1*(g*cos(Vj0-Vi0)+b*sin(Vj0-Vi0);sQji=-(Vj1*Vj1*sLl.B+Vj1*Vj1*b+Vi1*Vj1*(g*sin(Vj0-Vi0)-b*cos(Vj0-Vi0);else/*abnormal transformer ratio*/sPij=Vi1*Vi1*g/sLl.B/sLl.B-Vi1*Vj1*(g*cos(Vi0-Vj0)/sLl.B+b*sin(Vi0-Vj0)/sLl.B);sQij=-(Vi1*Vi1*b/sLl.B/sLl.B+Vi1*Vj1*(g*sin(V

25、i0-Vj0)/sLl.B-b*cos(Vi0-Vj0)/sLl.B);sPji=Vj1*Vj1*g-Vi1*Vj1*(g*cos(Vj0-Vi0)/sLl.B+b*sin(Vj0-Vi0)/sLl.B);sQji=-(Vi1*Vi1*b+Vi1*Vj1*(g*sin(Vj0-Vi0)/sLl.B-b*cos(Vj0-Vi0)/sLl.B);/*输电效率*/dsp=P4;dsq=Q4;sumgen=P4;.for(i=0;i<NBUS;i+)dsp+=sBi.GenP-sBi.LoadP;dsq+=sBi.GenQ-sBi.LoadQ;sumgen+=sBi.GenP;dp=dsp/su

26、mgen*100;/*输出功率情况*/if(fp=fopen("功率情况 .txt","w")=NULL)printf("Can not open the file named '功率情况 .txt' n");exit(0);fprintf(fp,"-功率情况 -n");fprintf(fp,"平衡节点功率S=%10.5f+ j%10.5fn",P4,Q4);for(i=1;i<NBUS+1;i+)for(j=1;j<NBUS+1;j+)if(fabs(sPij-0.

27、0)>0.000001)fprintf(fp,"S(%3d,%-3d)=(%10.5f,j%10.5f)n",i,j,sPij,sQij);fprintf(fp,"网损为 %10.5f+j%10.3f,输电效率为 %10.3fn",dsp,dsq,100-dp);if(fp!=NULL) fclose(fp);return 0;结果：.1导纳阵Y( 1,1 )=(1.01534, -8.19201)Y( 1,2 )=( -0.56148,2.30208)Y( 1,3 )=(0.00000,3.66667)Y( 1,4 )=( -0.45386,1.

28、89107)Y( 2,1 )=( -0.56148,2.30208)Y( 2,2 )=(1.04225, -4.67651)Y( 2,4 )=( -0.48077,2.40385)Y( 3,1 )=(0.00000,3.66667)Y( 3,3 )=(0.00000, -3.33333)Y( 4,1 )=( -0.45386,1.89107)Y( 4,2 )=( -0.48077,2.40385)Y( 4,4 )=(0.93463, -4.26159)2设定电压初值e1(0)jf 1( 0)1 j 0;e(20)jf2(0)1j3计算功率和电压偏移；4P( 0)PsP(0)Pse( 0)(Gj ej( 0)111111j1Q(0)QQ(0)Qf( 0)4G e( 0)(11s11s11jjj1同理可算出0;e(30)B1j f 1( 0) )B1j f 1( 0) )jf 3( 0)1. 1j 04f( 0)(Gj f j( 0)11j14e1(0)(G1j f j( 0)j1B1j e(10) )0. 27731B1j e(10) )0.05097P2( 0)P2