function [baseMVA,bus,gen,branch] = case33
baseMVA = 100;
%%bus data
% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin
bus = [
1 3 0 0 0 0 1 1 0 12.66 1 1.05 0.95;
2 1 0.100 0.060 0 0 1 1 0 12.66 1 1.05 0.95;
3 1 0.090 0.040 0 0 1 1 0 12.66 1 1.05 0.95;
4 1 0.120 0.080 0 0 1 1 0 12.66 1 1.05 0.95;
5 1 0.060 0.030 0 0 1 1 0 12.66 1 1.05 0.95;
6 1 0.060 0.020 0 0 1 1 0 12.66 1 1.05 0.95;
7 1 0.200 0.100 0 0 1 1 0 12.66 1 1.05 0.95;
8 1 0.200 0.100 0 0 1 1 0 12.66 1 1.05 0.95;
9 1 0.060 0.020 0 0 1 1 0 12.66 1 1.05 0.95;
10 1 0.060 0.020 0 0 1 1 0 12.66 1 1.05 0.95;
11 1 0.045 0.030 0 0 1 1 0 12.66 1 1.05 0.95;
12 1 0.060 0.035 0 0 1 1 0 12.66 1 1.05 0.95;
13 1 0.060 0.035 0 0 1 1 0 12.66 1 1.05 0.95;
14 1 0.120 0.080 0 0 1 1 0 12.66 1 1.05 0.95;
15 1 0.060 0.010 0 0 1 1 0 12.66 1 1.05 0.95;
16 1 0.060 0.020 0 0 1 1 0 12.66 1 1.05 0.95;
17 1 0.060 0.020 0 0 1 1 0 12.66 1 1.05 0.95;
18 1 0.090 0.040 0 0 1 1 0 12.66 1 1.05 0.95;
19 1 0.090 0.040 0 0 1 1 0 12.66 1 1.05 0.95;
20 1 0.090 0.040 0 0 1 1 0 12.66 1 1.05 0.95;
21 1 0.090 0.040 0 0 1 1 0 12.66 1 1.05 0.95;
22 1 0.090 0.040 0 0 1 1 0 12.66 1 1.05 0.95;
23 1 0.090 0.050 0 0 1 1 0 12.66 1 1.05 0.95;
24 1 0.420 0.200 0 0 1 1 0 12.66 1 1.05 0.95;
25 1 0.420 0.200 0 0 1 1 0 12.66 1 1.05 0.95;
26 1 0.060 0.025 0 0 1 1 0 12.66 1 1.05 0.95;
27 1 0.060 0.025 0 0 1 1 0 12.66 1 1.05 0.95;
28 1 0.060 0.020 0 0 1 1 0 12.66 1 1.05 0.95;
29 1 0.120 0.070 0 0 1 1 0 12.66 1 1.05 0.95;
30 1 0.200 0.100 0 0 1 1 0 12.66 1 1.05 0.95;
31 1 0.150 0.070 0 0 1 1 0 12.66 1 1.05 0.95;
32 1 0.210 0.100 0 0 1 1 0 12.66 1 1.05 0.95;
33 1 0.060 0.040 0 0 1 1 0 12.66 1 1.05 0.95;
];
%%generator data
% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin
gen = [
1 0 0 60 -60 1 100 1 100 0
% 17 p(1) p(1)*0.8 60 -60 1 100 1 100 0
% 21 p(2) p(2)*0.8 60 -60 1 100 1 100 0
% 24 p(3) p(3)*0.8 60 -60 1 100 1 100 0
% 28 p(4) p(4)*0.8 60 -60 1 100 1 100 0
% 31 p(5) p(5)*0.8 60 -60 1 100 1 100 0
];
%%branch data
% fbus tbus r x b rateA rateB rateC ratio angle status
branch = [
1 2 0.0922 0.0470 0 130 130 130 0 0 1;
2 3 0.4930 0.2512 0 130 130 130 0 0 1;
3 4 0.3661 0.1864 0 130 130 130 0 0 1;
4 5 0.3811 0.1941 0 130 130 130 0 0 1;
5 6 0.8190 0.7070 0 130 130 130 0 0 1;
6 7 0.1872 0.6188 0 130 130 130 0 0 1;
7 8 0.7115 0.2351 0 130 130 130 0 0 1;
8 9 1.0299 0.7400 0 130 130 130 0 0 1;
9 10 1.0440 0.7400 0 130 130 130 0 0 1;
10 11 0.1967 0.0651 0 130 130 130 0 0 1;
11 12 0.3744 0.1298 0 65 65 65 0 0 1;
12 13 1.4680 1.1549 0 32 32 32 0 0 1;
13 14 0.5416 0.7129 0 65 65 65 0 0 1;
14 15 0.5909 0.5260 0 65 65 65 0 0 1;
15 16 0.7462 0.5449 0 65 65 65 0 0 1;
16 17 1.2889 1.7210 0 65 65 65 0 0 1;
17 18 0.7320 0.5739 0 32 32 32 0 0 1;
2 19 0.1640 0.1565 0 32 32 32 0 0 1;
19 20 1.5042 1.3555 0 32 32 32 0 0 1;
20 21 0.4095 0.4784 0 16 16 16 0 0 1;
21 22 0.7089 0.9373 0 16 16 16 0 0 1;
3 23 0.4512 0.3084 0 16 16 16 0 0 1;
23 24 0.8980 0.7091 0 16 16 16 0 0 1;
24 25 0.8959 0.7071 0 32 32 32 0 0 1;
6 26 0.2031 0.1034 0 32 32 32 0 0 1;
26 27 0.2842 0.1447 0 32 32 32 0 0 1;
27 28 1.0589 0.9338 0 32 32 32 0 0 1;
28 29 0.8043 0.7006 0 32 32 32 0 0 1;
29 30 0.5074 0.2585 0 32 32 32 0 0 1;
30 31 0.9745 0.9629 0 16 16 16 0 0 1;
31 32 0.3105 0.3619 0 16 16 16 0 0 1;
32 33 0.3411 0.5302 0 16 16 16 0 0 1;
];
return;
评论2