前向回代法计算流程

前推回代法计算流程

要看懂前推回代法计算程序，汇报论述计算原理及计算流程。绘制计算流程框图。确定前推回代

支路次序(广度优先，或深度优先)，编写前推回代计算输入文献。进行时尚计算。

下列为节点配电网构造图及系统支路参数和系统负荷参数表C

1(

JkJkJ

76289

JOO

103H12

4')---------0

13

5(,：

图1-2节点P把电网构造图

表1系多、支路参数

支路

R(Q)X(Q.)B/2(TK)(S)

1~23.3673.6850.0

2~32.3562.M10.0

3~41.1451.280.0

4~54.5245.040.0

2~60.8561.140.0

6~72.7452.9650.0

2~83.7434.2510.0

8~92.2372.7560.0

3-104.144.6960.0

3~111.3281.7630.0

11-122.4362.8660.0

4~133.5213.9660.0

表2系统负荷参数

节点编号节点类型节点初始电压Pi(MVA)Qi(MVA)

(kV)

1根节点10.400

2中间节点10.00.03420.0301

3中间节点10.00.06930.0642

4中间节点10.00.08450.0763

5叶节点10.00.02950.0261

6中间节点10.00.04740.0409

7叶节点10.00.11760.0957

8中间节点10.00.09460.0857

9叶节点10.00.09160.0859

10叶节点10.00.02710.0229

11中间节点10.00.06960.0643

12叶节点10.00.06760.0579

13叶节点10.00.02980.0242

主程序清单:

[PQ,FT,RX]=casell4();%调用数据文献

NN=size(PQ,l);%节点数

NB=sizc(FTJ);%支路数数

V=PQ(:,1);%V初始电压相量

maxd=1

k=l

whilemaxd>0.0001

PQ2二PQ：%每一次迭代各节点H勺注入有功和无功相似

PL=O.O;

fori=l:NB

kf=FT(i,l);%前推始节点号

kt=FT(i,2);%前推终节点号

x=(PQ2(kf,2)A2+PQ2(kf,3)A2)/V(kD/V(kf);%计算沿线电流平方A

PQ1(i,l)=PQ2(kf,2)+RX(i,l)*x;%计算支路首端有功/MWRX(i,l)~R

PQ1(i,2)=PQ2(kf,3)+RX(i,2)*x;%计算沿支路的无功损耗/MvarRX(i,2)~X

PQ2(kt,2)=PQ2(kt,2)+PQl(i,l);%用PQ1去修正支路末端节点的有功P单位MW

PQ2(kt,3)=PQ2(kt,3)+PQKi,2);%用PQ1去修正支路末端节点的有功Q单位Mvar

PL=PL+RX(i,l)*x;

end

angle(l)=0.0;

for

kf=FT(i,2);%回代始节点号

k(=FT(i,l);%回代终节点号

dvl=(PQl(i』)*RX(i』)+PQl&2)*RX(i,2))/V(k。；%计算支路电压损耗的纵分量dvl

dv2=(PQl(i,l)*RX(i,2)-PQ1(i,2)*RX(i,l))/V(kf)；%计算支路电压损耗口勺横分量dv2

V2(kt)=sqrt((V(kf)-dvl)A2+dv2A2);%计算支路末端电压/kV

angle(kt)=angle(kf)+atand(dv2/(V(kf)-dv1));%计算支路

end

maxd=abs(V2(2)-V(2));

V2(1)=V(1);

fori=3:l:NN

ifabs(V2(i)-V(i))>maxd;

maxd=abs(V2(i)-V(i));

end

end

10.00.06760.0579

10.00.02980.0242

];

FT=[

%首端末端

54

134

43

103

1211

113

76

62

98

82

32

21

1；

RX=[

%RX

4.5245.04

3.5213.966

1.1451.28

4.144.696

2.4362.866

1.3281.763

2.7452.965

0.8561.14

2.2372.756

3.7434.251

2.3562.541

3.3673.685

］；

计算过程

maxd-

1

k=

1

maxd=

0.1780

k=

2

PQ1

0.02960.0262

0.02990.0243

0.14430.1272

0.02720.0230

0.06780.0581

0.13780.1230

0.11820.0964

0.16600.1378

0.09200.0863

0.18900.1748

0.38470.3439

0.80990.7260

Columns1through8

10.40009.88079.82209.96729.97349.97019.9390

9.8550

Columns9through13

9.95569.97809.96009.96689.9799

angle

Columns1through8

00.30110.39860.42110.43870.34210.3916

0.3878

Columns9through13

0.42250.41730.44440.47470.4400

PL=

0.0471

maxd=

0.1787

k=

3

PQ1=

0.02960.0262

0.02990.0243

0.14430.1272

0.02720.0230

0.06780.0581

0.13780.1230

0.11820.0964

0.16600.1378

0.09200.0863

0.18900.1749

0.38490.3442

0.81120.7274

V

Columns1through8

10.40009.87989.70049.78869.94059.85049.9089

9.7338

Columns9through13

9.81009.79969.78139.92679.9470

angle=

Columns1through8

00.30110.40110.42440.44210.34310.3929

0.3899

Columns9through13

0.42570.42040.44860.47910.4434

PL=

0.0484

maxd=

0.1793

k=

4

PQ1=

0.02960.0262

0.02990.0243

0.14440.1272

0.02720.0230

0.06780.0581

0.13790.1231

0.11820.0964

0.16610.1378

0.09200.0864

0.18910.1749

0.38510.3444

0.81150.7277

V=

Columns1through8

10.40009.87969.69949.66669.76149.84959.7884

9.7329

Columns9through13

9.68839.67779.65919.74749.7680

angle=

Columns1through8

00.30110.40110.42500.44330.34310.3942

0.3899

Columns9through13

0.42660.42090.44980.48140.4447

PL

0.0487

maxd=

0.1226

k=

5

PQ1=

0.02960.0262

0.02990.0243

0.14440.1272

0.02720.0230

0.06780.0581

0.13790.1231

0.11830.0964

0.16610.1379

0.09200.0864

0.18910.1750

0.38520.3444

0.81150.7278

V

Columns1through8

10.40009.87959.69919.66569.63919.84939.7875

9.7326

Columns9through13

9.68739.67679.65819.62489.6457

angle=

Columns1through8

00.30110.40110.42500.44380.34310.3942

0.3899

Columns9through13

0.42660.42090.44980.48220.4452

PL=

0.0487

maxd=

0.0010

k=

6

PQ1=

0.02960.0262

0.02990.0243

0.14440.1273

0.02720.0230

0.06780.0581

0.13790.1231

0.11830.0964

0.16610.1379

0.09200.0864

0.18910.1750

0.38520.3445

0.81160.7278

V=

Columns1through8

10.40009.87959.69919.66539.63819.84929.7873

9.7326

Columns9through13

9.68709.67649.65799.62389.6447

angle=

Columns1through8

00.30110.40110.42500.44380.34310.3942

0.3899

Columns9through13

0.42660.42090.44980.48230.4452

PL=

0.0488

maxd=

2.6C21e-004

k=

7

PQ1=

0.02960.0262

0.02990.0243

0.14440.1273

0.02720.0230

0.06780.0581

0.13790.1231

0.11830.0964

0.16610.1379

0.09200.0864

0.18910.1750

0.38520.3445

0.81160.7278

V=

Columns1through8

10.40009.87959.69919.66529.63789.84929.7872

9.7326

Columns9through13

9.68709.67649.65789.62359.6445

angle

Columns1through8

00.30110.40110.42500.44380.34310.3942

0.3899

Columns9through13

0.42660.42090.44980.48230.4452

PL

0.0488

maxd=

6.1C46e-005

k

8

PQ1

0.02960.0262

0.02990.0243

0.14440.1273

0.02720.0230

0.06780.0581

0.13790.1231

0.11830.0964

0.16610.1379

0.09200.0864

0.18910.1750

0.38520.3445

0.81160.7278

V=

Columns1through8

10.40009.87959.69919.66529.63779.84929.7872

9.7326

Columns9through13

9.68709.67649.65789.62359.6444

angle=