1、電力系統潮流上機課程設計報告院系：電氣與電子工程學院班級:電氣1108學號：1111550112學生姓名：龍日尚指導教師：劉寶柱設計周數：兩周成績：日期：2014年1月10日一、課程設計的目的與要求培養學生的電力系統潮流計算機編程能力，掌握計算機潮流計算的相關知識二、設計正文（詳細內容見附錄）1. 手算2. 計算機計算3. 思考題三、課程設計總結或結論（詳細內容見附錄）四、參考文獻1. 電力系統計算：電子數字計算機的應用，西安交通大學等合編。北京：水利電力出版社;2. 現代電力系統分析，王錫凡主編，科學出版社；3. 電力系統穩態分析，陳珩，中國電力出版社，1995年，第三版；附錄（設計流程圖、

2、程序、表格、數據等）4. 機算潮流程序及結果/ dierti.cpp :定義控制臺應用程序的入口點/#include "stdafx.h"struct Line /線路結構體int Num,Numl,NumJ; /線路號左節點名右節點名float R,X,B,K;/電阻電抗電納變比（K等于1為普通支路， 不等于1為變壓器支路的變比）;struct Bus/int Num ;節點結構體float Volt,Phase,GenP,GenQ,LoadP,LoadQ;int Type;；#include"stdio.h"#include"string.

3、h"#include"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+) JS

4、i=0;for(k=0;k<n;k+)d=0.0;for(j=k;j<n;j+)if(fabs(akj)>d) /*d=fabs(akj);JSk=j;if(fabs(d)<0.000001)/*L=0;else if(JSk!=k)for(i=0;i<n;i+)t=aik;aik=aiJSk;/*aiJSk=t;if(L=0)break;在一行中找到一個最大值賦值d,并用JSK記住這個最大值所在的列號*/如果d的數值太小，做為被除數將帶來很大的誤差*/進行列交換，讓最大值始終在對角元上*/for(j=k+1;j<n;j+)akj=akj/akk; /*對角

5、元上的元素消為1*/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;x*/if(fabs(an-1n-1)>0.00001)/用追趕法，解方程組，求未知數xn-1=bn-1;for(i=n-2;i>=0;i-)t=0.0;for(j=i+1;jvn;j+)t=t+a【ij*xj;x【i=(b【i-t);if(fp=fopen("gauss.txt","w")=NULL)/*將結果寫到 TXT文件中 */

6、printf("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;stru

7、ct Line sLNLINE; 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"

8、;,&nB,&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.

9、GenQ=d4;sBi.LoadP=d5,sBi.LoadQ=d6;sBi.T ype=i2;for(i=0;i<nL;i+)sLi.Num=sLi.Numl=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.Numl=i2;sLi.NumJ=i3;sLi.R=d1;sLi.X=d2;sLi.B=d3;sLi.K=

10、d4; if(fp!=NULL) fclose(fp);/*Make Y Matrix*/ for(i=1;i<nB+1;i+)for(j=1;j<nB+1;j+)YGij=O.O;YBij=O.O；;for(l=0; IvnL; l+)i=sLl.Numl;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=Y

11、Gij-g;YGji=YGji-g;YBij=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("C

12、an not open the 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(s

13、Bi-1.Type=0)ui0=0.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;*/e+jf:n",v);for(v=1;v+)/*迭代次數可以無限大/*節點電壓附初值*/printf(”迭代第d次賦予的電壓初值為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

15、-sBi-1.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;f

16、or(i=1;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(”迭代第4次的偏移量為：n",v);for(i=0;i<6;i+)printf("%f",P_PQi);printf("n");printf("n");printf("n");/*雅可比矩陣*/double H6 6,N 6 6,J66,L 6 6,R 6 6,S66,aa 6,bb 6;fo

17、r(i=1;i<5;i+)if(fabs(sBi-1.Type-2.0)<0.000001)continue;elsefor(j=1;j<5;j+)if(i!=j)H【ij=-YB【ij*u【i【1+YG【ij*u【i【0;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

18、+aai;J【i【i=-YG【i【i*u【i【1-YB【i【i*u【i【0+aa【i;L【i【i=YG【i【i*u【i【0-YB【i【i*u【i【1-bb【i;R【i【i=2*u【i【0;S【i【i=2*u【i【1;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+)ss【2*i-2【2*j-2=H【ij; 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

19、(j=1;j<4;j+)ss2*i-22*j-2=Hij;ss2*i-22*j-1=Nij;ss2*i-12*j-2=Rij;ss2*i-12*j-1=Sijprintf(”迭代第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+)u【i【0=u【i【0+mn【2*(

20、i-1);ui1=ui1+mn2*i-1;double max;max=fabs(P_PQ0);for(i=0;i<=5;i+)if (max<fabs(P_PQi)max=fabs(P_PQi);if(fabs(max)<0.0001)%dn",v);printf(”滿足精度要求，迭代終止，迭代次數為printf("n");printf("n");break;/*疊代循環的括號*/printf(”最終求得的節點電壓值為e+jf:n");for(i=1;i<nB+1;i+) printf("%lf,%

21、lfn",ui1,ui0);printf("n");printf("n");double uu5,Phase5;for(i=1;i<nB+1;i+) uu【i=sqrt(u【i【1*u【i【1+u【i【0*u【i【0); 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&l

22、t;nB+1;i+)for(j=1;j<nB+1;j+)tempP+=YGij*uj1-YBij*uj0; tempQ+=YGij*uj0+YBij*uj1;Pi=tempP*ui1+tempQ*ui0;Qi=tempP*ui0-tempQ*ui1;tempQ=0.0;tempP=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+)V【ij=u【ij;double sP55,sQ55;double

23、dsq,dsp,dp,sumgen;for(i=1;i<NBUS+1;i+)for(j=1;j<NBUS+1;j+)sPij=0.0;sQij=0.0;for(l=0; IvnL; l+)i=sLl.Numl;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.

24、B+Vi1*Vi1*b+Vi1*Vj1*(g*sin(Vi0-Vj0)-b*cos(Vi 【0-Vj【0)；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(Vj 0-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=-(V

25、i1*Vi1*b/sLl.B/sLl.B+Vi1*Vj1*(g*sin(Vi0-Vj0)/sLl.B-b*cos(Vi【0-Vj【0)/sL【l.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.G

26、enQ-sBi.LoadQ;sumgen+=sBi.GenP;dp=dsp/sumgen*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+

27、)for(j=1;j<NBUS+1;j+)if(fabs(sPij-0.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.00

28、000,3.66667)Y(1,4)=(-0.45386,1.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 .設定電壓初值e10)jff 1 jQ 醫 jf20)1 j0;e30)jf30)1.1 j03 計算功率和電壓偏移;P(0)P11sP(0)4Pse10)(Gje(0)j 14BJ®f10)j(朝0)1j,)0. 27731Q0)QsQ0)4Qsf10) (Gj*j 14耳牢)e10)j 1(GjfB1 e