%% load data

%% Step with wood pellets:

load('.\..\plant measurements and tests\xPC_Data_20071217_174202.mat') % good, model fitted after this. Sim. time: 8500
m_sf = [xpcdata6.time(5001:90000)-500 xpcdata6.Stok(5001:90000)];
m_a  = [xpcdata6.time(5001:90000)-500 xpcdata6.Blow(5001:90000)];
T_w  = [xpcdata6.time(5001:90000)-500 xpcdata6.Temp(5001:90000)];
T_iw = [xpcdata6.time(5001:90000)-500 xpcdata6.Retur(5001:90000)];
T_eg = [xpcdata6.time(5001:90000)-500 xpcdata6.Tr0xF8g(5001:90000)];
O2   = [xpcdata6.time(5001:90000)-500 xpcdata6.Lambda(5001:90000)];
fan  = [xpcdata6.time(5001:90000)-500 xpcdata6.Fan(5001:90000)];
u = [xpcdata6.time(5001:90000)-500 xpcdata6.Stok(5001:90000)*beta_motor xpcdata6.Blow(5001:90000)*beta_blower];
y_meas = [xpcdata6.time(5001:90000)-500 xpcdata6.Temp(5001:90000) xpcdata6.Lambda(5001:90000)];

% load('.\..\plant measurements and tests\xPC_Data_20071010_170227.mat') % not good, because lambda <0, data must be taken out to use this.
% m_sf = [xpcdata3.time xpcdata3.Stok];
% m_a  = [xpcdata3.time xpcdata3.Blow];
% T_w  = [xpcdata3.time xpcdata3.Temp];
% T_iw = [xpcdata3.time xpcdata3.Retur];
% T_eg = [xpcdata3.time xpcdata3.Tr0xF8g];
% O2   = [xpcdata3.time xpcdata3.Lambda];

% load('.\..\plant measurements and tests\xPC_Data_20071003_171127.mat') %Not good
% m_sf = [xpcdata3.time(101:25330)-10 xpcdata3.Stok(101:25330)];
% m_a  = [xpcdata3.time(101:25330)-10 xpcdata3.Blow(101:25330)];
% T_w  = [xpcdata3.time(101:25330)-10 xpcdata3.Temp(101:25330)];
% T_iw = [xpcdata3.time(101:25330)-10 xpcdata3.Retur(101:25330)];
% T_eg = [xpcdata3.time(101:25330)-10 xpcdata3.Tr0xF8g(101:25330)];
% O2   = [xpcdata3.time(101:25330)-10 xpcdata3.Lambda(101:25330)];

%% Steady state with wood pellets:

% load('.\..\plant measurements and tests\xPC_Data_20071217_121801.mat') %good. Sim. time: 4450
% m_sf = [xpcdata6.time(12001:56845)-xpcdata6.time(12001) xpcdata6.Stok(12001:56845)];
% m_a  = [xpcdata6.time(12001:56845)-xpcdata6.time(12001) xpcdata6.Blow(12001:56845)];
% T_w  = [xpcdata6.time(12001:56845)-xpcdata6.time(12001) xpcdata6.Temp(12001:56845)];
% T_iw = [xpcdata6.time(12001:56845)-xpcdata6.time(12001) xpcdata6.Retur(12001:56845)];
% T_eg = [xpcdata6.time(12001:56845)-xpcdata6.time(12001) xpcdata6.Tr0xF8g(12001:56845)];
% O2   = [xpcdata6.time(12001:56845)-xpcdata6.time(12001) xpcdata6.Lambda(12001:56845)];
% fan  = [xpcdata6.time(12001:56845)-xpcdata6.time(12001) xpcdata6.Fan(12001:56845)];

%% Step with corn:

% load('.\..\plant measurements and tests\xPC_Data_20080502_163308.mat') % sim time: 18946
% m_sf = [xpcdata6.time(102500:end)-xpcdata6.time(102500) xpcdata6.Stok(102500:end)];
% m_a  = [xpcdata6.time(102500:end)-xpcdata6.time(102500) xpcdata6.Blow(102500:end)];
% T_w  = [xpcdata6.time(102500:end)-xpcdata6.time(102500) xpcdata6.Temp(102500:end)];
% T_iw = [xpcdata6.time(102500:end)-xpcdata6.time(102500) xpcdata6.Retur(102500:end)];
% T_eg = [xpcdata6.time(102500:end)-xpcdata6.time(102500) xpcdata6.Tr0xF8g(102500:end)];
% O2   = [xpcdata6.time(102500:end)-xpcdata6.time(102500) xpcdata6.Lambda(102500:end)];
% fan  = [xpcdata6.time(102500:end)-xpcdata6.time(102500) xpcdata6.Fan(102500:end)];

%% Steady state with corn:

% load('.\..\plant measurements and tests\xPC_Data_20080205_153740.mat') % good.
% starttime = 90000;
% m_sf = [xpcdata6.time(starttime:end)-xpcdata6.time(starttime) xpcdata6.Stok(starttime:end)];
% m_a  = [xpcdata6.time(starttime:end)-xpcdata6.time(starttime) xpcdata6.Blow(starttime:end)];
% T_w  = [xpcdata6.time(starttime:end)-xpcdata6.time(starttime) xpcdata6.Temp(starttime:end)];
% T_iw = [xpcdata6.time(starttime:end)-xpcdata6.time(starttime) xpcdata6.Retur(starttime:end)];
% T_eg = [xpcdata6.time(starttime:end)-xpcdata6.time(starttime) xpcdata6.Tr0xF8g(starttime:end)];
% O2   = [xpcdata6.time(starttime:end)-xpcdata6.time(starttime) xpcdata6.Lambda(starttime:end)];
% fan  = [xpcdata6.time(starttime:end)-xpcdata6.time(starttime) xpcdata6.Fan(starttime:end)];

disp('Data loaded')