-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathcase5_risk_edited.m
More file actions
81 lines (74 loc) · 2.42 KB
/
Copy pathcase5_risk_edited.m
File metadata and controls
81 lines (74 loc) · 2.42 KB
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
% Case5 with loads only on generator buses.
% Solution should have no active branches, and buses 1,5 are inactive along with gens 1,5
function mpc = case5
mpc.version = '2';
mpc.baseMVA = 100.0;
% risk_weight
mpc.risk_weight = .1
%% area data
% area refbus
mpc.areas = [
1 4;
];
%% bus data
% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin
mpc.bus = [
1 2 0.0 0.0 0.0 0.0 1 1.00000 0.00000 230.0 1 1.10000 0.90000;
2 1 300.0 98.61 0.0 0.0 1 1.00000 0.00000 230.0 1 1.10000 0.90000;
3 2 300.0 98.61 0.0 0.0 1 1.00000 0.00000 230.0 1 1.10000 0.90000;
4 3 400.0 131.47 0.0 0.0 1 1.00000 0.00000 230.0 1 1.10000 0.90000;
5 2 0.0 0.0 0.0 0.0 1 1.00000 0.00000 230.0 1 1.10000 0.90000;
];
%column_names% power_risk base_risk
mpc.bus_risk = [
0.0 0.0; %1
1.0 0.0; %2
1.0 0.0; %3
0.0 0.0; %4
0.0 0.0; %5
];
%% generator data
% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin
mpc.gen = [
2 20.0 0.0 30.0 -30.0 1.0 100.0 1 40.0 0.0;
2 85.0 0.0 127.5 -127.5 1.0 100.0 1 300.0 0.0;
3 260.0 0.0 390.0 -390.0 1.0 100.0 1 300.0 0.0;
4 100.0 0.0 150.0 -150.0 1.0 100.0 1 400.0 0.0;
5 300.0 0.0 450.0 -450.0 1.0 100.0 1 600.0 0.0;
];
%column_names% power_risk base_risk
mpc.gen_risk = [
0.0 0.0;
0.0 0.0;
0.0 0.0;
0.0 0.0;
0.0 0.0;
];
%% generator cost data
% 2 startup shutdown n c(n-1) ... c0
mpc.gencost = [
2 0.0 0.0 3 0.000000 14.000000 0.000000;
2 0.0 0.0 3 0.000000 15.000000 0.000000;
2 0.0 0.0 3 0.000000 30.000000 0.000000;
2 0.0 0.0 3 0.000000 40.000000 0.000000;
2 0.0 0.0 3 0.000000 10.000000 0.000000;
];
%% branch data
% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax
mpc.branch = [
1 2 0.00281 0.0281 0.00712 400.0 400.0 400.0 0.0 0.0 1 -30.0 30.0;
1 4 0.00304 0.0304 0.00658 426 426 426 0.0 0.0 1 -30.0 30.0;
1 5 0.00064 0.0064 0.03126 426 426 426 0.0 0.0 1 -30.0 30.0;
2 3 0.00108 0.0108 0.01852 426 426 426 0.0 0.0 1 -30.0 30.0;
3 4 0.00297 0.0297 0.00674 426 426 426 0.0 0.0 1 -30.0 30.0;
4 5 0.00297 0.0297 0.00674 240.0 240.0 240.0 0.0 0.0 1 -30.0 30.0;
];
%column_names% power_risk base_risk
mpc.branch_risk = [
0.0 0.0;
0.0 0.0;
1.0 0.0;
1.0 0.0;
1.0 0.0;
0.0 0.0;
];