/
t_case3p_c.m
112 lines (97 loc) · 3.05 KB
/
t_case3p_c.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
function mpc = t_case3p_c
% t_case3p_c - Six bus hybrid test case, 2 single-phase buses, 4 3-phase buses.
%
% One bus is a hybrid PV bus (PV on single-phase side). Three phase bus
% solution should match t_case3p_a.
%% MATPOWER Case Format : Version 2
mpc.version = '2';
%%----- Power Flow Data -----%%
%% system MVA base
mpc.baseMVA = 100;
%% bus data
% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin
mpc.bus = [
1 3 0 0 0 0 1 1 0.20115203 12.47 1 1.1 0.9;
% 2 2 0 0 0 0 1 1 0 12.47 1 1.1 0.9;
% 3 2 0 0 0 0 1 1 0 12.47 1 1.1 0.9;
4 2 0 0 0 0 1 1 0 12.47 1 1.1 0.9;
% 5 1 90 30 0 0 1 1 0 12.47 1 1.1 0.9;
% 6 1 0 0 0 0 1 1 0 12.47 1 1.1 0.9;
% 7 1 100 35 0 0 1 1 0 12.47 1 1.1 0.9;
% 8 1 0 0 0 0 1 1 0 12.47 1 1.1 0.9;
% 9 1 125 50 0 0 1 1 0 12.47 1 1.1 0.9;
];
%% generator data
% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin Pc1 Pc2 Qc1min Qc1max Qc2min Qc2max ramp_agc ramp_10 ramp_30 ramp_q apf
mpc.gen = [
1 72.3 27.03 300 -300 1.0024291384 100 1 250 10 0 0 0 0 0 0 0 0 0 0 0;
4 0 0 300 -300 1 100 1 250 10 0 0 0 0 0 0 0 0 0 0 0;
% 2 163 6.54 300 -300 1.025 100 1 300 10 0 0 0 0 0 0 0 0 0 0 0;
% 3 85 -10.95 300 -300 1.025 100 1 270 10 0 0 0 0 0 0 0 0 0 0 0;
];
%% branch data
% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax
mpc.branch = [
1 4 0 0.0576 0 250 250 250 0 0 1 -360 360;
% 4 5 0.017 0.092 0.158 250 250 250 0 0 1 -360 360;
% 5 6 0.039 0.17 0.358 150 150 150 0 0 1 -360 360;
% 3 6 0 0.0586 0 300 300 300 0 0 1 -360 360;
% 6 7 0.0119 0.1008 0.209 150 150 150 0 0 1 -360 360;
% 7 8 0.0085 0.072 0.149 250 250 250 0 0 1 -360 360;
% 8 2 0 0.0625 0 250 250 250 0 0 1 -360 360;
% 8 9 0.032 0.161 0.306 250 250 250 0 0 1 -360 360;
% 9 4 0.01 0.085 0.176 250 250 250 0 0 1 -360 360;
];
%%----- OPF Data -----%%
%% generator cost data
% 1 startup shutdown n x1 y1 ... xn yn
% 2 startup shutdown n c(n-1) ... c0
mpc.gencost = [
2 1500 0 3 0.11 5 150;
2 1500 0 3 0.11 5 150;
% 2 2000 0 3 0.085 1.2 600;
% 2 3000 0 3 0.1225 1 335;
];
%%----- 3 Phase Model Data -----%%
%% system data
mpc.freq = 60; %% frequency, Hz
mpc.basekVA = 1000; %% system kVA base
%% bus data
% busid type basekV Vm1 Vm2 Vm3 Va1 Va2 Va3
mpc.bus3p = [
1 1 12.47 1 1 1 0 -120 120;
2 1 12.47 1 1 1 0 -120 120;
3 1 4.16 1 1 1 0 -120 120;
4 1 4.16 1 1 1 0 -120 120;
];
%% buslink data
% linkid busid bus3pid status
mpc.buslink = [
1 4 1 1;
];
%% branch data
% brid fbus tbus status lcid len
mpc.line3p = [
1 1 2 1 1 2000/5280;
2 3 4 1 1 2500/5280;
];
%% transformer
% xfid fbus tbus status R X basekVA basekV
mpc.xfmr3p = [
1 2 3 1 0.01 0.06 6000 12.47;
];
%% load
% ldid ldbus status Pd1 Pd2 Pd3 ldpf1 ldpf2 ldpf3
mpc.load3p = [
1 4 1 1275 1800 2375 0.85 0.9 0.95;
];
%% gen
% genid gbus status Vg1 Vg2 Vg3 Pg1 Pg2 Pg3 Qg1 Qg2 Qg3
mpc.gen3p = [
% 1 1 1 1 1 1 2000 2000 2000 0 0 0;
];
%% line construction
% lcid R11 R21 R31 R22 R32 R33 X11 X21 X31 X22 X32 X33 C11 C21 C31 C22 C32 C33
mpc.lc = [
1 0.457541 0.15594 0.153474 0.466617 0.157996 0.461462 1.078 0.501648 0.384909 1.04813 0.423624 1.06502 15.0671 -4.86241 -1.85323 15.875 -3.09098 14.3254
];