% Simulação de Equações de Balanço para Reator Descontínuo (Batelada)

clear all
global mixmax ks alfa yxs yps kp
% parâmetros cinéticos empregados nos modelos
mixmax=0.30;%0.23;
ks=0.2;
alfa=130;
kp=5;
% fatores de conversão empregados nos modelos
yxs=0.3;%0.04903475;
yps=0.50;%0.13513514;
% Condições iniciais das variáveis dependentes
x0=0.015;
s0=155.5;
p0=0;
% intervalo da variável independente, tempo
tspan=[0,21];
y0=[x0;s0;p0];
% Chamando a rotina Runge-Kutta 4/5 ordem para integrar y0 num tempo span
[t,xcn]=ode45('fbatelada',tspan,y0);
% Como a minha saída é a matriz xcn, preciso identificar as colunas
x=xcn(:,1);
s=xcn(:,2);
p=xcn(:,3);
xcn(length(xcn),1);
produto=xcn(length(xcn),3);
tempo = t(length(t));
Prod=produto/tempo;
tx=[t,xcn];
%save ('C:\Users\rpiccoli\Desktop\Aula 05 maio 2020\batelada_2.txt','tx','-ASCII');
% Caso tenhamos dados experimentais podemos colocar tudo num mesmo gráfico
% dados experimnentais representado por uma matriz a qualquer(t,s,p,x):
a=[0	155.5	0	0.015;
1.5	151.1	0.52	0.015;
3	150	2	0.017;
4.5	144.3	2.71	0.025;
6	137.2	1.93	0.034;
7.5	120.3	6.57	0.05;
9	113.6	10.35	0.058;
10.5	105.4	14.81	0.08;
12	89.4	20.9	0.106;
13.5	80	27.28	0.2;
15	65.8	33.91	0.238;
16.5	44.5	44.05	0.383;
18	25.6	52.66	0.6;
19.5	5.3	60.87	0.667;
21	0.1	64.26	0.777];
te=a(:,1);
xe=a(:,4);
se=a(:,2);
pe=a(:,3);

% montando os gráficos (o subscrito e = experimental)
subplot(221), plot(t,x,te,xe,'og')
title('X(g/L)');
subplot(222), plot(t,s,te,se,'ok')
title('S(g/L)');
subplot(223), plot(t,p,te,pe,'ob')
title('P(g/L)');