电网潮流计算程序C版本

1、细心整理看大家在人人上发了各种电力系统作业，刚编完我也来得瑟一下吧，根本上可以解决开式单直网络和树状网络计算。树状网络计算时要自己先设定好支路起始节点和终止节点标号以及计算依次，原来想做一个树状图深度遍历来进展自动排序，不过从网上找到一个深度遍历代码太长，在此就偷个懒把依次自己输进去了。测试后3-4,3-3,以及例3-2答案根本正确。另外fun3是不参与横向电压重量计算函数，在例3-2计算用是这个函数。另外就是迭代终止条件问题上我只是设定固定迭代次数为3，也可以用精度来限定，大致就是如此吧。请大家多多指教了。源代码：#include #include #include #include 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

2、_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

3、_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); n

4、ode *nod; nod=new nodenum_n;/节点数目 line *lin; lin=new linenum_l;/线路数目 nod0.delta=0; double *u;/节点额定电压 u=new doublenum_n; for(int i=0;iui; ; double *p;/节点有功功率 p=new doublenum_n; for(int i=0;ipi; ; double *q;/节点无功功率 q=new doublenum_n; for(int i=0;iqi; ; for(int i=0;inum_n;i+)/设定节点标号 参数 nodi.i=i; nodi.P=pi; nodi.Q=qi; nodi.U=ui; ; double *r;/线路电阻 r=new doublenum_l; for(int i=0;iri; ; double *x;/线路电抗 x=new doublenum_l; for(int i=0;ixi; ; double *b;/线路电纳 b=new doublenum_l; for(int i=0;ibi; ; for(int i

5、=0;inum_l;i+)/设定线路参数 lini.R=ri; lini.X=xi; lini.B=bi; ; for(int i=0;ilini.f_node.i; finlini.s_node.i; ; for(int i=0;inum_l;i+)/计算节点运算负荷 nodlini.f_node.i.Q-=lini.B*nodlini.f_node.i.U*nodlini.f_node.i.U; nodlini.s_node.i.Q-=lini.B*nodlini.s_node.i.U*nodlini.s_node.i.U; ; double *P_c,*Q_c;/保存运算负荷数据 P_c=new doublenum_n; Q_c=new doublenum_n; for(int tc=0;tcnum_n;tc+) P_ctc=nodtc.P; Q_ctc=nodtc.Q; ; for(int i=0;inum_l;i+)/设定线路节点数据 lini.Set_node(nodlini.f_node.i,nodlini.s_node.i); ; for(int re=0;re3;re

6、+)/迭代运算起先 fout第re+1次迭代endl; double to_P(0),to_dP(0); for(int i=1;inum_n;i+) to_P+=nodi.P; ; for(int tc=0;tcnum_n;tc+)/重置运算负荷 nodtc.P=P_ctc; nodtc.Q=Q_ctc; ; for(int ts=0;tsnum_l;ts+)/置各线路初始输出功率为子节点运算负荷 lints.P_out=lints.s_node.P; lints.Q_out=lints.s_node.Q; ; for(int i=0;i0;j-) fun1(linj); nodlinj.f_node.i.P+=linj.P_in; nodlinj.f_node.i.Q+=linj.Q_in; for(int i=0;inum_l;i+)/设定线路节点数据 lini.Set_node(nodlini.f_node.i,nodlini.s_node.i); ; for(int ts=0;tsnum_l;ts+) lints.P_out=lints.s_node.P; lints.Q_o

7、ut=lints.s_node.Q; ; ; fun1(linj); for(int i=0;inum_l;i+)/设定线路节点数据 lini.Set_node(nodlini.f_node.i,nodlini.s_node.i); ; int t=0;/求线路各点电压 for(t=0;tnum_l-1;t+) fun2(lint); nodlint.s_node.i.U=lint.s_node.U; nodlint.s_node.i.delta=lint.s_node.delta; for(int i=0;inum_l;i+)/设定线路节点数据 lini.Set_node(nodlini.f_node.i,nodlini.s_node.i); ; ; fun2(lint); nodlint.s_node.i.U=lint.s_node.U; nodlint.s_node.i.delta=lint.s_node.delta; fout支路信息:endl; for(int i=0;inum_l;i+)/输出线路信息 fout支路lini.f_node.i+1-lini.s_node.i+1:endl; fout始端功率:lini.P_in+jlini.Q_inendl; fout末端功率:lini.P_out+jlini.Q_outendl; fout功率损耗:lini.d_P+jlini.d_Qendl; fout电压损耗lini.f_node.U-lini.s_node.Uendl; ; fout节点信息:endl; for(int i=0;inum_n;i+)/输出节点信息 fout节点i+1endl; fout电压:nodi.U相角:nodi.delta*180/3.14endl; ; double *lu;/求最低电压及最低电压点 lu=new doublenum_n; int *lua; lua=new intnum_n;

《电网潮流计算程序C版本》由会员公****分享，可在线阅读，更多相关《电网潮流计算程序C版本》请在金锄头文库上搜索。