原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
IEEE33节点matlab潮流程序
Bus [1,0,0 ;
2,100, 60;
3,90,40;
4,120,80;
5,60,30;
6,60,20;
7, 200, 100 ;
8,200, 100 ;
9, 60, 20;
10,60, 20;
11,45,30;
12,60, 35;
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 ;
while s 0
t find Branch1 :,2 Branch1 s,3 ;
if isempty t
T1 n,: Branch1 s,: ;
n n+1;
else
T2 m,: Branch1 s,: ;
m m+1;
end;
s s-1;
end;
Branch1 T2;
T2 [];
end;
while e 1.0e-05
P zeros busnum,1 ;
Q zeros busnum,1 ;
for s 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 ;
end;
for s branchnum:-1:1
i T1 s,2 ;
j T1 s,3 ;
R T1 s,4 ;
X T1 s,5 ;
Vbu
文档评论（0）