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