%Este archivo toma los datos de cada workspace resultado de las simulaciones del
%proceso Eastman, para graficar y desplegarlos en orden similar a los
%resultados presentados por Luyben para su sistema de Control.
%Autor: Cinda Luz Sandoval Torres
%CENIDET. Ingenieria Mecatronica.
%2007-2009
clc;
figure1 = figure('Name','Pruebas del sistema de control en el Proceso Eastman (1 de 2)');
d = XMEAS7;
s=size(d);
FM=1/abs(d(1,1)-d(2,1));
vector=[0:(1/FM):size(d)*(1/FM)];
subplot(4,4,1);
plot(vector(1:7201),[d(1:s(1),2)]);
xlim([0 20]);
% El comando plot grafica los vectores (2)
% contra el vector de tiempo (1)
title('Presión del reactor');
xlabel('Tiempo(h)');
ylabel('Presión (kPa)');
grid

d = ScopeData4;
s=size(d);
FM=1/abs(d(1,1)-d(2,1));
vector=[0:(1/FM):size(d)*(1/FM)];
subplot(4,4,2);
plot(vector(1:7201),[d(1:s(1),2)]);
xlim([0 20]);
% El comando plot grafica los vectores (2)
% contra el vector de tiempo (1)
title('Flujo de alim. A/C');
xlabel('Tiempo (h)');
ylabel('% apertura XMV4');
grid

d = XMEAS27;
s=size(d);
FM=1/abs(d(1,1)-d(2,1));
vector=[0:(1/FM):size(d)*(1/FM)];
subplot(4,4,3);
plot(vector(1:7201),[d(1:s(1),2)]);
xlim([0 20]);
% El comando plot grafica los vectores (2)
% contra el vector de tiempo (1)
title('A (Corriente 9 o de purga)');
xlabel('Tiempo (h)');
ylabel('yA (Mol %))');
grid

d = ScopeData3;
s=size(d);
FM=1/abs(d(1,1)-d(2,1));
vector=[0:(1/FM):size(d)*(1/FM)];
subplot(4,4,4);
plot(vector(1:7201),[d(1:s(1),2)]);
xlim([0 20]);
% El comando plot grafica los vectores (2)
% contra el vector de tiempo (1)
title('Flujo de alim. A');
xlabel('Tiempo (h)');
ylabel('% apertura XMV3');
grid

d = XMEAS9;
s=size(d);
FM=1/abs(d(1,1)-d(2,1));
vector=[0:(1/FM):size(d)*(1/FM)];
subplot(4,4,5);
plot(vector(1:7201),[d(1:s(1),2)]);
xlim([0 20]);
% El comando plot grafica los vectores (2)
% contra el vector de tiempo (1)
title('Temp. del reactor');
xlabel('Tiempo (h)');
ylabel('TR (°C)');
grid

d = ScopeData9;
s=size(d);
FM=1/abs(d(1,1)-d(2,1));
vector=[0:(1/FM):size(d)*(1/FM)];
subplot(4,4,6);
plot(vector(1:7201),[d(1:s(1),2)]);
xlim([0 20]);
% El comando plot grafica los vectores (2)
% contra el vector de tiempo (1)
title('Flujo agua enfriamiento reactor');
xlabel('Tiempo (h)');
ylabel('% apertura XMV10');
grid

d = XMEAS28;
s=size(d);
FM=1/abs(d(1,1)-d(2,1));
vector=[0:(1/FM):size(d)*(1/FM)];
subplot(4,4,7);
plot(vector(1:7201),[d(1:s(1),2)]);
xlim([0 20]);
% El comando plot grafica los vectores (2)
% contra el vector de tiempo (1)
title('B (Corriente 9 o de purga');
xlabel('Tiempo (h)');
ylabel('yB (Mol %)');
grid

d = XMEAS10;
s=size(d);
FM=1/abs(d(1,1)-d(2,1));
vector=[0:(1/FM):size(d)*(1/FM)];
subplot(4,4,8);
plot(vector(1:7201),[d(1:s(1),2)]);
xlim([0 20]);
ylim([0 1]);
% El comando plot grafica los vectores (2)
% contra el vector de tiempo (1)
title('Vel. de purga');
xlabel('Tiempo (h)');
ylabel('Purga(kscmh)');
grid

d = XMEAS11;
s=size(d);
FM=1/abs(d(1,1)-d(2,1));
vector=[0:(1/FM):size(d)*(1/FM)];
subplot(4,4,9);
plot(vector(1:7201),[d(1:s(1),2)]);
xlim([0 20]);
% El comando plot grafica los vectores (2)
% contra el vector de tiempo (1)
title('Temp. del separador');
xlabel('Tiempo (h)');
ylabel('Tsep (°C)');
grid

d = XMEAS1;
s=size(d);
FM=1/abs(d(1,1)-d(2,1));
vector=[0:(1/FM):size(d)*(1/FM)];
subplot(4,4,10);
plot(vector(1:7201),[d(1:s(1),2)]);
xlim([0 20]);
% El comando plot grafica los vectores (2)
% contra el vector de tiempo (1)
title('Flujo inferior columna (Producto)');
xlabel('Tiempo (h)');
ylabel('B (m3/h)');
grid

d = XMEAS29;
s=size(d);
FM=1/abs(d(1,1)-d(2,1));
vector=[0:(1/FM):size(d)*(1/FM)];
subplot(4,4,11);
plot(vector(1:7201),[d(1:s(1),2)]);
xlim([0 20]);
% El comando plot grafica los vectores (2)
% contra el vector de tiempo (1)
title('C (Corriente 9 o de purga)');
xlabel('Tiempo (h)');
ylabel('yC (Mol %)');
grid

d = XMEAS22;
s=size(d);
FM=1/abs(d(1,1)-d(2,1));
vector=[0:(1/FM):size(d)*(1/FM)];
subplot(4,4,12);
plot(vector(1:7201),[d(1:s(1),2)]);
xlim([0 20]);
% El comando plot grafica los vectores (2)
% contra el vector de tiempo (1)
title('Flujo de reciclado');
xlabel('Tiempo (h)');
ylabel('Reciclado(Kscmh)');
grid

d = XMEAS5;
s=size(d);
FM=1/abs(d(1,1)-d(2,1));
vector=[0:(1/FM):size(d)*(1/FM)];
subplot(4,4,13);
plot(vector(1:7201),[d(1:s(1),2)]);
xlim([0 20]);
% El comando plot grafica los vectores (2)
% contra el vector de tiempo (1)
title('Temp. de columna');
xlabel('Tiempo (h)');
ylabel('Tstrip (°C)');
grid

d = ScopeData8;
s=size(d);
FM=1/abs(d(1,1)-d(2,1));
vector=[0:(1/FM):size(d)*(1/FM)];
subplot(4,4,14);
plot(vector(1:7201),[d(1:s(1),2)]);
xlim([0 20]);
% El comando plot grafica los vectores (2)
% contra el vector de tiempo (1)
title('Flujo de vapor a la columna');
xlabel('Tiempo (h)');
ylabel('% apertura XMV9');
grid

d = XMEAS37;
s=size(d);
FM=1/abs(d(1,1)-d(2,1));
vector=[0:(1/FM):size(d)*(1/FM)];
subplot(4,4,15);
plot(vector(1:7201),[d(1:s(1),2)]);
xlim([0 20]);
% El comando plot grafica los vectores (2)
% contra el vector de tiempo (1)
title('E (Corriente 11 o de producto)');
xlabel('Tiempo (h)');
ylabel('xBE (Mol %)');
grid

d = ScopeData8;
s=size(d);
FM=1/abs(d(1,1)-d(2,1));
vector=[0:(1/FM):size(d)*(1/FM)];
subplot(4,4,16);
plot(vector(1:7201),[d(1:s(1),2)]);
xlim([0 20]);
% El comando plot grafica los vectores (2)
% contra el vector de tiempo (1)
title('Flujo de vapor a la columna');
xlabel('Tiempo (h)');
ylabel('% apertura XMV9');
grid

figure2 = figure('Name','Pruebas del sistema de control en el Proceso Eastman (2 de 2)');

d = XMEAS8;
s=size(d);
FM=1/abs(d(1,1)-d(2,1));
vector=[0:(1/FM):size(d)*(1/FM)];
subplot(3,4,1);
plot(vector(1:7201),[d(1:s(1),2)]);
xlim([0 20]);
% El comando plot grafica los vectores (2)
% contra el vector de tiempo (1)
title('Nivel en reactor');
xlabel('Tiempo (h)');
ylabel('HR (%)');
grid

d = ScopeData;
s=size(d);
FM=1/abs(d(1,1)-d(2,1));
vector=[0:(1/FM):size(d)*(1/FM)];
subplot(3,4,2);
plot(vector(1:7201),[d(1:s(1),2)]);
xlim([0 20]);
% El comando plot grafica los vectores (2)
% contra el vector de tiempo (1)
title('Flujo de alim. D (Corriente 2)');
xlabel('Tiempo (h)');
ylabel('Fod (kg/h)');
grid

d = XMEAS12;
s=size(d);
FM=1/abs(d(1,1)-d(2,1));
vector=[0:(1/FM):size(d)*(1/FM)];
subplot(3,4,5);
plot(vector(1:7201),[d(1:s(1),2)]);
xlim([0 20]);
% El comando plot grafica los vectores (2)
% contra el vector de tiempo (1)
title('Nivel en separador');
xlabel('Tiempo (h)');
ylabel('Hsep (%)');
grid

d = ScopeData10;
s=size(d);
FM=1/abs(d(1,1)-d(2,1));
vector=[0:(1/FM):size(d)*(1/FM)];
subplot(3,4,6);
plot(vector(1:7201),[d(1:s(1),2)]);
xlim([0 20]);
% El comando plot grafica los vectores (2)
% contra el vector de tiempo (1)
title('Flujo agua de enfriam. del condensador');
xlabel('Tiempo (h)');
ylabel('% apertura XMV11');
grid

d = XMEAS39;
s=size(d);
FM=1/abs(d(1,1)-d(2,1));
vector=[0:(1/FM):size(d)*(1/FM)];
subplot(3,4,7);
plot(vector(1:7201),[d(1:s(1),2)]);
xlim([0 20]);
% El comando plot grafica los vectores (2)
% contra el vector de tiempo (1)
title('G (Corriente 11 o producto)');
xlabel('Tiempo (h)');
ylabel('xBG (Mol %)');
grid

d = ScopeData;
s=size(d);
FM=1/abs(d(1,1)-d(2,1));
vector=[0:(1/FM):size(d)*(1/FM)];
subplot(3,4,8);
plot(vector(1:7201),[d(1:s(1),2)]);
xlim([0 20]);
% El comando plot grafica los vectores (2)
% contra el vector de tiempo (1)
title('Flujo de alim. D (Corriente 2)');
xlabel('Tiempo (h)');
ylabel('% apertura XMV1');
grid

d = XMEAS3;
s=size(d);
FM=1/abs(d(1,1)-d(2,1));
vector=[0:(1/FM):size(d)*(1/FM)];
subplot(3,4,9);
plot(vector(1:7201),[d(1:s(1),2)]);
xlim([0 20]);
% El comando plot grafica los vectores (2)
% contra el vector de tiempo (1)
title('Nivel en columna');
xlabel('Tiempo (h)');
ylabel('Hstrip (%)');
grid

d = ScopeData6;
s=size(d);
FM=1/abs(d(1,1)-d(2,1));
vector=[0:(1/FM):size(d)*(1/FM)];
subplot(3,4,10);
plot(vector(1:7201),[d(1:s(1),2)]);
xlim([0 20]);
% El comando plot grafica los vectores (2)
% contra el vector de tiempo (1)
title('Flujo de líq. del separador (Corriente 10)');
xlabel('Tiempo (h)');
ylabel('% apertura XMV7');
grid

d = XMEAS40;
s=size(d);
FM=1/abs(d(1,1)-d(2,1));
vector=[0:(1/FM):size(d)*(1/FM)];
subplot(3,4,11);
plot(vector(1:7201),[d(1:s(1),2)]);
xlim([0 20]);
% El comando plot grafica los vectores (2)
% contra el vector de tiempo (1)
title('H (Corriente 11 o producto)');
xlabel('Tiempo (h)');
ylabel('xBH (Mol %)');
grid

d = ScopeData1;
s=size(d);
FM=1/abs(d(1,1)-d(2,1));
vector=[0:(1/FM):size(d)*(1/FM)];
subplot(3,4,12);
plot(vector(1:7201),[d(1:s(1),2)]);
xlim([0 20]);
% El comando plot grafica los vectores (2)
% contra el vector de tiempo (1)
title('Flujo de alim. E (Corriente 3)');
xlabel('Tiempo (h)');
ylabel('% apertura XMV2');
grid
