完整版潮流計算C++程序

1、程序為計算書3-4的過程程序可以解決開式單直網(wǎng)絡和樹狀網(wǎng)絡的計算.樹狀網(wǎng)絡計算時要自己先設定好支路的起始節(jié)點和終止節(jié)點標號以及計算順序源代碼：#include <iostream>#include <fstream>#include<iomanip>#include<math.h>using namespace std;struct node節(jié)點類int i;/節(jié)點編號double U,P,Q,delta;/額定電壓計算負荷電壓相角;struct line線路類連接父節(jié)點子節(jié)點node f_node,s_node;/ 父節(jié)點子節(jié)點double R

2、,X,B;/線路參數(shù) R X B/2double 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

3、_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);/子節(jié)點電壓lin.s_node.delta=lin.f_node.delta-atan(lin.d_U/

4、(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;/ 子節(jié)點電壓lin.s_node.delta=lin.f_node.delta-atan(lin.d_U/(lin.f_node.U-lin.D_U);void main()i

5、nt num_l;int num_n;/支路數(shù)節(jié)點數(shù)ifstream fin;fin.open("E:data.txt");fin>>num_n>>num_l;輸入節(jié)點數(shù)支路數(shù)ofstream fout;fout.open("E:databak.txt");node *nod;nod=new nodenum_n;/ 節(jié)點數(shù)目line *lin;lin=new linenum_l;/ 線路數(shù)目nod0.delta=0;double *u;/節(jié)點額定電壓u=new doublenum_n;for(int i=0;i<num_n;

6、i+)fin>>ui;double *p;/節(jié)點有功功率p=new doublenum_n;for(int i=0;i<num_n;i+)fin>>pi;double *q;/節(jié)點無功功率q=new doublenum_n;for(int i=0;i<num_n;i+)fin>>qi;for(int i=0;i<num_n;i+)/ 設定節(jié)點標號參數(shù)nodi.i=i;nodi.P=pi;nodi.Q=qi;nodi.U=ui;double *r;/線路電阻r=new doublenum_l;for(int i=0;i<num_l;i+

7、)fin>>ri;)；double *x;/線路電抗x=new doublenum_l;for(int i=0;i<num_l;i+)fin>>xi;);double *b;/線路電納b=new doublenum_l;for(int i=0;i<num_l;i+)fin>>bi;);for(int i=0;i<num_l;i+)/ 設定線路參數(shù)lini.R=ri;lini.X=xi;lini.B=bi;);for(int i=0;i<num_l;i+)/ 確定線路節(jié)點關系fin>>lini.f_node.i;fin>

8、;>lini.s_node.i;);for(int i=0;i<num_l;i+)/ 計算節(jié)點運算負荷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;/保存運算負荷數(shù)據(jù)P_c=new doublenum_n;Q_c=new doublenum_n;for(int tc=0;tc<num_n;tc+)P_ctc=nodtc.P;Q_

9、ctc=nodtc.Q;);for(int i=0;i<num_l;i+)/ 設定線路節(jié)點數(shù)據(jù)lini.Set_node(nodlini.f_node.i,nodlini.s_node.i););for(int re=0;re<3;re+)/ 迭代運算開始fout<<"第''<<re+1<<"次迭代"<<endl;double to_P(0),to_dP(0);for(int i=1;i<num_n;i+)to_P+=nodi.P;);for(int tc=0;tc<num_n

10、;tc+)/ 重置運算負荷nodtc.P=P_ctc;nodtc.Q=Q_ctc;);for(int ts=0;ts<num_l;ts+)/置各線路初始輸出功率為子節(jié)點運算負荷 lints.P_out=lints.s_node.P;lints.Q_out=lints.s_node.Q;);for(int i=0;i<num_l;i+)/ 設定線路節(jié)點數(shù)據(jù)lini.Set_node(nodlini.f_node.i,nodlini.s_node.i););int j=num_l-1;/反向求各支路功率損耗和功率分布for(j=num_l-1;j>0;j-)fun1(linj);

11、nodlinj.f_node.i.P+=linj.P_in;nodlinj.f_node.i.Q+=linj.Q_in;for(int i=0;i<num_l;i+)/ 設定線路節(jié)點數(shù)據(jù)lini.Set_node(nodlini.f_node.i,nodlini.s_node.i););for(int ts=0;ts<num_l;ts+)lints.P_out=lints.s_node.P;lints.Q_out=lints.s_node.Q;););fun1(linj);for(int i=0;i<num_l;i+)/ 設定線路節(jié)點數(shù)據(jù)lini.Set_node(nodli

12、ni.f_node.i,nodlini.s_node.i););int t=0;/求線路各點電壓for(t=0;t<num_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;i<num_l;i+)/ 設定線路節(jié)點數(shù)據(jù)lini.Set_node(nodlini.f_node.i,nodlini.s_node.i);););fun2(lint);nodlint.s_node.i.U=lint.s_node.U;nodlint.s_

13、node.i.delta=lint.s_node.delta;fout<<"支路信息:"<<endl;for(int i=0;i<num_l;i+)/ 輸出線路信息fout<<"支路"<<lini.f_node.i+1<<"-'<<lini.s_node.i+1<<":"<<endl;fout<<"始端功率:"<<lini.P_in<<"+j"

14、;<<lini.Q_in<<endl;fout<<"末端功率:"<<lini.P_out<<"+j"<<lini.Q_out<<endl;fout«"功率損耗:"«lini.d_P«"+j"«lini,d_Q«endl;fout<<"電壓損耗"«lini.f_node.U-lini.s_node.U«endl;；foutvv&quo

15、t;節(jié)點信息:"«endl;for(int i=0;i<num_n;i+)/ 輸出節(jié)點信息fout«"節(jié)點"«i+1«endl;fout«"電壓:"«nodi.U«"相角:"«nodi.delta*180/3.14«endl;；double *lu;求最低電壓及最低電壓點lu=new doublenum_n;int *lua;lua=new intnum_n;for(int i=0;i<num_n;i+)lui=nodi.U

16、;luai=i;；for(int i=0;i<num_n-1;i+)if(lui<lui+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;i<num_l;i+)to_dP+=lini.d_P;)；foutvv"全網(wǎng)信息："«endl;foutvv"總電源有功："«lin0.P_in«endl;foutvv"總負荷有功："«to_P«endl;fou

17、tvv"總有功損耗："«to_dP«endl;fout«"網(wǎng)損率:"«to_dP/(to_P+lin0.P_in)«endl;foutvv"最低電壓："«lunum_n-1«"最低電壓點："vvluanum_n-1+1vendl;deletelu;deletelua;；deletenod;deletelin;deleteu;deletep;deleteq;deleter;deletex;deleteb;)附：3-4的data文件3 2113 11

18、0 1000 0,17 200 1,7 158.5 1.2220,5 20,20.000282 00 11 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.202138+j0.151875末端功率:0.2+j0.15功率損耗:0.0009375+j0.001875電壓損耗0.07396

19、43節(jié)點信息：節(jié)點1電壓:10.5相角:0節(jié)點2電壓:10.2259 相角:-0,864932節(jié)點3電壓:10.1169 相角:-1,25281節(jié)點4電壓:10,152 相角:-1,072全網(wǎng)信息：總電源有功:1.02165總負荷有功：1 總有功損耗：0.0216481網(wǎng)損率: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末端功率:0.5

20、+j0.3功率損耗:0.00332186+j0.00664371電壓損耗0.108957支路2-4:始端功率:0.20211+j0.151819末端功率:0.2+j0.15功率損耗:0.000909642+j0.00181928電壓損耗0.0739403節(jié)點信息：節(jié)點1電壓:10.5相角:0節(jié)點2電壓:10.2264 相角:-0.86489節(jié)點3電壓:10.1175 相角:-1.25273節(jié)點4電壓:10.1525 相角:-1.07194全網(wǎng)信息：總電源有功:1.02078總負荷有功：1.70434總有功損耗:0.0207799網(wǎng)損率:0.00762533最低電壓:10.1175最低電壓點:3

21、第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.20211+j0.151819末端功率:0.2+j0.15功率損耗:0.000909549+j0.0018191電壓損耗0.0739402節(jié)點信息：節(jié)點1電壓:10.5相角:0節(jié)點2電壓:10.2264 相角:-0.86

22、489節(jié)點3電壓:10.1175 相角:-1.25273節(jié)點4電壓:10.1525 相角:-1.07194全網(wǎng)信息：總電源有功:1.02078總負荷有功:1.70423總有功損耗:0.0207778網(wǎng)損率:0.00762487最低電壓:10.1175最低電壓點:3例 3-2 的 data :4 310.5 10 10 100 0.3 0.5 0.20 0.2 0.3 0.151.2 1 1.52.4 2 30 0 00 11 21 3輸出：第1次迭代支路信息：支路1-2:始端功率:1.02165+j0.693296末端功率：1.00434+j0.658675功率損耗:0.0173106+j0.

23、0346211電壓損耗0.275227支路2-3:始端功率:0.5034+j0.3068末端功率:0.5+j0.3功率損耗:0.0034+j0.0068電壓損耗0.109244支路2-4:始端功率:0.202138+j0.151875末端功率:0.2+j0.15功率損耗:0.0009375+j0.001875電壓損耗0.0740389節(jié)點信息：節(jié)點1電壓:10.5相角:0節(jié)點2電壓:10.2248 相角:-0.864932節(jié)點3電壓:10.1155 相角:-1.2529節(jié)點4電壓:10.1507 相角:-1.07205全網(wǎng)信息：總電源有功：1.02165總負荷有功:1總有功損耗：0.0216481網(wǎng)損率:0.0107081最低電壓:10.1155最低電壓點:3第2次迭代支路信息：支路1-2:始端功率:1.02078+j0.69157末端功率：1.00423+j0.658465功率損耗:0.0165523+j0.0331045電