DISEO EXPERIMENTAL OCTOGONAL

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 continuacin:

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 clculo 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 matemtico

([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 matemtico para el diseo octogonal propuesto a escala codificada es la siguiente

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-(b11) - (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*1) =

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.5345829

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

a2Z2= 4.3714221234

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 interaccin 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

Anlisis de varianza para diseos rotables y compuestos

. Para los datos del problema del diseo octogonal para doce pruebas experimentales su diagonal siempre ser:
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.0291666667
desviacion estandar del error experimental es: se= 0.17078
Ahora podemos calcular la Desviacin estndar 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.1183614692

En la siguiente tabla se muestran el resumen para los dems 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

Evaluacin del modelo matemtico

 [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

Calculo delos residuales:

Y (Y- ) (Y- )^2

92.5953 92.2 -0.3953 0.1563
93.4837 94 0.5163 0.2665
94.1911 93.6 -0.5911 0.3493
92.9241 93.5 0.5759 0.3316
91.3797 90.9 -0.4797 0.2301
91.8413 92.2 0.3587 0.1287
93.0839 92.8 -0.2839 0.0806
93.0009 93.3 0.2991 0.0895
94.0750 94.1 0.0250 0.0006
94.0750 94.3 0.2250 0.0506
94.0750 93.9 -0.1750 0.0306
94.0750 94 -0.0750 0.0056
SSR= 1.7202

La media de cuadrados del residual

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

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

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 matemtico ajusta o representa adecuadamente los datos

El siguiente cuadro se resume el anlisis de Varianza del modelo matemtico.

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

Anlisis de Mximos y Mnimos

a^2. Y/a.x1^2= -4.175 a^2/ax1^2=0.875 a^2.Y/ax1.x2 = -0.3

Remplazando estos valores en la matriz Hessiana tenemos

H11= -4.175

H22= -4.175 -0.3

-0.3 -0.875

Luego H22 = [(-4.175)(-0.875)] - [(0.3)(0.3) ] = +3.563

Como podemos observar, el signo de las determinantes es alternativo (H 11 = -4.175) (H22 = +3.563).

Soluciones ptimas Sin Restricciones

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

empleamos la ecuacin siguiente para el clculo de las variables X 1 y X2 y lo resolveremos matricialmente:

[A] [B]
-4.175 -0.3 -0.6078
-0.3 -0.875 -0.5536
[A]-1 [B] [X]
-0.245570952 0.084195755 * -0.6078 = X1
0.084195755 -1.171724259 -0.5536 X2

hallamos el centro y radio del diseo Octogonal

Z1 Z2
Zj 0.025 27.5
Zj 0.01 2.5
El valor de las variables en escala natural se halla por medio de

Zj = Zj + XjZj

X1 = 0.102647255 X2 =0.597492370
Z1= 0.025+(0.102647255*0.01) = 0.02602647
Z2=27.5+(0.597492370*2.5) =28.9937309

Las condiciones ptimas son:

Dosificacin de colector 0.026 (Lb/tn)
% de slidos 28.9936

Remplazando estos valores en el modelo matemtico a escala natural o a escala cod

 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

[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

es la siguiente

^2 - 0.3X1.X2

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

3458
ntales su diagonal siempre ser:

significancia
SIGNIF.
SIGNIF.
SIGNIF. t (0.975,3)=3.18
SIGNIF.
SIGNIF.
NO SIGNIF.
 [B] []
94.0750 92.5953
0.6078 93.4837
0.5536 94.1911
-2.0875 92.9241
* -0.4375 = 91.3797
-0.3000 91.8413
93.0839
93.0009
94.0750
94.0750
94.0750
94.0750

o se determina la significancia de los coeficientes, su valor es (SS E = 0.0875) y su media de cuadrados MSE = 0.0291667.
enos los grados de libertad del error experimental (f 2), luego: f3=6-3=3
presenta adecuadamente los datos experimentales.

media de cuadrados Fo aceptar modelo

0.544233 18.6594 29.46 SI
0.029167
0.286700

175) (H22 = +3.563).

^2 - 0.3X1.X2

solveremos matricialmente:

[X]=[A]-1*[B]
= ###
###
scala natural o a escala codificado se obtiene la mxima recuperacin que corresponde a 94.27%
 [X]T*[Y]
1118.8
2.43136
= 2.2142
369.6
372.9
-0.3
SE = 0.0291667.
onde a 94.27%
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