电网潮流计算程序C++版本Word文件下载.docx

电网潮流计算程序C++版本Word文件下载.docx

《电网潮流计算程序C++版本Word文件下载.docx》由会员分享，可在线阅读，更多相关《电网潮流计算程序C++版本Word文件下载.docx（10页珍藏版）》请在冰豆网上搜索。

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;

voidfun2（line&

lin）{//由前往后推电压

doublep=lin.P_in;

doubleq=lin.Q_in;

doubleu=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））;

voidfun3（line&

lin）{//由前往后推电压不计横向分量

lin.s_node.U=lin.f_node.U-lin.D_U,2;

voidmain（）{

intnum_l;

intnum_n;

//支路数节点数

ifstreamfin;

fin.open（"

E:

\\data.txt"

）;

fin>

>

num_n>

num_l;

//输入节点数支路数

ofstreamfout;

fout.open（"

\\databak.txt"

node*nod;

nod=newnode[num_n];

//节点数目

line*lin;

lin=newline[num_l];

//线路数目

nod[0].delta=0;

double*u;

//节点额定电压

u=newdouble[num_n];

for（inti=0;

i<

num_n;

i++）{

u[i];

double*p;

//节点有功功率

p=newdouble[num_n];

p[i];

double*q;

//节点无功功率

q=newdouble[num_n];

q[i];

i++）{//设定节点标号参数

nod[i].i=i;

nod[i].P=p[i];

nod[i].Q=q[i];

nod[i].U=u[i];

double*r;

//线路电阻

r=newdouble[num_l];

r[i];

double*x;

//线路电抗

x=newdouble[num_l];

x[i];

double*b;

//线路电纳

b=newdouble[num_l];

b[i];

i++）{//设定线路参数

lin[i].R=r[i];

lin[i].X=x[i];

lin[i].B=b[i];

i++）{//确定线路节点关系

lin[i].f_node.i;

lin[i].s_node.i;

i++）{//计算节点运算负荷

nod[lin[i].f_node.i].Q-=lin[i].B*nod[lin[i].f_node.i].U*nod[lin[i].f_node.i].U;

nod[lin[i].s_node.i].Q-=lin[i].B*nod[lin[i].s_node.i].U*nod[lin[i].s_node.i].U;

double*P_c,*Q_c;

//保存运算负荷数据

P_c=newdouble[num_n];

Q_c=newdouble[num_n];

for（inttc=0;

tc<

tc++）{

P_c[tc]=nod[tc].P;

Q_c[tc]=nod[tc].Q;

i++）{//设定线路节点数据

lin[i].Set_node（nod[lin[i].f_node.i],nod[lin[i].s_node.i]）;

for（intre=0;

re<

3;

re++）{//迭代运算开始

fout<

<

"

第"

re+1<

次迭代"

endl;

doubleto_P（0）,to_dP（0）;

for（inti=1;

to_P+=nod[i].P;

tc++）{//重置运算负荷

nod[tc].P=P_c[tc];

nod[tc].Q=Q_c[tc];

for（intts=0;

ts<

ts++）{//置各线路初始输出功率为子节点运算负荷

lin[ts].P_out=lin[ts].s_node.P;

lin[ts].Q_out=lin[ts].s_node.Q;

intj=num_l-1;

//反向求各支路功率损耗和功率分布

for（j=num_l-1;

j>

0;

j--）{

fun1（lin[j]）;

nod[lin[j].f_node.i].P+=lin[j].P_in;

nod[lin[j].f_node.i].Q+=lin[j].Q_in;

ts++）{

intt=0;

//求线路各点电压

for（t=0;

t<

num_l-1;

t++）{

fun2（lin[t]）;

nod[lin[t].s_node.i].U=lin[t].s_node.U;

nod[lin[t].s_node.i].delta=lin[t].s_node.delta;

支路信息:

i++）{//输出线路信息

支路"

lin[i].f_node.i+1<

-"

lin[i].s_node.i+1<

:

始端功率:

lin[i].P_in<

+j"

lin[i].Q_in<

末端功率:

lin[i].P_out<

lin[i].Q_out<

功率损耗:

lin[i].d_P<

lin[i].d_Q<

电压损耗"

lin[i].f_node.U-lin[i].s_node.U<

节点信息:

i++）{//输出节点信息

节点"

i+1<

电压:

nod[i].U<

相角:

nod[i].delta*180/3.14<

double*lu;

//求最低电压及最低电压点

lu=newdouble[num_n];

int*lua;

lua=newint[num_n];

lu[i]=nod[i].U;

lua[i]=i;

num_n-1;

if（lu[i]<

lu[i+1]）{

doublest;

inta;

st=lu[i];

a=lua[i];

lu[i]=lu[i+1];

lua[i]=lua[i+1];

lu[i+1]=st;

lua[i+1]=a;

to_dP+=lin[i].d_P;

全网信息:

总电源有功:

lin[0].P_in<

总负荷有功:

to_P<

总有功损耗:

to_dP<

网损率:

to_dP/（to_P+lin[0].P_in）<

最低电压:

lu[num_n-1]<

最低电压点:

lua[num_n-1]+1<

delete[]lu;

delete[]lua;

delete[]nod;

delete[]lin;

delete[]u;

delete[]p;

delete[]q;

delete[]r;

delete[]x;

delete[]b;

}

附：

3-4的data文件

32

113110100

00.1720

01.715

8.51.22

20.520.2

0.0002820

01

12输出为：

第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

节点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.02078+j0.69156

1.00423+j0.658463

0.0165484+j0.0330969

电压损耗0.273567

0.503322+j0.306644

0.00332186+j0.00664371

电压损耗0.108957

0.20091+j0.151819

0.000909642+j0.00181928

电压损耗0.0739403

10.2264相角:

-0.86489

10.1175相角:

-1.25273

10.1525相角:

-1.07194

1.02078

1.70434

0.0207799

0.00762533

10.1175最低电压点:

第3次迭代

1.02078+j0.691556

1.00423+j0.658462

0.0165468+j0.0330936

电压损耗0.273566

0.503322+j0.306643

0.0033215+j0.00664301

0.000909549+j0.0018191

电压损耗0.0739402

1.70423

0.0207778

0.00762487

3

例3-2的data：

43

10.5101010

00.30.50.2

00.20.30.15

1.211.5

2.423

000

12

13

输出：

第1次迭代

1.02165+j0.693296

1.00434+j0.658675

0.0173106+j0.0346211

电压损耗0.275227

0.5034+j0.3068

0.5+j0.3

0.0034+j0.0068

电压损耗0.109244

0.200938+j0.151875

0.2+j0.15

0.0009375+j0.001875

电压损耗0.0740389

节点1

节点2

10.2248相角:

-0.864932

节点3

10.1155相角:

-1.2529

节点4

10.1507相角:

-1.07205

1.02165

1

0.0216481

0.0107081

10.1155最低电压点:

第2次迭代

1.02078+j0.69157

1.00423+j0.658465

0.0165523+j0.0331045

电压损耗0.274734

0.503323+j0.306646

0.00332278+j0.00664556

电压损耗0.109201

0.20091+j0.15182

0.000909864+j0.00181973

电压损耗0.0740151

10.2253相角:

-0.86489

10.1161相角:

-1.25282

10.1513相角:

-1.07199

1.02078

1.70434

0.0207849

0.00762714

10.1161最低电压点:

第3次迭代

1.02078+j0.691566

1.00423+j0.658464

0.0165506+j0.0331013

电压损耗0.274733

0.503322+j0.306645

0.00332243+j0.00664486

0.000909771+j0.00181954

电压损耗0.074015

1.70423

0.0207828

0.00762669