前推回代法计算流程Word下载.docx
《前推回代法计算流程Word下载.docx》由会员分享，可在线阅读，更多相关《前推回代法计算流程Word下载.docx（15页珍藏版）》请在冰豆网上搜索。
2.965
2~8
3.743
4.251
8~9
2.237
2.756
3~10
4.14
4.696
3~11
1.328
1.763
11~12
2.436
2.866
4~13
3.521
3.966
表2系统负荷参数
节点编号
节点类型
节点初始电压（kV）
Pi（MVA）
Qi（MVA）
1
根节点
10.4
2
中间节点
10.0
0.0342
0.0301
3
0.0693
0.0642
4
0.0845
0.0763
5
叶节点
0.0295
0.0261
6
0.0474
0.0409
7
0.1176
0.0957
8
0.0946
0.0857
9
0.0916
0.0859
10
0.0271
0.0229
11
0.0696
0.0643
12
0.0676
0.0579
13
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:
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））;
%计算支路
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
maxd
k=k+1
PQ1%潮流分布即支路首端潮流MVA
V=V2%节点电压模计算结果kV
angle%节点电压角度计算结果单位度
PL%网损单位MW
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=
0.1780
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=
00.30110.39860.42110.43870.34210.39160.3878
0.42250.41730.44440.47470.4400
PL=
0.0471
0.1787
3
0.18900.1749
0.38490.3442
0.81120.7274
10.40009.87989.70049.78869.94059.85049.90899.7338
9.81009.79969.78139.92679.9470
00.30110.40110.42440.44210.34310.39290.3899
0.42570.42040.44860.47910.4434
0.0484
0.1793
4
0.14440.1272
0.13790.1231
0.16610.1378
0.09200.0864
0.18910.1749
0.38510.3444
0.81150.7277
10.40009.87969.69949.66669.76149.84959.78849.7329
9.68839.67779.65919.74749.7680
00.30110.40110.42500.44330.34310.39420.3899
0.42660.42090.44980.48140.4447
0.0487
0.1226
5
0.11830.0964
0.16610.1379
0.18910.1750
0.38520.3444
0.81150.7278
10.40009.87959.69919.66569.63919.84939.78759.7326
9.68739.67679.65819.62489.6457
00.30110.40110.42500.44380.34310.39420.3899
0.42660.42090.44980.48220.4452
0.0010
6
0.14440.1273
0.38520.3445
0.81160.7278
10.40009.87959.69919.66539.63819.84929.78739.7326
9.68709.67649.65799.62389.6447
0.42660.42090.44980.48230.4452
0.0488
2.6021e-004
7
10.40009.87959.69919.66529.63789.84929.78729.7326
9.68709.67649.65789.62359.6445
6.1046e-005
8
10.40009.87959.69919.66529.63779.84929.78729.7326
9.68709.67649.65789.62359.6444
计算结果清单：
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
10.40009.87959.69919.66529.63779.84929.78729.73269.68709.67649.65789.62359.6444
00.30110.40110.42500.44380.34310.39420.38990.42660.42090.44980.48230.4452
PL=0.0488
参考文献
[1]何仰赞温增银《电力系统分析》．华中科技大学出版社.
[2]李维波.《MATLAB在电气工程中应用》．中国电力出版社.2007
升级会员 