% Fifteen-Units Test System with Ramp-Rate Limits, Prohibited-
% Operating Zones, and Valve-Point Loading Effects
% This data is entered by: Ali R. Al-Roomi in the 7th of May, 2016
% Refer to http://www.al-roomi.org for more details about this system

% The system data 
n=15; % number of units
Load=2630; % Load Demand (in MW)
Pmin=[150;150;20;20;150;135;135;60;25;25;20;20;25;15;15];% for P1, P2, ..., Pn (in MW)
Pmax=[455;455;130;130;470;460;465;300;162;160;80;80;85;55;55]; % for P1, P2, ..., Pn (in MW)

% Fuel-Cost Curves Coefficients
% Curves constants
a=[671;574;374;374;461;630;548;227;173;175;186;230;225;309;323]; % for unit 1, 2, ..., n
b=[10.1;10.2;8.8;8.8;10.4;10.1;9.8;11.2;11.2;10.7;10.2;9.9;13.1;12.1;12.4]; % for unit 1, 2, ..., n
c=[0.000299;0.000183;0.001126;0.001126;0.000205;0.000301;0.000364;0.000338;0.000807;0.001203;0.003586;0.005513;0.000371;0.001929;0.004447]; % for unit 1, 2, ..., n

% Valve-Point Loading Effects Coefficients
d=[100;100;100;150;120;100;200;200;200;200;200;200;300;300;300]; % for unit 1, 2, ..., n
e=[0.084;0.084;0.084;0.063;0.077;0.084;0.042;0.042;0.042;0.042;0.042;0.042;0.035;0.035;0.035]; % for unit 1, 2, ..., n

% Ramp-Rate Limits
Pnow=[400;300;105;100;90;400;350;95;105;110;60;40;30;20;20]; % for unit 1, 2, ..., n
UR=[80;80;130;130;80;80;80;65;60;60;80;80;80;55;55]; % for unit 1, 2, ..., n
DR=[120;120;130;130;120;120;120;100;100;100;80;80;80;55;55]; % for unit 1, 2, ..., n

% Prohibited Operating Zones
M=NaN(15,6); % create m by n empty matrix
M(2,1)=185;M(2,2)=225;M(2,3)=305;M(2,4)=335;M(2,5)=420;M(2,6)=450;
M(5,1)=180;M(5,2)=200;M(5,3)=305;M(5,4)=335;M(5,5)=390;M(5,6)=420;
M(6,1)=230;M(6,2)=225;M(6,3)=365;M(6,4)=395;M(6,5)=430;M(6,6)=455;
M(12,1)=30;M(12,2)=40;M(12,3)=55;M(12,4)=65;
PZL1=M(:,1); % 1st Lower Limit of Prohibited Zone 1
PZU1=M(:,2); % 1st Upper Limit of Prohibited Zone 1
PZL2=M(:,3); % 2nd Lower Limit of Prohibited Zone 2
PZU2=M(:,4); % 2nd Upper Limit of Prohibited Zone 2
PZL3=M(:,5); % 3rd Lower Limit of Prohibited Zone 3
PZU3=M(:,6); % 3rd Upper Limit of Prohibited Zone 3

% B-Coefficients
Bij=[0.0017,0.0012,0.0007,-0.0001,-0.0005,-0.0002;...
0.0012,0.0014,0.0009,0.0001,-0.0006,-0.0001;...
0.0007,0.0009,0.0031,0.0000,-0.0010,-0.0006;...
-0.0001,0.0001,0.0000,0.0024,-0.0006,-0.0008;...
-0.0005,-0.0006,-0.0010,-0.0006,0.0129,-0.0002;...
-0.0002,-0.0001,-0.0006,-0.0008,-0.0002,0.0150];
B0i=0.001*[-0.3908,-0.1297,0.7047,0.0591,0.2161,-0.6635];
B00=0.056;