Taller 3

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 10

TALLER #3

Informática para ingeniería química

AJUSTE E INTERPOLACIÓN DE CURVAS

Presentado por:

Presentado a:

PROFESOR FRANCISCO MUÑOZ PABA M.SC

UNIVERSIDAD DEL ATLÁNTICO


BARRANQUILLA, ATLÁNTICO
19/07/19

1. Tres organismos portadores de enfermedades se descomponen en agua de mar según el


siguiente modelo:
−1.5 t −0.3t −0.05 t
y= A e +Be +C e
Estimar la concentración inicial de cada organismo (A, B y C) dadas las siguientes
medidas.
t 0.5 1 2 3 4 5 6 7 9
y 6 4.4 3.2 2.7 2 1.9 1.7 1.4 1.1
⮚ Entrada en el editor
1- Clear all clf clc
2- t=[0.5 1 2 3 4 5 6 7 9]’;
3- y=[6 4.4 3.2 2.7 2 1.9 1.7 1.4 1.1]’;
4- A=[ones(size(t)) t exp(-1.5*t) exp(-0.3*t) exp(-0.05*t)];
5- coef=A\y;
6- a=coef(1)
7- a1=coef(2)
8- a2=coef(3)
9- T=(0.5:0.1:9)’;
10- Y=[ones(size(T)) T exp(-1.5*T) exp(-0.3*T) exp(-0.05*T)]*coef;
11- plot(t,y,’o’,T,Y,’-’)

⮚ Resultados

a = -20.8824
a1 = 0.6767
a2 = 5.1087

⮚ Vista del editor


⮚ Gráfica Matlab

⮚ Vista general desde Matlab


2. Un investigador ha informado de los datos tabulados a continuación. Se sabe que estos
datos pueden ser modelados por la siguiente ecuación:
( y−b)
a
x=e
Donde a y b son parámetros. Desarrolle un programa en Matlab, usando regresión no lineal
para determinar a y b. En función de su ajuste halle y para x=2.6.
x 1 2 3 4 5
y 0.5 2 2.9 3.5 4

⮚ Linealización de la ecuación
( y−b)
a
x=e
( y−b )
ln ln x=
a
a ln ln x= y −b
y=b +a ln ln x

⮚ Entrada en el editor
1- Clear all clf clc
2- x=[1 2 3 4 5]’;
3- y=[0.5 2 2.9 3.5 4]’;
4- X=log(x);
5- A=[ones(size(X)) X];
6- coef=A\y;
7- a=coef(2)
8- b=coef(1)
9- XX=(log(1.8):0.1:log(10))’;
10- Y=[ones(size(XX)) XX]*coef;
11- plot(X,y,’o’,XX,Y,’-’)

⮚ Resultados

a = 2.1729
b = 0.4994

⮚ El valor de y para cuando x = 2.6:

y=b +a ln ln x

y=0.4994 +(2.1729)2.6 ¿

y=2.5756

⮚ Vista del editor


⮚ Gráfica Matlab
⮚ Vista general desde Matlab

3. Se reunieron los siguientes datos para determinar la relación entre presión y temperatura
de un volumen fijo de 1 kg de nitrógeno. El volumen es de 10 m3.
t, °C -20 0 20 40 50 70 100 120
P, N/m2 7500 8104 8700 9300 9620 10200 11200 11700

Emplee la Ley de los gases ideales para determinar R con base a estos datos.
PV = nRT
Donde: T = 273.15 + t

⮚ Nota
Se usó para la presión (P) unidades en pascales (Pa).
Se usó para el volumen (V) unidades en metros cúbicos (m3).
El valor de n = 35.71, el cual corresponde al número de moles que hay en 1kg de nitrógeno
molecular.
Se expresó los valores de la temperatura (T) en kelvin (K).
⮚ Entrada en el editor
1- Clear all clf clc
2- T=[253.15 273.15 293.15 313.15 323.15 343.15 373.15 393.15]’;
3- P=[7500 8104 8700 9300 9620 10200 11200 11700]’;
4- A=[ones(size(T)) T];
5- coef=A\P;
6- a=(coef(2)*10)/35.71
7- TT=(253.15:0.1:393.15)’;
8- PP=[ones(size(TT)) TT]*coef;
9- plot(T,P,’o’,TT,PP,’-’)

⮚ Resultados

3
Pa∗m
a = 8.4896
K∗mol

Este resultado corresponde al valor de la constante de los gases ideales (R).

⮚ Vista del editor


⮚ Gráfica Matlab

⮚ Vista general desde Matlab

También podría gustarte