Scribd
Cargar
262 vistas10 páginas

Ajuste de Curvas en Ingeniería Química

Este documento presenta tres ejemplos de ajuste de curvas y regresión utilizando Matlab. El primer ejemplo ajusta una curva exponencial a datos para estimar las concentraciones iniciales de tres organismos. El segundo ejemplo usa regresión no lineal para determinar los parámetros de una ecuación que modela datos tabulados. El tercer ejemplo emplea la ley de los gases ideales para determinar la constante de los gases R a partir de datos de presión y temperatura de nitrógeno.

Cargado por

JORGE LUIS HERRERA ARROYO
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
262 vistas10 páginas

Ajuste de Curvas en Ingeniería Química

Este documento presenta tres ejemplos de ajuste de curvas y regresión utilizando Matlab. El primer ejemplo ajusta una curva exponencial a datos para estimar las concentraciones iniciales de tres organismos. El segundo ejemplo usa regresión no lineal para determinar los parámetros de una ecuación que modela datos tabulados. El tercer ejemplo emplea la ley de los gases ideales para determinar la constante de los gases R a partir de datos de presión y temperatura de nitrógeno.

Cargado por

JORGE LUIS HERRERA ARROYO
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd

TALLER #3

Informática para ingeniería química

AJUSTE E INTERPOLACIÓN DE CURVAS

Presentado por:

Presentado a:

PROFESOR FRANCISCO MUÑOZ PABA [Link]

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