前推回代法计算流程.docx
《前推回代法计算流程.docx》由会员分享，可在线阅读，更多相关《前推回代法计算流程.docx（17页珍藏版）》请在冰豆网上搜索。
前推回代法计算流程
前推回代法计算流程
要瞧懂前推回代法计算程序，报告叙述计算原理及计算流程。
绘制计算流程框图。
确定前推回代支路次序（广度优先,或深度优先），编写前推回代计算输入文件。
进行潮流计算。
下列为节点配电网结构图及系统支路参数与系统负荷参数表。
图1-2节点配电网结构图
表１ 系统支路参数
支路
Ｒ（Ω）
Ｘ（Ω、）
B/２（ＴK）（S）
1~2
３、36７
３、６８5
0、0
2~３
2、３56
2、５41
0、0
３~4
1、145
1、2８
0、0
４~5
4、524
5、0４
0、０
2~6
0、8５6
1、14
0、0
6~7
2、74５
2、９6５
0、0
2～8
3、74３
4、25１
0、0
8~9
2、23７
２、756
0、0
３~10
4、14
4、６96
０、0
３~１１
1、3２8
1、763
0、０
１1~12
2、4３6
2、866
0、0
４~13
3、521
３、９６6
0、0
表2 系统负荷参数
节点编号
节点类型
节点初始电压（kＶ）
Pi（MVA）
Qｉ（MVA）
１
根节点
10、4
0
０
2
中间节点
10、0
0、０342
0、0３01
３
中间节点
１0、0
0、０693
0、０642
4
中间节点
１0、0
0、0845
0、０７63
5
叶节点
1０、0
0、０2９5
0、0261
6
中间节点
１0、0
０、047４
0、０４09
7
叶节点
１0、0
０、1１76
0、0957
8
中间节点
1０、0
0、0９46
0、0８57
9
叶节点
10、0
0、0９１6
０、0859
10
叶节点
1０、0
0、０271
0、0229
11
中间节点
１0、０
0、0696
0、0643
１２
叶节点
1０、0
0、0676
0、0579
1３
叶节点
１0、0
0、０29８
０、0242
主程序清单：
[PQ,FT,RX]＝case114（）；％调用数据文件
NN=sｉｚｅ（PQ，1）; %节点数
ＮB=size（ＦT,1）;％支路数数
Ｖ=ＰQ（:
1）;%V初始电压相量
maxｄ=１
ｋ=１
wｈilｅmaxd>0、0001
PQ２=PQ;%每一次迭代各节点得注入有功与无功相同
PＬ=0、0;
for ｉ=1：
NB
ｋf＝FT（i,1）;%前推始节点号
kｔ=ＦT（i,2）;%前推终节点号
x=（PQ2（kf，2）＾2+PQ2（ｋf,３）^2）/V（kｆ）/V（kｆ）; %计算沿线电流平方A
PQ1（i,１）=PQ2（kf,2）+RX（i，1）*x; %计算支路首端有功/MＷ RX（i,1）~R
PQ1（i,2）=PＱ2（kf,3）＋RX（i，２）*x; %计算沿支路得无功损耗/Mｖaｒ RX（ｉ，２）~Ｘ
PQ２（ｋｔ,2）=PＱ2（ｋt,2）+ＰQ1（i，1）； %用PQ1去修正支路末端节点得有功P单位MW
PQ２（kｔ,３）=ＰQ2（kt,3）+PQ1（i，２）; %用PQ１去修正支路末端节点得有功Ｑ单位Mｖar
ＰL＝PL+RＸ（i,1）*ｘ;
ｅnd
anｇｌe（１）＝０、0;
foｒi=ＮＢ：
-1:
1
kf=FT（i,2）; ％回代始节点号
kt=FT（i,1）; %回代终节点号
dv1=（PQ1（i,1）*ＲX（i，１）+PQ1（i,2）*ＲX（i,２））/V（kf）; %计算支路电压损耗得纵分量dv1
dv2=（ＰＱ1（i，1）＊RX（i,2）-PQ1（ｉ,2）＊ＲX（ｉ，1））/V（ｋf）；%计算支路电压损耗得横分量dｖ2
Ｖ２（kt）=sqrt（（Ｖ（ｋｆ）－dv1）^2＋dｖ2＾2）; %计算支路末端电压/kV
ａngle（ｋt）=anｇｌe（ｋｆ）＋atand（dv2/（Ｖ（kf）-dｖ1））； %计算支路
enｄ
mａｘｄ=ａbs（Ｖ２
（2）－V（２））;
V2
（1）=V（１）;
fｏｒi=3:
1:
NN
if ａｂｓ（Ｖ2（i）-V（i））>ｍaxd;
ｍaxd＝abs（V2（i）-Ｖ（i））;
enｄ
eｎd
ｍaxｄ
ｋ=k+1
PＱ1 %潮流分布即支路首端潮流MＶA
V=Ｖ2 ％节点电压模计算结果kV
anｇｌe ％节点电压角度计算结果单位度
PL %网损单位MW
ｅnd
ｃleａr
输入文件清单：
functｉｏn[ＰQ,ＦT,ＲX]=cａse1１４（）
PQ＝[
%节点电压 有功无功
１０、４ 0 ０
10、0 0、０３４20、0301
10、0 0、0６９3 ０、０642
1０、0 ０、0８45０、0763
10、０ ０、02950、0261
10、0 0、０47４0、0409
10、00、１176０、09５7
10、0 0、09460、085７
1０、0 ０、091６0、0859
１０、０ 0、0２71 0、0229
10、0 ０、0６96０、０６4３
１0、０ 0、067６ 0、０579
10、０ 0、０29８0、0242
]；
ＦT=［
%首端末端
5 4
13 4
４ ３
10 3
12 1１
11 3
76
６ 2
9８
8 2
3 2
２ 1
];
RX=［
％ＲＸ
4、５2４ 5、0４
3、52１ 3、９66
1、14５ 1、2８
4、１44、696
2、43６ 2、866
1、3２81、7６３
2、74５２、965
0、8561、14
２、２37 2、7５6
３、7４34、25１
2、356２、54１
3、367 3、６８5
];
计算过程
ｍaｘd＝
１
k =
1
mａxd=
0、178０
k =
2
PQ1=
0、029６ 0、02６2
0、０299 ０、02４３
0、14430、12７2
0、０272 ０、02３0
0、０６78 0、０５81
0、１37８ ０、123０
0、１１82 0、０964
0、16６0 0、1３78
0、092０ ０、0８63
0、1890 ０、1748
0、３８47 ０、343９
0、8099０、7260
V=
Coｌｕmns1ｔhｒough８
１0、40009、８8079、82209、9672 ９、9734９、９7019、9３90９、8550
Cｏlumnｓ9ｔhrough１3
9、9556 9、9７８０ 9、9600 9、9668 9、9７99
ａｎgｌe =
Coｌumns1through8
00、３０1１0、３９860、421１０、43８7 0、34210、391６0、3８78
Colｕmns ９ tｈrough 13
０、４22５ ０、4１73 ０、４4４４０、４７４７ 0、４400
PL=
0、0４７１
mａxd＝
0、1787
k =
3
PQ1=
０、０296 ０、0262
０、02990、024３
0、14430、1272
0、02720、023０
0、0678 0、０581
0、1378 0、1230
0、1182 0、09６4
0、16６0 0、13７8
0、09200、08６３
0、1890 0、１749
0、3849０、３442
0、８1１20、７２74
Ｖ＝
Coluｍｎs１ ｔhrough8
10、４0０0 9、８798 ９、7004 9、78８6 ９、９４０59、８５04 9、９０899、733８
Columnｓ９throuｇh 13
9、8100 9、79９69、7８139、9267 9、９４7０
angle =
Columns1through８
00、30110、４011 0、4２44０、44210、3431 ０、3929 0、3899
Coｌumｎｓ9 ｔhrough １3
0、4257０、4２0４ 0、４4８６ ０、4791０、4434
PL=
0、0484
maxｄ＝
0、1793
k =
4
PQ1 =
0、0２960、０２６２
0、０2990、０2４３
0、14４4 0、１27２
0、027２ 0、０230
0、06７80、0５8１
0、13790、123１
0、11８2 0、０９64
0、１６61 0、13７8
0、09200、０864
０、1891 0、17４9
0、３851 ０、3４44
0、8115 0、727７
Ｖ=
Columnｓ１through 8
10、40０0９、879６ ９、6９949、６666 9、７614 9、８495 9、78849、7329
Columns9throｕgh１３
９、６88３９、６7779、65９1 9、７4７4 ９、7６8０
angle＝
Ｃolumｎｓ1ｔｈrｏugh8
0 ０、3０11 0、40１10、42５０ 0、4４3３ 0、3431０、3９42 0、3899
Coｌumnｓ9 thrｏugh１3
0、426６0、420９ 0、４498 0、４8１40、444７
PL=
0、0487
ｍaｘd=
0、1２2６
k＝
5
PQ1=
０、02960、０262
0、0299 0、0243
0、１444 0、1２72
０、0272 0、0230
0、0678 0、0５81
０、1３79 0、12３1
0、118３ ０、09６4
0、１6610、1379
0、0920０、0864
０、１８９10、1７50
０、３８５２ ０、344４
0、81150、７2７８
V ＝
Ｃｏｌｕmnｓ 1ｔhrｏuｇh8
10、４000 9、87９5 ９、6991 9、6656 9、639１ ９、8493 9、7875 9、7326
Ｃｏlｕmnｓ 9through13
9、６87３ 9、676７ 9、658１ 9、62４8 9、6457
anglｅ ＝
Cｏlumns1throuｇh 8
０ 0、3011 ０、4011 0、４2500、4438 0、34310、39４2 0、3８99
Columｎs9ｔhrouｇh 13
0、4266 ０、４２０9 0、44９80、４8２20、4４52
ＰL=
0、0487
ｍaxd=
0、0010
k =
6
ＰQ1 ＝
0、02960、0262
０、0２99 ０、0243
0、1４440、1273
0、０27２0、023０
０、0678 0、0581
０、1３７９ ０、1２31
0、1１８3 0、0964
0、16610、１37９
０、0９2０ 0、0８64
0、1891 0、17５０
0、3852 0、3４４5
０、811６0、7278
V =
Ｃolumns１tｈrｏugｈ8
1０、40０0 ９、879５９、6991９、665３ 9、６38１ 9、849２9、78７3 9、732６
Cｏluｍnｓ９tｈrｏｕgh１3
９、687０ ９、676４9、6５79 9、6２38９、644７
anｇlｅ =
Ｃｏｌｕmns1tｈrouｇh8
00、3０110、401１0、４25０0、443８ 0、3431 0、3９42 0、3899
Colｕｍnｓ 9ｔhrougｈ １3
0、４2６6 0、4209 0、4498 0、４823 0、44５2
PL=
０、04８8
maxｄ =
２、6021e－０04
k =
7
PQ1=
0、0296 0、02６2
0、０29９0、０２43
0、14４４ 0、1２7３
0、0272 0、０２30
0、０6７8 0、058１
0、1３７9 0、1２31
0、1１8３0、0９64
0、1661 0、1379
０、0９2０ ０、08６4
0、18910、175０
0、3８５2 ０、3445
０、８11６ 0、7278
V=
Coｌumns1thrｏugh８
1０、400０ ９、8795 ９、6９９1 9、66５2９、６378 9、８492 9、78729、７326
Columns ９tｈｒoｕgｈ１3
9、68７0 9、6764 ９、6５7８９、6235 9、644５
ａnｇle =
Colｕmns 1 ｔｈrｏuｇh８
0 ０、3011 0、40１1 ０、４25０ 0、4438 0、３43１ 0、３9420、38９9
Coｌumｎs9throuｇh 13
0、42660、42０9 0、4４９8 0、4823 0、445２
PL=
0、0488
maxd=
6、1０46e－005
k=
8
PＱ１＝
0、0２96 0、0２62
0、02990、０２43
0、1４４４ 0、12７３
０、027２ 0、02３０
0、06７80、０581
0、1３7９ 0、1231
0、11830、0964
０、166１ 0、１379
0、0920 ０、０８64
0、１８91 0、17５0
0、３852 0、3445
0、８116０、72７8
V=
Columｎs1 ｔhrougｈ8
10、400０9、8７9５ 9、６991 9、6652 9、６3７79、８492 ９、7８７2 9、7３26
Cｏlumｎs 9 thrｏｕｇh13
9、68７０ 9、6764 9、6578９、6２3５ 9、6444
angｌe =
Cｏlumns1tｈｒoｕgh 8
０ ０、3011 0、4011 0、4２50０、4438 0、３431 0、3942 0、38９9
Coluｍns9tｈrｏugh1３
0、４２66 ０、42090、449８ 0、4８２３ ０、445２
ＰL=
0、0４８８
计算结果清单:
maxd=1
ｋ=1
ｍaxd= ０、１７80
k =2
ｍａｘｄ=0、1787
k＝3
mａxｄ=0、1793
k =4
mａxd＝0、122６
k=5
ｍaｘd=0、0010
k= 6
mａｘd=2、6０21e－０4
k =７
mａxd=６、104６e－05
k=８
PQ1=
0、0２96 0、０262
0、０２9９ 0、0２43
0、14４４ ０、1273
0、０2７2０、0230
0、067８ 0、0５８１
0、13790、1231
0、11８30、096４
０、１6610、１379
０、０920 0、0864
0、18910、1７50
0、38520、３4４５
0、8116 0、7２78
V =
1０、4000 9、８7９５ ９、69９1 ９、6６52 ９、6377 9、８４92 ９、78７29、7326９、68７0 9、6７64 9、65７8 ９、6235 ９、６444
ａngle＝
0 0、3011 ０、4０1１ 0、42５00、443８ 0、3431 0、３942 ０、３8990、４２６60、４２0９ 0、4４９8 0、4８2３ 0、4452
ＰＬ=0、0488
参考文献
[１] 何仰赞 温增银《电力系统分析》.华中科技大学出版社、
ﻩ[2]李维波、《MATLAB在电气工程中应用》.中国电力出版社、20０7
升级会员 