IEEE 33节点配电网

本次上传的IEEE 33节点配电网是根据M. A. Kashem等发表的《A novel method for loss minimization in distribution networks》进行搭建MATLAB模型，此模型适应于进行配电网潮流计算，不适用于配电网故障定位，若是想要研究故障定位方向的，请找我之前发的一个资源。

clear all;
Bus=[100   ;
2100 60;
39040;
412080;
56030;
66020;
7 200 100 ;
8200  100 ;
9 60 20;
1060 20;
114530;
1260 35;
1360 35;
1412080;
1560 10;
1660 20;
1760 20;
1890 40;
1990 40;
2090 40;
2190 40;
2290 40; 
2390 40; 
24420 200;
25420 200 ;
2660 25;
2760 25;
2860 20;
29120 70;
30200600 ;
31150 70;
32210 100 ;
336040;
];
Branch=[1 1 2 0.09220.0407;
  2 2 3 0.49300.2511;
  3 3 4 0.36600.1864;
4 4 5 0.38110.1941;
5 5 6 0.81900.7070;
6 6 7 0.18720.6188;
7 7 8 0.71440.2351;
8 8 9 1.03000.7400;
9 9 101.04400.7400;
1010110.19660.065;
1111120.37440.1238;
1212131.46801.1550;
1313140.54160.7129;
1414150.59100.5260;
1515160.74630.5450;
1616171.28901.7210;
1717180.73200.5740;
182 190.16400.1565;
1919201.50421.3554;
2020210.40950.4784;
2121220.70890.9373;
223 230.45120.3083;
2323240.89800.7091;
2424250.89600.7011;
256 260.20300.1034;
2626270.28420.1447;
2727281.05900.9337;
2828290.80420.7006;
2929300.50750.2585;
3030310.97440.9630;
3131320.31050.3619;
3232330.34100.5302;
 ];
[busnumrow]=size（Bus）;
[branchnumrow]=size（Branch）;
soubus=Branch（:2）;
mobus=Branch（:3）;
Vbus=ones（busnum1）;
Vbus（:1）=12.66;
Vbus1=Vbus;
Ploss=zeros（busnum1）;
Qloss=zeros（busnum1）;
e=1;
k=0;
Branch1=Branch;
n=1;
while ~isempty（Branch1）
m=1;
[srow]=size（Branch1）;
while s>0
t=find（Branch1（:2）== Branch1（s3））;
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（busnum1）;
Q=zeros（busnum1）;
for s=1:branchnum
i=T1（s2）;
j=T1（s3）;
R=T1（s4）;
X=T1（s5）;
Pload=Bus（j2）;
Qload=Bus（j3）;
II=（（Pload+P（j））^2+（Qload+Q（j））^2）/（Vbus（j）^2*1000）;
Ploss（ij）=II*R;
Qloss（ij）=II*X;
P（ij）=Pload+Ploss（ij）+P（j）;
Q（ij）=Qload+Ploss（ij）+Q（j）;
P（i）=P（i）+P（ij）;
Q（i）=Q（i）+Q（ij）;
end;
for s=branchnum:-1:1
i=T1（s2）;
j=T1（s3）;
R=T1（s4）;
X=T1（s5）;
Vbus（j）=（Vbus（i）-（P（ij）*R+Q（ij）*X）/（Vbus（i）*1000））^2+（（P（ij）*X-Q（ij）*R）/（Vbus（i）*1000））^2;
Vbus（j）=sqrt（Vbus（j））;
end;
e=max（abs（Vbus1-Vbus））;
Vbus1=Vbus;
k=k+1;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      28504  2019-08-22 16:37  IEEE 33节点配电网\IEEE 33节点配电网.rar

     文件      32050  2019-11-26 15:38  IEEE 33节点配电网\ieee33.slx

     文件      27788  2019-08-21 17:20  IEEE 33节点配电网\ieee33.slx.r2012b

     文件     650066  2019-11-26 15:33  IEEE 33节点配电网\IEEE33配电网系统.pdf

     文件       2315  2019-08-21 17:27  IEEE 33节点配电网\mieee33.m

     目录          0  2019-11-26 15:40  IEEE 33节点配电网

----------- ---------  ---------- -----  ----

               740723                    6