前推回代法计算流程Word文档下载推荐.docx
《前推回代法计算流程Word文档下载推荐.docx》由会员分享，可在线阅读，更多相关《前推回代法计算流程Word文档下载推荐.docx（17页珍藏版）》请在冰点文库上搜索。
1、3２8
1、763
１1~12
2、4３6
2、866
４~13
3、521
３、９６6
表2 系统负荷参数
节点编号
节点类型
节点初始电压（kＶ）
Pi（MVA）
Qｉ（MVA）
１
根节点
10、4
０
2
中间节点
10、0
0、０342
0、0３01
３
１0、0
0、０693
0、０642
4
0、0845
0、０７63
5
叶节点
1０、0
0、０2９5
0、0261
6
０、047４
0、０４09
7
０、1１76
0、0957
8
0、0９46
0、0８57
9
0、0９１6
０、0859
10
0、０271
0、0229
11
１0、０
0、0696
0、0643
１２
0、0676
0、0579
1３
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
ｃ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０
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
3
０、０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
0、0484
maxｄ＝
0、1793
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７
0、0487
ｍaｘd=
0、1２2６
k＝
5
０、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=
ｍaxd=
0、0010
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
０、04８8
maxｄ =
２、6021e－０04
7
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
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２
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
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２
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=８
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
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
升级会员 