-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy patht_case9mod_opf.m
executable file
·65 lines (58 loc) · 1.91 KB
/
t_case9mod_opf.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
function [baseMVA, bus, gen, branch, areas, gencost] = t_case9mod_opf
%T_CASE9MOD_OPF Power flow data for 9 bus, 3 generator case, with OPF data.
% This is the same as t_case9_opf except for a line-flow limit of 100 MVA
% rather than 40 MVA for the line between buses 6 and 7. This less strict
% line-flow limit results in an exact semidefinite relaxation for this
% problem. Also all lines have minimum resistances of 1e-4 per unit.
% Please see CASEFORMAT for details on the case file format.
%% MATPOWER Case Format : Version 1
%%----- Power Flow Data -----%%
%% system MVA base
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 345 1 1.1 0.9;
2 2 0 0 0 0 1 1 0 345 1 1.1 0.9;
30 2 0 0 0 0 1 1 0 345 1 1.1 0.9;
4 1 0 0 0 0 1 1 0 345 1 1.1 0.9;
5 1 90 30 0 0 1 1 0 345 1 1.1 0.9;
6 1 0 0 0 0 1 1 0 345 1 1.1 0.9;
7 1 100 35 0 0 1 1 0 345 1 1.1 0.9;
8 1 0 0 0 0 1 1 0 345 1 1.1 0.9;
9 1 125 50 0 0 1 1 0 345 1 1.1 0.9;
];
%% generator data
% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin
gen = [
1 0 0 300 -300 1 100 1 250 90;
2 163 0 300 -300 1 100 1 300 10;
30 85 0 300 -300 1 100 1 270 10;
];
%% branch data
% fbus tbus r x b rateA rateB rateC ratio angle status
branch = [
1 4 0.0001 0.0576 0 0 250 250 0 0 1;
4 5 0.0170 0.0920 0.158 0 250 250 0 0 1;
5 6 0.0390 0.1700 0.358 150 150 150 0 0 1;
30 6 0.0001 0.0586 0 0 300 300 0 0 1;
6 7 0.0119 0.1008 0.209 100 150 150 0 0 1;
7 8 0.0085 0.0720 0.149 250 250 250 0 0 1;
8 2 0.0001 0.0625 0 250 250 250 0 0 1;
8 9 0.0320 0.1610 0.306 250 250 250 0 0 1;
9 4 0.0100 0.0850 0.176 250 250 250 0 0 1;
];
%%----- OPF Data -----%%
%% area data
% area refbus
areas = [
1 5;
];
%% generator cost data
% 1 startup shutdown n x1 y1 ... xn yn
% 2 startup shutdown n c(n-1) ... c0
gencost = [
1 0 0 4 0 0 100 2500 200 5500 250 7250;
2 0 0 2 24.035 -403.5 0 0 0 0 0 0;
1 0 0 3 0 0 200 3000 300 5000 0 0;
];