2022年潮流程序 .pdf

《2022年潮流程序 .pdf》由会员分享，可在线阅读，更多相关《2022年潮流程序 .pdf（18页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、1 电力系统分析潮流计算潮流计算框图：输入电网节点注入数据计算节点导纳矩阵参数设置节点电压初值设置 k=0 及最大迭代次数kmax 计算 fi(x(k)以及雅各比矩阵x(k+1)=x(k)+dx(k) 求解修正方程获得dx max|fi(x(k)|kmax? 输出潮流不收敛信息结束N Y N Y 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 18 页 - - - - - - - - - 2 程序部分：#include #include #include FuShu.h

2、#include #include using namespace std; double Function1(PowerGrid &y,double *e,double *f,double *w,double *P,double *Q,double *U,int n,int m); void Function2(int n,double *pa); void Function3(int n,int m,double *Ya,PowerGrid &Y,double *e,double *f); void Function4(int &n); void Function5(int k,doubl

3、e *ya,int n,PowerGrid &y,FuShu *z); void Function6(double *e,double *f,double *P,double *Q,PowerGrid &y,int n); void Function7(double *e,double *f,double *P,double *Q,PowerGrid &y,int n,int m); void Function8(PowerGrid &y,double *e,double *P,double *Q,double *U,int n,int m); void Function9(FuShu &If

4、,int &k,double *U,FuShu *z,FuShu zf,double *e,double *f,int n); void main() int n(0),m; int shuliang(0); double MAX(0.0); int i,j; int kmax(4); double *e,*de,*f,*df,*P,*dP,*Q,*dQ,*U,*dU2; double h=1e-5; double *w,*ef,*Ya; int tell(1); FuShu *Z; int k(4); FuShu If(0,0),Zf(0,0); char nomeaning; Functi

5、on4(n); P=new double n+1; Q=new double n+1; U=new double n+1; Z=new FuShu n+1; PowerGrid Y(n); Y.FromFile(Data1.txt); Y.FromFile(Data4.txt); coutnn节点导纳矩阵实部G 为:endl; Y.PrintG(); coutnn节点导纳矩阵虚部B 为:endl; Y.PrintB(); 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 1

6、8 页 - - - - - - - - - 3 coutiPjQjUj; if(Uj=0) m+; j+; fin2.close(); e=new double n+1; f=new double n+1; de=new double n+1; df=new double n+1; dU2=new double n+1; dP=new double n+1; dQ=new double n+1; ifstream fin3; fin3.open(Data3.txt); while(!fin3.eof() for(i=1;iei; dU2i=dPi=dQi=0.0; for(i=1;ifi; f

7、or(i=1;in+1;i+) dei=dfi=0.0; fin3.close(); w=new double 2*n-1; ef=new double 2*n-1; Ya=new double*2*n+1; for(i=1;i2*n+1;i+) *(Ya+i)=new double 2*n+1; if(ih) shuliang+; fout.open(out.txt); Function3(n,m,Ya,Y ,e,f); for(i=1;i2*n-1;i+) for(j=1;j2*n-1;j+) foutYaij ; coutendl; fout.close(); i=1; fin4.ope

8、n(out.txt); while(i2*n-1) for(j=1;jYaij; i+; fin4.close(); Function2(2*n-2,Ya); for(i=1;i2*n-1;i+) for(j=1;j2*n-1;j+) efi+=Yaij*wj; for(i=1;i2*n-1;i+) efi=-efi; for(i=2;in+1;i+) dei=efi-1; dfi=efi-1+n-1; for(i=1;in+1;i+) ei=ei+dei; fi=fi+dfi; MAX=Function1(Y,e,f,w,P ,Q,U,n,m); for(i=1;ikmax-1) break

9、; cout迭代次数为 :shuliangendl; cout各节点电压为:endl; for(i=1;in+1;i+) cout 节点 i: e= setw(10)eisetw(10)f= fiendl; Function6(e,f,P,Q,Y ,n); Function7(e,f,P,Q,Y ,n,m); cout极坐标形式为 :endl; for(i=1;in+1;i+) ei=sqrt(ei*ei+fi*fi); fi=atan(fi/ei); for(i=1;in+1;i+) cout节点i: 幅值: setw(10)ei 相角: setw(10)fiendl; cout各节点功率为

10、：endl; for(i=1;in+1;i+) cout 节 点 i: 有 功 功 率 :setw(10)Pi 无 功 功 率 : setw(10)Qiendl; Function8(Y,e,P,Q,U,n,m); Function5(k,Ya,n,Y ,Z); cout节电阻抗矩阵的第k列为 :endl; for(i=1;in+1;i+) coutsetw(5); Zi.Print(); coutendl; coutendl; Function9(If,k,U,Z,Zf,e,f,n); coutnn短路电流为： ; If.Print(); coutendl; cout短路后各节点电压为:n;

11、 for(i=1;in+1;i+) cout节点i: 幅值: setw(12)ei 相角: setw(10)fiendl; cout短路点为 :kendl; delete P; delete dP; 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 18 页 - - - - - - - - - 6 delete Q; delete dQ; delete e; delete de; delete f; delete df; delete U; delete dU2; dele

12、te w; delete ef; delete Z; Y.Deletep(); for(i=1;inomeaning; double Function1(PowerGrid &y,double *e,double *f,double *w,double *P,double *Q,double *U,int n,int m) int i,j; double temp2(0),temp3(0.0),temp(0); for(i=2;in+1;i+) for(j=1;jn+1;j+) temp2+=(y.pij.shibu*ej-y.pij.xubu*fj); temp3+=(y.pij.shibu

13、*fj+y.pij.xubu*ej); temp=Pi-ei*temp2-fi*temp3; wi-1=temp; temp3=temp2=0.0; temp2=temp3=0.0; for(i=n-m+1;in+1;i+) for(j=1;jn+1;j+) temp2+=(y.pij.shibu*ej-y.pij.xubu*fj); temp3+=(y.pij.shibu*fj+y.pij.xubu*ej); temp=Qi-fi*temp2+ei*temp3; wi+m-1=temp; temp3=temp2=0.0; 名师资料总结 - - -精品资料欢迎下载 - - - - - - -

14、- - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 18 页 - - - - - - - - - 7 temp2=temp3=0.0; for(i=2;in-m+1;i+) temp=Ui*Ui-ei*ei-fi*fi; wn+m+i-2=temp; temp2=0; for(i=1;itemp2) temp2=fabs(wi); return temp2; void Function2(int n,double *pa) int i; double *pb; pb=new double*n+1; for(i=1;in+1;i+) *(pb+i

15、)=new doublen+1; double *pe; pe=new double*n+1; for(i=1;in+1;i+) *(pe+i)=new doublen+1; int j,k,l; double m(1.0),h(1.0); for(i=1;in+1;i+) for(j=1;jn+1;j+) if(j!=i) peij=0; else peij=1; for(i=1;in+1;i+) for(j=1;jn+1;j+) if(paij!=0.000) m=paij; k=j; break; 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - -

16、- - - - - - 名师精心整理 - - - - - - - 第 7 页，共 18 页 - - - - - - - - - 8 if(j=(n+1) k=n+1; break; for(j=1;jn+1;j+) paij/=m; peij/=m; for(l=1;ln+1;l+) h=palk; if(l!=i) for(j=1;jn+1;j+) palj-=(h*paij); pelj-=(peij*h); if(kn+1) for(i=1;in+1;i+) for(j=1;jn+1;j+) if(paij!=0.0) k=j; break; for(j=1;jn+1;j+) pbkj=

17、peij; else cout矩阵不可逆！endl; for(i=1;in+1;i+) for(j=1;jn+1;j+) 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 8 页，共 18 页 - - - - - - - - - 9 paij=pbij; for(i=1;in+1;i+) delete *(pb+i); delete *(pe+i); delete pb; delete pe; void Function3(int n,int m,double *Ya,PowerGri

18、d &Y,double *e,double *f) int i,j,k; double temp1=0,temp2=0; double *H,*N,*M,*L,*R,*S; H=new double*n; N=new double*n; for(i=1;in;i+) *(H+i)=new doublen; *(N+i)=new doublen; M=new double*m+1; L=new double*m+1; for(i=1;im+1;i+) *(M+i)=new doublen; *(L+i)=new doublen; R=new double*n-m; S=new double*n-

19、m; for(i=1;in-m;i+) *(R+i)=new doublen; *(S+i)=new doublen; temp1=temp2=0.0; for(i=1;in;i+) for(j=1;jn;j+) if(j+1)!=(i+1) Hij=-(Y.pi+1j+1.shibu*ei+1+Y.pi+1j+1.xubu*fi+1); Nij=(Y.pi+1j+1.xubu*ei+1-Y.pi+1j+1.shibu*fi+1); 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 9

20、 页，共 18 页 - - - - - - - - - 10 else for(k=1;kn+1;k+) temp1+=(Y.pi+1k.shibu*ek-Y.pi+1k.xubu*fk); temp2+=(Y.pi+1k.shibu*fk+Y.pi+1k.xubu*ek); Hij=-temp1-(Y.pi+1i+1.shibu*ei+1+Y.pi+1i+1.xubu*fi+1); Nij=-temp2+(Y.pi+1i+1.xubu*ei+1-Y.pi+1i+1.shibu*fi+1); temp1=temp2=0.0; temp1=temp2=0.0; temp1=temp2=0.0;

21、for(i=1;im+1;i+) for(j=1;jn;j+) if(j+1)!=(n-m+i) Mij=Y.pn-m+ij+1.xubu*en-m+i-Y.pn-m+ij+1.shibu*fn-m+i; Lij=Y.pn-m+ij+1.shibu*en-m+i+Y.pn-m+ij+1.xubu*fn-m+i; else for(k=1;kn+1;k+) temp1+=(Y.pn-m+ik.shibu*fk+Y.pn-m+ik.xubu*ek); temp2+=(Y.pn-m+ik.shibu*ek-Y.pn-m+ik.xubu*fk); Mij=temp1+Y.pn-m+in-m+i.xub

22、u*en-m+i-Y.pn-m+in-m+i.shibu*fn-m+i; Lij=-temp2+Y.pn-m+in-m+i.shibu*en-m+i+Y.pn-m+in-m+i.xubu*fn-m+i; temp1=temp2=0.0; temp1=temp2=0.0; temp1=temp2=0.0; for(i=1;in-m;i+) for(j=1;jn;j+) if(j!=i) Rij=0.0; 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 10 页，共 18 页 - - -

23、 - - - - - - 11 Sij=0.0; else Rii=-2*ei+1; Sii=-2*fi+1; for(i=1;in;i+) for(j=1;j2*n-1;j+) if(jn) Yaij=Hij; else Yaij=Nij-n+1; for(i=n;in+m;i+) for(j=1;j2*n-1;j+) if(jn) Yaij=Mi-n+1j; else Yaij=Li-n+1j-n+1; for(i=n+m;in*2-1;i+) for(j=1;j2*n-1;j+) if(jn) Yaij=Ri-n-m+1j; else Yaij=Si-n-m+1j-n+1; for(i=

24、1;in;i+) delete *(H+i); delete *(N+i); delete H; delete N; for(i=1;im+1;i+) delete *(M+i); delete *(L+i); delete M; delete L; for(i=1;istendY.shibuY.xubuB2k; if(stn) n=st; if(endn) n=end; while(!fin2.eof() fin2stendY.shibuY.xubuB2k; if(stn) n=st; if(endn) n=end; fin1.close(); fin2.close(); void Func

25、tion5(int k,double *ya,int n,PowerGrid &y,FuShu *z) int i,j; for(i=1;i2*n+1;i+) for(j=1;j2*n+1;j+) if(i=n&j=n) yaij=y.pij.shibu; if(in) yaij=-y.pij-n.xubu; if(in&jn&jn) 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 12 页，共 18 页 - - - - - - - - - 13 yaij=y.pi-nj-n.shi

26、bu; Function2(2*n,ya); for(i=1;in+1;i+) zi.shibu=yaik; zi.xubu=yai+nk; void Function6(double *e,double *f,double *P,double *Q,PowerGrid &y,int n) double temp1(0),temp2(0); int j; for(j=1;jn+1;j+) temp1+=(y.p1j.shibu*ej-y.p1j.xubu*fj); temp2+=(y.p1j.shibu*fj+y.p1j.xubu*ej); P1=e1*temp1+f1*temp2; Q1=f

27、1*temp1-e1*temp2; temp1=temp2=0.0; void Function7(double *e,double *f,double *P,double *Q,PowerGrid &y,int n,int m) int i,j; double temp1(0),temp2(0); for(i=2;in-m+1;i+) for(j=1;jn+1;j+) temp1+=(y.pij.shibu*ej-y.pij.xubu*fj); temp2+=(y.pij.shibu*fj+y.pij.xubu*ej); Qi=fi*temp1-ei*temp2; temp1=temp2=0

28、.0; void Function8(PowerGrid &y,double *e,double *P,double *Q,double *U,int n,int m) y.FromFile(Data5.txt); int num(0),i; for(i=1;in+1;i+) if(Pi=0&Ui!=0) 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 13 页，共 18 页 - - - - - - - - - 14 num=i; break; if(num!=(n+1) y.pnu

29、mnum.shibu+=fabs(Pnum)/(enum*enum); y.pnumnum.xubu+=(Qnum)/(enum*enum); for(i=n-m+1;in+1;i+) y.pii.shibu+=fabs(Pi)/(ei*ei); y.pii.xubu+=(Qi)/(ei*ei); void Function9(FuShu &If,int &k,double *U,FuShu *z,FuShu zf,double *e,double *f,int n) int i; FuShu t(0,0); zk=zk+zf; zk.ToRad(); If.shibu=ek/(zk.shib

30、u); If.xubu=fk-zk.xubu; for(i=1;in+1;i+) zi.ToRad(); t.shibu=zi.shibu*If.shibu; t.xubu=zi.xubu+If.xubu; t.d=1; t.ToShiXu(); ei=ei*cos(fi); fi=ei*sin(fi); ei=ei-t.shibu; fi=fi-t.xubu; ei=sqrt(ei*ei+fi*fi); fi=atan(fi/ei); 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第

31、 14 页，共 18 页 - - - - - - - - - 15 程序附加部分：定义复数运算：#include #include #include #include using namespace std; class FuShu public: double shibu; double xubu; int d; FuShu() shibu=0; xubu=0; d=0; FuShu Set(double s,double x) return FuShu(s,x); FuShu(double s,double x) shibu=s; xubu=x; d=0; FuShu() void Pri

32、nt() cout=0) cout+xubuj; else coutxubud!=1) this-shibu=sqrt(shibu*shibu+xubu*xubu); this-xubu=atan(xubu/shibu); this-d=!d; void ToShiXu() if(this-d!=0) this-shibu=shibu*cos(xubu); this-xubu=shibu*sin(xubu); this-d=!d; ; class PowerGrid public: int NodeNumber; FuShu *p; PowerGrid(int nn) NodeNumber=n

33、n; p=new FuShu*NodeNumber+1; for(int i=0;iNodeNumber+1;i+) *(p+i)=new FuShuNodeNumber+1; 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 16 页，共 18 页 - - - - - - - - - 17 void FromFile(char s20) ifstream fin(s); int startpoint=1,endpoint=1; double B2=0,k=1; FuShu Y(0,0

34、),YB; int i=0; while(istartpoint; finendpoint; finY.shibuY.xubuB2k; YB=YB.Set(0.0,B2); if(k!=0&startpoint!=0&endpoint!=0) pstartpointstartpoint=pstartpointstartpoint+Y.DaoShu(); pendpointendpoint=pendpointendpoint+(Y.Kbei(k*k).DaoShu(); pstartpointendpoint=pstartpointendpoint-(Y.Kbei(k).DaoShu(); pe

35、ndpointstartpoint=pendpointstartpoint-(Y.Kbei(k).DaoShu(); if(k=0&startpoint!=0&endpoint!=0) pstartpointstartpoint=pstartpointstartpoint+Y.DaoShu()+YB; pendpointendpoint=pendpointendpoint+Y.DaoShu()+YB; pstartpointendpoint=pstartpointendpoint-Y.DaoShu(); pendpointstartpoint=pendpointstartpoint-Y.Dao

36、Shu(); if(startpoint!=0&endpoint=0) pstartpointstartpoint=pstartpointstartpoint+Y.DaoShu(); fin.close(); void Deletep() for(int i=1;iNodeNumber+1;i+) delete *(p+i); delete p; void Print() 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 17 页，共 18 页 - - - - - - - - - 18

37、 for(int i=1;i=NodeNumber;i+) for(int j=1;j=NodeNumber;j+) pij.Print(); cout ; coutendl; void PrintG() for(int i=1;i=NodeNumber;i+) for(int j=1;j=NodeNumber;j+) coutsetw(10)pij.shibu; /cout ; coutendl; void PrintB() for(int i=1;i=NodeNumber;i+) for(int j=1;j=NodeNumber;j+) coutsetw(10)pij.xubu; /cout ; coutendl; ; 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 18 页，共 18 页 - - - - - - - - -