Minimos Cuadrados Multivariable Ejercicios Resueltos

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

TAREA 8.

MÍNIMOS CUADRADOS MULTIVARIABLE


Problema 7

Diversos tipos de hojas de acero procesadas en frío tienen diferentes

composiciones de cobre y temperaturas de templado. Al medir su dureza

resultante se obtuvieron los siguientes valores:

Dureza Rockell 30-T Contenido de Cobre Temperatura de


[%] templado [°F]
78.9 0.02 1000
65.1 0.02 1100
55.2 0.02 1200
56.4 0.02 1300
80.9 0.10 1000
69.7 0.10 1100
57.4 0.10 1200
55.4 0.10 1300
85.3 0.18 1000
71.8 0.18 1100
60.7 0.18 1200
58.9 0.18 1300

Si se sabe que la dureza depende en forma lineal del contenido de Cobre ( u) y

de la temperatura de templado (v) de la siguiente manera:

y = a0 + a1 u + a2 v

Calcule los parámetros a0, a1 y a2 siguiendo el criterio de los Mínimos

Cuadrados.

Asignamos la variable “y” a la dureza, “u” al contenido de cobre, y “v” a la


temperatura de templado asignamos estos valores y hacemos las operaciones
correspondientes para obtener las siguientes tablas de valores

Dureza Contenid temperatur


Rockell o de a de
30-t cobre templado
Puntos y u v
1 78.9 0.02 1000
2 65.1 0.02 1100
3 55.2 0.02 1200
4 56.4 0.02 1300
5 80.9 0.1 1000
6 69.7 0.1 1100
7 57.4 0.1 1200
8 55.4 0.1 1300
9 85.3 0.18 1000
10 71.8 0.18 1100
11 60.7 0.18 1200
12 58.9 0.18 1300
Sumatori
795.7 1.2 13800
a

Operaciones
2
Puntos u u*v u*y v2 v*y
1 0.0004 20 1.578 1000000 78900
2 0.0004 22 1.302 1210000 71610
3 0.0004 24 1.104 1440000 66240
4 0.0004 26 1.128 1690000 73320
5 0.01 100 8.09 1000000 80900
6 0.01 110 6.97 1210000 76670
7 0.01 120 5.74 1440000 68880
8 0.01 130 5.54 1690000 72020
9 0.0324 180 15.354 1000000 85300
10 0.0324 198 12.924 1210000 78980
11 0.0324 216 10.926 1440000 72840
12 0.0324 234 10.602 1690000 76570
Sumatori
0.1712 1380 81.258 16020000 902230
a

Con las sumatorias pasamos a armar el sistema de ecuaciones para calcular el

valor de a0, a1, y a2

m a0 + a1u + a2v = y

a0u + a1u2 + a2uv = uy


a0v + a1uv + a2v2 = vy

Resolviendo el sistema de ecuaciones por matrices tenemos que

Matriz A Matriz B Matriz inversa de A


12 1.2 13800 795.7
9.0953125 -1.953125 -0.00766667
1.2 0.1712 1380 81.258-1.953125 19.53125 -4.1028E-17
1380 -
1380 16020000 902230 -3.3249E-17 6.6667E-06
0 0.00766667

Obtenemos los valores de a0, a1, y a2

a0 =161.3365

a1 =32.9688

a2 = -0.0855
dando la ecuación resultante

y = 161.3365 + 32.9688 u - 0.0855v

Problema 8
Veinte tipos de hojas de acero procesadas en frío tienen diferentes

composiciones de Cobre a diferentes temperaturas de templado. Al medir su

dureza, se obtuvieron los siguientes valores:

Contenido de Temperatura de templado [°F]


Cobre [%] 1000 1100 1200 1300
0.02 78.9 65.1 55.2 56.4
0.10 80.9 69.7 57.95 55.4
0.18 85.3 71.8 60.7 58.9

Realice un ajuste lineal por mínimos cuadrados para calcular la dureza cuando el

contenido de Cobre es de 0.05% a una temperatura de templado de 1125°F.

Solución:

Debido a que son los mismos valores para la dureza contenido de cobre y

temperatura respectivamente

El valor de a0, a1, y a2 no cambia por lo que se sustituye en la ecuación los

valores de 0.05% y 1125°F para conocer la dureza del material

y = 161.3365 + 32.9688 u - 0.0855v

y = 161.3365 + 32.9688(0.05) - 0.0855(1125)

Y = 67.7865

La dureza del material cuando tiene un contenido de cobre de 0.05% y esta

templado a una temperatura de 1125°F es de 67.7865

Problema 9
En un estudio veterinario se midió el porcentaje de sobrevivencia de

cierto tipo de semen animal en varias combinaciones de concentraciones de

tres materiales que se utilizan para aumentar su oportunidad de sobrevivencia.

Los datos son los siguientes:


 
y(% sobrevivencia) x1(peso %) x2(peso %) x3(peso %)

25.5 1.74 5.30 10.80

31.2 6.32 5.42 9.40

25.9 6.22 8.41 7.20

38.4 10.52 4.63 8.50

18.4 1.19 11.60 9.40

26.7 1.22 5.85 9.90

26.4 4.10 6.62 8

25.9 6.32 8.72 9.10


32 4.08 4.42 8.70

25.2 4.15 7.60 9.20

39.7 10.15 4.83 9.40

35.7 1.72 3.12 7.60

26.5 1.70 5.30 8.20

Calcule los coeficientes de regresión del siguiente modelo, usando la técnica de

Mínimos Cuadrados:

y = b0 + b1 x1 + b2 x2 + b3 x3

cambiando el nombre de las variables para acomodar de manera mas fácil a la

ecuación de mínimos cuadrados multivariables tenemos que:

x1=u, x2=v, x3=z


Al observar la ecuación de mínimos cuadrados nos notamos que nos hacen falta

calcular las sumatorias de u,v,z,y,u2, v2, u*v, u*y, u*z….. etc

m a0 + a1u + a2v + a3z = y

a0u + a1u2 + a2uv + a3uz = uy

a0v + a1uv + a2v2 + a3vz = vy

a0z + a1uz + a2vz + a3z2 = zy

Obteniendo los valores de sumatorias y productos tenemos las tablas

Y (%
sobrevivencia x1 (peso %) x2 (peso %) x3 (peso %)
Punto )
y u v z
1 25.5 1.74 5.3 10.8
2 31.2 6.32 5.42 9.4
3 25.9 6.22 8.41 7.2
4 38.4 10.52 4.63 8.5
5 18.4 1.19 11.6 9.4
6 26.7 1.22 5.85 9.9
7 26.4 4.1 6.62 8
8 25.9 6.32 8.72 9.1
9 32 4.08 4.42 8.7
10 25.2 4.15 7.6 9.2
11 39.7 10.15 4.83 9.4
12 35.7 1.72 3.12 7.6
13 26.5 1.7 5.3 8.2
Sumatoria 377.5 59.43 81.82 115.4

P. u2 u*v u*z U*y v2 v*z v*y z2 z*y


1 3.0276 9.222 18.792 44.37 28.09 57.24 135.15 116.64 275.4
2 39.9424 34.2544 59.408 197.184 29.3764 50.948 169.104 88.36 293.28
3 38.6884 52.3102 44.784 161.098 70.7281 60.552 217.819 51.84 186.48
4 110.6704 48.7076 89.42 403.968 21.4369 39.355 177.792 72.25 326.4
5 1.4161 13.804 11.186 21.896 134.56 109.04 213.44 88.36 172.96
6 1.4884 7.137 12.078 32.574 34.2225 57.915 156.195 98.01 264.33
7 16.81 27.142 32.8 108.24 43.8244 52.96 174.768 64 211.2
8 39.9424 55.1104 57.512 163.688 76.0384 79.352 225.848 82.81 235.69
9 16.6464 18.0336 35.496 130.56 19.5364 38.454 141.44 75.69 278.4
10 17.2225 31.54 38.18 104.58 57.76 69.92 191.52 84.64 231.84
11 103.0225 49.0245 95.41 402.955 23.3289 45.402 191.751 88.36 373.18
12 2.9584 5.3664 13.072 61.404 9.7344 23.712 111.384 57.76 271.32
13 2.89 9.01 13.94 45.05 28.09 43.46 140.45 67.24 217.3
∑ 394.7255 360.6621 522.078 1877.567 576.7264 728.31 2246.661 1035.96 3337.78

m a0 + a1u + a2v + a3z = y

a0u + a1u2 + a2uv + a3uz = uy

a0v + a1uv + a2v2 + a3vz = vy

a0z + a1uz + a2vz + a3z2 = zy

Sustituyendo en la ecuación los valores obtenidos tenemos la siguiente matriz

Matriz A Matriz B
13 59.43 81.82 115.4 377.5
59.43 394.7255 360.6621 522.078 1877.567
81.82 360.6621 576.7264 728.31 2246.661
115.4 522.078 728.31 1035.96 3337.78

Resolviendo la matriz obtenemos los valores de

a0 39.15735
a1 1.01610044
a2 -1.8616492
-
a3 0.34326049

Sustituyendo en la ecuación tenemos

y = 39.15735 + 1.0161 x1 – 1.8616 x2 + 0.3432 x3

También podría gustarte