-
Notifications
You must be signed in to change notification settings - Fork 0
/
WB5.m
executable file
·61 lines (51 loc) · 1.82 KB
/
WB5.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
function mpc = WB5
%% 5 bus case in MATPOWER format
% W. A. Bukhsh, Feb 2013
%% References:
% [1] W. A. Bukhsh, Andreas Grothey, Ken McKinnon, Paul trodden, "Local Solutions of Optimal Power Flow Problem"
% submitted to IEEE Transactions on Power Systems, 2013
% [2] W. A. Bukhsh, Andreas Grothey, Ken McKinnon, Paul trodden, "Local Solutions of Optimal Power Flow Problem"
% Technical Report ERGO, 2011
%% 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 0 345 1 1.05 0.95;
2 1 130 20 0 0 1 1.0 -10 345 1 1.05 0.95;
3 1 130 20 0 0 1 1.0 -20 345 1 1.05 0.95;
4 1 65 10 0 0 1 1.0 -135 345 1 1.05 0.95;
5 2 0 0 0 0 1 1.0 -140 345 1 1.05 0.95;
];
%% 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 500 50 1800 -30 1 100 1 5000 0 0 0 0 0 0 0 0 0 0 0 0;
5 0 0 1800 -30 1 100 1 5000 0 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 2 0.04 0.09 0.0 2500 2500 2500 0 0 1 -360 360;
1 3 0.05 0.10 0.0 2500 2500 2500 0 0 1 -360 360;
2 4 0.55 0.90 0.45 2500 2500 2500 0 0 1 -360 360;
3 5 0.55 0.90 0.45 2500 2500 2500 0 0 1 -360 360;
4 5 0.06 0.1 0.0 2500 2500 2500 0 0 1 -360 360;
2 3 0.07 0.09 0.0 2500 2500 2500 0 0 1 -360 360;
];
%%----- OPF Data -----%%
%% area data
% area refbus
mpc.areas = [
1 5;
];
%% generator cost data
% 1 startup shutdown n x1 y1 ... xn yn
% 2 startup shutdown n c(n-1) ... c0
mpc.gencost = [
2 2 0 3 0 4.00 0;
2 2 0 3 0 1.00 0;
];