潮流计算C++程序.doc

程序为计算书3-4的过程 程序可以解决开式单直网络和树状网络的计算。树状网络计算时要自己先设定好支路的起始节点和终止节点标号以及计算顺序源代码： #include iostream #include fstream #includeiomanip #includemath.h using namespace std; struct node{//节点类 int i;//节点编号 double U,P,Q,delta;//额定电压 计算负荷 电压相角 }; struct line{//线路类 连接父节点子节点 node f_node,s_node;//父节点 子节点 double R,X,B;//线路参数R X B/2 double P_in,Q_in,P_out,Q_out,d_P,d_Q,D_U,d_U;//线路输入输出功率以及线路消耗功率 void Set_node(node nod1,node nod2){ f_node=nod1; s_node=nod2; } }; void fun1(line lin){//由后往前递推功率 double p=lin.P_out; double q=lin.Q_out; double u=lin.s_node.U; lin.d_P=(p*p+q*q)/u/u*lin.R; lin.d_Q=(p*p+q*q)/u/u*lin.X; lin.P_in=lin.d_P+lin.P_out; lin.Q_in=lin.d_Q+lin.Q_out; }; void fun2(line lin){//由前往后推电压 double p=lin.P_in; double q=lin.Q_in; double u=lin.f_node.U; lin.D_U=(p*lin.R+q*lin.X)/u; lin.d_U=(p*lin.X-q*lin.R)/u; lin.s_node.U=sqrt(pow(lin.f_node.U-lin.D_U,2)+pow(lin.d_U,2));//子节点电压 lin.s_node.delta=lin.f_node.delta-atan(lin.d_U/(lin.f_node.U-lin.D_U)); }; void fun3(line lin){//由前往后推电压不计横向分量 double p=lin.P_in; double q=lin.Q_in; double u=lin.f_node.U; lin.D_U=(p*lin.R+q*lin.X)/u; lin.d_U=(p*lin.X-q*lin.R)/u;? lin.s_node.U=lin.f_node.U-lin.D_U,2;//子节点电压 lin.s_node.delta=lin.f_node.delta-atan(lin.d_U/(lin.f_node.U-lin.D_U));? }; void main(){ int num_l; int num_n;//支路数 节点数 ifstream fin; fin.open(E:\\data.txt); finnum_nnum_l;//输入节点数支路数 ofstream fout; fout.open(E:\\databak.txt); node *nod; nod=new node[num_n];//节点数目 line *lin; lin=new line[num_l];//线路数目 nod[0].delta=0; double *u;//节点额定电压 u=new double[num_n]; for(int i=0;inum_n;i++){ finu[i]; }; double *p;//节点有功功率 p=new double[num_n]; for(int i=0;inum_n;i++){ finp[i]; }; double *q;//节点无功功率 q=new double[num_n]; for(int i=0;inum_n;i++){ finq[i]; }; for(int i=0;inum_n;i++){//设定节点标号 参数 nod[i].i=i; nod[i].P=p[i]; nod[i].Q=q[i]; nod[i].U=u[i]; }; double *r;//线路电阻 r=new double[num_l]; for(int i=0;inum_l;i++){ finr[i]; }; double *x;//线路电抗 x=new double[num_l]; for(int