前推回代法计算流程.docx
《前推回代法计算流程.docx》由会员分享，可在线阅读，更多相关《前推回代法计算流程.docx（17页珍藏版）》请在冰点文库上搜索。
前推回代法计算流程
前推回代法计算流程
要看懂前推回代法计算程序，报告叙述计算原理及计算流程。
绘制计算流程框图。
确定前推回代支路次序（广度优先，或深度优先），编写前推回代计算输入文件。
进行潮流计算。
下列为节点配电网结构图及系统支路参数和系统负荷参数表。
图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
升级会员 