Minimos Cuadrados

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

TAREA 6.

MÍNIMOS CUADRADOS UNIVARIABLE

Problema 1

En un reactor CSTR (Tipo Tanque Perfectamente Agitado) se obtuvieron

los siguientes datos de conversión ( xA ) y velocidad de reacción ( rA ):

xA 0.4 0.5 0.7 0.8 0.85

rA 0.0040 0.0033 0.0018 0.00175 0.0010

a).- Genere un polinomio de segundo orden que ajuste a estos


datos para conocer el valor de rA cuando xA = 0.6 mediante la
técnica de Mínimos Cuadrados.
b).- Calcule el volumen del reactor para xA = 0.6 empleando para

ello la siguiente ecuación:

𝟏
V = xA( )FA0
𝒓𝑨

Donde FA0 = 4 (Flujo molar de alimentación al reactor)


NOTA.- Utilice cuatro cifras decimales redondeadas en su resultado.
Para resolver el inciso a) acomodamos los valores de xA y rA en columnas
rA XA
0.004 0.4
0.0033 0.5
0.0018 0.7
0.00175 0.8
0.001 0.85
Utilizando el siguiente sistema de ecuaciones para mínimos cuadrados
Ma0 + a1∑x1 + a2∑(xi2) = ∑yi

a0∑(xi) + a1∑(xi2) + a2∑(xi3) = ∑(xi*yi)

a0∑(xi2) + a1∑(xi3) + a2∑(xi4) = ∑(xi2*yi)


Calculando los valores tenemos la siguiente tabla y calculando la sumatoria para
de cada valor

m rA XA X2 X3 X4 X*Y X2*Y
1 0.004 0.4 0.16 0.064 0.0256 0.0016 0.0006
2 0.0033 0.5 0.25 0.125 0.0625 0.0017 0.0008
3 0.0018 0.7 0.49 0.343 0.2401 0.0013 0.0009
4 0.00175 0.8 0.64 0.512 0.4096 0.0014 0.0011
5 0.001 0.85 0.7225 0.6141 0.5220 0.00085 0.0007
Sumatoria 0.01185 3.25 2.2625 1.6581 1.2598 0.0068 0.0042

Sustituyendo en la ecuación los valores obtenidos tenemos

5a0 + 3.25a1 + 2.2625a2 =0.01185


3.25a0 + 2.2625a1 + 1.6581a2 =0.0068
2.2625a0+ 1.6581a1+ 1.2598a2 = 0.0042
Resolviendo el sistema de ecuaciones por Matrices tenemos las siguientes
operaciones
Matriz A Matriz B Matriz inversa de A
5 3.25 2.2625 0.01185 81.5621 -277.0606 218.1818
3.25 2.2625 1.6581 0.0068 -277.0606 953.6364 -757.5758
2.2625 1.6581 1.2598 0.0042 218.1818 -757.5758 606.0606

Multiplicando la matriz inversa de A por la Matriz B tenemos los resultados para


a0, a1 y a2
a0 0.0077
a1 -0.01045
a2 0.0033

Tenemos ahora la ecuación P(X)

P(xi)= a0+a0x+a2x2

Sustituyendo los valores


P(xi)= 0.0077-0.01045x+0.0033x2
Sustituyendo x=0.6 en la ecuación obtenemos el valor de r

P(xi)= 0.0077-0.01045(0.6) + 0.0033(0.6)2

ra=0.0026

Para el inciso b) x=0.6 ra=0.0638 y F=4

𝟏
V = x( )FA0
𝒓𝑨

𝟏
𝑽 = 𝟎. 𝟔 ∗ ( ) ∗ 𝟒 = 𝟗𝟐𝟑. 𝟎𝟕𝟔𝟗𝟐𝟑𝟏~𝟗𝟐𝟑. 𝟎𝟕𝟔𝟖
𝟎.𝟎𝟎𝟐𝟔

V=923.0769

Reactor CSTR
0.0045
0.004
Velocidad de Reacción

0.0035
0.003
0.0025
0.002
0.0015
y = 0.0033x2 - 0.0104x + 0.0077
0.001
0.0005
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9
Conversión

Solución: a) rA=0.0077-0.0105x+0.0033x2
b) V=923.0769
Problema 2

Las siguientes cifras corresponden al tiempo de secado de cierto barniz (y) y a

la cantidad de un aditivo con el que se intenta reducir el tiempo (x):

cantidad
de 0 1 2 3 4 5 6 7 8
aditivo
(gr)
tiempo
de 12.0 10.5 10.0 8.0 7.0 8.0 7.5 8.5 9.0
secado
(hr)

Calcule las constantes del polinomio y = a0 + a1x + a2x2

Tabulando los valores obteniendo y las sumatorias obtenemos la siguiente tabla


de valores

m Y X X2 X3 X4 X*Y X2*Y
1 12 0 0 0 0 0 0
2 10.5 1 1 1 1 10.5 11
3 10 2 4 8 16 20 40
4 8 3 9 27 81 24 72
5 7 4 16 64 256 28 112
6 8 5 25 125 625 40 200
7 7.5 6 36 216 1296 45 270
8 8.5 7 49 343 2401 59.5 417
9 9 8 64 512 4096 72 576
Sumatoria 80.5 36 204 1296 8772 299 1697

Calculamos el polinomio con las ecuaciones sustituyendo los valores y


obteniendo la matriz correspondiente
Ma0 + a1∑x1 + a2∑(xi2) = ∑yi

a0∑(xi) + a1∑(xi2) + a2∑(xi3) = ∑(xi*yi)


a0∑(xi2) + a1∑(xi3) + a2∑(xi4) = ∑(xi2*yi)

Matriz A Matriz B Matriz inversa de A


9 36 204 80.5 0.6606 -0.3091 0.0303
36 204 1296 299 -0.3091 0.2245 -0.0260
204 1296 8772 1697 0.0303 -0.0260 0.0032

Los valores de las constantes del polinomio resultan

a0 12.1848485
a1 -1.8465368
a2 0.18290043

Sustituyendo las contantes en el polinomio y = a0 + a1x + a2x2

Solución: Y= 12.1848 -1.8465x + 0.1829x2

Problema 2
14

12
y = 0.1829x2 - 1.8465x + 12.185
10
Tiempo de secado

0
0 1 2 3 4 5 6 7 8 9
Cantidad de aditivo
Problema 3

Calcule el pH de una solución de suero cuyo potencial (E) es de 20.3 mV de

acuerdo a los datos de la siguiente tabla. Utilice el método de Mínimos

Cuadrados, sabiendo que existe una relación lineal entre el pH y el potencial E.

pH E [mV]

5 -53.8

4 -27.7

3 2.7

2 31.9

1 65.1

Tomando E=x y pH=f(x) y sabiendo que para mínimos cuadrados en su

forma lineal o de primer grado se tiene la forma f(x)=mx+b se calcula el

valor de m (la pendiente) y el valor de b utilizando la siguiente formula

na0 + a1∑x1 = ∑yi

a0∑(xi)+a1∑(xi2) = ∑(xi*yi)

se utilizan los siguientes valores en la formula y se resuelve el


sistema de ecuaciones con dos incógnitas

n E pH X2 x*y
1 -53.8 5 2894.44 -269
2 -27.7 4 767.29 -110.8
3 2.7 3 7.29 8.1
4 31.9 2 1017.61 63.8
5 65.1 1 4238.01 65.1
∑ 18.2 15 8924.64 -242.8
Los valores obtenidos para a0 y a1 son

a0=3.1222

a1=-0.0336

Sustituyendo se obtiene la ecuación lineal

f(x)=-0.0336x + 3.1222

sustituimos el pH = 20.3

f(x)=-0.0336(20.3) + 3.1222

pH=2.4401

Problema 3
6

pH 3

1
y = -0.0336x + 3.1222
0
-60 -40 -20 0 20 40 60 80
Potencial

Solución: pH(20.3) = 2.4

También podría gustarte