IEEE节点matlab潮流程序Word文档格式.doc
《IEEE节点matlab潮流程序Word文档格式.doc》由会员分享，可在线阅读，更多相关《IEEE节点matlab潮流程序Word文档格式.doc（3页珍藏版）》请在冰豆网上搜索。
13,60,35;
14,120,80;
15,60,10;
16,60,20;
17,60,20;
18,90,40;
19,90,40;
20,90,40;
21,90,40;
22,90,40;
23,90,40;
24,420,200;
25,420,200;
26,60,25;
27,60,25;
28,60,20;
29,120,70;
30,200,600;
31,150,70;
32,210,100;
33,60,40;
];
Branch=[1,1,2,0.0922,0.0407;
2,2,3,0.4930,0.2511;
3,3,4,0.3660,0.1864;
4,4,5,0.3811,0.1941;
5,5,6,0.8190,0.7070;
6,6,7,0.1872,0.6188;
7,7,8,0.7144,0.2351;
8,8,9,1.0300,0.7400;
9,9,10,1.0440,0.7400;
10,10,11,0.1966,0.065;
11,11,12,0.3744,0.1238;
12,12,13,1.4680,1.1550;
13,13,14,0.5416,0.7129;
14,14,15,0.5910,0.5260;
15,15,16,0.7463,0.5450;
16,16,17,1.2890,1.7210;
17,17,18,0.7320,0.5740;
18,2,19,0.1640,0.1565;
19,19,20,1.5042,1.3554;
20,20,21,0.4095,0.4784;
21,21,22,0.7089,0.9373;
22,3,23,0.4512,0.3083;
23,23,24,0.8980,0.7091;
24,24,25,0.8960,0.7011;
25,6,26,0.2030,0.1034;
26,26,27,0.2842,0.1447;
27,27,28,1.0590,0.9337;
28,28,29,0.8042,0.7006;
29,29,30,0.5075,0.2585;
30,30,31,0.9744,0.9630;
31,31,32,0.3105,0.3619;
32,32,33,0.3410,0.5302;
[busnum,row]=size（Bus）;
[branchnum,row]=size（Branch）;
soubus=Branch（:
2）;
mobus=Branch（:
3）;
Vbus=ones（busnum,1）;
Vbus（:
1）=12.66;
Vbus1=Vbus;
Ploss=zeros（busnum,1）;
Qloss=zeros（busnum,1）;
e=1;
k=0;
Branch1=Branch;
n=1;
while~isempty（Branch1）
m=1;
[s,row]=size（Branch1）;
whiles>
t=find（Branch1（:
2）==Branch1（s,3））;
ifisempty（t）
T1（n,:
）=Branch1（s,:
）;
n=n+1;
else
T2（m,:
m=m+1;
end;
s=s-1;
Branch1=T2;
T2=[];
whilee>
1.0e-05
P=zeros（busnum,1）;
Q=zeros（busnum,1）;
fors=1:
branchnum
i=T1（s,2）;
j=T1（s,3）;
R=T1（s,4）;
X=T1（s,5）;
Pload=Bus（j,2）;
Qload=Bus（j,3）;
II=（（Pload+P（j））^2+（Qload+Q（j））^2）/（Vbus（j）^2*1000）;
Ploss（i,j）=II*R;
Qloss（i,j）=II*X;
P（i,j）=Pload+Ploss（i,j）+P（j）;
Q（i,j）=Qload+Ploss（i,j）+Q（j）;
P（i）=P（i）+P（i,j）;
Q（i）=Q（i）+Q（i,j）;
fors=branchnum:
-1:
1
Vbus（j）=（Vbus（i）-（P（i,j）*R+Q（i,j）*X）/（Vbus（i）*1000））^2+（（P（i,j）*X-Q（i,j）*R）/（Vbus（i）*1000））^2;
Vbus（j）=sqrt（Vbus（j））;
e=max（abs（Vbus1-Vbus））;
k=k+1;
end
升级会员 