前推回代法计算流程.docx
《前推回代法计算流程.docx》由会员分享，可在线阅读，更多相关《前推回代法计算流程.docx（16页珍藏版）》请在冰豆网上搜索。
前推回代法计算流程
前推回代法计算流程
要看懂前推回代法计算程序,报告叙述计算原理及计算流程。
绘制计算流程框图.确定前推回代支路次序（广度优先，或深度优先），编写前推回代计算输入文件。
进行潮流计算.
下列为节点配电网结构图及系统支路参数和系统负荷参数表。
图1-2节点配电网结构图
表1系统支路参数
支路
R（Ω）
X（Ω.）
B/2（TK）（S）
1～2
3.367
3。
685
0.0
2~3
2。
356
2。
541
0。
0
3~4
1。
145
1.28
0。
0
4~5
4。
524
5。
04
0.0
2～6
0.856
1。
14
0。
0
6～7
2.745
2.965
0。
0
2~8
3。
743
4。
251
0。
0
8~9
2.237
2。
756
0.0
3~10
4.14
4.696
0.0
3~11
1.328
1。
763
0.0
11～12
2。
436
2。
866
0。
0
4~13
3.521
3。
966
0.0
表2系统负荷参数
节点编号
节点类型
节点初始电压（kV）
Pi（MVA）
Qi（MVA）
1
根节点
10.4
0
0
2
中间节点
10。
0
0.0342
0。
0301
3
中间节点
10。
0
0.0693
0。
0642
4
中间节点
10.0
0.0845
0。
0763
5
叶节点
10.0
0.0295
0.0261
6
中间节点
10。
0
0.0474
0.0409
7
叶节点
10。
0
0。
1176
0.0957
8
中间节点
10。
0
0.0946
0。
0857
9
叶节点
10。
0
0。
0916
0。
0859
10
叶节点
10.0
0.0271
0.0229
11
中间节点
10.0
0。
0696
0。
0643
12
叶节点
10。
0
0。
0676
0。
0579
13
叶节点
10。
0
0.0298
0。
0242
主程序清单：
[PQ，FT,RX]=case114（）;％调用数据文件
NN=size（PQ,1）；％节点数
NB=size（FT，1）;％支路数数
V=PQ（:
1）；％V初始电压相量
maxd=1
k=1
whilemaxd>0.0001
PQ2=PQ；%每一次迭代各节点的注入有功和无功相同
PL=0.0；
fori=1：
NB
kf=FT（i,1）；%前推始节点号
kt=FT（i，2）；%前推终节点号
x=（PQ2（kf,2）^2+PQ2（kf,3）^2）/V（kf）/V（kf）;％计算沿线电流平方A
PQ1（i，1）=PQ2（kf，2）+RX（i,1）*x;％计算支路首端有功/MWRX（i，1）~R
PQ1（i,2）=PQ2（kf，3）+RX（i,2）＊x；%计算沿支路的无功损耗/MvarRX（i,2）～X
PQ2（kt,2）=PQ2（kt，2）+PQ1（i，1）；％用PQ1去修正支路末端节点的有功P单位MW
PQ2（kt,3）=PQ2（kt,3）+PQ1（i,2）；%用PQ1去修正支路末端节点的有功Q单位Mvar
PL=PL+RX（i,1）＊x；
end
angle
（1）=0。
0;
fori=NB:
—1：
1
kf=FT（i，2）；％回代始节点号
kt=FT（i,1）；%回代终节点号
dv1=（PQ1（i,1）*RX（i,1）+PQ1（i,2）*RX（i,2））/V（kf）；％计算支路电压损耗的纵分量dv1
dv2=（PQ1（i，1）＊RX（i，2）-PQ1（i，2）＊RX（i，1））/V（kf）；％计算支路电压损耗的横分量dv2
V2（kt）=sqrt（（V（kf）—dv1）^2+dv2^2）；％计算支路末端电压/kV
angle（kt）=angle（kf）+atand（dv2/（V（kf）-dv1））;%计算支路
end
maxd=abs（V2
（2）-V
（2））;
V2
（1）=V
（1）；
fori=3:
1:
NN
ifabs（V2（i）-V（i））〉maxd;
maxd=abs（V2（i）—V（i））;
end
end
maxd
k=k+1
PQ1％潮流分布即支路首端潮流MVA
V=V2％节点电压模计算结果kV
angle%节点电压角度计算结果单位度
PL％网损单位MW
end
clear
输入文件清单:
function［PQ，FT，RX]=case114（）
PQ=［
％节点电压有功无功
10.400
10.00。
03420.0301
10.00.06930。
0642
10.00。
08450。
0763
10.00。
02950.0261
10。
00。
04740。
0409
10.00.11760.0957
10。
00。
09460.0857
10.00.09160。
0859
10。
00.02710。
0229
10。
00。
06960.0643
10.00。
06760.0579
10.00。
02980.0242
］；
FT=［
%首端末端
54
134
43
103
1211
113
76
62
98
82
32
21
]；
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
V=
Columns1through8
10.40009.88079.82209.96729。
97349.97019。
93909.8550
Columns9through13
9。
95569.97809。
96009。
96689.9799
angle=
Columns1through8
00.30110。
39860.42110.43870.34210.39160.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。
90899。
7338
Columns9through13
9。
81009。
79969。
78139.92679.9470
angle=
Columns1through8
00。
30110.40110.42440.44210。
34310.39290.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。
78849.7329
Columns9through13
9.68839。
67779。
65919。
74749.7680
angle=
Columns1through8
00.30110。
40110。
42500。
44330。
34310.39420.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.78759。
7326
Columns9through13
9.68739.67679.65819。
62489。
6457
angle=
Columns1through8
00。
30110.40110。
42500.44380.34310。
39420。
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。
78739.7326
Columns9through13
9.68709。
67649。
65799.62389。
6447
angle=
Columns1through8
00。
30110.40110.42500。
44380。
34310.39420.3899
Columns9through13
0.42660。
42090.44980。
48230.4452
PL=
0.0488
maxd=
2.6021e-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.78729。
7326
Columns9through13
9。
68709。
67649.65789。
62359。
6445
angle=
Columns1through8
00。
30110。
40110.42500。
44380。
34310。
39420.3899
Columns9through13
0。
42660。
42090。
44980。
48230。
4452
PL=
0。
0488
maxd=
6.1046e—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。
78729。
7326
Columns9through13
9.68709.67649。
65789.62359.6444
angle=
Columns1through8
00.30110。
40110.42500。
44380.34310.39420.3899
Columns9through13
0.42660。
42090。
44980.48230.4452
PL=
0.0488
计算结果清单:
maxd=1
k=1
maxd=0。
1780
k=2
maxd=0。
1787
k=3
maxd=0.1793
k=4
maxd=0.1226
k=5
maxd=0.0010
k=6
maxd=2。
6021e-04
k=7
maxd=6.1046e-05
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=
10。
40009.87959。
69919。
66529。
63779。
84929。
78729。
73269.68709.67649.65789.62359。
6444
angle=
00。
30110.40110。
42500.44380。
34310。
39420.38990.42660.42090.44980。
48230.4452
PL=0.0488
参考文献
[1]何仰赞温增银《电力系统分析》．华中科技大学出版社.
[2］李维波.《MATLAB在电气工程中应用》．中国电力出版社.2007
升级会员 