《潮流计算c++程序》由会员分享，可在线阅读，更多相关《潮流计算c++程序（11页珍藏版）》请在金锄头文库上搜索。
1、程序为计算书 3-4 的过程程序可以解决开式单直网络和树状网络的计算。树状网络计算时要自己先设定好支路的起始节点和终止节点标号以及计算顺序源代码：#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_
2、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 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 double q=l
3、in.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 double q=lin.Q_in; double u=lin.f_node.U; lin.D_U=(p*l
4、in.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:databa
5、k.txt“); node *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;iri; ; double *x;/线路电抗 x=ne
6、w doublenum_l; for(int i=0;ixi; ; double *b;/线路电纳 b=new doublenum_l; for(int i=0;ibi; ; for(int i=0;ilini.f_node.i; finlini.s_node.i; ; 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.
7、s_node.i); ; for(int ts=0;tsnum_l;ts+) lints.P_out=lints.s_node.P; lints.Q_out=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.delt
8、a=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“始端功
9、率:“lini.P_in“+j“lini.Q_inendl; fout“末端功率:“lini.P_out“+j“lini.Q_outendl; fout“功率损耗:“lini.d_P“+j“lini.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 doub
10、lenum_n; int *lua; lua=new intnum_n; for(int i=0;inum_n;i+) lui=nodi.U; luai=i; ; for(int i=0;inum_n-1;i+) if(luilui+1) double st; int a; st=lui; a=luai; lui=lui+1; luai=luai+1; lui+1=st; luai+1=a; ; ; for(int i=0;inum_l;i+) to_dP+=lini.d_P; ; fout“全网信息:“endl; fout“总电源有功:“lin0.P_inendl; fout“总负荷有功:“
11、to_Pendl; fout“总有功损耗:“to_dPendl; fout“网损率:“to_dP/(to_P+lin0.P_in)endl; fout“最低电压:“lunum_n-1“最低电压点:“luanum_n-1+1endl; deletelu; deletelua; ; deletenod; deletelin; deleteu; deletep; deleteq; deleter; deletex; deleteb; 附：3-4 的 data 文件3 2 113 110 100 0 0.17 20 0 1.7 15 8.5 1.22 20.5 20.2 0.000282 0 0 1
12、1 2 输出为：第 1 次迭代 支路信息: 支路 1-2: 始端功率:1.02165+j0.693296 末端功率:1.00434+j0.658675 功率损耗:0.0173106+j0.0346211 电压损耗 0.274063 支路 2-3: 始端功率:0.5034+j0.3068 末端功率:0.5+j0.3 功率损耗:0.0034+j0.0068 电压损耗 0.109 支路 2-4: 始端功率:0.200938+j0.151875 末端功率:0.2+j0.15 功率损耗:0.0009375+j0.001875 电压损耗 0.0739643 节点信息: 节点 1 电压:10.5 相角:0
13、节点 2 电压:10.2259 相角:-0.864932 节点 3 电压:10.1169 相角:-1.25281 节点 4 电压:10.152 相角:-1.072 全网信息: 总电源有功:1.02165 总负荷有功:1 总有功损耗:0.0216481 网损率:0.0107081 最低电压:10.1169 最低电压点:3 第 2 次迭代 支路信息: 支路 1-2: 始端功率:1.02078+j0.69156 末端功率:1.00423+j0.658463 功率损耗:0.0165484+j0.0330969 电压损耗 0.273567 支路 2-3: 始端功率:0.503322+j0.306644
14、末端功率:0.5+j0.3 功率损耗:0.00332186+j0.00664371 电压损耗 0.108957 支路 2-4: 始端功率:0.20091+j0.151819 末端功率:0.2+j0.15 功率损耗:0.000909642+j0.00181928 电压损耗 0.0739403 节点信息: 节点 1 电压:10.5 相角:0 节点 2 电压:10.2264 相角:-0.86489 节点 3 电压:10.1175 相角:-1.25273 节点 4 电压:10.1525 相角:-1.07194 全网信息: 总电源有功:1.02078 总负荷有功:1.70434 总有功损耗:0.0207
15、799 网损率:0.00762533 最低电压:10.1175 最低电压点:3 第 3 次迭代 支路信息: 支路 1-2: 始端功率:1.02078+j0.691556 末端功率:1.00423+j0.658462 功率损耗:0.0165468+j0.0330936 电压损耗 0.273566 支路 2-3: 始端功率:0.503322+j0.306643 末端功率:0.5+j0.3 功率损耗:0.0033215+j0.00664301 电压损耗 0.108957 支路 2-4: 始端功率:0.20091+j0.151819 末端功率:0.2+j0.15 功率损耗:0.000909549+j0.0018191 电压损耗 0.0739402 节点信息: 节点 1 电压:10.5 相角:0 节点 2 电压:10.2264 相角:-0.8648
