通用潮流程序.doc

电力系统分析——潮流计算 潮流计算框图： 事业第一,但别把钱看得太重；学会体谅父母；交上好朋友；不要沉迷于任何东西；年轻没有失败；不要轻易崇拜或者鄙视一个人；要有责任心；外貌并不重要；学会保护身体；别觉得一事无成；认真工作；认真对待感情；留一点童心。努力做到这些，年轻没有失败， #include iostream #include fstream #include FuShu.h #include math.h #include iomanip 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,double **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,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; Function4(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); cout\n\n节点导纳矩阵实部G为:endl; Y.PrintG(); cout\n\n节点导纳矩阵虚部B为:endl; Y.PrintB(); coutendl; ifstream fin2(Data2.txt); i=1;j=1; m=0; while(!fin2.eof()) { fin2iP[j]Q[j]U[j]; if(U[j]==0) m++; j++; } fin2.close(); e=new double [n+1]; f=new double [n+1]; de=new