Estudio_Cinetico

Esta es una vista previa del archivo. Inicie sesión para ver el archivo original

clc;clear all;close all
%% Figura 1
% Datos experimentales
t=[0,10,20,30,40,50,60];
Fenol=[0,0.025,0.05,0.075,0.1,0.125,0.15];
pcl_fenol=[0,0.05,0.08,0.14,0.17,0.24,0.26];
p_cresol=[0,0.053,0.099,0.15,0.20,0.26,0.30];
% Cálculo de los parámetros de regresión lineal
p1=polyfit(t,Fenol,1); r2_1=corrcoef(t,Fenol); m_1=p1(1); b_1=p1(2);
p2=polyfit(t,pcl_fenol,1); r2_2=corrcoef(t,pcl_fenol); m_2=p2(1); b_2=p2(2);
p3=polyfit(t,p_cresol,1); r2_3=corrcoef(t,p_cresol); m_3=p3(1); b_3=p3(2);
% Datos a mostrar (Fenol)
fprintf('Fenol\nPendiente (m) =%5.3f\n',m_1);
fprintf('Intercepto (b) =%5.3f\n',b_1);
fprintf('Coef. de regresión =%5.3f\n\n',r2_1(2,1));
% Cálculo de los nuevos valores según la recta de ajuste
Fenol_C=polyval(p1,t);
pcl_fenol_C=polyval(p2,t);
p_cresol_C=polyval(p3,t);
% Tabla de resultados (Fenol)
fprintf('t(min) Ln_A0_A Ln_A0_A_C \n\n');
for i=1:length(t)
 fprintf('%2.0f',t(i));
 fprintf('%12.3f',Fenol(i));
 fprintf('%12.3f\n',Fenol_C(i));
end
% Datos a mostrar (p-cl Fenol)
fprintf('\np-cl Fenol\nPendiente (m) =%5.3f\n',m_2);
fprintf('Intercepto (b) =%5.3f\n',b_2);
fprintf('Coef. de regresión =%5.3f\n\n',r2_2(2,1));
%Tabla de resultados (p-cl Fenol)
fprintf('t(min) Ln_A0_A Ln_A0_A_C \n\n');
for i=1:length(t)
 fprintf('%2.0f',t(i));
 fprintf('%12.3f',pcl_fenol(i));
 fprintf('%12.3f\n',pcl_fenol_C(i));
end
% Datos a mostrar (p_cresol)
fprintf('\np Cresol\nPendiente (m) =%5.3f\n',m_3);
fprintf('Intercepto (b) =%5.3f\n',b_3);
fprintf('Coef. de regresión =%5.3f\n\n',r2_3(2,1));
%Tabla de resultados (p_cresol)
fprintf('t(min) Ln_A0_A Ln_A0_A_C \n\n');
for i=1:length(t)
 fprintf('%2.0f',t(i));
 fprintf('%12.3f',p_cresol(i));
 fprintf('%12.3f\n',p_cresol_C(i));
end
% Gráfico de los datos experimentales y de la recta de ajuste
f1=(plot(t,Fenol,'d'));
title('Ln[(A_\infty-A_0)/(A_\infty-A_t)] vs Tiempo')
ylabel('Ln[(A_\infty-A_0)/(A_\infty-A_t)]')
xlabel('Tiempo (min)')
set(figure(1),'Color','w')
hold on 
plot(t,Fenol_C)
f2=(plot(t,pcl_fenol,'^')); plot(t,pcl_fenol_C);
f3=(plot(t,p_cresol,'s')); plot(t,p_cresol_C);
legend([f1 f2 f3],{'Fenol y=0.002X R^2=1','p-cl Fenol y=0.004X R^2=0.995','p-Cresol y=0.005X R^2=0.999'})
hold off
%% Figura 2
% Datos experimentales
t=[0,10,20,30,40,50,60];
Fenol=[0,0.02,0.04,0.062,0.08,0.095,0.125];
pcl_fenol=[0,0.032,0.064,0.083,0.12,0.14,0.19];
p_cresol=[0,0.14,0.24,0.35,0.47,0.60,0.75];
% Cálculo de los parámetros de regresión lineal
p1=polyfit(t,Fenol,1); r2_1=corrcoef(t,Fenol); m_1=p1(1); b_1=p1(2);
p2=polyfit(t,pcl_fenol,1); r2_2=corrcoef(t,pcl_fenol); m_2=p2(1); b_2=p2(2);
p3=polyfit(t,p_cresol,1); r2_3=corrcoef(t,p_cresol); m_3=p3(1); b_3=p3(2);
% Datos a mostrar (Fenol)
fprintf('Fenol\nPendiente (m) =%5.3f\n',m_1);
fprintf('Intercepto (b) =%5.3f\n',b_1);
fprintf('Coef. de regresión =%5.3f\n\n',r2_1(2,1));
% Cálculo de los nuevos valores según la recta de ajuste
Fenol_C=polyval(p1,t);
pcl_fenol_C=polyval(p2,t);
p_cresol_C=polyval(p3,t);
% Tabla de resultados (Fenol)
fprintf('t(min) Ln_A0_A Ln_A0_A_C \n\n');
for i=1:length(t)
 fprintf('%2.0f',t(i));
 fprintf('%12.3f',Fenol(i));
 fprintf('%12.3f\n',Fenol_C(i));
end
% Datos a mostrar (p-cl Fenol)
fprintf('\np-cl Fenol\nPendiente (m) =%5.3f\n',m_2);
fprintf('Intercepto (b) =%5.3f\n',b_2);
fprintf('Coef. de regresión =%5.3f\n\n',r2_2(2,1));
%Tabla de resultados (p-cl Fenol)
fprintf('t(min) Ln_A0_A Ln_A0_A_C \n\n');
for i=1:length(t)
 fprintf('%2.0f',t(i));
 fprintf('%12.3f',pcl_fenol(i));
 fprintf('%12.3f\n',pcl_fenol_C(i));
end
% Datos a mostrar (p_cresol)
fprintf('\np Cresol\nPendiente (m) =%5.3f\n',m_3);
fprintf('Intercepto (b) =%5.3f\n',b_3);
fprintf('Coef. de regresión =%5.3f\n\n',r2_3(2,1));
%Tabla de resultados (p_cresol)
fprintf('t(min) Ln_A0_A Ln_A0_A_C \n\n');
for i=1:length(t)
 fprintf('%2.0f',t(i));
 fprintf('%12.3f',p_cresol(i));
 fprintf('%12.3f\n',p_cresol_C(i));
end
% Gráfico de los datos experimentales y de la recta de ajuste
figure(2)
f1=(plot(t,Fenol,'d'));
title('Ln[(A_\infty-A_0)/(A_\infty-A_t)] vs Tiempo')
ylabel('Ln[(A_\infty-A_0)/(A_\infty-A_t)]')
xlabel('Tiempo (min)')
set(figure(2),'Color','w')
hold on
plot(t,Fenol_C)
f2=(plot(t,pcl_fenol,'^')); plot(t,pcl_fenol_C);
f3=(plot(t,p_cresol,'s')); plot(t,p_cresol_C);
legend([f1 f2 f3],{'Fenol y=0.002X R^2=0.998','p-cl Fenol y=0.003X R^2=0.995','p-Cresol y=0.012X R^2=0.998'})
hold off
%% Figura 3
% Datos experimentales
Sub=[0.005,0.01,0.015,0.02,0.025,0.05];
Fenol=[0.33,0.53,0.9,1.2,1.46,2.8];
p_Cresol=[0.55,1.11,1.6,2,2.8,5.6];
pBr_Fenol=[1.03,1.7,2.57,3.7,4.45,8.86];
% Cálculo de los parámetros de regresión lineal
p1=polyfit(Sub,Fenol,1); r2_1=corrcoef(Sub,Fenol); m_1=p1(1); b_1=p1(2);
p2=polyfit(Sub,p_Cresol,1); r2_2=corrcoef(Sub,p_Cresol); m_2=p2(1); b_2=p2(2);
p3=polyfit(Sub,pBr_Fenol,1); r2_3=corrcoef(Sub,pBr_Fenol); m_3=p3(1); b_3=p3(2);
% Datos a mostrar (Fenol)
fprintf('Fenol\nPendiente (m) =%5.3f\n',m_1);
fprintf('Intercepto (b) =%5.3f\n',b_1);
fprintf('Coef. de regresión =%5.3f\n\n',r2_1(2,1));
% Cálculo de los nuevos valores según la recta de ajuste
Fenol_C=polyval(p1,Sub);
pcl_fenol_C=polyval(p2,Sub);
p_cresol_C=polyval(p3,Sub);
% Tabla de resultados (Fenol)
fprintf('[Sustrato] k k_C \n\n');
for i=1:length(Sub)
 fprintf('%2.3f',Sub(i));
 fprintf('%12.3f',Fenol(i));
 fprintf('%12.3f\n',Fenol_C(i));
end
% Datos a mostrar (p-Cresol)
fprintf('\np-Cresol\nPendiente (m) =%5.3f\n',m_2);
fprintf('Intercepto (b) =%5.3f\n',b_2);
fprintf('Coef. de regresión =%5.3f\n\n',r2_2(2,1));
%Tabla de resultados (p-Cresol)
fprintf('[Sustrato] k k_C \n\n');
for i=1:length(Sub)
 fprintf('%2.3f',Sub(i));
 fprintf('%12.3f',p_Cresol(i));
 fprintf('%12.3f\n',pcl_fenol_C(i));
end
% Datos a mostrar (p-Br Fenol)
fprintf('\np-Br Fenol\nPendiente (m) =%5.3f\n',m_3);
fprintf('Intercepto (b) =%5.3f\n',b_3);
fprintf('Coef. de regresión =%5.3f\n\n',r2_3(2,1));
%Tabla de resultados (p-Br Fenol)
fprintf('[Sustrato] k k_C \n\n');
for i=1:length(Sub)
 fprintf('%2.3f',Sub(i));
 fprintf('%12.3f',pBr_Fenol(i));
 fprintf('%12.3f\n',p_cresol_C(i));
end
% Gráfico de los datos experimentales y de la recta de ajuste
figure(3)
f1=(plot(Sub,Fenol,'d'));
title('k vs [Sustrato]')
xlabel('[Sustrato](mol/dm^3)')
ylabel('10^2k/min')
set(figure(3),'Color','w')
hold on
plot(Sub,Fenol_C)
f2=(plot(Sub,p_Cresol,'s')); plot(Sub,pcl_fenol_C);
f3=(plot(Sub,pBr_Fenol,'^')); plot(Sub,p_cresol_C);
legend([f1 f2 f3],{'Fenol y=55.502X+0.047 R^2=0.999','p-Cresol y=112.734X-0.072 R^2=0.999','p-Br Fenol y=176.387X+0.044 R^2=0.999'})
hold off
%% Figura 4
% Datos experimentales
Sub=[0.005,0.01,0.015,0.02,0.025,0.05];
Fenol=[0.2,0.4,0.56,0.72,0.8,1.6];
p_Cresol=[1.4,2.3,3.5,4.7,5.5,12];
% Cálculo de los parámetros de regresión lineal
p1=polyfit(Sub,Fenol,1); r2_1=corrcoef(Sub,Fenol); m_1=p1(1); b_1=p1(2);
p2=polyfit(Sub,p_Cresol,1); r2_2=corrcoef(Sub,p_Cresol); m_2=p2(1); b_2=p2(2);
% Datos a mostrar (Fenol)
fprintf('Fenol\nPendiente (m) =%5.3f\n',m_1);
fprintf('Intercepto (b) =%5.3f\n',b_1);
fprintf('Coef. de regresión =%5.3f\n\n',r2_1(2,1));
% Cálculo de los nuevos valores según la recta de ajuste
Fenol_C=polyval(p1,Sub);
pcl_fenol_C=polyval(p2,Sub);
% Tabla de resultados (Fenol)
fprintf('[Sustrato] k k_C \n\n');
for i=1:length(Sub)
 fprintf('%2.3f',Sub(i));
 fprintf('%12.3f',Fenol(i));
 fprintf('%12.3f\n',Fenol_C(i));
end
% Datos a mostrar (p-Cresol)
fprintf('\np-Cresol\nPendiente (m) =%5.3f\n',m_2);
fprintf('Intercepto (b) =%5.3f\n',b_2);
fprintf('Coef. de regresión =%5.3f\n\n',r2_2(2,1));
%Tabla de resultados (p-Cresol)
fprintf('[Sustrato] k k_C \n\n');
for i=1:length(Sub)
 fprintf('%2.3f',Sub(i));
fprintf('%12.3f',p_Cresol(i));
 fprintf('%12.3f\n',pcl_fenol_C(i));
end
% Gráfico de los datos experimentales y de la recta de ajuste
figure(4)
f1=(plot(Sub,Fenol,'d'));
title('k vs [Sustrato]')
xlabel('[Sustrato](mol/dm^3)')
ylabel('10^2k/min')
set(figure(4),'Color','w')
hold on
plot(Sub,Fenol_C)
f2=(plot(Sub,p_Cresol,'s')); plot(Sub,pcl_fenol_C);
legend([f1 f2],{'Fenol y=30.400X+0.080 R^2=0.998','p-Cresol y=237.246X-0.043 R^2=0.998'})
hold off
%% Figura 5
% Datos experimentales
t=[0,60,122,185,260];
Fenol=[6.6,7.7,9,9.9,11.2];
% Cálculo de los parámetros de regresión lineal
p1=polyfit(t,Fenol,1); r2_1=corrcoef(t,Fenol); m_1=p1(1); b_1=p1(2);
% Datos a mostrar (Fenol)
fprintf('Fenol\nPendiente (m) =%5.3f\n',m_1);
fprintf('Intercepto (b) =%5.3f\n',b_1);
fprintf('Coef. de regresión =%5.3f\n\n',r2_1(2,1));
% Cálculo de los nuevos valores según la recta de ajuste
Fenol_C=polyval(p1,t);
% Tabla de resultados (Fenol)
fprintf('Tiempo Fenol Fenol_C \n\n');
for i=1:length(t)
 fprintf('%2.0f',t(i));
 fprintf('%12.3f',Fenol(i));
 fprintf('%12.3f\n',Fenol_C(i));
end
% Gráfico de los datos experimentales y de la recta de ajuste
figure(5)
f1=(plot(t,Fenol,'d'));
title('1/(A_\infty-A_t) vs Time')
xlabel('Time(min)')
ylabel('1/(A_\infty-A_t)')
set(figure(5),'Color','w')
hold on
plot(t,Fenol_C)
legend([f1],{'Fenol y=0.018X+6.667 R^2=0.998'})
hold off
%% Figura 6
% Datos experimentales
t=[0,5,10,15,22,26,32,37,40,46,51,56,60];
p_Cresol=[0.78,0.85,0.92,0.98,1.05,1.12,1.2,1.23,1.3,1.32,1.4,1.5,1.55];
% Cálculo de los parámetros de regresión lineal
p1=polyfit(t,p_Cresol,1); r2_1=corrcoef(t,p_Cresol); m_1=p1(1); b_1=p1(2);
% Datos a mostrar (p-Cresol)
fprintf('\np-Cresol\nPendiente (m) =%5.3f\n',m_1);
fprintf('Intercepto (b) =%5.3f\n',b_1);
fprintf('Coef. de regresión =%5.3f\n\n',r2_1(2,1));
% Cálculo de los nuevos valores según la recta de ajuste
p_Cresol_C=polyval(p1,t);
%Tabla de resultados (p-Cresol)
fprintf('Tiempo p-Cresol p-Cresol_C \n\n');
for i=1:length(t)
 fprintf('%2.0f',t(i));
 fprintf('%12.3f',p_Cresol(i));
 fprintf('%12.3f\n',p_Cresol_C(i));
end
% Gráfico de los datos experimentales y de la recta de ajuste
figure(6)
f1=(plot(t,p_Cresol,'d'));
title('1/(A_\infty-A_t) vs Time')
xlabel('Time(min)')
ylabel('1/(A_\infty-A_t)')
set(figure(6),'Color','w')
hold on
plot(t,p_Cresol_C)
legend([f1],{'p-Cresol y=0.012X+0.787 R^2=0.997'})
hold off
%% Figura 7
% Datos experimentales
t=[0,5.11,9.76,14.42,19.76,24.42,30,34.88,40,45.3,50,54.65,60];
p_Cresol=[0.87,0.94,0.98,1.05,1.11,1.15,1.2,1.35,1.45,1.55,1.65,1.75,1.85];
% Cálculo de los parámetros de regresión lineal
p1=polyfit(t,p_Cresol,1); r2_1=corrcoef(t,p_Cresol); m_1=p1(1); b_1=p1(2);
% Datos a mostrar (p-Cresol)
fprintf('\np-Cresol\nPendiente (m) =%5.3f\n',m_1);
fprintf('Intercepto (b) =%5.3f\n',b_1);
fprintf('Coef. de regresión =%5.3f\n\n',r2_1(2,1));
% Cálculo de los nuevos valores según la recta de ajuste
p_Cresol_C=polyval(p1,t);
%Tabla de resultados (p-Cresol)
fprintf('Tiempo p-Cresol p-Cresol_C \n\n');
for i=1:length(t)
 fprintf('%2.0f',t(i));
 fprintf('%12.3f',p_Cresol(i));
 fprintf('%12.3f\n',p_Cresol_C(i));
end
% Gráfico de los datos experimentales y de la recta de ajuste
figure(7)
f1=(plot(t,p_Cresol,'d'));
title('1/(A_\infty-A_t) vs Time')
xlabel('Time(min)')
ylabel('1/(A_\infty-A_t)')
set(figure(7),'Color','w')
hold on
plot(t,p_Cresol_C)
legend([f1],{'p-Cresol y=0.016X+0.809 R^2=0.989'})
hold off
%% Figura 8
% Datos experimentales
t=[0,60,122,182,241];
P_cFENOL=[4.3,4.53,4.75,5,5.3];
% Cálculo de los parámetros de regresión lineal
p1=polyfit(t,P_cFENOL,1); r2_1=corrcoef(t,P_cFENOL); m_1=p1(1); b_1=p1(2);
% Datos a mostrar (p–cloro Fenol)
fprintf('\np–cloro Fenol\nPendiente (m) =%5.3f\n',m_1);
fprintf('Intercepto (b) =%5.3f\n',b_1);
fprintf('Coef. de regresión =%5.3f\n\n',r2_1(2,1));
% Cálculo de los nuevos valores según la recta de ajuste
P_cFENOL_C=polyval(p1,t);
%Tabla de resultados (p–cloro Fenol)
fprintf('Tiempo p–cloro Fenol p–cloro Fenol_C \n\n');
for i=1:length(t)
 fprintf('%2.0f',t(i));
 fprintf('%12.3f',P_cFENOL(i));
 fprintf('%12.3f\n',P_cFENOL_C(i));
end
% Gráfico de los datos experimentales y de la recta de ajuste
figure(8)
f1=(plot(t,P_cFENOL,'d'));
title('1/(A_\infty-A_t) vs Time')
xlabel('Time(min)')
ylabel('1/(A_\infty-A_t)')
set(figure(8),'Color','w')
hold on
plot(t,P_cFENOL_C)
legend([f1],{'p–cloro Fenol y=0.004X+4.281 R^2=0.998'})
hold off
%% Figura 9
% Datos experimentales
T=[303,308,313,318,329];
Fenol=[75650,76200,76800,77800,78600,];
p_clFenol=[75960,77200,78600,79800,81700];
% Cálculo de los parámetros de regresión lineal
p1=polyfit(T,Fenol,1); r2_1=corrcoef(T,Fenol); m_1=p1(1); b_1=p1(2);
p2=polyfit(T,p_clFenol,1); r2_2=corrcoef(T,p_clFenol); m_2=p2(1); b_2=p2(2);
% Datos a mostrar (Fenol)
fprintf('Fenol\nPendiente (m) =%5.3f\n',m_1);
fprintf('Intercepto (b) =%5.3f\n',b_1);
fprintf('Coef. de regresión =%5.3f\n\n',r2_1(2,1));
% Cálculo de los nuevos valores según la recta de ajuste
Fenol_C=polyval(p1,T);
p_clFenol_C=polyval(p2,T);
% Tabla de resultados (Fenol)
fprintf('k G G_C \n\n');
for i=1:length(T)
 fprintf('%2.0f',T(i));
 fprintf('%12.0f',Fenol(i));
 fprintf('%12.0f\n',Fenol_C(i));
end
% Datos a mostrar (p-cl Fenol)
fprintf('\np-cl Fenol\nPendiente (m) =%5.3f\n',m_2);
fprintf('Intercepto (b) =%5.3f\n',b_2);
fprintf('Coef. de regresión =%5.3f\n\n',r2_2(2,1));
%Tabla de resultados (p-cl Fenol)
fprintf('k G G_C \n\n');
for i=1:length(T)
 fprintf('%2.0f',T(i));
 fprintf('%12.0f',p_clFenol(i));
 fprintf('%12.0f\n',p_clFenol_C(i));
end
% Gráfico de los datos experimentales y de la recta de ajuste
figure(9)
f1=(plot(T,Fenol,'d'));
title('Gibbs-Helmholtz')
xlabel('Temperatura (K)')
ylabel('\DeltaG(J/mol)')
set(figure(9),'Color','w')
hold on
plot(T,Fenol_C)
f2=(plot(T,p_clFenol,'s')); plot(T,p_clFenol_C);
legend([f1 f2],{'Fenol y=117.954X+39948.897 R^2=0.986','p-cl Fenol y=222.387X+8777.954 R^2=0.994'})
hold off
%% Figura 10
% Datos experimentales
T=[303,308,313,318,329];
Fenol=[77000,77400,78000,78600,79300,];
p_clFenol=[78600,78800,79000,79200,79400];
% Cálculo de los parámetros de regresión lineal
p1=polyfit(T,Fenol,1); r2_1=corrcoef(T,Fenol); m_1=p1(1); b_1=p1(2);
p2=polyfit(T,p_clFenol,1); r2_2=corrcoef(T,p_clFenol); m_2=p2(1); b_2=p2(2);
% Datos a mostrar (Fenol)
fprintf('Fenol\nPendiente (m) =%5.3f\n',m_1);
fprintf('Intercepto (b) =%5.3f\n',b_1);
fprintf('Coef. de regresión =%5.3f\n\n',r2_1(2,1));
% Cálculo de los nuevos valores según la recta de ajuste
Fenol_C=polyval(p1,T);
p_clFenol_C=polyval(p2,T);
% Tabla de resultados (Fenol)
fprintf('k G G_C \n\n');
for i=1:length(T)
 fprintf('%2.3f',T(i));
 fprintf('%12.3f',Fenol(i));
 fprintf('%12.3f\n',Fenol_C(i));
end
% Datos a mostrar (p-Cresol)
fprintf('\np-Cresol\nPendiente (m) =%5.3f\n',m_2);
fprintf('Intercepto (b) =%5.3f\n',b_2);
fprintf('Coef. de regresión =%5.3f\n\n',r2_2(2,1));
%Tabla de resultados (p-Cresol)
fprintf('k G G_C \n\n');
for i=1:length(T)
 fprintf('%2.3f',T(i));
 fprintf('%12.3f',p_clFenol(i));
 fprintf('%12.3f\n',p_clFenol_C(i));
end
% Gráfico de los datos experimentales y de la recta de ajuste
figure(10)
f1=(plot(T,Fenol,'d'));
title('Gibbs-Helmholtz')
xlabel('Temperatura (K)')
ylabel('\DeltaG(J/mol)')
set(figure(10),'Color','w')
hold on
plot(T,Fenol_C)
f2=(plot(T,p_clFenol,'s')); plot(T,p_clFenol_C);
legend([f1 f2],{'Fenol y=91.374X+49350.251 R^2=0.991','p-cl Fenol y=31.093X+69230.491 R^2=0.982'})
hold off