1*c4762a1bSJed Brownfunction mpc = case9 2*c4762a1bSJed Brown% CASE9 Power flow data for 9 bus, 3 generator case. 3*c4762a1bSJed Brown% MATPOWER case9 data 4*c4762a1bSJed Brown% http://www.pserc.cornell.edu//matpower/ 5*c4762a1bSJed Brown 6*c4762a1bSJed Brown% R. D. Zimmerman, C. E. Murillo-Sánchez, and R. J. Thomas, 7*c4762a1bSJed Brown% "MATPOWER: Steady-State Operations, Planning and Analysis Tools 8*c4762a1bSJed Brown% for Power Systems Research and Education," Power Systems, IEEE 9*c4762a1bSJed Brown% Transactions on, vol. 26, no. 1, pp. 12-19, Feb. 2011 10*c4762a1bSJed Brown 11*c4762a1bSJed Brown%%----- Power Flow Data -----%% 12*c4762a1bSJed Brown%% system MVA base 13*c4762a1bSJed Brownmpc.baseMVA = 100; 14*c4762a1bSJed Brown 15*c4762a1bSJed Brown%% bus data 16*c4762a1bSJed Brown% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin 17*c4762a1bSJed Brownmpc.bus = [ 18*c4762a1bSJed Brown 1 3 0 0 0 0 1 1 0 345 1 1.1 0.9; 19*c4762a1bSJed Brown 2 2 0 0 0 0 1 1 0 345 1 1.1 0.9; 20*c4762a1bSJed Brown 3 2 0 0 0 0 1 1 0 345 1 1.1 0.9; 21*c4762a1bSJed Brown 4 1 0 0 0 0 1 1 0 345 1 1.1 0.9; 22*c4762a1bSJed Brown 5 1 90 30 0 0 1 1 0 345 1 1.1 0.9; 23*c4762a1bSJed Brown 6 1 0 0 0 0 1 1 0 345 1 1.1 0.9; 24*c4762a1bSJed Brown 7 1 100 35 0 0 1 1 0 345 1 1.1 0.9; 25*c4762a1bSJed Brown 8 1 0 0 0 0 1 1 0 345 1 1.1 0.9; 26*c4762a1bSJed Brown 9 1 125 50 0 0 1 1 0 345 1 1.1 0.9; 27*c4762a1bSJed Brown]; 28*c4762a1bSJed Brown 29*c4762a1bSJed Brown%% generator data 30*c4762a1bSJed Brown% 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 31*c4762a1bSJed Brownmpc.gen = [ 32*c4762a1bSJed Brown 1 0 0 300 -300 1 100 1 250 10 0 0 0 0 0 0 0 0 0 0 0; 33*c4762a1bSJed Brown 2 163 0 300 -300 1 100 1 300 10 0 0 0 0 0 0 0 0 0 0 0; 34*c4762a1bSJed Brown 3 85 0 300 -300 1 100 1 270 10 0 0 0 0 0 0 0 0 0 0 0; 35*c4762a1bSJed Brown]; 36*c4762a1bSJed Brown 37*c4762a1bSJed Brown%% branch data 38*c4762a1bSJed Brown% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax 39*c4762a1bSJed Brownmpc.branch = [ 40*c4762a1bSJed Brown 1 4 0 0.0576 0 250 250 250 0 0 1 -360 360; 41*c4762a1bSJed Brown 4 5 0.017 0.092 0.158 250 250 250 0 0 1 -360 360; 42*c4762a1bSJed Brown 5 6 0.039 0.17 0.358 150 150 150 0 0 1 -360 360; 43*c4762a1bSJed Brown 3 6 0 0.0586 0 300 300 300 0 0 1 -360 360; 44*c4762a1bSJed Brown 6 7 0.0119 0.1008 0.209 150 150 150 0 0 1 -360 360; 45*c4762a1bSJed Brown 7 8 0.0085 0.072 0.149 250 250 250 0 0 1 -360 360; 46*c4762a1bSJed Brown 8 2 0 0.0625 0 250 250 250 0 0 1 -360 360; 47*c4762a1bSJed Brown 8 9 0.032 0.161 0.306 250 250 250 0 0 1 -360 360; 48*c4762a1bSJed Brown 9 4 0.01 0.085 0.176 250 250 250 0 0 1 -360 360; 49*c4762a1bSJed Brown]; 50*c4762a1bSJed Brown 51*c4762a1bSJed Brown%%----- OPF Data -----%% 52*c4762a1bSJed Brown%% area data 53*c4762a1bSJed Brown% area refbus 54*c4762a1bSJed Brownmpc.areas = [ 55*c4762a1bSJed Brown 1 5; 56*c4762a1bSJed Brown]; 57*c4762a1bSJed Brown 58*c4762a1bSJed Brown%% generator cost data 59*c4762a1bSJed Brown% 1 startup shutdown n x1 y1 ... xn yn 60*c4762a1bSJed Brown% 2 startup shutdown n c(n-1) ... c0 61*c4762a1bSJed Brownmpc.gencost = [ 62*c4762a1bSJed Brown 2 1500 0 3 0.11 5 150; 63*c4762a1bSJed Brown 2 2000 0 3 0.085 1.2 600; 64*c4762a1bSJed Brown 2 3000 0 3 0.1225 1 335; 65*c4762a1bSJed Brown]; 66