0001 function mpc = case33bw
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013 mpc.version = '2';
0014
0015
0016
0017 mpc.baseMVA = 100;
0018
0019
0020
0021 mpc.bus = [
0022 1 3 0 0 0 0 1 1 0 12.66 1 1 1;
0023 2 1 100 60 0 0 1 1 0 12.66 1 1.1 0.9;
0024 3 1 90 40 0 0 1 1 0 12.66 1 1.1 0.9;
0025 4 1 120 80 0 0 1 1 0 12.66 1 1.1 0.9;
0026 5 1 60 30 0 0 1 1 0 12.66 1 1.1 0.9;
0027 6 1 60 20 0 0 1 1 0 12.66 1 1.1 0.9;
0028 7 1 200 100 0 0 1 1 0 12.66 1 1.1 0.9;
0029 8 1 200 100 0 0 1 1 0 12.66 1 1.1 0.9;
0030 9 1 60 20 0 0 1 1 0 12.66 1 1.1 0.9;
0031 10 1 60 20 0 0 1 1 0 12.66 1 1.1 0.9;
0032 11 1 45 30 0 0 1 1 0 12.66 1 1.1 0.9;
0033 12 1 60 35 0 0 1 1 0 12.66 1 1.1 0.9;
0034 13 1 60 35 0 0 1 1 0 12.66 1 1.1 0.9;
0035 14 1 120 80 0 0 1 1 0 12.66 1 1.1 0.9;
0036 15 1 60 10 0 0 1 1 0 12.66 1 1.1 0.9;
0037 16 1 60 20 0 0 1 1 0 12.66 1 1.1 0.9;
0038 17 1 60 20 0 0 1 1 0 12.66 1 1.1 0.9;
0039 18 1 90 40 0 0 1 1 0 12.66 1 1.1 0.9;
0040 19 1 90 40 0 0 1 1 0 12.66 1 1.1 0.9;
0041 20 1 90 40 0 0 1 1 0 12.66 1 1.1 0.9;
0042 21 1 90 40 0 0 1 1 0 12.66 1 1.1 0.9;
0043 22 1 90 40 0 0 1 1 0 12.66 1 1.1 0.9;
0044 23 1 90 50 0 0 1 1 0 12.66 1 1.1 0.9;
0045 24 1 420 200 0 0 1 1 0 12.66 1 1.1 0.9;
0046 25 1 420 200 0 0 1 1 0 12.66 1 1.1 0.9;
0047 26 1 60 25 0 0 1 1 0 12.66 1 1.1 0.9;
0048 27 1 60 25 0 0 1 1 0 12.66 1 1.1 0.9;
0049 28 1 60 20 0 0 1 1 0 12.66 1 1.1 0.9;
0050 29 1 120 70 0 0 1 1 0 12.66 1 1.1 0.9;
0051 30 1 200 600 0 0 1 1 0 12.66 1 1.1 0.9;
0052 31 1 150 70 0 0 1 1 0 12.66 1 1.1 0.9;
0053 32 1 210 100 0 0 1 1 0 12.66 1 1.1 0.9;
0054 33 1 60 40 0 0 1 1 0 12.66 1 1.1 0.9;
0055 ];
0056
0057
0058
0059 mpc.gen = [
0060 1 0 0 10 -10 1 100 1 10 0 0 0 0 0 0 0 0 0 0 0 0;
0061 ];
0062
0063
0064
0065 mpc.branch = [
0066 1 2 0.0922 0.0470 0 0 0 0 0 0 1 -360 360;
0067 2 3 0.4930 0.2511 0 0 0 0 0 0 1 -360 360;
0068 3 4 0.3660 0.1864 0 0 0 0 0 0 1 -360 360;
0069 4 5 0.3811 0.1941 0 0 0 0 0 0 1 -360 360;
0070 5 6 0.8190 0.7070 0 0 0 0 0 0 1 -360 360;
0071 6 7 0.1872 0.6188 0 0 0 0 0 0 1 -360 360;
0072 7 8 0.7114 0.2351 0 0 0 0 0 0 1 -360 360;
0073 8 9 1.0300 0.7400 0 0 0 0 0 0 1 -360 360;
0074 9 10 1.0440 0.7400 0 0 0 0 0 0 1 -360 360;
0075 10 11 0.1966 0.0650 0 0 0 0 0 0 1 -360 360;
0076 11 12 0.3744 0.1238 0 0 0 0 0 0 1 -360 360;
0077 12 13 1.4680 1.1550 0 0 0 0 0 0 1 -360 360;
0078 13 14 0.5416 0.7129 0 0 0 0 0 0 1 -360 360;
0079 14 15 0.5910 0.5260 0 0 0 0 0 0 1 -360 360;
0080 15 16 0.7463 0.5450 0 0 0 0 0 0 1 -360 360;
0081 16 17 1.2890 1.7210 0 0 0 0 0 0 1 -360 360;
0082 17 18 0.7320 0.5740 0 0 0 0 0 0 1 -360 360;
0083 2 19 0.1640 0.1565 0 0 0 0 0 0 1 -360 360;
0084 19 20 1.5042 1.3554 0 0 0 0 0 0 1 -360 360;
0085 20 21 0.4095 0.4784 0 0 0 0 0 0 1 -360 360;
0086 21 22 0.7089 0.9373 0 0 0 0 0 0 1 -360 360;
0087 3 23 0.4512 0.3083 0 0 0 0 0 0 1 -360 360;
0088 23 24 0.8980 0.7091 0 0 0 0 0 0 1 -360 360;
0089 24 25 0.8960 0.7011 0 0 0 0 0 0 1 -360 360;
0090 6 26 0.2030 0.1034 0 0 0 0 0 0 1 -360 360;
0091 26 27 0.2842 0.1447 0 0 0 0 0 0 1 -360 360;
0092 27 28 1.0590 0.9337 0 0 0 0 0 0 1 -360 360;
0093 28 29 0.8042 0.7006 0 0 0 0 0 0 1 -360 360;
0094 29 30 0.5075 0.2585 0 0 0 0 0 0 1 -360 360;
0095 30 31 0.9744 0.9630 0 0 0 0 0 0 1 -360 360;
0096 31 32 0.3105 0.3619 0 0 0 0 0 0 1 -360 360;
0097 32 33 0.3410 0.5302 0 0 0 0 0 0 1 -360 360;
0098 21 8 2.0000 2.0000 0 0 0 0 0 0 0 -360 360;
0099 9 15 2.0000 2.0000 0 0 0 0 0 0 0 -360 360;
0100 12 22 2.0000 2.0000 0 0 0 0 0 0 0 -360 360;
0101 18 33 0.5000 0.5000 0 0 0 0 0 0 0 -360 360;
0102 25 29 0.5000 0.5000 0 0 0 0 0 0 0 -360 360;
0103 ];
0104
0105
0106
0107
0108
0109 mpc.gencost = [
0110 2 0 0 3 0 20 0;
0111 ];
0112
0113
0114
0115 [PQ, PV, REF, NONE, BUS_I, BUS_TYPE, PD, QD, GS, BS, BUS_AREA, VM, ...
0116 VA, BASE_KV, ZONE, VMAX, VMIN, LAM_P, LAM_Q, MU_VMAX, MU_VMIN] = idx_bus;
0117 [F_BUS, T_BUS, BR_R, BR_X, BR_B, RATE_A, RATE_B, RATE_C, ...
0118 TAP, SHIFT, BR_STATUS, PF, QF, PT, QT, MU_SF, MU_ST, ...
0119 ANGMIN, ANGMAX, MU_ANGMIN, MU_ANGMAX] = idx_brch;
0120 Vbase = mpc.bus(1, BASE_KV) * 1e3;
0121 Sbase = mpc.baseMVA * 1e6;
0122 mpc.branch(:, [BR_R BR_X]) = mpc.branch(:, [BR_R BR_X]) / (Vbase^2 / Sbase);
0123
0124
0125 mpc.bus(:, [PD, QD]) = mpc.bus(:, [PD, QD]) / 1e3;