Ejercicios Matematica financiera
Enviado por Danna Medina • 24 de Septiembre de 2023 • Práctica o problema • 1.573 Palabras (7 Páginas) • 101 Visitas
Soluciones de EDO
[pic 1]
Taller N°3 |
Angie Hernández, Andrés Martínez, Danna Medina, Camila Molina. |
Facultad de Ingeniería, Ingenieria Quimica, Universidad del Atlántico, Barranquilla, Colombia. |
- Una ley de velocidad de enésimo orden a menudo se usa para modelar reacciones químicas que dependen únicamente de la concentración de un solo reactivo:
[pic 2]
Donde:
c = Concentración (moles)
t =Tiempo(min)
n=Orden de la reacción
k= Constante de velocidad (min – 1mol1 – n).
El método diferencial se puede usar para evaluar los parámetros k y n. Esto implica aplicar una transformación logarítmica a la ley de velocidad para producir:
[pic 3]
Por lo tanto, si la ley de enésimo orden tiene una gráfica de log (– dc/dt) versus log(c) debería producir una línea recta con una pendiente n y una intersección de log(k). Use el método de diferencias finitas divididas y regresión lineal para determinar k y n dados los siguientes datos para la conversión de cianato de amonio a urea.
[pic 4]
- Solucion.
Por medio del método de diferencias divididas se ajusta una serie de datos a un polinomio, la cual está representada por su ecuación general:
[pic 5]
Para lo cual es necesario realizar una tabla en la cual se determinará una aproximación de la razón de cambio de los datos, haciendo uso de la siguiente ecuación:
Donde f representa la razón de cambio de los datos. Los cálculos se muestran a continuación:
[pic 6]
Tabla 1. Datos y aproximaciones de las derivadas.
x | C | f’ | f’’ | f’’’ | f4’ |
0 | 0.75 | -- | -- | -- | -- |
5 | 0.594 | -0.0312 | -- | -- | -- |
15 | 0.42 | -0.0174 | 0.00092 | -- | -- |
30 | 0.291 | -0.0086 | 0.00035 | [pic 7] | -- |
45 | 0.223 | -0.0045 | 0.00013 | [pic 8] | [pic 9] |
A partir de esta tabla se determina que:
𝑎0 = 0,75
𝑎1 = −0.0312
𝑎2 = 0,00092
𝑎3 = −1,8933x10−5
𝑎4 = 3,0049𝑥10−7
Por tanto, al reemplazar los datos obtenidos en la ecuación 1, agrupar términos semejantes y simplificar obtenemos la siguiente ecuación:
[pic 10]
[pic 11]
Obteniendo los resultados:
[pic 12]
La línea roja en la gráfica representa a la curva del polinomio de Newton, mientras que las cruces azules representan los datos experimentales observándose que la curva corta cada uno de estos puntos, además de que los datos arrojados por el software son casi idénticos.
[pic 13]
[pic 14]
Ahora se procede a determinar el orden (n) de la reacción y la constante de velocidad cinética (k), por medio del siguiente algoritmo:
Clear All Clc
disp('Este programa calcula la constante de velocidad y el orden de la reacción')
disp('con los siguientes datos suministrados de tiempo y concentración')
disp(' ')
t=[0 5 15 30 45]';
C=[0.75 0.594 0.42 0.291 0.223]';
M=5;N=M-1;
disp(' Tiempo Concentración ')
disp([t C]);
for i=1:N
c(i,1)=(C(i+1)-C(i))/(t(i+1)-t(i));
end
c; % Calcula Dc/dt
disp( ' Dc/dt Cpromedio Cprom^1.5 ')
for i=1:M-1
Cprom(i)=(C(i+1)+C(i))/2;
CCprom(i)=Cprom(i)^1.5;
end
Cprom=Cprom';CCprom=CCprom';
disp([c Cprom CCprom])
cc=log(abs(c));CC=log(Cprom);Cn=log(CCprom);
A=[ones(size(CC)) CC];
coef=A\cc;
k =exp(coef(1));
n =coef(2);
ccc=(log(0.75):-0.01:log(0.223))';
Dcdt=[ones(size(ccc)) ccc]*coef;
loglog(CC,cc,'o',ccc,Dcdt,'-m');
title('Grafica de concentración vs velocidad de reacción');
xlabel('Concentración','Fontsize',10);
ylabel('Velocidad de reacción','Fontsize',10);
text(-10^0.04,-10^0.7,'-dC/dt=K*C^n','Fontsize',10);
...