Algoritmo Simplex
Algoritmo Simplex
Prueba de Optimalidad:
FO
Max Min
Todo (Cj-Zj)<=0 Todo (Cj-Zj)>=0
CASOS DE INFACTIBILIDAD:
1. Cuando existe una VARIABLE ARTIFICIAL como VARIABLE BASICA
2. Cuando existe una VARIABLE BÁSICA NEGATIVA
SIGNO DEL PRECIO DUAL
Signo
FO Restricción PD
MAX ó <= +
MIN >= -
El uso de una variable artificial permite crear una soluci[on factible básica para
iniciar el algoritmo simplex (genera matriz identidad).
Una variable artificial no tiene significado físico en términos de un problema de
programación lineal del mundo real, por lo tanto no se le permite aparecer en la
solución final del problema, en caso contrario estamos en un caso de solución no
factible.
El uso de una variable artificial permite crear una soluci[on factible básica para
iniciar el algoritmo simplex (genera matriz identidad).
Una variable artificial no tiene significado físico en términos de un problema de
programación lineal del mundo real, por lo tanto no se le permite aparecer en la
solución final del problema, en caso contrario estamos en un caso de solución no
factible.
ndo. Se construye la tabla inicial del simplex Tercero. Se obtiene la PRIMERA IT
F25*(-2)
4. Se hace la PRUEBA DE
5. De no haber llegado a
ra o precio dual:
Fila Cj-Zj de la Columna S1
Fila Cj-Zj de la Columna S2
Fila Cj-Zj de la Columna S3
Fila Cj-Zj de la Columna S4
Fila Cj-Zj de la Columna S5
RELACIÓN PRIMAL DUAL
MODELO PRIMAL: MODELO DUAL:
MIN W = 430Y1 + 460Y2
Problema 2.- Para el modelo matemático siguiente, el cual representa un problema de producción: Y1+ 3Y2 + Y3 >= 3
Y1 + 4Y3 >= 2
MAX Z= 3X1+2X2+5X3 ganancia en soles MODELO ESTANDARIZADO: Y1+ 2Y2 >= 5
ST MAX Z= 3X1+2X2+5X3+0S1+0S2+0S3
X1+X2+X3<=430 horas-hombre ST
3X1+2X3<=460 horas-máquina X1+X2+X3+S1=430
X1+4X2<=420 Kg. de Materia Prima 3X1+2X3+S2=460
X1,X2,X3>=0 X1+4X2+S3=420
Método Simplex
Cj ®
Cj ® 3 2 5 0 0 0 ¯
¯ VB x1 x2 x3 s1 s2 s3 b q 1000
0 s1 1 1 1 1 0 0 430 430 1000
0 s2 3 0 2 0 1 0 460 230 1000
0 s3 1 4 0 0 0 1 420 M
Zj 0 0 0 0 0 0 0
Cj-Zj 3 2 5 0 0 0
Determinar:
a) El plan de producción óptimo. La ganancia máxima.
Plan de producción óptimo:
0 unidades de X1
105 unidades de X2
230 unidades de X3
Ganancia máxima = 1360 soles
b) Los precios sombra de cada Lado Derecho
PS(LD1)=0
PS(LD2)=+2.5
PS(LD3)=+0.5
c) La holgura o excedente de cada LD
Holgura de LD1=95 horas hombre
Holgura de LD2=0 horas máquina
Holgura de LD3=0 Kg. De materia prima
d) ¿Cuánto estaría Ud. dispuesto a pagar como máximo por una hora-hombre extra?,
¿Por una hora-máquina extra?, ¿Por un Kg. adicional de Materia Prima?
Por una hora hombre extra = 0 soles
Por una hora máquina extra = 2.5 soles
Por un kg. adicional de materia prima = 0.5 soles
e) ¿Cuánto debería de aumentar la utilidad unitaria de x1 para que justifique su producción?
Debería de aumentar una cantidad >=5 soles.
MODELO DUAL:
MIN W = 430Y1 + 460Y2 + 420Y3
Y1+ 3Y2 + Y3 >= 3
Y1 + 4Y3 >= 2
Y1+ 2Y2 >= 5
Modelo Estandarizado
MIN W = 430Y1 + 460Y2 + 420Y3 + 0S1 + 0S2 + 0S3 + MA1 + MA2 + MA3
ST
Y1+ 3Y2 + Y3 -S1 +1000A1 = 3
Y1 + 4Y3 - S2 + 1000A2 = 2
Y1 + 2Y2 -S3 + 1000A3 = 5
M al ser un valor sumamente grande, le damos el valor 1000
VB Y1 Y2 Y3 S1 S2 S3 A1 A2 A3 b q
A1 0.75 3 0 -1 0.25 0 1 -0.25 0 2.5 0.8333 Elemento Nuevo (
Y3 0.25 0 1 0 -0.25 0 0 0.25 0 0.5 inf
A3 1 2 0 0 0 -1 0 0 1 5 2.5
Zj --> 1855 5000 420 -1000 145 -1000 1000 -145 1000 7710 Para encontra
Cj - Zj -- -1425 -4540 0 1000 -145 1000 0 1145 0 los elemntos
partir de la 2d
VB Y1 Y2 Y3 S1 S2 S3 A1 A2 A3 b q tabla, excepto
Y2 0.25 1 0 -0.3333 0.0833 0 0.3333 -0.0833 0 0.8333 -2.5 de la fila pivo
Y3 0.25 0 1 0 -0.25 0 0 0.25 0 0.5 inf puede aplicar
A3 0.5 0 0 0.6667 -0.1667 -1 -0.6667 0.1667 1 3.3333 5 siguiente
Zj --> 720 460 420 513.33 -233.33 -1000 -513.33 233.33 1000 3926.7 relación:
Cj - Zj -- -290 0 0 -513.33 233.33 1000 1513.3 766.67 0
EN=ET-SPF*SP
VB Y1 Y2 Y3 S1 S2 S3 A1 A2 A3 b
Y2 1 1 0 0 0 -1 0 0 1 3
Y3 0.25 0 1 0 -0.25 0 0 0.25 0 0.5
S1 0.75 0 0 1 -0.25 -1.5 -1 0.25 1.5 5
Zj --> 335 460 420 0 -105 -230 0 105 230 1360
Cj - Zj -- 95 0 0 0 105 230 1000 895 770
Costos reducidos Precios duales
Solución del Modelo DUAL: Solución del Modelo PRIMAL:
Y1 0 X1 0
Y2 2 X2 105
Y3 0.5 X3 230
W 1360 Z 1360
S1 5 S1 95
S2 0 S2 0
S3 0 S3 0
PS1 0 PS1 0
PS2 -105 PS2 2
PS3 -230 PS3 0.5
CR Y1 95 CR X1 5
CR Y2 0 CR X2 0
CR Y3 0 CR X3 0
Semipivot de Columna (SPC)
Elemento a Transformar
(ET)
Pivot (P)
Para encontrar
los elemntos a
partir de la 2da
tabla, excepto los
de la fila pivot, se
puede aplicar la
siguiente
relación:
EN=ET-SPF*SPC/P
Problema 3.- Para el modelo matemático siguiente, el cual representa un problema de minimización de costos:
MIN Z = 4X1 + 3X2
ST Modelo Estandarizado:
X1+3X2>=15 MIN Z = 4X1 + 3X2+0s1+0s2+0s3+Ma1+Ma2+Ma3
2X1+ X2>=10 ST
X1+ X2>=8 X1+3X2-s1+a1=15
X1,X2>=0 2X1+ X2-s2+a2=10
X1+ X2-s3+a3=8
Se pide:
a) El gráfico respectivo y la solución óptima (valor de las variables de decisión y valor de la función objetivo)
Z=26
X1=2
X2=6
Método Simplex
Cj ® 4 3 0 0 0 M M M
¯ VB x1 x2 s1 s2 s3 a1 a2 a3 b q
M a1 1 3 -1 0 0 1 0 0 15 5
I M a2 2 1 0 -1 0 0 1 0 10 10
M a3 1 1 0 0 -1 0 0 1 8 8
Zj 4M 5M -M -M -M M M M 33M
Cj-Zj 4-4M 3-5M M M M 0 0 0
3 x2 0 1 0 1 -2 0 -1 2 6
IV 4 x1 1 0 0 -1 1 0 1 -1 2
0 s1 0 0 1 2 -5 -1 -2 5 5
Zj 4 3 0 -1 -2 0 1 2 26
Cj-Zj 0 0 0 1 2 M M-1 M-2
Precios Sombra
Cj ® 15 10 8 0 0
¯ VB y1 y2 y3 s1 s2 b q
0 s1 1 2 1 1 0 4 4
I 0 s2 3 1 1 0 1 3 1
Zj 0 0 0 0 0 0
Cj-Zj 15 10 8 0 0
10 y2 -2 1 0 1 -1 1
IV 8 y3 5 0 1 -1 2 2
Zj 20 10 8 2 6 26
Cj-Zj -5 0 0 -2 -6
Precios Sombra
a) Muestre la gráfica respectiva y la tabla de evaluación de los puntos intersección de la región factible.
Tabla:
Punto x1 x2 Z
P1 3 5 21
P2 6,5 1,5 17,5
P3 5 1 13 Min
Z’ = 2(5.4284)+3(1.1428) = 14.26
Þ PD(RHS2)= 13 - 14.2852= -1.285
d) Muestre la salida del WinQsb:
Método Simplex
Cj ® 2 3 0 0 0 M
¯ VB x1 x2 s1 s2 s3 a1 b q
0 s1 1 -3 1 0 0 0 2 2
I M a1 2 1 0 -1 0 1 11 5.5
0 s3 1 1 0 0 1 0 8 8
Zj 2M M 0 -M 0 M 11M
Cj-Zj 2-2M 3-M 0 M 0 0
2 x1 1 -3 1 0 0 0 2 -0.67
II M a1 0 7 -2 -1 0 1 7 1
0 s3 0 4 -1 0 1 0 6 1.5
Zj 2 -6+7M 2-2M -M 0 M 4+7M
Cj-Zj 0 9-7M -2+2M M 0 0
Solución óptima:
x1=5
x2=1
Z=13
Estado de los lados derechos:
s1=0 (agotado)
s2=0 (agotado)
s3=2 (abundante)
Precios sombra:
PS1=0.571
PS2=-1.29
Ps3=0
Costos reducidos:
CR1=0
CR2=0
CASO ESPECIAL: SOLUCIONES OPTIMAS MULTIPLES
Modelo Estandarizado:
Max 6x1+2x2+0s1+0s2+0s3-Ma1
st
3x1+x2+s1=48
3x+4x2+s2=120
3x1+x2-s3+a1=36
Cj ® 6 2 0 0 0 -M
¯ VB x1 x2 s1 s2 s3 a1 b q
0 s1 3 1 1 0 0 0 48 16
I 0 s2 3 4 0 1 0 0 120 40
-M a1 3 1 0 0 -1 1 36 12
Zj -3M -M 0 0 M -M -36M
Cj-Zj 6+3M 2+M -2 0 -M ###
0 s1 0 0 1 0 1 -1 12 12
II 0 s2 0 3 0 1 1 -1 84 84
6 x1 1 0.333 0 0 -0.333 0.333 12 -36
Zj 6 2 0 0 -2 2 72
Cj-Zj 0 0 0 0 2 -M-2
0 s3 0 0 1 0 1 -1 12 INF
III 0 s2 0 3 -1 1 0 0 72 24
6 x1 1 0.333 0.333 0 -3E-05 3E-05 16 48
Zj 6 2 2 0 0 0 96
Cj-Zj 0 0 -2 0 0 -M (Cj-Zj)=0 para x2 que es una VNB
0 S3 0 0 1 0 1 -1 12
IV 2 X2 0 1 -0.333 0.333 0 0 24
6 X1 1 3E-05 0.444 -0.111 -3E-05 3E-05 8.0004
Zj 6 2 2 7E-05 0 0 96
Cj-Zj 0 0 -2 0 0 -M (Cj-Zj)=0 para x1 que es una VNB
0
0
II 0
40
-M
0
0
III 30
40
-M
0
0
IV 30
40
-M
Método Simplex
40 30 0 0 0 0 0 -M -M
VB x1 x2 s1 s2 s3 s4 s5 a1 a2 b q
s1 0.4 0.2 1 0 0 0 0 0 0 20 50
s2 0 0.2 0 1 0 0 0 0 0 5 INF
s3 0.6 0.3 0 0 1 0 0 0 0 21 35
a1 1 0 0 0 0 -1 0 1 0 30 30
a2 0 1 0 0 0 0 -1 0 1 15 INF
Zj -M -M 0 0 0 M M -M -M -45M
Cj-Zj 40+M 30+M 0 0 0 -M -M 0 0
s1 0 0 1 0 -0.666667 0 0 0 0 6 INF
s2 0 0 0 1 -0.666667 -0.4 0 0.4 0 3 7.5
x2 0 1 0 0 3.333333 2 0 -2 0 10 -5
x1 1 0 0 0 0 -1 0 1 0 30 30
a2 0 0 0 0 -3.333333 -2 -1 2 1 5 2.5
Zj 40 30 0 0 100+3.3M 20+2M M -20-2M -M 1500-5M
Cj-Zj 0 0 0 0 -100-3.3M -20-2M -M 20+M 0
s1 0 0 1 0 -0.666667 0 0 0 0 6
s2 0 0 0 1 0 0 0.2 0 -0.2 2
x2 0 1 0 0 0 0 -1 0 1 15
x1 1 0 0 0 1.666667 0 0.5 0 -0.5 27.5
a1 0 0 0 0 -1.666667 -1 -0.5 1 0.5 2.5
Zj 40 30 0 0 66.6+1.66M M -10+0.5M -M 10-0.5M 1550
Cj-Zj 0 0 0 0 -66.6-1.66M -M 10-0.5M 0 -10-0.5M <=0, Solución óptima, PERO EXISTE una Variable Arti
tima, PERO EXISTE una Variable Artificial como VB
CASO ESPECIAL: SOLUCIÓN INFINITA Método Simplex
Modelo estandarizado: Cj ® 1 2
Max x1+2x2+0s1+0s2 ¯ VB x1 x2
st 0 s1 -4 3
-4x1+3x2+s1=3 I 0 s2 1 -1
x1-x2+s2=3 Zj 0 0
Cj-Zj 1 2
2 x2 -1.333 1
II 0 s2 -0.333 0
Zj -2.667 2
Cj-Zj 3.667 0
0 0
s1 s2 b q
1 0 3 1
0 1 3 -3
0 0 0
0 0
0.333 0 1 -0.75
0.333 1 4 -12
0.667 0 2
-0.667 0 NO EXISTE UN <q+