Diseño Octogonal

Descargar como xlsx, pdf o txt
Descargar como xlsx, pdf o txt
Está en la página 1de 16

VARIABLES Nivel (-) Nivel (+)

Z1:dosificacion de colector (Lb/Tn) 0.015 0.035


Z2:porcentaje de solidos (%S) 25 30

Z1 Z2
Zj° 0.025 27.5
ΔZj 0.01 2.5

CODIFICADO NATURAL
N° X1 X2 Z1 Z2 Y
1 1.0000 0.0000 0.0350 27.5000 92.20
2 0.7071 0.7071 0.0321 29.2678 94.00
3 0.0000 1.0000 0.0250 30.0000 93.60
4 -0.7071 0.7071 0.0179 29.2678 93.50
5 -1.0000 0.0000 0.0150 27.5000 90.90
6 -0.7071 -0.7071 0.0179 25.7323 92.20
7 0.0000 -1.0000 0.0250 25.0000 92.80
8 0.7071 0.7071 0.0321 29.2678 93.30
9 0.0000 0.0000 0.0250 27.5000 94.10
10 0.0000 0.0000 0.0250 27.5000 94.30
11 0.0000 0.0000 0.0250 27.5000 93.90
12 0.0000 0.0000 0.0250 27.5000 94.00

la matriz X y el vector y estan definidas como:


[X]
X0 X1 X2 X1^2 X2^2 X1.X2
1 1 0 1 0 0
1 0.7071 0.7071 0.5 0.5 0.5
1 0 1 0 1 0
1 -0.7071 0.7071 0.5 0.5 -0.5
1 -1 0 1 0 0
1 -0.7071 -0.7071 0.5 0.5 0.5
1 0 -1 0 1 0
1 0.7071 -0.7071 0.5 0.5 -0.5
1 0 0 0 0 0
1 0 0 0 0 0
1 0 0 0 0 0
1 0 0 0 0 0
la matriz transpuesta [X]T esta dada por:

1 1 1 1 1 1
1 0.7071 0 -0.7071 -1 -0.7071
0 0.7071 1 0.7071 0 -0.7071
1 0.5 0 0.5 1 0.5
0 0.5 1 0.5 0 0.5
0 0.5 0 -0.5 0 0.5

El producto de [X]T•[X] se muestra a continuación:

12 0 0 4 4 0
0 4 0 0 0 0
0 0 4 0 0 0
4 0 0 3 1 0
4 0 0 1 3 0
0 0 0 0 0 1

El cálculo de la matriz inversa ([X]T* [X])-1

0.250 0.000 0.000 -0.250 -0.250 0.000


0.000 0.250 0.000 0.000 0.000 0.000
0.000 0.000 0.250 0.000 0.000 0.000
-0.250 0.000 0.000 0.625 0.125 0.000
-0.250 0.000 0.000 0.125 0.625 0.000
0.000 0.000 0.000 0.000 0.000 1.000

El producto de la matriz transpuesta de [X]T por el vector de respuestas [Y]

[X]T
1 1 1 1 1 1
1 0.7071 0 -0.7071 -1 -0.7071
0 0.7071 1 0.7071 0 -0.7071
1 0.5 0 0.5 1 0.5
0 0.5 1 0.5 0 0.5
0 0.5 0 -0.5 0 0.5
Luego aplicando la (Ec.70) encontramos los coeficientes del modelo matemático
([X]T* [X])-1
0.250 0.000 0.000 -0.250 -0.250 0.000
0.000 0.250 0.000 0.000 0.000 0.000
[B]= 0.000 0.000 0.250 0.000 0.000 0.000
-0.250 0.000 0.000 0.625 0.125 0.000
-0.250 0.000 0.000 0.125 0.625 0.000
0.000 0.000 0.000 0.000 0.000 1.000

El modelo matemático para el diseño octogonal propuesto a escala codificada es la sigu

y= 94.0750 + 0.6078X1 + 0.5536X2 - 2.0875X1^2 - 0.4375X2^2 - 0.3X1.X2

aplicando ecucaciones para el modelo:

N° Z1 Z2
(-) 0.015 25
(+) 0.035 30
Z° 0.025 27.5
ΔZ 0.01 2.5
ε 2.5 11

entonces para a0 tenemos:


a0= b0-(b1ε1) - (b2.ε2) + (b11.ε1^2) + (b22.ε2^2) + (b12.ε1.ε2)
a0=
94.750 - (0.6878*2.5) - 0.5536*11) + (-0.20875*(2.5)^2) + (-0.4375*(11)^2) + (-0.3*2.5

para los terminos lineales:

a1Z1 = (0.6078/0.01)-2((-2.0875*2.5)/0.01)-((-0.3*11)/0.01) = 1434.53458


a1Z1 = 1434.5346

a2Z2= (0.5536/2.5)-2((-0.4375*11)/2.5)-((-0.3*2.5)/2.5) = 4.37142


a2Z2= 4.3714221

para los terminos cuadraticos:

a11Z1^2= (-2.0875/(0.01)^2) = -20875


a11Z1^2= -20875

a22Z2^2= (-0.4375/(2.5)^2) = -0.07


a22Z2^2= -0.07
para los terminos de interacción empleamos:

a12Z1Z2= (-0.3/(0.01*2.5)) = -12


a12Z1Z2= -12

el modelo matematico a escala natural es:

Y = 12.2319-1434.53458.Z1 + 4.37142.Z2 - 20875.Z1^2 - 0.07.Z2^2 - 12.Z1Z2

Análisis de varianza para diseños rotables y compuestos


. Para los datos del problema del diseño octogonal para doce pruebas experimentales su
Cjj
0.250
0.250
0.250
0.625
0.625
1.000

calculando elerror experimental:

N° X1 X2 Z1 Z2 Y
9 0 0 0.025 27.5 94.1
10 0 0 0.025 27.5 94.3
11 0 0 0.025 27.5 93.9
12 0 0 0.025 27.5 94

y°= 94.075
SSE= 0.0875
G° de libertad = 4 - 1 =3
MSE= 0.0291667
desviacion estandar del error experimental es: se= 0.17078
Ahora podemos calcular la Desviación estándar asociada a los coeficientes:
Sb1 = 0.085391
Sb11 = 0.135015
Sb12 = 0.170783
Calculo de la significancia de los coeficientes
tj = (b1/Sb1)
tj = 7.1183615

En la siguiente tabla se muestran el resumen para los demás coeficientes:


Se Cjj Sbj bj tj
X0 0.17078 0.250 0.085391 94.0750 1101.694
X1 0.17078 0.250 0.085391 0.6078 7.118
X2 0.17078 0.250 0.085391 0.5536 6.483
X1^2 0.17078 0.625 0.135015 -2.0875 -15.461
X2^2 0.17078 0.625 0.135015 -0.4375 -3.240
X1.X2 0.17078 1.000 0.170783 -0.3000 -1.757

Evaluación del modelo matemático

[X]
X0 X1 X2 X1^2 X2^2 X1.X2
1 1 0 1 0 0
1 0.7071 0.7071 0.5 0.5 0.5
1 0 1 0 1 0
1 -0.7071 0.7071 0.5 0.5 -0.5
1 -1 0 1 0 0
1 -0.7071 -0.7071 0.5 0.5 0.5
1 0 -1 0 1 0
1 0.7071 -0.7071 0.5 0.5 -0.5
1 0 0 0 0 0
1 0 0 0 0 0
1 0 0 0 0 0
1 0 0 0 0 0

Recordando que la suma de cuadrados del error experimental se calcula cuando se dete
grados de libertad f3 es la diferencia de los grados de libertad del residual (f 1) menos los

MSFA= 1.6327/3 = 0.5442

Por ultimo calculamos el Fo del modelo

Fo = 18.658
entonces tenemos:
F(3,3)99% = 29.46
Fo=18.659 < F=29.46 el modelo matemático ajusta o representa adecuadamen
El siguiente cuadro se resume el análisis de Varianza del modelo matemático.

fuente de variacion suma de cuadrados grados de libertad


error falta de ajuste 1.6327 3
error experimental 0.0875 3
error residual 1.7202 6
DISEÑO EXPERIMENTAL OCTOGONAL
1 1 1 1 1 1
0 0.7071 0 0 0 0
-1 -0.7071 0 0 0 0
0 0.5 0 0 0 0
1 0.5 0 0 0 0
0 -0.5 0 0 0 0

respuestas [Y]

[Y]
1 1 1 1 1 1 92.2
0 0.7071 0 0 0 0 94
-1 -0.7071 0 0 0 0 * 93.6
0 0.5 0 0 0 0 93.5
1 0.5 0 0 0 0 90.9
0 -0.5 0 0 0 0 92.2
92.8
93.3
94.1
94.3
93.9
94
[X]T*[Y] [bj] [bj]
1118.8 b0 94.0750
2.43136 b1 0.6078
* 2.2142 = b2 = 0.5536
369.6 b11 -2.0875
372.9 b22 -0.4375
-0.3 b12 -0.3000

escala codificada es la siguiente

75X2^2 - 0.3X1.X2

+ (-0.4375*(11)^2) + (-0.3*2.5*1)122319
=

1434.53458
2^2 - 12.Z1Z2

e pruebas experimentales su diagonal siempre será:

os coeficientes:

significancia
SIGNIF.
SIGNIF.
SIGNIF. t (0.975,3)=3.18
SIGNIF.
SIGNIF.
NO SIGNIF.

Calculo delos residuales:

[B] [Ȳ] Ȳ Y (Y-Ȳ )


94.0750 92.5953 92.5953 92.2 -0.3953
0.6078 93.4837 93.4837 94 0.5163
0.5536 94.1911 94.1911 93.6 -0.5911
-2.0875 92.9241 92.9241 93.5 0.5759
* -0.4375 = 91.3797 91.3797 90.9 -0.4797
-0.3000 91.8413 91.8413 92.2 0.3587
93.0839 93.0839 92.8 -0.2839
93.0009 93.0009 93.3 0.2991
94.0750 94.0750 94.1 0.0250
94.0750 94.0750 94.3 0.2250
94.0750 94.0750 93.9 -0.1750
94.0750 94.0750 94 -0.0750
SSR=

La media de cuadrados del residual


tanto los grados de libertad son (12 - 6 = 6)
MSR= 0.286696

al se calcula cuando se determina la significancia de los coeficientes, su valor es (SS E = 0.0875) y su media d
d del residual (f 1) menos los grados de libertad del error experimental (f 2), luego: f3=6-3=3

a o representa adecuadamente los datos experimentales.


media de cuadrados Fo aceptar modelo
0.544233 18.6594 29.46 SI
0.029167
0.286700
[X]T*[Y]
1118.8
2.43136
= 2.2142
369.6
372.9
-0.3
(Y-Ȳ )^2
0.1563
0.2665
0.3493
0.3316
0.2301
0.1287
0.0806
0.0895
0.0006
0.0506
0.0306
0.0056
1.7202

del residual
rtad son (12 - 6 = 6)

S E = 0.0875) y su media de cuadrados MSE = 0.0291667.


1 1 0 1 0 0
1 0.5 0.866 0.25 0.74995 0.433
1 -0.5 0.866 0.25 0.74995 -0.433
1 -1 0 1 0 0
1 -0.5 -0.866 0.25 0.74995 0.433
1 0.5 -0.866 0.25 0.74995 -0.433
1 0 0 0 0 0
1 0 0 0 0 0
1 0 0 0 0 0

1 1 1 1 1 1 1
1 0.5 -0.5 -1 -0.5 0.5 0
0 0.866 0.866 0 -0.866 -0.866 0
1 0.25 0.25 1 0.25 0.25 0
0 0.74995 0.74995 0 0.74995 0.74995 0
0 0.433 -0.433 0 0.433 -0.433 0
1 1
0 0
0 0
0 0
0 0
0 0

También podría gustarte