Libro de Investigacion Operativa I
Libro de Investigacion Operativa I
Libro de Investigacion Operativa I
DE
OPERACIONES
NDICE
INTRODUCCIN
.
CAPTULO 1:
1.1 EJEMPLOS
1.2 PROBLEMAS RESUELTOS
1.3 PROBLEMAS DE P.L PREPARADOS CON LINGO.
1.4 ASPECTOS DEL ALGEBRA LINEAL Y ANLISIS CONVEXO
1.4.1 VECTORES
1.4.2 OPERACIONES CON VECTORES.
1.4.3 MATRICES
1.4.4 ECUACIONES LINEALES SIMULTNEAS.
1.4.5 CONJUNTOS CONVEXOS.
7
19
71
92
92
92
94
100
103
CAPITULO 2:
PROGRAMACIN LINEAL: TABLERO SIMPLEX
2.1 MTODO GRFICO.
2.2 MTODO SIMPLEX..
2.3 MTODO DE PENALIZACIN...
2.4 MTODO DE LAS DOS FASES...
105
107
113
115
CAPITULO 3:
DUALIDAD
3.1 DUALIDAD: UN ENFOQUE CONCEPTUAL.
3.2 RELACIONES PRIMAL DUAL.
3.3 HOLGURA COMPLEMENTARIA...
3.4 MTODO DUAL SIMPLEX.
3.5 MTODO PRIMAL DUAL
3.6 PROBLEMAS RESUELTOS..
118
122
125
129
131
136
CAPITULO 4:
ANLISIS DE SENSIBILIDAD
4.1 ANLISIS GRFICO DE SENSIBILIDAD
4.2 CAMBIOS EN LOS COEFICIENTES DE LA FUNCIN OBJETIVO
4.3 CAMBIOS EN LA DISPONIBILIDAD DE RECURSOS..
4.4 PRECIO DUAL
4.5 CAMBIOS EN LA MATRIZ DE COEFICIENTES TECNOLGICOS
4.6 ADICIN DE UNA VARIABLE
4.7 ADICIN DE UNA RESTRICCIN.
4.8 REGLA DEL 100% .
4.9 INTERPRETACIN DEL PROGRAMA LINDO...
4.10 INTERPRETACIN DEL PROGRAMA LINGO.
153
157
163
165
166
168
170
172
177
195
CAPITULO 5:
PROGRAMACIN ENTERA
229
246
290
INTRODUCCIN
En el mundo real, las organizaciones de diferentes naturalezas tienen problemas de
decisin en el uso de sus recursos escasos. Como por ejemplo: Un empresario
dedicado al servicio de mantenimiento y reparacin de computadoras tiene cinco
tcnicos que atienden pedidos de diversas empresas en Lima y provincias, est
interesado en determinar el lugar ms apropiado para su sede central. Recursos
escasos: tiempo no productivo, pasajes, etc. Otro ejemplo: Un empresario
propietario de 5 automviles dedicados al servicio de taxi en la ciudad de Lima est
interesado en determinar el grifo que debe abastecer a sus vehculos. En este caso,
los recursos escasos son: las llantas, el tiempo dedicado para abastecerse de
gasolina, el mismo combustible, etc.
El proceso para alcanzar este objetivo consiste ms en formular el problema que en
construir y resolver modelos matemticos. En forma especfica, los problemas de
decisin a menudo incluyen importantes factores que muchas veces no se pueden
incluir en el modelo matemtico. El factor principal es el hombre y su
comportamiento. El modelo puede ser muy bueno, pero si la influencia de las
personas es muy fuerte, la solucin ptima del modelo es impracticable.
El vector columna del lado derecho representa los requerimientos mnimos que
deben satisfacer.
Las restricciones X 1 , X 2 ,........., X n 0 son las condiciones de no negatividad
de cada variable.
El mtodo simplex est diseado para resolver programas lineales donde las
variables de decisin son no negativas.
A continuacin se presenta una serie de problemas con sus respectivos programas
lineales, el objetivo que se persigue es mostrar la mayor cantidad posible de
mecanismos necesarios para formular cualquier problema lineal.
1.1 EJEMPLOS
CASO: PRODUCCIN
1. Una compaa elabora dos productos P1 y P2 cada uno requiere de
componentes c1 y c 2 la disponibilidad de componentes y precio venta se
muestra en el siguiente cuadro.
Componentes
c1
c2
Producto
P1
P2
Dispone
(Unid.)
1
3
2
1
15000
10000
Precio Venta
(S/. / Unidad)
4
3
Max Z = 4 X 1 + 3 X
Sujeto a:
X1 + 3 X
15,000
7
10,000
X1 , X 2 0
2 X1 + X
X 3 = Unidades producidas de P3
X 31 = Unidades producidas de P3 que se venden.
X 32 = Unidades producidas de P3 que se destruyen.
Max z = 4 X 1 + 3 X 2 + 2 X 31 - 0.5 X 32
Sujeto a:
X 1 + 3 X 2 15,000
2 X1 +
X 2 10,000
X 3 = 3 X 31
X 31
500
X 31 + X 32 = X 3
X 1 , X 2 , X 3 , X 31 , X 32 0
Componente c1
Rg.
De
Componente c 2
S/. /
Unid
Rg.
De
S/. /
Unid
8
1
2
3
1
5,000
5,001 12,000
12,001 15,000
0.3
0.4
0.5
1
2
1
8,000
8,001 10,000
0.2
0.4
X 1 + 3 X 2 15,000
2 X1 +
X 2 10,000
X c11 + X c1 2 + X c1 3 = X 1 + 3 X 2
X c2 1 + X c2 2 = 2 X 1 + X 2
X c11 5,000
X c1 2 7,000
X c1 3 3,000
X c2 1 8,000
X c2 2 2,000
X 1 , X 2 , X c11 ,, X c2 2 0
CASO: METAS DE TRABAJO
P1
P2
Disponibilidad
(Unid.)
Componentes
c1
c2
1
3
2
1
15000
10000
Solucin:
X i = Unidades del producto i (i = 1, 2, 3)
Max z = X 3
Sujeto a:
X 1 + 3 X 2 15,000
2 X 1 + X 2 10,000
X1 X 3
X2 X3
X1 , X 2 , X 3 0
Mes
1
2
3
Costo de
Produccin
100
150
200
Venta
(Unidades)
300
350
400
10
FINAL
MES 1
X 1 - 300 = Y1
X 1 700
MES 2
11
Y1 + X 2 - 350 = Y2
X 2 700
MES 3
Y2 + X 3 - 400 = 0
X3
700
700
X 1 , X 2 , X 3 , Y1 , Y2 , Y3 0
Mes
1
2
3
Costo de
Produccin
300
200
100
Venta Total
Venta
(Unidades)
300
350
400
1050
Si GAMMA no cumple con la venta mensual a BETA tendr que pagar una
multa de S/. 30.00 por unidad mensual faltante.
GAMMA est obligada a cumplir con la entrega de las 1.050 unidades al final
del tercer mes.
Solucin:
X i = Produccin en el mes i (i = 1, 2, 3)
Yi = Excedente o dficit en el mes i (i = 1, 2, 3)
Wi = Costo mensual de almacenamiento o multa en el mes i (i =1, 2, 3)
PRODUCCIN MENSUAL
INV. INICIAL + PRODUCCIN - VENTA = INV. FINAL (DFICIT)
MES 1
X 1 - 300 = Y1
X 1 800
MES 2
Y1 + X 2 - 350 = Y2
X 2 800
MES 3
Y2 + X 3 - 400 = 0
X3
800
X 1 + X 2 + X 3 = 1,050
COSTOS DE ALMACENAMIENTO O MULTA
MES 1
Si: Y1 0:
10 Y1
-30 Y1
W1
W1
13
Si: Y1 < 0:
-30 Y1 W1
10 Y1 W1
MES 2
W1
W1
10 Y2 W2
-30 Y2 W2
800
0
0
10 ( Y21 - Y22 ) - W1 0
-30 ( Y21 - Y22 ) - W2 0
10 ( Y11 - Y12 ) - W1
-30 ( Y11 - Y12 ) - W1
14
Excedente, Faltante)
MES 1
X 1 - 300 = Y1e Y1 f
X 1 800
MES 2
Y1e Y1 f + X 2 - 350 = Y2 e Y2 f
X 2 800
MES 3
Y2 e Y2 f + X 3 400 = 0
X 3 800
X 1 + X 2 + X 3 = 1,050
15
Fabrica
A
B
Demanda
(Unid)
Costo de Transporte
c1
c2
5
10
12
3
250
Produccin
(Unidades)
300
400
350
1,2)
Min z 5 X 11 10 X 12 12 X 21 3 X 22
Sujeto a:
PRODUCCIN
DEMANDA
X 11 X 12 300
X 21 X 22 400
X 11 X 21 250
X 12 X 22 350
16
X 11 , X 12 , X 21 , X 22 0
Si se cambia por en la restriccin de la demanda, entonces cuando se
resuelva el problema el valor de la funcin objetiva es igual a cero; porque no
se transporta nada y eso no es lo queremos.
CASO: PROCESOS DE MEZCLA
9. Un Kg de P es el resultado de mezclar A, B y C cuyas caractersticas son las
siguientes:
Producto
A
B
C
Elemento 1
(%)
20
30
10
Elemento 2
(%)
40
15
30
Precio (S/. /
Kg)
70
40
60
Obtenga la mezcla ptima si se desea que un kg. P tenga al menos 25% y 30%
de los elementos 1 y 2 respectivamente
Solucin
X i Cantidad del producto i (i = A, B, C) a utilizar en un Kg de P.
Min z 70 X A 40 X B 60 X C
sujeto a :
0.2 X A 0.3 X B 0.1X C 0.25 1Kg
0.4 X A 0.15 X B 0.3 X C 0.30 1Kg
X A X B X C 1Kg
X A, X B, XC 0
Cedulas de Servicio
Turno
Horas
Entrada
1
2
3
4
5
6
0
4
8
12
16
20
Salida
8
12
16
20
0
Intervalo de Tiempo
Turno
00 04 04 08 08 12 12 16 16 20 20 00
X1
X1
1
X2
X2
2
X3
X3
3
X4
X4
4
X5
X5
5
X6
X6
6
Requerimiento
5
9
12
10
6
10
Min z X 1 X 2 X 3 X 4 X 5 X 6
18
sujeto a :
X1 X 6 5
X1 X 2 9
X 2 X 3 12
X 3 X 4 10
X4 X6 6
X 5 X 6 10
X1, X 2 , X 3 , X 4 , X 5 , X 6 0
Venta por
Barril
5
2
Composicin
Cebada
Lpulo
5
2
2
1
Solucin
X i Unidades producidas i (i = 1, 2)
Max z 5 X 1 2 X 2
sa :
5 X 1 2 X 2 60
2 X 1 X 2 25
X1, X 2 0
19
Horas de
Trabajo
2
3
60
Materia
Prima
(lb.)
1
2
40
Producto (oz.)
A
2
3
3
2
Solucin
X i Nmero de procesos de tipo i (i =1,2)
Y j Cantidad producida de j (j = A, B)
YBk Cantidad del proceso B con k (k = V, D)
Producto
1
2
Precio de
Venta
(Por Unidad)
450
700
Costo
Variable
Costo Fijo
240
360
150,000
240,000
Solucin
X i Unidades del producto i (i = 1,2)
Viga
Mquina
B
C
21
Pequea
Mediana
Larga
Extra Larga
300
250
200
100
600
400
350
200
800
700
600
300
Supongamos que cada mquina se puede usar hasta 50 horas por semana y que
los costos de operacin por hora de estas mquinas son $ 30, $ 50 y $ 80
respectivamente. Supngase adems, que semanalmente se requiere 10 000, 8
000, 6 000 y 6 000 pies de los distintos tamaos de las vigas I.
Formular el problema de programacin de mquinas como un programa lineal.
Solucin
X ij Cantidad de horas para producir la viga i (i = pequea, mediana, larga
Min Z 30 X 11 X 21 X 31 X 41 50 X 12 X 22 X 32 X 42
80 X 13 X 23 X 33 X 43
22
Costo por
Hora ($)
1
2
3
32
39
46
151
147
155
72
61
57
118
126
121
89
81
84
Taller
(tiempo
disponible,
Hr.)
160
160
160
Solucin
X ij Articulo producido en el taller i (i = 1, 2, 3) y del tipo de artculo j
(j=A, B, C, D)
Min Z 89 X 1 A X 1B X 1C X 1D 81 X 2 A X 2 B X 2C X 2 D
84 X 3 A X 3 B X 3C X 3 D
Sujeto a:
X 1 A , X 2 A , X 3 A , X 1B , X 2 B , X 3 B , X 1C , X 2C , X 3C , X 1D , X 2 D , X 3 D 0
24
Mquina A
12
15
Mquina B
06
12
25
Acero 2
Tiempo
Costo ($)
(min)
11
22
9
18
10
30
Solucin
Xi j = Cantidad de acero tipo j (j = 1, 2) producido en la acera i (i = 1, 2, 3)
Min Z = 10X11 + 12X21 + 14X31 + 11X12 + 9X22 + 10X32
Sujeto a:
20X11 + 22X12 12000
24X21 + 18X22 12000
28X31 + 30X32 12000
X11 + X21 + X31 500
X12 + X22 + X32 600
De los
ngeles
De Chicago
A Houston
A Nueva
York
20000
15000
18000
17000
Solucin
Xij = Cantidad de barriles anuales provenientes de i con destino j
Yij = Cantidad de barriles (x milln) provenientes de la ampliacin en i con
destino j.
Max z = 20000X11 + 15000X12 + 18000X21 + 17000X22 - 120000 (Y11 + Y12) 150000 (Y21 + Y22)
Sujeto a:
X11 + X21 + Y11 + Y21 5
X12 + X22 + Y12 + Y22 5
X11 + X12 2
X21 + X22 3
10. Para realizar una encuesta por telfono, un grupo de investigacin de mercado
necesita comunicar por lo menos a 150 esposas, 120 maridos, 100 varones
adultos solteros y 110 mujeres adultas solteras. Cuestan 2 dlares realizar una
llamada telefnica durante el da, y 5 dlares durante la noche (debido a
mayores costos laborales). Estos resultados se muestran la tabla sgte. Se
pueden realizar a lo ms la mitad de estas llamadas en la noche, por disponer
de un nmero limitado de empleados. Formule un PL que minimice los costos
para completar la encuesta.
Persona que
Contesto
Esposa
% de llamadas % de llamadas
diurnas
nocturnas
30
30
27
Marido
Soltero
Soltera
Nadie
10
10
10
40
30
15
20
05
Solucin
Xi = Cantidad de llamadas realizadas en el da o en la noche i (i = 1, 2)
Min z = 2X1 + 5X2
Sujeto a:
0.30X1 + 0.30X2 150
0.10X1 + 0.30X2 120
0.10X1 + 0.15X2 100
0.10X1 + 0.20X2 110
0.4X1 + 0.05X2 0
2X2 X1
11. CSL es una cadena de tiendas de servicio para computadoras. El nmero de
horas de reparacin especializada que requiere CSL durante los prximos
cinco meses, se dan a continuacin:
Mes 1 (enero) = 6000 horas
Mes 2 (febrero) = 7000 horas
Mes 3 (marzo) = 8000 horas
Mes 4 (abril)
= 9500 horas
Solucin
Roble (p.t.)
Mesas
Sillas
Disponibilidad
17
05
150
Pino (p.t.)
30
13
210
Precio
de
Venta (US$)
40
15
30
Madera
tipo 1
Comprar a
externos
Madera
tipo 2
Madera
Procesada
Producir el
mismo
Tiempos de secado
2 seg.
0.8 seg.
1.3 seg.
Tipo de madera
Tipo1
Tipo2
Tronco
Solucin
X1 = madera tipo 1 costo (3 + 4 dlares/pie3) = 07 dlar/pie3
X2 = madera tipo 2 costo (7 + 4 dlares/pie3) = 11 dlar/pie3
31
Grupo 2 de
Clientes
10
12
15
Valor de A
para(dlares)
Valor de B
para(dlares)
Solucin
0.3
0.5
Taladr
o
--0.1
Fresador
a
--0.3
24,000
7,200
14,400
Producto
Forja
Torno
Esmeril
Barra
Bastidor
Horas
Disponibles
0.5
0.8
0.2
--
48,000
12,000
Max. Z = X3
Sujeto a:
16. Con rubes y zafiros zales Jewelers producen dos tipos de anillos. Un anillo
tipo 1 requiere 2 rubes, 3 zafiros, y 1 h de trabajo de un joyero. Un anillo tipo
2 requiere 3 rubes, 2 zafiros, y 2 h de trabajo de un joyero. Cada anillo tipo 1
se vende a 400 dlares, y cada anillo tipo 2, a 500 dlares. Se pueden vender
todos los anillos producidos por zales. Actualmente zales dispone de 100
rubes, 120 zafiros y 70 horas de trabajo de un joyero. Se puede comprar ms
rubes a un costo de 100 dlares el rub. La demanda del mercado requiere una
34
17. Suponga que la planta en San Luis fabrica al producto 1, que sirve como
componente (insumo) para la fabricacin de un producto final 2, en Monterrey
y otro producto final 3 en Monclova. As mismo el producto 3 requiere como
insumo adicional el producto 2. La siguiente figura muestra el flujo de
fabricacin.
Planta San
Luis
Producto 1
Planta
Monterey
Producto 2
Ventas
Planta
Monclova
Producto 3
Ventas
35
Fabrica
San Luis
Monterrey
Monclova
Capacidad de
Produccin
(miles de unidades )
200
120
100
Insumo
Producto 1
Producto 2
4
2
1
Venta Nacional por mes
Mnima
Mxima
10000
30000
25000
50000
40000
60000
X1 200,000
X2 120,000
X3 100,000
VENTA DE PRODUCTO 1
La venta es el resultado de la diferencia entre la produccin y el requerimiento
de unidades que participan como insumo para la produccin de otros
productos.
X1 = Y11 + Y12
Y11 10,000
Y11 30,000
Y12 = 0.10 * Y11
Y11 = 4X2 + X3 + Y12
Costo por
Onza
0.05
0.02
Protena (%)
Grasa (%)
40
15
15
18
Solucin
X1 = Cantidad de onzas del producto de soya tipo 1.
X2 = Cantidad de onzas del producto de soya tipo 2.
Min Z = 0.05X1 + 0.02X2
Sujeto a:
0.40X1 + 0.15X2 5
0.15X1 + 0.18X2 2
X 1, X 2 0
30
20
40
20
20
60
15
30
40
30
25
20
20
40
30
15
38
El nuevo producto consiste del 20% del elemento A, al menos 30% del
elemento B y al menos 20% del elemento C. Debido a los efectos laterales de
los compuestos 1 y 2, estos no deben de exceder del 30% y 40% del contenido
del nuevo producto.
Formular como programa lineal el problema de encontrar la forma menos
costosa de obtener un nuevo producto.
Solucin
Xi = Cantidad del compuesto qumico i (i = 1, 2, 3, 4)
Min Z = 20X1 + 30X2 + 20X3 + 15X4
Un kilogramo del nuevo producto tiene las siguientes caractersticas:
0.3X1 + 0.20X2 + 0.40X3 + 0.2X4 = 0.2
0.2X1 + 0.60X2 + 0.30X3 + 0.4X4 0.3
0.4X1 + 0.15X2 + 0.25X3 + 0.3X4 0.2
X1 0.3
X2 0.4
X1+ X2 + X3 + X4 = 1
X1, X2, X3, X4 0
3. Una compaa produce dos salsas para carne, la aromtica Diablo y la suave
Barn Rojo. Estas salsas se obtienen mezclando dos ingredientes A y B. Se
permite cierto nivel de flexibilidad en la frmula de estos productos. De hecho
las restricciones son:
La Barn debe contener un mximo del 75% del ingrediente A;
La Diablo debe contener por lo menos 25% de A y por lo menos 50% de B.
Se pueden vender ms de 40 cuartos de A y 30 cuartos de B. La compaa
puede vender la salsa que produzca al precio por cuarto de $ 3.35 La Diablo y
$ 2.85 la Barn Rojo.
A y B cuestan $ 1.60 y $ 2.95 por cuarto respectivamente se desea maximizar
el ingreso neto por venta de las salsas.
Formule el problema como programa lineal.
Solucin
X1 = Produccin en cuartos de salsa Diablo
39
Otras restricciones:
X2 0.75X21
X1 0.25X11
X1 0.50X12
X11 + X21 40
X12 + X22 30
X11 + X12 = X1
X21 + X22 = X2
X1, X11........, X22 0
Fertilizantes
Precio ($/lb.)
40
(lb.)
25
10
5
I
II
III
(lb.)
10
5
10
(lb.)
5
10
5
10
8
7
Solucin
Xi= Cantidad de fertilizantes i (i = 1, 2, 3) dado en fraccin de unidad.
Min Z = 10X1 + 8X2 + 7X3
Sujeto a:
25X1 + 10X2 + 5X3 10
10X1 + 5X2 + 10X3 7
5X1 + 10X2 + 5X3 5
X1 + X2 + X3 = 1
X1, X2, X3 0
5. Un vinatero desea mezclar vinos de 5 aos diferentes i= (1,.., 5) para hacer tres
tipos de vinos mezclados. La oferta disponible (en galones) de vino del ao i
es Si, i = 1,2,....,5. La mezcla 1 se considera especial, por lo que no se
producirn ms de 100 galones. En la figura se dan las restricciones de cada
una de las mezclas. Se pide formular un programa lineal.
Mezcla
1
2
3
C3
Solucin
41
ESPECIFICACIONE
PRECI
42
CA
O
DE
VENTA
S
No menos del 60% de
ESCOCS A
KILT
680
No ms del 20% de C
No ms del 60% de C
No menos del 15% de
570
No ms del 50% de C
450
A.
TARTAN
LITROS
DISPONIBLES
2000
2500
1200
PRECIO DE COSTO
$/LITRO
700
500
400
X11 0.60X1
X13 0.20X1
X11 + X21 + X13 = X1
43
KILT
X23 0.60X2
X51 0.15X2
X21 + X22 + X23 = X2
TARTAN
X33 0.5X3
X31 + X32 + X33 = X3
Disponibilidad de los licores A, B, C.
X11 + X21 + X31 2,000
X12 + X22 + X32 2,500
X13 + X23 + X33 1,200
X1, X11,, X33 0
7. Una compaa petrolera produce dos tipos de gasolina que vende a 18 y 21
centavos de dlar por galn. La refinera puede comprar cuatro diferentes
crudos con los siguientes anlisis y costos:
Crudo
1
2
3
4
A
0.80
0.30
0.70
0.40
B
0.10
0.30
0.10
0.50
C
0.10
0.40
0.20
0.10
D
0.14
0.10
0.15
0.12
La gasolina cuyo precio de venta es 21 centavos de dlar por galn debe tener
cuando menos 60% de A y no ms de 35% de B. La de 18 centavos de dlar
por galn no debe tener ms de 30% de C. En el proceso de mezclado se
pierde, por evaporacin 2% de A y 1% de B y C.
Demustrese como se determinan las cantidades relativas de crudo que se
deben utilizar.
44
Solucin
= Cantidad de crudo i (i = 1, 2, 3) que intervienen en la gasolina j (j = 1, 2)
La funcin objetivo es:
Min. Z = 0.14(X11 + X12)+ 0.10(X21 + X22) + 0.15(X31 + X32) + 0.12(X41+X42)
Para determinar las cantidades de crudo a utilizar se parte de la produccin de un
galn de gasolina de cada tipo.
Como en el proceso de mezclado se pierde por evaporacin parte de los elementos
A, B, C; se registra a continuacin los porcentajes que quedan de cada elemento y
la suma total de estos componentes.
Crudo
1
2
3
4
A
0.784
0.294
0.686
0.392
B
0.099
0.297
0.099
0.495
C
0.099
0.396
0.198
0.099
D
0.982
0.987
0.983
0.986
Por ejemplo: el crudo 1 antes del proceso tiene el 80& del elemento A, en el
proceso de mezclado pierde el 2% de A por evaporacin, entonces queda slo: 0.80
x 0.98 = 0.784% de A.
Finalmente sumando los porcentajes da como resultado 0.982.
0.982X11 + 0.987X21 + 0.983X31 + 0.986X41 = 1
0.982X12 + 0.987X22 + 0.983X32 + 0.986X42 = 1
Caractersticas de la gasolina tipo 2:
0.784X12 + 0.294X22 + 0.686X32 + 0.392X42
0.983X32 + 0.986X42)
0.099X12 + 0.297X22 + 0.099X32 + 0.495X42
0.983X32 + 0.986X42)
45
0.983X31 + 0.986X41)
8. Una fbrica de vidrio produce dos tipos de vidrio para uso industrial que se
hacen a base de Borosilicato de Plomo y, la mayor parte de las veces, a base de
sustitutos. La empresa tiene almacenado Slice, Plomo, Brax y pedecera de
vidrio, y dispone de dos mezcladoras y dos hornos para preparar sus productos,
cada tipo de vidrio se procesa en cualquiera de las mezcladoras y en cualquier
horno. Todo el vidrio plano se lamina en la misma mquina de modo que no es
necesario considerar esta operacin. Los productos y los factores de
produccin estn relacionados como se muestra en las siguientes tablas:
Materia
Brax (A)
Plomo (B)
Silice (C)
Pedecera
(D)
Mquina
Mezcladora
L
Mezcladora
M
Horno X
Horno Y
Composicin (Tn)
Vidrio 1
Vidrio 2
0.1
0.2
0.1
0.2
0.8
0.5
0.0
0.1
Composicin (Tn)
Vidrio 1
Vidrio 2
Abastecimiento
Costo (Ton)
(Ton)
25000
100
35000
300
50000
60
15000
30
0.4
0.2
2000
Costo
Variable
(Ton)
30
0.1
0.2
1000
50
0.2
0.5
0.4
0.2
2000
1800
40
30
Capacidad
(Hor)
Los tipos de vidrio no se pueden sustituir uno con otro, por lo que es necesario
producir cuando menos 100 toneladas de cada tipo para pedidos especiales. Si el
precio de venta del vidrio 1 es de $ 200 la tonelada y el vidrio 2 es de $ 300 la
46
47
Condiciones de Produccin:
X1 100
X2 100
X1, X2,......, X2LM, X2MY 0
INGREDIENTES
Maz
Cal
Harina
25
15
25
15
30
20
48
Maiz
Carne
Cal
Pollo
Harina
18
25 X 11 15 X 21 25 X 31
22
X 11 X 21 X 31
20
15 X 11 30 X 21 20 X 31
23
X 11 X 21 X 31
20
25 X 11 15 X 21 25 X 31
X 11 X 21 X 31
49
20
15 X 11 30 X 21 20 X 31
25
X 11 X 21 X 31
la
Aleacin 1
190
2%
1%
3%
Aleacin 2
200
2.5 %
1.5 %
4%
42000
50000
Solucin
= cantidad de aleacin i, (i = 1,2)
Componentes
Silicio
Nquel
Carbono
Costo ($)
Ton
Tipo de Aleacin
Aleacin 1
Aleacin 2
2
2.5
1
1.5
3
4
190
Especificacin
%
1.8 2.5
0.9 1.2
3.2 3.5
200
X1 + X2 = 1
0.01 X1 + 0.015 X2 0.009
0.03 X1 + 0.04 X2 0.035
0.03 X1 + 0.04 X2 0.032
[X1/(X1 + X2)] (42000) + [X2/(X1 + X2)] (50000) 450000 (*)
Simplificando (*)
3000 X1 - 5000 X2 = 0
X1, X2 0
11. Feedco produce dos tipos de alimentos para ganado. Ambos productos estn
hechos completamente de trigo y de alfalfa. El alimento 1 debe contener por lo
menos 80% de trigo, y el alimento 2 por lo menos 60% de alfalfa. El alimento
1 se vende a 1.50 U$ / lb, y el alimento 2 a 1.30 U$ / lb. Feedco puede comprar
hasta 1000 lb de trigo a 0.50 U$ / lb y hasta 80 lb de alfalfa, a 0.40 U$ / lb. La
demanda de ambos tipos de alimento no tiene lmite. Formule un P.L. para
maximizar las ganancias de Feedco.
Solucin
Insumos
Trigo
Alfalfa
Precio ($ /
lb)
Alimento
1
>= 80%
<= 20%
Alimento
2
>= 40%
<= 60%
1.50
1.30
Compra max
(lb)
1000
800
Precio
($ / lb)
0.50
0.40
Sujeto a:
X11 + X12 1000
X21 + X22 800
0.2X11 - 0.8X21 0
0.4X22 - 0.6X12 0
12. Feedco decidi otorgar a su cliente (supngase que hay solamente un cliente)
un descuento, dependiente de la cantidad comprada. Si el cliente compra ms
de 300 lb del producto 1, se le vender cada libra que rebase las primeras 300
lb, a solo 1,25 dlares. Similarmente, si el cliente compra ms de 300 lb del
producto 2, se le vender cada libra que rebase las primeras 300 lb, a slo 1,00
dlar. Modifique el PL del problema 11 para tomar en cuenta los descuentos
por la cantidad comprada. (Sugerencia: defina variables para el alimento
vendido a cada precio).
Solucin
Si la compra es mayor de 300 lb
Alimento 1
Alimento 2
T 1000
A 800
X1T + X2T = T
X1A + X2A = A
Y11 + Y12 = Y1
Y21 + Y22 = Y2
Y11 300
Y21 300
X1T, X1A,......., Y11, Y22 0
13. Feedco decidi otorgar a su cliente (supngase que hay solo un cliente) un
descuento, dependiente de la cantidad comprada. Si el cliente compra ms de
300 lb del producto 1, se le vender cada libra que rebase las primeras 300 lb, a
solo 1.25 dlares. Similarmente, si el cliente compra ms de 300 lb del
producto 2, se le vender cada libra que rebase las primeras 300 lb, a solo 1.00
dlar. Modifique el PL del problema 11 para tomar en cuenta los descuentos
por la cantidad comprada. (Sugerencia: defina variables para el alimento
vendido a cada precio.
Solucin (Verificar con el Problema anterior)
Sean:
i: alimento1, alimento2
j: trigo, alfalfa
Y21 300
Y12 300
Y22 300
X11 - 0.8X1 0
X22 - 0.6X2 0
X11 + X21 1000
X12 + X22 800
14. Un fabricante de gasolina para aviacin vende dos clases de combustibles: A
y B. El combustible A tiene 25% de gasolina de grado 1, 25% de gasolina de
grado 2 y 50% de grado 3. El combustible B tiene 50% de gasolina de grado 2
y 50% de grado 3.. Hay 500 gln/hr. De grado 1 y 200 gln./hr de los grados 2
y3, disponible para su produccin. Los costos son de 30 ctvs. ($0.30) por gln
de grado 1, $0.60 por gln de grado 2 y $0.50 por gln. de grado 3. La clase A
puede venderse a $0.75 por gln., mientras que la clase B alcanza $0.90/gln.
Qu cantidad puede producirse de cada combustible?
Solucin
La informacin se resume en el siguiente cuadro:
Gasolina
Grado 1
Grado 2
Grado 3
Precio ($ /
gl)
Combustible
A
B
0.25
0.25
0.50
0.50
0.50
0.75
Costo
($ / gl.)
0.30
0.60
0.50
Disponibilidad
(gl. / hr.)
500
200
200
0.90
Sea:
X1 =La cantidad de galones a producirse del combustible A
X2 =La cantidad de galones a producirse del combustible B
La cantidad de gasolina de cada grado a usarse ser:
55
Nmero mnimo
de mozos
56
1
2
3
4
5
6
02 10
06 14
10 18
14 22
18 02
22 06
04
08
10
07
12
04
Solucin
= Nmero de mozos contratados en el turno i, (i= 1,..,6)
Trminos
1
2
3
4
5
6
10
14
18
22
X1
X2
X3
X4
X6
X5
Min Z = X1 + X2 + X3 + X4 + X5 + X6
Sujeto a:
Restricciones de personal de mozos en el turno
X1 + X6 4
X1 + X2 8
X2 + X3 10
X3+X4 7
X4 + X5 12
X5 + X6 4
Restricciones de signo:
X1, X2, X3, X4, X5, X6
2. Una aerolnea desea asignar dos tipos de aviones a tres rutas. Cada avin puede
hacer a lo ms dos vueltas diarias. Adems, se dispone de tres aviones del tipo
57
Tipo de
avin
A
B
Solucin
XAi = Cantidad de vuelos por da en la ruta i (i = 1, 2, 3) de los aviones tipo A.
XBi = Cantidad de vuelos por da en la ruta i (i = 1, 2, 3) de los aviones tipo B.
Min. Z = 3000XA1 + 2500XA2 + 200XA3 + 2400XB1 + 2000XB2 + 1800XB3
Sujeto a:
XA1 + XA2 + XA3 6
XB1 + XB2 + XB3 8
140XA1 + 100XB1 300
140XA2 + 100XB2 700
140XA3 + 100XB3 220
XA1, XA2,., XB3 0
3. El Ghotam City National Bank abre de lunes a viernes, de las 9 a.m. hasta las 5
p.m. De experiencias anteriores, el banco sabe que necesita el nmero de
cajeras, indicado en la tabla A. El banco contrata dos tipos de cajeras. Las
cajeras de tiempo completo trabajan de 9 a 5, los cinco das de la semana, y
tienen 1 hora de descanso para comer. ( El banco determina cuando una
empleada de tiempo completo puede comer, pero cada cajera tiene que comer
entre medioda y la 1 p.m. o entre la 1 y las 2 p.m.) Se les paga 8 dlares
58
Turnos
Tipo de cajera
T1 9 -10
Tiempo
Completo
X1
X1
X1
Y1
Y1
Y2
Y1
Y2
Y3
Tiempo
Parcial
Requerimiento
X4
Y2
Y3
Y4
6
X5
Y3
Y4
Y5
5
T8 4 - 5
X1
X1
X1
Y4
Y5
Y6
6
Y5
Y6
8
Y6
8
59
0.5X1 10
0.6X2 39
60
0.6X2 24
Tiempo Total de Visitas = Tiempo Disponible:
1.5X1 + 0.5X2 80
X1, X2 0
Se pide optimizar el nmero de visitas:
Max Z = 0.1 [150(0.5) X1 + 70(0.6) X2] - 1.5 [1,5X1 + 0,5X2]
2. Alden Enterprises produce dos productos. Se puede fabricar cada producto en
cualquiera de dos mquinas. En la tabla A, se dan los tiempos necesarios (en
horas) para producir cada producto en cada mquina. Cada mes los clientes
estn dispuestos a comprar los productos hasta las cantidades y a los precios
indicados en la tabla B. La meta de la compaa es maximizar los ingresos
obtenidos mediante la venta de los productos durante los prximos dos meses.
Formule un PL para ayudar alcanzar esta meta.
TABLA A
MAQUINA
1
MAQUINA
2
Producto 1
Producto 2
TABLA B
DEMANDAS PRECIO(Dlares)
Mes 1 Mes 2 Mes 1
Mes 2
Producto 1
100
190
55
12
Producto 2
140
130
65
32
61
Solucin
Sea:
Xijk = Cantidad de producto i, fabricado en maquina j, en el mes k (i, j, k =1, 2)
Mes 1:
4X111
3X121
+ 7X211 500
+ 4X221 500
4X112
3X122
+ 7X212 500
+ 4X222 500
Mes 2:
Sea:
Cik = La cantidad de producto i, vendida en el mes k (i, k = 1, 2)
Nik = La cantidad de producto i, que no se vende en el mes k (i, k = 1, 2)
Est sujeta a las siguientes restricciones:
X111 + X121 = C1 +N11
C11 100
X211 + X22 = C21 + N21
C21 140
X112 + X122 + N11 = C12 + N12
C12 190
X212 + X222 + N21 = C22 + N22
C22 130
Luego la funcin objetivo viene a ser la maximizacin de la venta de los
productos durante los prximos dos meses.
Max z = 55C11 + 65C21 + 12C12 + 32C22
CASO: PROBLEMAS DE DISTRIBUCIN DE TIERRAS
62
1. Una empresa agrcola explota una finca de 200 Ha., de regado, que puede
dedicarse en principio a dos cultivos C 1 y C2. Los ingresos y costos variables
por hectreas para cada cultivo figuran en la siguiente tabla:
Costos
Variables (S/.
El cultivo C1
puede repetirse
/ Ha)
Cultivo C1
14.000
6.000
Cultivo C2
15.000
6.000
indefinidamente todos los aos en la misma parcela; en cambio el cultivo C 2 ha de
implementarse en parcelas que el ao anterior llevaron otro cultivo; pues sino se
sigue esta norma tcnica (rotacin de cosechas), disminuirn apreciablemente los
rendimientos.
El agua para riego es de 1 lt/seg. y por hectrea, es decir 610,000 m2 al mes para
toda la finca.
Las necesidades de agua de los cultivos en el mes prximo:
Cultivo
Ingresos
(S/. / Ha)
X 1, X 2 0
2. La Canadian Parks Comission vigila dos terrenos. El terreno 1 est formado
de 300 acres y el terreno 2 por 100 acres. Se puede utilizar cada acre del
terreno 1 para abetos, la caza o para ambas cosas. Se puede utilizar cada acre
del terreno 2 para abetos, para acampar o para ambas cosas. En la tabla, se da
el capital, (en cientos de dlares) Y la mano de obra (das hombre) que se
necesitan para mantener un acre de cada terreno, y la ganancia (miles de
dlares) por acre, para cada uso posible del suelo. Se dispone un capital de
150000 y 200 das-hombre de trabajo. Cmo se tiene que asignar el suelo a
los usos diferentes, para maximizar la ganancia recibida de los dos terrenos?
Terreno 1
Abetos
Terreno 1
Caza
Terreno 1
Ambas cosas
Terreno 2
Abetos
Terreno 2
Acampar
Terreno 2
Ambas cosas
Capital
Mano de
obra
Ganancia
0.1
0.2
0.2
0.4
0.2
0.5
0.05
0.06
30
0.09
10
1.01
1.1
Solucin
Xij = # de acres del terreno i (1 ,2) para la actividad j (1, 2, 3)
Max Z = 0.2X11 + 0.4X12 + 0.5X13 + 0.06X21 + 0.09X22 + 1.1X23
Sujeto a:
X11 + X12 + X13 = 300
X21 + X22 + X23 = 100
300X11 + 300X12 + 400X13 + 100X21 + 3000X22 + 1000X23 150000
0.1X11 + 0.2X12 + 0.2X13 + 0.05X21 + 5X22 + 1.01X23 200
64
4,000
6,000
2,000
10,000
8,000
COSTOS DE DISTRIBUCION
De la tienda
1
2
1
0.06
0.15
Al vendedor mayorista
2
3
4
0.04
0.12
0.12
0.09
0.05
0.08
5
0.05
0.08
Solucin
Xij = Cantidad de paquetes entregados por la tienda i al vendedor j (i = 1,2) (j
= 1,2,3,4,5 ).
Se debe minimizar el costo del paquete y distribucin de las tiendas a los
vendedores.
Min Z = 5.31X1 1 + 5.29X12 + 5.37X13 + 5.37X14 + 5.3X15 + 5.85X21 + 5.79X22 +
5.75X23 + 5.78X24 + 5.78X25
Sujeto a:
LONGITUD DEL
PEDIDO
(m)
18,000
9,000
9,000
ANCHO
(cm)
64
60
35
68
Se aclara que los cortes deben efectuarse en sentido longitudinal y que los
mismos no necesitan estar formados por una sola tira.
Solucin
Se debe establecer los posibles patrones de corte, o sea las distintas maneras
que se ha de cortar la bobina a fin de satisfacer los pedidos.
= Longitud de la tira en metros del patrn i
ANCHO (cm)
LONGITUD (m)
X1
X2
X3
X4
X5
X6
X7
X8
X9
64
60
35
Ancho Del
Recorte
3
2
2
1
1
1
-
1
2
1
2
3
-
2
2
4
2
1
6
23
27
17
31
21
11
25
-5
70
Una compaa elabora dos productos P1 y P2, cada uno requiere e componentes
C1 y C2, la disponibilidad de los componentes y precio de venta de muestra en el
siguiente cuadro:
Producto
P1
P2
Dispone
Componentes
C1
1
3
15000
C2
2
1
10000
Precio de Venta
(S/./Unidad)
4
3
71
Costo de Produccin
(S/. / unidades)
100
150
200
Venta (Unidades)
300
350
400
INGREDIENTE
CANTIDAD DISPONIBLE
(kg)
COSTOS
(S/./ kg)
4.000
1.300
B
C
6.000
2.000
1.500
1.000
2.0 S/./ kg
3.0 S/./ kg
1.5 S/./ kg
Con todos los datos que FLORANID S.A. nos ha facilitado, nos piden que
determinemos: Cunta cantidad de cada tipo de abono hay que producir de forma
que se maximice el beneficio de la compaa?
Solucin:
Con los datos podemos construir un primer esquema que nos permitir desarrollar
el modelo de programacin lineal para la resolucin del problema:
INGREDIENTES
A
B
C
1
X11
X21
X31
ABONOS
2
X12
X22
X32
3
X13
X23
X33
CANTIDAD
DISPONIBLE (kg)
4000
6000
2000
COSTOS
(S/. /kg)
1300
1500
1000
CACAHUATE
80%
50%
NUEZ
20%
50%
Un vigilante debe ser entrenado durante cinco das antes de estar disponible para
asignarlo a los clientes.
73
Existe 65 das de trabajo en cada trimestre y al inicio del primer trimestre existen
120 vigilantes calificados en la nmina. Los vigilantes son pagados por la empresa
y no por el cliente; ellos ganan un salario de S/.500 al mes. Durante cada trimestre
la empresa pierde el 15) de su personal (incluyendo vigilantes entrenados en el
trimestre anterior). Formular la PL.
B1
B2
B3
Azcar
Chocolate
Ganancia
unitaria
1
2
3
1
3
7
1
1
5
Cantidad
disponible
50
100
Solucin:
!D=CANTIDAD DISPONIBLE G=GANANCIA UNITARIA;
!IN=MATERIA PRIMA B=TIPO DE BARRA DE CHOCOLATE;
!P=REQUERIMIENTOS PARA CADA PRODUCTO ;
sets:
in/1..2/:d;
b/1..3/:p,g;
ca(in,b):uso;
endsets
data:
g=3 7 5;
d=50 100;
uso= 1 1 1
2 3 1;
enddata
max=@sum(b:p*g);
@for(in(i):@sum(b(j):uso(i,j)*p(j))<=d(i));
end
74
MAX
3 P( 1) + 7 P( 2) + 5 P( 3)
SUBJECT TO
2] P( 1) + P( 2) + P( 3) <=
50
3] 2 P( 1) + 3 P( 2) + P( 3) <=
100
END
7. Las fabricas F1 y F2 tienen una capacidad de produccin de 30 y 20 unidades
respectivamente, se tiene adems 3 centros de demanda C1, C2 y C3, con
capacidades de 10, 25 y 15 unidades respectivamente; finalmente el costo
unitario de transporte de las Fabricas a los Centros es como sigue:
Fabrica / Centro
F1
F2
C1
2
7
C2
4
10
C3
6
1
sets:
f/1..2/:cp;
!Fabricas con su respectivo costo de
produccin
c/1..3/:d;
!Centros de demanda y su respectivo
valor de demanda
rutas(f,c):ct,x;
endsets
data:
cp=30,20;
d=10,25,15;
ct=2 4 6,
7 10 1;
enddata
75
min = @sum(rutas:ct*x);
Funcin Objetivo
@for(c(j):@sum(f(i):x(i,j))>=d(j));
@for(f(i):@sum(c(j):x(i,j))<=cp(i));
end
Formulacin:
MIN
2 X(
1) + 10 X( 2,
+ X( 2,
SUBJECT TO
2] X( 1, 1)
3] X( 1, 2)
4] X( 1, 3)
5] X( 1, 1)
6] X( 2, 1)
END
1, 1) + 4 X( 1, 2) + 6 X( 1, 3) + 7 X( 2,
2)
3)
+
+
+
+
+
X(
X(
X(
X(
X(
2,
2,
2,
1,
2,
1)
2)
3)
2)
2)
>=
>=
>=
+ X(
+ X(
10
25
15
1, 3) <=
2, 3) <=
30
20
8. Steelco produce dos tipos de acero en tres diferentes aceras. Durante un mes
dado, cada acera dispone de 200 horas de alto horno. El tiempo y el costo de
produccin de una tonelada de acero, difiere de una fbrica a otra, debido a las
diferencia en los hornos de cada fbrica. En el cuadro siguiente se muestran el
tiempo y el costo de produccin para cada fbrica. Cada mes, Steelco tiene que
producir por lo menos 500 toneladas de acero 1 y 600 toneladas de acero2.
formule un PL, para minimizar los costos para producir el acero deseado.
ACERIA
ACERIA 1
ACERIA 2
ACERIA 3
COSTO
10
12
14
ACERO 1
TIEMPO
20
24
28
COSTO
11
9
10
ACERO 2
TIEMPO
22
18
30
Solucion:
!X=CANTIDAD PRODUCIDA;
SETS:
aceria/1..3/:horas;
76
acero/1..2/:cantidad;
rutas(aceria,acero):costo,tiempo,x;
ENDSETS
DATA:
horas=12000,12000,12000;
cantidad=500,600;
costo=10 11,
12 9,
14 10;
tiempo=20 22,
24 18,
28 30;
ENDDATA
min=@sum(rutas:costo*x);
@for(acero(j):@sum(aceria(i):x(i,j))>=cantidad(j));
@for(aceria(i):@sum(acero(j):tiempo(i,j)*x(i,j))<=horas(i));
END
Protenas
(Unid /Oz)
20
30
40
40
45
30
Carbohidratos
(Unid /Oz)
50
30
20
25
50
20
Grasa
(Unid / Oz)
4
9
11
10
9
10
Costo
(Oz)
2
3
5
6
8
8
Solucin:
!TIPO=TIPO DE ALIMENTO
CO=COSTO DEL ALIMENTO POR ONZA;
!MACRO=TIPO DE MACRONUTRIENTE UNID=UNIDADES QUE DEBE RECIBIR
EL HAMSTER;
77
X=CANTIDAD A MEZCLAR;
SETS:
TIPO/1..6/:CO,X;
MACRO/1..3/:UNID;
MATRIZ1(TIPO,MACRO):CANT;
ENDSETS
DATA:
CO=2,3,5,6,8,8;
UNID=70,100,20;
CANT=20,50,4,
30,30,9,
40,20,11
40,25,10,
45,50,9,
30,20,10;
ENDDATA
MIN=@SUM(TIPO:CO*X);
@FOR(MACRO(J):@SUM(TIPO(I):CANT(I,J)*X(I))>=UNID(J));
END
10. Las capacidades de produccin del producto P de las fbricas A y B, los costos
por unidad transportada a los centros de consumo C 1 y C2 y las demandas de
estos son como sigue:
Fabrica
A
B
Demanda
Costos de Transporte
C1
C2
5
10
12
3
250
350
Produccin
300
400
Solucin:
!FAB=FABRICA(1,2) PROD=PRODUCCION DE CADA FABRICA;
!CEN=CENTRO (1,2) DEM=DEMANDA;
!CO= COSTO DE TRANSPORTE X=UNIDADES TRANSPORTADAS;
SETS:
78
FAB/1..2/:PROD;
CEN/1..2/:DEM;
MATRIZ1(FAB,CEN):CO,X;
ENDSETS
DATA:
PROD=300,400;
DEM=250,350;
CO= 5,10,
12,3;
ENDDATA
MIN=@SUM(MATRIZ1:CO*X);
@FOR(FAB(I):@SUM(CEN(J):X(I,J))<=PROD(I));
@FOR(CEN(J):@SUM(FAB(I):X(I,J))>=DEM(J));
END
Capacidad
(hora)
500
380
!MAQ=MAQUINA(1,2)
CAP=CAPACIDAD
DE
MAQUINA(HORAS);
!PROD=PRODUCTO(1,2,3,4) PV=PRECIO DE VENTA;
!TM= TIEMPO X=UNIDADES PRODUCIDAS UTI= UTILIDAD;
CADA
SETS:
MAQ/1..2/:CAP;
PROD/1..4/:PV;
MATRIZ1(MAQ,PROD):TM,X,uti;
ENDSETS
DATA:
79
CAP=500,380;
PV=65,70,55,45;
TM= 2,3,4,2,
3,2,1,2;
UTI=45,40,15,25
50,60,50,35;
ENDDATA
MAX=@SUM(MATRIZ1:UTI*X);
@FOR(MAQ(I):@SUM(PROD(J):TM(I,J)*X(I,J))<=CAP(I));
END
12. Para una jornada de 24 horas un hospital esta requiriendo el siguiente personal
para el rea de enfermera, se define 6 turnos de 4 horas cada uno.
Turno
2:00 - 6:00
6:00 - 10:00
10:00 - 14:00
14:00 - 18:00
18:00 - 20:00
20:00 - 24:00
Nmero mnimo
de personal
4
8
10
7
12
4
Horas
2:00 6:00
X1
20:00 24:00
80
X3
X3
X4
X4
X5
X6
Personal
10
12
X5
X6
4
MIN Z = X1 + X2 + X3 + X4 + X4 + X5 + X6
Sujeto a:
Turno 1:
Turno 2:
Turno 3:
Turno 4:
Turno 5:
Turno 6:
X1 + X6 >= 4
X1 + X2 >=8
X2 + X3 >=10
X3 + X4 >=7
X4 + X5 >=12
X5 + X6 >=4
EN
EL
SETS:
HORAS/1..6/:X;
PERS/1..6/:MIN;
MATRIZ1(HORAS,PERS):CANT;
ENDSETS
DATA:
MIN=4,8,10,7,12,4;
CANT=1,1,0,0,0,0,
0,1,1,0,0,0,
0,0,1,1,0,0,
0,0,0,1,1,0,
0,0,0,0,1,1,
1,0,0,0,0,1;
ENDDATA
MIN=@SUM(HORAS:X);
81
@FOR(PERS(J):@SUM(HORAS(I):CANT(I,J)*X(I))>=MIN(J));
END
MIN
X( 1) +
SUBJECT TO
2] X( 1) + X(
3] X( 1) + X(
4] X( 2) + X(
5] X( 3) + X(
6] X( 4) + X(
7] X( 5) + X(
END
X( 2) + X( 3) + X( 4) + X( 5) + X( 6)
6)
2)
3)
4)
5)
6)
>=
>=
>=
>=
>=
>=
4
8
10
7
12
4
Tipo
de Interes Anual Factor
Inversion
(%)
Riesgo
1
2
3
4
5
6
8.5
9
8.5
14.3
6.7
13
0.02
0.01
0.38
0.45
0.07
0.35
de
Plazo
promedio de
inversion
8
2
5
6
2
4
!DAT= CARACTERISTICAS;
SETS:
TIPO/1..4/:INV;
CAR/1..6/:INT,X;
MATRIZ1(TIPO,CAR):DAT;
ENDSETS
DATA:
INT=8.5,9,8.5,14.3,6.7,13;
INV= 2000,400,10000,500;
DAT= 1,1,1,1,1,1,
0.02,0.01,0.38,0.45,0.07,0.35,
8,2,5,6,2,4,
0,0,0,1,0,1;
ENDDATA
MIN=@SUM(CAR:0.01*INT*X);
@FOR(TIPO(I):@SUM(CAR(J):DAT(I,J)*X(J))>=INV(I));
END
14. Salvaje Oeste produce dos clases de sombrero vaquero. Un sombrero de la
clase 1 requiere el doble de mano de obra que uno de la clase 2. Si toda la
mano de obra se dedicara solo a la clase 2, la empresa podra producir
diariamente 400 de estos sombreros. Los lmites de mercado respectivos son
150 y 200 sombreros diarios para esas clases. La utilidad es $8 por cada
sombrero de la clase 1, y $5 por cada uno de la clase 2.
Solucin:
!Rhs=Recursos U=Utilidad de cada sombrero;
!Aij=Coeficientes de las variables
X= Cantidad de
sombreros a producir;
SETS:
VARI/1..3/:Rhs;
VARJ/1..2/:U,X;
ConsVar(VARI,VARJ):Aij;
ENDSETS
83
DATA:
Rhs=400,150,200;
U=8,5;
Aij= 2,1,
1,0,
0,1;
ENDDATA
MAX=@SUM(VARJ:U*X);
@FOR(VARI(I):@SUM(VARJ(J):Aij(I,J)*X(J))<=Rhs(I));
END
15. Blubber Maid, Inc. Fabrica tres productos de caucho: Airtex (material
esponjoso), Extendex (material elstico) y Resistex (material rgido). Los tres
productos requieren los mismos tres polmeros qumicos y una base. La
cantidad de cada ingrediente usado por libra del producto final se muestra en la
siguiente tabla.
Producto
Airtex
Extendex
Resistex
Inventario
Base
6
9
2
1100
Cada producto tiene una utilidad de 7, 7 y 6 S/. ; mientras que la demanda de cada
uno es de 1000, 500 y 400 unidades respectivamente.
Solucin:
!PROD=PRODUCTO GAN=UTILIDAD DEM=DEMANDA;
!X= CANTIDAD A PRODUCIR INV=INVENTARIO CANT=CANTIDAD DE
INGREDIENTES;
SETS:
PROD/1..3/:GAN,DEM,X;
ING/1..4/:INV;
VECTOR(PROD,ING):CANT;
ENDSETS
DATA:
GAN=7,7,6;
84
DEM=1000,500,400;
INV=500,425,650,1100;
CANT=4,2,4,6,
3,2,2,9,
6,3,5,2;
ENDDATA
MAX=@SUM(PROD:GAN*X);
@FOR(ING(J):@SUM(PROD(I):CANT(I,J)*X(I))<=16*INV(J));
@FOR(PROD(I):X>=DEM(I));
END
16. Walnut Orchard tiene dos granjas que cultivan trigo y maz. Debido a las
diferentes condiciones el suelo, existen diferencias en la produccin y en los
costos e produccin de las dos granjas. En la tabla se encuentran los costos y la
produccin para las dos granjas. Cada granja dispone de 100 acres para los
cultivos. Hay que producir 11000 busheles de trigo y 7000 busheles de maz.
Determinar un plan de siembra que minimice los costos para satisfacer estas
demandas.
Granja
Granja 1
Granja 2
Costo
100
120
Maz
Produccin
500
650
Costo
90
80
Trigo
Produccin
400
350
Solucin:
!TAM=TAMAO DE ACRES DISPONIBLES PARA CADA GRANJA;
!UNID=PRODUCCION DE CADA GRANJA X=CANTIDAD DE LA
GRANJA(1,2);
SETS:
GRANJA/1..2/:TAM;
PRODUCTO/1..2/:DEMANDA;
RUTAS(GRANJA,PRODUCTO):UNID,COSTO,X;
ENDSETS
DATA:
TAM=100,100;
DEMANDA=11000,7000;
UNID=500,650,
85
400,350;
COSTO=100,120,
90,80;
ENDDATA
MIN=@SUM(RUTAS:COSTO*X);
@FOR(GRANJA(I):@SUM(PRODUCTO(J):X(I,J))<=TAM(I));
@FOR(PRODUCTO(J):@SUM(GRANJA(I):UNID(I,J)*X(I,J))>=DEMANDA(J)
);
Planta
1
2
3
4
Disponibilidad
Solucin:
86
18. Una empresa que fabrica un producto nico, tiene 3 fbricas y 4 clientes. Las 3
fabricas producen 3 000, 5 000 y 5 000 unidades respectivamente, durante el
siguiente periodo. La empresa se comprometi a vender 4 000 unidades al
87
1
(dlares)
65
68
63
AL CLIENTE
2
3
4
(dlares) (dlares) (dlares)
63
62
64
67
65
62
60
59
60
SETS:
FABR/FAB1,FAB2,FAB3/: CAPAC;
CLIEN/CEN1,CEN2,CEN3,CEN4/:DEMAN;
VIAS(FABR,CLIEN):GANAN,UNID;
ENDSETS
DATA:
CAPAC=3000 5000 5000;
DEMAN=4000 3000 3000 6000;
GANAN=65,63,62,64,
68,67,65,62,
63,60,59,60;
ENDDATA
MAUNID=@SUM(VIAS:GANAN*UNID);
@FOR(CLIEN(J):@SUM(FABR(I):UNID(I,J))>=DEMAN(J));
@FOR(FABR(I):@SUM(CLIEN(J):UNID(I,J))<=CAPAC(I));
END
88
19. OILCO tiene campos petrolferos es San Diego y en Los ngeles. El campo de
San Diego puede producir diariamente hasta 500 000 barriles por da.
Se manda el petrleo desde los campos hacia una refinera en Dallas o en
Houston (suponga que cada refinera tiene capacidad ilimitada).
Cuesta 700 dlares refinar 100 000 barriles de petrleo en Dallas y 900 dlares
100 000 barriles en Houston. Se enva petrleo refinado hacia clientes en
Chicago y en New York. Los clientes en Chicago necesitan diariamente 400 000
barriles de petrleo refinado y los clientes de Nueva York slo 300 000 barriles
de petrleo refinado. En la tabla se muestran los costos de envo de 100 000
barriles de petrleo (refinado o no) entre las ciudades. Formule un modelo de
transporte balanceado para esta situacin.
DESDE
L.A.
San Diego
Dallas
Houston
Dallas
300
420
-
Houston
110
100
-
HACIA
New York
450
470
Chicago
550
530
Solucin:
!S=COSTO DE ENVIAR DE LOS CAMPOS A LAS REFINERIAS;
!C=COSTO DE ENVIAR DE LAS REFIENRIAS A LOS CLIENTES;
!Y=CANTIDAD ENVIADA A LAS REFINERIAS X=CANTIDAD ENVIADA A LOS
CLIENTES;
SETS:
CAMPOS/C1,C2/:PRODUCCION;
CLIENTES/Z1,Z2,Z3/:DEMANDA;
!se ha creado un cliente ficticio al cual llamamos z3;
RUTA1(CAMPOS,CAMPOS):S,Y;
RUTAS(CAMPOS,CLIENTES):C,X;
ENDSETS
MIN=@SUM(RUTAS:C*X)+@SUM(RUTA1:S*Y);
@FOR(CAMPOS(I):@SUM(CLIENTES(J):X(I,J))>=PRODUCCION(I))
;
@FOR(CLIENTES(J):@SUM(CAMPOS(I):X(I,J))>=DEMANDA(J));
89
Z2)
Z2)
Z1)
Z2)
Z3)
+ X(
+ X(
>=
>=
>=
500
400
de 1
10
20
4800
30
6
8
3600
50
16
12
6000
30
12
10
6000
50
90
Producto 1
Producto 2
4
7
Maquina 2
3
4
Tabla 2
Demanda mes Demanda mes Precio mes 1
1
2
100
90
55
140
70
65
Precio mes 2
42
62
UNIDADES
840
760
670
1030
n-1
Ejemplos:
a.- (1, 3, -1, 5) es un vector fila de dimensin
4
b.- es un vector columna de dimensin
8
n = 4.
n=2
Suma De Vectores
Espacio Euclidiano
92
Combinacin Lineal
b=
j 1
reales.
R a
j 1
Ejemplo:
a1 = (3, 5) y a2 = (1, 7), estos vectores son linealmente Independientes puesto que:
R1 (3, 5) + R2 (1, 7) = (0, 0)
(3R1 + R2, 5R1 + 7R2) = (0, 0)
3R1 + R2 = 0
5R1 + 7R2 = 0
La solucin es R1 = R2 = 0
Si para alguna R j = R1, R2,. . ., R k donde no todos son ceros se dice que los
vectores son linealmente dependientes.
Ejemplo:
a1 = (3, 5) y a2 = (6, 10)
R1 (3,5) + R2 (6, 10) = (0, 0)
93
1.4.3 MATRICES
Una matriz es un arreglo rectangular de nmeros denotados por A =[a i j] mxn donde
m = # de filas y n = # de columnas.
Matriz Cero
0
A
0
0
0
0
0
Suma de Matrices
94
Si A = [ a
ij
]yB[b
ij
Ejemplo:
3
A
4
1
3
5
B
6
8
8
8 9
C A B
10 11
Sea A = [ a
ij
4 3
K 2
6
2A
8
2
6
95
Multiplicacin de Matrices
d i j = a i k * b k j Para: i = 1, 2,. . ., m y
j = 1, 2,. . ., n
Ejemplo:
1
1
A 4
2
2
0
5
1
5
B 3
1
0
0
1
1x0 1x 0 1x1
4 x0 2 x 0 5 x1
2 x0 0 x0 1x1
1
5
1
Matriz Transpuesta
i j
de A, es igual al
Ejemplo:
96
2
A 3
4
5
6
2
AT
5
3
6
4
7
Matriz Identidad
Sea A = [a i j ] una matriz nxn, se dice que es una matriz identidad, denotada por I,
si todos los elementos de la diagonal son iguales a uno y todos los dems
elementos son iguales a cero.
Ejemplo:
1
I 0
0
0
1
0
0
1
Inversin de Matrices
Sea A = [a i j] una matriz cuadrada nxn. Si B = [b i j] es una matriz nxn tal que AB =
I y BA = I, entonces B se llama inversa de A. La matriz inversa, si existe, es nica
y se denota por A-1.
Si A tiene una inversa, entonces A se llama no singular; en caso contrario se llama
singular. Una matriz dada A = [a i j] nxn tiene inversa, si y solo si, las filas de A son
linealmente independientes o, de manera equivalente, si las columnas de A son
linealmente independientes.
97
3 1 X 1 9
5 2 X 4
2
La solucin de X y la inversa de la matriz base pueden obtenerse directamente
considerando:
(A | I)(x) = b y omitiendo (x) (A | I | b)
Multiplicando por A-1 (A-1) (A | I | b)
Obteniendo finalmente (I | A-1 | A-1b)
Por consiguiente, aplicando una operacin de transformacin de filas, se obtiene
las siguientes iteraciones:
Para el sistema A x = b, le damos la forma (A | I | b)
98
1/ 3
1/ 3
11 / 3
5/3
11
0 15 / 33
2 / 11
1 / 11
3 / 11
1 / 11
3 / 11
99
2
3
1
1
2
13
15
2
0
1 8
1 13
0 6
100
Rango de (A, b) = 3 y
AX = b no tiene solucin
Caso 2:
X1 + X2 = 8
2X1 + X2 = 13
1
1 8
1 13
0 5
1 3
Por consiguiente X1 = 5 Y X2 = 3
Caso 3:
X1 + X2 + X3 = 8
2X1 + X2
= 13
1
1 1 8
1 0 13
0 1 5
1 2 3
101
n!
m! n m !
3!
2! 3 2 !
102
CONJUNTO CONVEXO
CONJUNTO NO CONVEXO
Puntos Extremos
Sea un conjunto X en E n , se dice que los puntos extremos son aquellos que no
pueden ser representados como una combinacin convexa estricta de dos puntos
distintos en X.
Ejemplo:
Hiperplano
103
Conjunto Polidrico
104
X1, X2 0
Si por el momento se considera a estas desigualdades como igualdades, se obtienen
puntos que luego los llevaremos a una grfica, que se muestra en la siguiente
pgina.
2X1 + X2 = 6
Si X1 = 0 X2 = 6
Si X2 = 0 X1 = 3
X1 - X2 = 1
X1
0
3
X2
6
0
X1
0
1
X2
-1
0
105
Z = CX
Sujeto a:
AX = b
X0
CB: Coeficientes de las variables bsicas
CN: Coeficientes de las variables no bsicas
XB: Variables bsicas
106
XB
XN
Z (C B C N )
Sujeto a:
XB
b
( B N )
XN
Luego:
Z = CB XB + CN XN
Sujeto a:
BXB + NXN = b
A continuacin multiplicamos por B-1 tanto en la parte derecha como en la parte
izquierda de la restriccin:
B-1(BXB + NXN) = B-1b
Resolviendo el siguiente producto:
B-1BXB + B-1NXN = B-1b
XB = B-1b - B-1NXN
Sabemos que: XN = 0
Reemplazamos 2
en 1
Z = CB (B-1b - B-1NXN) + CN XN
Z = CB B-1b (CBB-1NXN - CN XN)
107
. 3
:
Z = CB B-1b (CBB-1N - CN) XN
1Z + 0 XB + (CBB-1N - CN) XN = CB B-1b 4
De 2 :
XB = B-1b - B-1NXN
0 Z + 1XB + (B-1N)XN = B-1b
De
4 y
. 5
5 :
1 0 C B B 1 N C N
0 1
B 1 N
Z
XB
X
N
C B B 1b
B
b
XB
XN
LD
CBB-1N - CN
CB B-1b
XB
B-1N
B-1b
108
VARIABLES DE HOLGURA
Es una variable positiva que representa la diferencia entre los dos lados de una
restriccin.
VARIABLES ARTIFICIALES
Despus de introducir las variables de holgura y observar que no existe una
submatriz identidad para tener una solucin bsica factible inicial, entonces se
introducir variables denominadas como Variables Artificiales para obtener la
submatriz identidad.
Se va ilustrar con un ejemplo los pasos a dar para la resolucin de un problema.
Max Z = 3X1 + 10X2
Sujeto a:
2X1 + 3X2 8
8X1 + 3X2 20
X 1, X2 0
Z
X3
Z
1
0
X1
-3
2
X2
-10
3
X3
0
1
X4
0
0
LD
0
8
109
X4
20
En el ejemplo, la variable X 2 es la que tiene el coeficiente ms negativo (10), por lo tanto se convertira en la variable de entrada.
c) Identificacin de la variable de salida de la base: Se denomina variable de
salida a aquella variable, cuyo valor se aproxime ms rpidamente o cero a
medida que el valor de la variable de entrada vaya creciendo, esto se hace
mediante el siguiente procedimiento algebraico.
br
X K min
, YrK 0
YrK
En el ejemplo:
K
por la variable X .
Z
X2
X4
Z
1
0
0
X1
-3
2
8
X2
-10
3
3
X3
0
1
0
X4
0
0
1
LD
0
8
20
110
Se ha de convertir en
X2
0
1
0
8)
x 1/3
2/3
1
-3
-10
11/3
0
1/3
0
10/3
0
0
0
8/3 )
0)
80/3
x 10
2/3
1/3
8/3 )
8
6
3
0
0
-1
1
1
20 )
12
x (3)
111
Z
X2
X4
Z
1
0
0
X1
11/3
2/3
6
X2
0
1
0
X3
10/3
1/3
-1
X4
0
0
1
LD
80/3
8/3
12
minimizacin.
Se procede a solucionar el problema.
Ejemplo:
Min Z = 3X1 + 8X2
Sujeto a:
X1 + X2 = 200
X1 800
X2 60
112
Z
X3
X4
X6
Z
1
0
0
0
X1
-3
1
1
0
X2
-8
1
0
1
X3
-M
1
0
0
X4
0
0
1
0
X5
0
0
0
-1
X6
-M
0
0
1
LD
0
200
80
60
M-3
X3
X4
X6
0
0
0
1
1
0
2M8
1
0
1
-M
260M
1
0
0
0
1
0
0
0
-1
0
0
1
200
80
60
M-3
M-8
X3
X4
X2
0
0
0
1
1
0
0
0
1
1
0
0
0
1
0
1
0
-1
3-M
M-8
X3
X1
0
0
0
1
0
0
1
0
-1
1
1
0
82M
-1
0
1
82M
-1
0
140M+480
140
80
60
60M+720
60
80
113
X2
-1
60
Z
X5
X1
X2
1
0
0
0
0
0
1
0
0
0
0
1
8-M
1
0
1
-5
-1
1
-1
0
1
0
0
-M
-1
0
0
1200
60
80
120
X2 = 120,
ZMIN = 1200
Ejemplo:
Con el ejemplo utilizado en el Mtodo de Penalizacin.
FASE I: Se tiene que la funcin objetivo para la primera fase es:
Min = X3 + X6
114
Z
X3
X4
X6
Z
1
0
0
0
X1
0
1
1
0
X2
0
1
0
1
X3
-1
1
0
0
X4
0
0
1
0
X5
0
0
0
-1
X6
-1
0
0
1
LD
0
200
80
60
Z
X3
X4
X6
1
0
0
0
1
1
1
0
2
1
0
1
0
1
0
0
0
0
1
0
-1
0
0
-1
0
0
0
1
260
200
80
60
Z
X3
X4
X2
1
0
0
0
1
1
1
0
0
0
0
1
0
1
0
0
0
0
1
0
1
1
0
-1
-2
-1
0
1
140
140
80
60
-1
-2
60
115
X3
X1
X2
0
0
0
0
1
0
0
0
1
1
0
0
-1
1
0
1
0
-1
-1
0
1
60
80
60
Z
X5
X1
X2
1
0
0
0
0
0
1
0
0
0
0
1
-1
1
0
1
0
-1
1
-1
0
1
0
0
-1
-1
0
0
0
60
80
120
Como se observa en ltimo tablero las variables artificiales tienen valor cero lo
cual significa que el problema tiene solucin.
FASE II: Por consiguiente la segunda fase comprende de la funcin objetivo
inicial y la informacin de las variables bsicas del ltimo tablero de la primera
fase donde, si se desea, se puede omitir la informacin referente a las variables
artificiales.
Z
X5
X1
X2
Z
1
0
0
0
X1
-3
0
1
0
X2
-8
0
0
1
Z
X5
X1
X2
1
0
0
0
0
0
1
0
0
0
0
1
X3
X4
0
-1
1
-1
X5
0
1
0
0
-5
-1
1
-1
0
1
0
0
X6
LD
0
60
80
120
1200
60
80
120
X2 = 120,
ZMIN = 1200
116
117
x0
Entonces el programa lineal dual est definido por:
D: Maximizar wb
Sujeto a:
wA c
w0
Ntese que existe exactamente una variable dual por cada restriccin primal, y
exactamente una restriccin dual por cada variable primal. Despus se dir ms
sobre esto.
Considrese el siguiente programa lineal y su dual:
P: Minimizar 6x1 + 8x2
Sujeto a:
3x1 + x2 4
5x1 + 2x2 7
x1 , x2 0
Su dual ser:
D: Maximizar 4w1 + 7w2
Sujeto a:
3w1 + 5w2 6
w1 + 2w2 8
w1, w2 0
En teora para aplicar la definicin cannica de dualidad primero se debe convertir
el programa lineal primal al formato anterior. Sin embargo, en la prctica es posible
escribir inmediatamente el dual de cualquier programa lineal.
Forma estndar de dualidad
.
Otra definicin equivalente se aplica cuando las restricciones son igualdades.
Supngase que el programa lineal primal est dado en la forma:
P: Minimizar Cx
118
Sujeto a:
Ax = b
x=0
Entonces el programa lineal dual est definido por:
D: Maximizar Wb
Sujeto a:
wA = c
w no restringida
Considrese el siguiente programa lineal y su dual:
P: Minimizar 6x1 + 8x2
Sujeto a:
3x1 + x2 - x3 = 4
5x1 + 2x2 - x4 = 7
x1, x2, x3, x4 = 0
D: Maximizar 4w1 + 7w2
Sujeto a:
3w1 + 5w2 = 6
w1 + 2w2 = 8
-w1
=0
-w2 = 0
w1, w2 no restringidas
Observacin 1: El dual del dual es el primal
Este lema indica que las definiciones se pueden aplicar al revs. Los trminos
primal y dual son relativos al marco de referencia que se seleccione.
119
teora, esto no presenta problema alguno porque se pueden aplicar las tcnicas de
transformacin para convertir cualquier problema mixto a una de las formas primal
o dual.
Considere el siguiente programa lineal.
PASO 1: Max z = 2x1 + x2
Sujeto a:
x1 + x2 2
2x1 - x2 3
x1 - x2 = 1 => x1 - x2 1
=> x1 - x2 1
x1 0, x2 nrs
PASO 2 : Max z = 2x1 + x2 - x2
Sujeto a :
x1 + x2 - x2 2 => -x1 - x2 + x2 -2
2x1 x2 + x2 3 => - 2x1 + x2 - x2 -3
x1 x2 + x2 1 => x1 x2 + x2 1
x1 x2 + x2 1 => -x1 + x2 - x2 -1
x1, x2, x2 0
PASO 3 : Min w = -2y1 - 3y2 + y3 - y3
Sujeto a:
-y1 - 2y2 + y3 - y4 2
-y1 + y2 - y3 + y4 1
2y1 - y2 + y3 y4 -1
y1, y2, y3, y4 0
De este ejemplo se ve que las restricciones de la forma mayor o igual que en el
problema de minimizacin dan origen a variables = 0 en el problema de
maximizacin dan origen a variables no restringidas.
Tabla: Relaciones entre problemas primario y dual.
120
Variables
Restric_
ciones
MINIMIZACION
DE PROBLEMA
=0
=0
No restringido
MAXIMIZACION
DE PROBLEMA
=
=
=
Variables
=
=
=
=0
=0
No restringido
Restricciones
121
122
-w1 + w2 = -1
w1 , w2 no restringidas
Graficando ambos se observa que ninguno tiene solucin factible
x2
w2
x1
w1
(c - w*A) x* = 0
123
ptimo
ptimo
124
no acotado
no factible
no acotado
no factible
no factible
no acotado o no factible
no factible
no acotado o no factible
I = 1,, m
Este es un teorema muy importante que relaciona los problemas primal y dual.
obviamente indica que al menos uno de los factores en cada una de las expresiones
anteriores debe ser cero. En particular,
x*j > 0 w*aj = cj
w*aj < cj x*i = 0
w*i > 0 aix* = bi
125
126
Variable
xj
cj
z
bi
aij
Significado
Nivel de actividad j ( j = 1, 2,n)
Utilidad unitaria de la actividad j
Utilidad total
Cantidad de recurso i disponible
por cada unidad de actividad j
127
variables bsicas para obtener la solucin primal; es decir, y* (la solucin optima)
representa el precio sombra del recurso i -valor marginal de este recurso.
Como cada unidad de actividad j en el problema primal consume unidades del
recurso i
m
a y
ij
i 1
a y crecursos debe ser por lo menos tanto como si fuese usados por una
ij
i 1
Minimizar y =
b y
i 1
z = cx
Ax = b
x=0
Minimizar
Sujeto a
yo = yb
Ay = c
y =0
Dual
128
mueve de una solucin bsica factible del problema dual a una solucin bsica
factible mejorada, hasta alcanzar la optimalidad del dual (y tambin del primal), o
bien hasta concluir que el dual es no acotado y que el primal es no factible.
Resumen del mtodo dual simplex (problema de minimizacin)
PASO INICIAL:
Encuntrese una base B del primal tal que zj - cj = cBB-1aj - cj = 0 para todo j.
PASO PRINCIPAL:
1. Si b = B-1 b 0, el proceso termina; la solucin presente es ptima. En
caso contrario, seleccinese el rengln pivote r con b r < 0, digamos
b r = Mnimo { b i}.
2. Si yrj 0 para todo j, el proceso termina; el dual es no acotado y el
primal es no factible. En caso contrario, seleccinese la columna pivote
k mediante la siguiente prueba de la razn mnima:
z k ck
3. Pivotese en
Minimo
rk
rk
z c :
y
y
rj
rj
y regrese al paso 1.
129
Z
x4
x1
Z
x2
x1
Z
1
0
0
x1
-2
-1
-2
z
1
0
0
z
1
0
0
x1
0
0
1
x1
0
0
1
x2
-3
-2
1
x2
-4
-5/2
-1/2
x2
0
1
0
x3
-4
-1
-3
x4
0
1
0
x5
0
0
1
LD
0
-3
-4
x3
-1
1/2
3/2
x4
0
1
0
x5
-1
-1/2
-1/2
LD
4
-1
2
x3
-9/5
-1/5
7/5
x4
-8/5
-2/5
-1/5
x5
-1/5
1/5
-2/5
LD
28/5
2/5
11/5
130
pivoteos, a travs de una serie de soluciones bsicas factibles duales hasta que la
solucin bsica primal complementaria asociada es factible, satisfaciendo as todas
las condiciones para optimalidad de Kuhn - Tucker.
En esta seccin se describir un mtodo, llamado el algoritmo primal - dual similar
al mtodo dual simplex, el cual empieza con factibilidad dual y proceda a obtener
factibilidad primaria, manteniendo durante el proceso holgura complementaria.
Una diferencia importante entre el mtodo dual simplex y el mtodo primal - dual
es que este ltimo no requiere que una solucin factible dual sea bsica. Dada una
solucin factible dual, se determina las variables primales que corresponden a
restricciones duales ligantes o activas (de tal manera que la holgura
complementaria se satisface). Usando la fase I del mtodo simplex, se trata de
alcanzar la factibilidad primal, se cambia la solucin factible dual en tal forma que
se admita al menos una nueva variable en el problema de la fase I. Esto se contina
hasta que, o bien la solucin primal se hace factible, o bien, la solucin dual se
hace no acotada.
Resumen del algoritmo primal - dual (problema de minimizacin)
PASO INICIAL:
Seleccinese un vector w tal que waj - cj = 0 para todo j.
PASO PRINCIPAL:
1. Sea Q = { j : waj - cj = 0} y resuelvase el siguiente problema restringido:
Minimizar 0 x j 1 x o
j Q
Sujeta a
a x x
j Q
xj 0 para j Q
xa 0
131
( w a j c j )
Minimo
v *a j
: v * a j 0 0
132
Una solucin factible dual inicial est dada por w = (w 1 , w2 ) = (0,0). Sustituyendo
w en cada restriccin dual, se encuentra que las dos ltimas restricciones duales
son estrictas de manera que Q = {6,7}. Denotando las variables artificiales por X 8
y X9, el problema primal restringido resulta ser el siguiente:
Minimizar x8 + x9
Sujeto a:
- x6 + x8 = 6
-x7 + x9 = 3
x6, x7, x8, x9 0
Es claro que la solucin ptima primal restringido es (x 6, x7, x8, x9 ) = (0,0,6,3) y el
objeto ptimo es xo = 9. El dual de este problema primal restringido es el siguiente:
Maximizar 6v1 + 3v2
Sujeto a:
-v1 = 0
-v2 = 0
v1 = 1
v2 = 1
v1, v2 no restringida
Utilizando holgura complementaria, se ve que, puesto que x 8 y x9 son bsicas, las
dos ltimas restricciones duales deben ser holgura y v* = (v* 1 , v*2) = (1,1).
Calculando v*aj para cada columna j, se obtiene v*a1 = 3, v*a2 = 0, v*a3 = 0, v*a4 =
7, y v*a5 = -3. Por lo tanto, se determina como sigue:
3 6 7
Minimo , , 1
3 3 7
y w1 = (0,0) +1(1,1) = (1,1).
133
Con la nueva solucin dual w1, se calcula de nuevo Q y se obtiene Q ={1,4} esto
da el siguiente problema primal restringido:
Minimizar x8 + x9
Sujeto a:
2x1 + 6x4 + x8 = 6
x1 + x4 +x9 = 3
x1 , x4 , x8 , x9 0
Esta vez una solucin al problema restringido est dada por:
(x1, x4 , x8 , x9 ) = (3 , 0, 0 , 0)
Con Xo = 0. As pues, se tiene una solucin ptima del problema original con
soluciones ptimas primal y dual dadas por:
(x*1, x*2 , x*3 , x*4, x*5, x*6 , x*7 ) = (3 , 0, 0 , 0 , 0, 0 , 0) y,
(w*1, w*2 ) = (1,1)
134
PELOTA TIPO
CUERO
GANANCIA
(por tipo)
(pie /tipo)
(S/./unid.)
1
2
0.3
0.2
2
1
200
150
Disponibilidad
18
200
(1)
Cuero
2X1 +
(2)
X2 < 200
X1, X2 > 0
La solucin del problema es:
Z = 13500;
X1 = 0;
X2 = 90
Y1
135
Cuero
Y2
Observe que se ha utilizado todas las horas mquina (restriccin limitante) y que
hay un excedente de 110 pies2 de cuero (restriccin no limitante).
Esto implica que para aumentar la ganancia e necesario aumentar la disponibilidad
del recurso 1.
Ahora procederemos a hallar el dual para realizar algunas discusiones:
El programa dual esta dado por.
Min W = 18Y1 + 200Y2
Sujeto a:
0.3Y1 + 2Y2 > 200
(1)
(2)
Y1, Y2 > 0
Si se aumenta, por ejemplo, en 2 unidades las horas mquina (recurso limitante) del
programa primal, se tiene:
Max Z = 200X1 + 150X2
Sujeto a:
0.3X1 + 0.2X2 < 20
2X1 + X2 < 200
La solucin del problema es:
Z = 15000;
X1 = 0;
X2 = 100
750 Y1
b
20 18
136
Este valor representa el precio dual del primer recurso, es decir que Z aumentar en
750 unidades si las horas mquina aumentan en 1 hora.
Si ahora aumentamos, por ejemplo, en 2 unidades el recurso cuero, se tiene:
Max Z = 200X1 + 150X2
Sujeto a:
0.3X1 + 0.2X2 < 18
2X1 + X2 < 202
La solucin del problema es:
Z = 13500;
X1 = 0;
X2 = 90
0 Y2
b
202 200
Como se sabe el recurso cuero esta en exceso y por consiguiente incrementar este
recurso no aporta en nada a la funcin objetivo.
La formulacin estndar y el ltimo tablero del primal es como sigue:
Max Z = 200X1 + 150X2 + 0X3 + 0X4
Sujeto a:
0.3X1 + 0.2X2 + X3 = 18
2X1 + X2 + X4 = 200
X1
X2
X3
X4
LD
25
75
1350
137
X2
1.5
0.5
X4
25
75
191
Como se observa (Z3 C3) = 75 y (Z4 C4) = 0 [debajo de X3 y X4] son los valores
de las variables U1 y U2 respectivamente, esto se explica por lo siguiente:
El tablero puede ser escrito de la siguiente forma:
Z
XB
Z
1
0
XB
0
I
XN1
CBB N1 CN1
B-1N1
-1
XH
CBB-1I 0
B-1
LD
CBB-1b
B-1b
Donde:
XN = (XN1 XH)
XN1 = (Xm+1 Xm+2 ... Xn)
XH = (Xn+1 Xn+2 ...
Xn+m)
Nutriente
A
Ingredientes
1
2
3
4
3
2
Requerimientos
u./comida.
20
138
B
C
D
E
Costo $/g
5
6
3
1
2
1
2
1
2
2
3
1
200 300 250
30
10
5
10
Sujeto a:
4X1 + 3X2 + 2X3 20
Nutriente A
Nutriente B
Nutriente C
Nutriente D
Nutriente E
Con Xi
0, i =1, 2,3.
139
Un gramo del alimento 1 cuesta $200 y aporta cuatro unidades del nutriente A,
cinco del nutriente B, uno del C, dos del D y dos del E. Por lo tanto, por esas
cantidades de los nutrientes puede pagarse en total un mximo de $200.
Similarmente, como un gramo del alimento 2 cuesta $300 y aporta tres unidades
del nutriente A, seis del B, dos del C, uno del D y tres del E, lo mximo que
podemos pagar conjuntamente por esas cantidades de los nutrientes es $300.
Si denotamos respectivamente con las variables YA, YB, YC, YD, YE, los precios
mximos que se pueden pagar por la pastilla con una unidad de cada uno de los
nutrientes, y efectuamos un anlisis para todos los ingredientes, obtenemos el
siguiente modelo de programacin lineal.
Maximizar Ventas ZD = 20YA + 30YB + 10YC + 5YD + 10YE
Sujeto a:
4YA + 5YB + 1YC + 2YD + 2YE 200 ingrediente 1
3YA + 6YB + 2YC + 1YD + 3YE 300 ingrediente 2
2YA + 3YB + 1YC + 2YD + 1YE 250 ingrediente 3
Con YA, YB, YC, YD, YE 0
Este segundo modelo representa el enfoque dual del primero y de nuevo podemos
verificar que se presentan ciertas relaciones estructurales, a saber
1.
140
4.3 Las variables de ambos problemas solo pueden tomar valores mayores o
iguales que cero.
Pero las relaciones de forma no son las ms importantes para nuestro estudio de la
dualidad en Programacin lineal, como si lo son las relaciones lgicas existentes
entre sus soluciones ptimas y el significado econmico de las variables del
modelo dual.
Resolviendo ambos modelos obtenemos, para el PRIMAL:
PRIMAL
141
Nutriente
A
B
C
D
E
cantidad
cantidad
contenida
20
34
10
8
16
requerida
20
30
10
5
10
exceso
0
4
0
3
6
Podemos decir que si nos aumentaran, por ejemplo en una unidad, las exigencias
de alguno de los nutrientes B, D o E, el costo (y la mezcla actual de ingredientes)
no cambiara ya que las unidades de esos nutrientes con que efectivamente queda
la comida, son superiores a las exigencias, as estas suban en una unidad.
142
143
144
En cambio si, por ejemplo, nos piden que la comida debe contener una unidad mas
del nutriente A,( la exigencia ser de 21 unidades y no de 20) la solucin actual
( X1=2 y X2= 4,), no cumplira esta nueva condicin y por ello no sera ms la
solucin ptima. Ser necesario encontrar una mezcla diferente de los
ingredientes, de tal forma que el contenido del nutriente A sea de 21 unidades y
esta mayor exigencia elevar el costo de la comida.
Escribamos el modelo modificado al incluir esta nueva exigencia para el contenido
del nutriente A y resolvmoslo.
Minimizar Costo:
145
Sujeto a:
4X1 +
5X1 +
1X1 +
2X1 +
2X1 +
Con Xi
3X2
6X2
2X2
1X2
3X2
+ 2X3
+ 3X3
+ 1X3
+ 2X3
+ 1X3
21
30
10
5
10
Nutriente A
Nutriente B
Nutriente C
Nutriente D
Nutriente E
0, i =1, 2,3.
146
147
148
149
150
2.
L a demanda flucta.
3.
antiguas.
4.
5.
Los gerentes y ejecutivos desean determinar la forma en que estos cambios afectan
a la solucin ptima del problema primitivo de programacin lineal.
EJEMPLO
Una compaa textil, incursiona en el mercado de bolsas de tela por despacho de
mercaderas y fabrica dos clases de bolsas: el modelo estndar por tiendas y
bodegas y el de lujo para supermercados y grandes almacenes. El proceso de
fabricacin es corte y teido, costura, terminado e inspeccin y embalaje, cuya
programacin lineal es la siguiente:
151
Sujeto a:
Corte y Teido
Costura
Terminado
Inspeccin y Emb.
La solucin optima X1= 540 bolsas estndar, X2 = 252 bolsas de lujo y Z = 7688,
donde X1 da $10 de utilidades y X2 da $9.0.
Supongamos que posteriormente debido a una reduccin en el precio, la
contribucin a las utilidades de las bolsas estndar se reduce a $7.0, puede
utilizarse el anlisis de sensibilidad para determinar si el programa de produccin
de 540 bolsas estndar y 252 bolsas de lujo sigue siendo la mejor solucin, si lo es,
no habr necesidad de resolver un programa lineal modificado que tenga 7X 1 + 9X2
como funcin objetivo
152
Despejando
La ecuacin de la recta B:
Despejando
153
154
Para calcular el intervalo de optimidad para la contribucin a las utilidades por las
bolsas estndares, se mantiene fija la contribucin a las utilidades por las bolsas de
lujo, en su valor inicial
= 9, luego:
155
CB
Caso I: Maximizacin
En la K-sima columna (K = m+1, m+2,... n), se tiene:
ZK CK 0
CBB-1aK CK 0
156
(C B1 , C B 2 ,..., C Br C Br ,..., C Bm )
Y1k
Y
2k
.
Ck 0
Yrk
.
.
Ymk
Operando resulta:
; K = m+1, m+2,, n
(1)
(2)
EJEMPLO
Se desea optimizar la produccin de mesas y sillas de una fbrica industrial, para lo
cual se dispone del nmero de horas mquina en cada una de las secciones
siguientes:
MAQ.
A
B
UTILIDAD
MESAS
2
4
10
SILLAS
2
2
8
HORAS - MAQUINA
20
28
157
Sujeto a:
Z
X2
X1
Variacin de
Z
1
0
0
X1
0
0
1
X2
0
1
0
X4
1
-1/2
1/2
LD
88
6
4
X3
3
1
-1/2
k = 3, 4;
r=2
Para k = 3
Entonces:
C1 6
158
Para k = 4
Entonces:
Por tanto:
Variacin de
C1 -2
-2 C1 6
k = 3, 4;
r=1
Para k = 3
Entonces:
C2 -3
159
Para k = 4
Entonces:
Por tanto: -3
C2 -2
C2 2
puesto que:
= CBj - CBj
Zj - CBj = 0
160
Por tanto:
161
Por tanto:
Caso II: Minimizacin
Zr - Cr 0
Cr Zr Cr
4.3 CAMBIOS EN LA DISPONIBILIDAD DE RECURSOS
Aqu se calcula la variacin de cada uno de los recursos, manteniendo la misma
base.
Se sabe que:
b1
b1
b2
b2
.
-1
-1
XB = B b y que:
b = br , B b = br
.
.
bm
bm
Si se altera el recurso
162
b1
b2
.
b =
br br
bm
v11
v 21
.....
.....
.....
vm1
.....
v1r
.....
.....
.....
.....
.....
.....
v2r
.....
.....
.....
vmr
.....
.....
.....
.....
.....
v1n
v 2 n
.....
.....
.....
vmn
.....
.....
.....
vm1
.....
v1r
.....
.....
.....
.....
.....
.....
v2r
.....
.....
.....
vmr
.....
.....
.....
.....
.....
b1 v1r
v1n
v2n
.....
.....
.....
vmn
b1
b2
br br
bm
br 0
b 2 v 2 r br 0
.....................
-1
B b = br vrr br 0
.....................
bm vmr br 0
En el problema anterior:
163
20 b1
b
2 8 b2
Entonces:
164
Entonces:
Restricci
n
<b
>b
Variacin en el
recurso b
Si el recurso b
aumenta una
unidad
Si el recurso b
disminuye en
una unidad
Si el recurso b
aumenta en una
unidad
Consecuencia en
funcin Objetivo
la F.O. mejora en
una tasa igual al
precio dual
la F.O. empeora en
una tasa igual al
precio dual
la F.O. empeora en
una tasa igual al
precio dual
165
si el recurso b
disminuye en
una unidad
MAX
=b
MIN
si el recurso b
aumenta en una
unidad
si el recurso b
disminuye en
una unidad
si el recurso b
aumenta en una
unidad
si el recurso b
disminuye en
una unidad
la F.O. mejora en
una tasa igual al
precio dual
la F.O. mejora en
una tasa igual al
precio dual
la F.O. empeora en
una tasa igual al
precio dual
la F.O. empeora en
una tasa igual al
precio dual
la F.O. mejora en
una tasa igual al
precio dual
Sujeto a:
AX b
X0
El vector columna aj se cambia por aj, entonces:
Sujeto a:
AX b
X0
-1
Este cambio afecta al producto B (aj). Luego tambin a Zj - Cj entonces se tiene:
166
-1
Zj - Cj = CBB a'j - Cj
Que deber cumplir con la funcin de optimizacin (Z j - Cj 0), en caso contrario
se pivotea para encontrar el ptimo.
EJEMPLO
Winco vende 3 productos 1, 2 y 3. En la tabla se dan los recursos requeridos para
producir una unidad de cada producto y los precios de venta de cada producto.
Materia Prima
Horas de Trabajo
Precio Venta
Producto 1
1
2
6
Producto 2
2
1
10
Producto 3
4
3
8
Z
X2
X1
Z
1
0
0
X1
0
0
1
X2
0
1
0
X3
38/3
5/3
2/3
X4
14/3
2/3
-1/3
X5
2/3
-1/3
2/3
LD
136/3
4/3
16/3
167
Se desea cambiar a3 =
-1
Y3 = B a'3
2/3
Luego:
1 / 3 2
1 / 3 2 / 3 2
Z3 - C3 =
10
2 / 3
2 / 3
2 / 3
- 8 = 8/3
2 / 3
EJEMPLO
En el problema de las mesas y sillas se desea incorporar la lnea de produccin de
repisas, si se necesitan una hora por cada mquina para elaborar una repisa siendo
la utilidad de S/.8 por unidad Cmo cambia la solucin? Recordemos que este es
el ltimo tablero de la funcin objetivo:
Z
X2
Z
1
0
X1
0
0
X2
0
1
X3
3
1
X4
1
-1/2
LD
88
6
168
X1
-1/2
1/2
Solucin:
La informacin para este problema ser:
Mquina
A
B
Utilidad
Mesas
2
4
10
Sillas
2
2
8
Repisas
1
1
8
Horas Mquina
20
28
El nuevo PL es:
Sujeto a:
2X1 + 2X2 + X3 20
4X1 + 2X2 + X3 28
X1, X2, X3 0
Hallando el valor de Z3 - C3, conociendo los valores de C3 = 8 (utilidad)
-1
1
y a 3 = Y 3 = B a3 =
1
1
1/ 2
1 / 2
1/ 2
1
1
1 / 2
Entonces:
Z3 - C3 =
10
1 / 2
- 8= -4
0
Z
X2
X1
Z
1
0
0
X1
0
0
1
X2
0
1
0
X3
3
1
-1/2
X4
1
-1/2
1/2
X5
-4
1/2
0
LD
88
6
4
169
Z
1
0
0
X1
6
2
2
X2
8
2
0
X3
8
1
-1
X4
0
0
1
X5
0
1
0
LD
160
20
8
Se concluye que se debe producir slo repisas, dado que las Variables bsicas (de
decisin) son ahora X5 y X4 (no se considera X4 porque es una variable de
holgura).
4.7 ADICIN DE UNA RESTRICCIN
Se necesita analizar si esta hace variar el vector solucin, si ocurre esto se pivotea,
si para un Programa lineal al introducir la restriccin m +1 el sistema tiene m + 1
filas y m + n + 1 columnas.
La adicin de una nueva restriccin puede dar origen a una de dos condiciones:
1. La restriccin satisface la solucin actual y en este caso la restriccin es
redundante, y, por lo tanto, su adicin no alterara la solucin.
2. La solucin actual no satisface la restriccin. En este caso, la nueva solucin se
obtiene utilizando el Mtodo Simplex Dual.
EJEMPLO
Para el problema de mesas y sillas se ha adicionado la mquina C la que produce
2 mesas y 1 silla, para lo cual dispone de 28 horas - mquina. Se desea saber si la
solucin ptima vara con el cambio.
Maquina
A
B
C
Mesas
2
4
2
Sillas
2
2
1
Horas - Mquina
20
28
28
170
Utilidad
10
El nuevo PL es:
Sujeto a:
2X1 + 2X2 20
4X1 + 2X2 28
2X1 + X2 28
X1, X2 0
Incorporando esta informacin en el tablero resulta:
Z
X2
X1
X5
Z
1
0
0
0
X1
0
0
1
2
X2
0
1
0
1
X3
3
1
-1/2
0
X4
1
-1/2
1/2
0
X5
0
0
0
1
LD
88
6
4
28
Por teora de matrices podemos ingresar toda una fila sin tener problemas
Z
X2
X1
X5
Z
1
0
0
0
X1
0
0
1
0
X2
0
1
0
0
X3
3
1
-1/2
0
X4
1
-1/2
1/2
-1/2
X5
0
0
0
1
LD
88
6
4
14
171
Caloras
Bizcocho
Helado de chocolate
(1 bola)
Refresco de Cola
(1 botella)
Pastel de Queso con pia
400
Chocolate
(onzas)
3
Azcar
(onzas)
2
Grasa
(onzas)
2
200
150
500
172
Sujeto a:
400X1 + 200X2 + 150X3 + 500X4 500 (Rest. de caloras)
3X1 + 2X2 6
(Rest. de chocolate)
(Rest. de grasa)
SUBJECT TO
2)
3)
4)
5)
END
OBJECTIVE FUNCTION VALUE
1)
VARIABLE
X1
90.000000
VALUE
0.000000
REDUCED COST
27.500000
173
3.000000
1.000000
0.000000
X2
X3
X4
ROW
2)
3)
4)
5)
SLACK OR SURPLUS
250.000000
0.000000
0.000000
5.000000
0.000000
0.000000
50.000000
DUAL PRICES
0.500000
-2.000000
-7.000000
0.000000
N ITERATIONS = 5
RANGES IN WHICH THE BASIC IS UNCHANGED:
OBJ COEFFICIENT RANGES
VARIABLE
CURRENT COEF
X1
X2
X3
X4
50.000000
20.000000
30.000000
80.000000
ALLOWABLE
INCREASE
INFINITY
18.333334
10.000000
INFINITY
ALLOWABLE
DECREASE
27.500000
5.000000
30.000000
50.000000
CURRENT RHS
2)
3)
4)
5)
500.000000
6.000000
10.000000
8.000000
ALLOWABLE
INCREASE
250.000000
4.000000
INFINITY
5.000000
ALLOWABLE
DECREASE
INFINITY
2.857143
4.000000
INFINITY
174
Solucin:
Como los bizcochos
y el pastel de queso con pia
tienen costos reducido
diferente de cero, se presenta el caso I.
De la salida del LINDO se ve que la base actual es ptima, si y slo si:
22.5 = 50 - 27.5
30 = 80 - 50
50 + =
80 + =
175
CASO I
A. Suponga que las caloras necesarias disminuyen hasta 400 y que el
requerimiento de grasa aumenta hasta 10 onzas. Permanecer ptima la base
actual?, Cul ser la nueva solucin optima?
Solucin:
Como ambas restricciones no son obligatorias, se presenta el caso I. De la
corrida del lindo, observamos que los intervalos permisibles para las
restricciones de las caloras de la grasa son:
500 calorias necesarias 500 250 750
8 requerimento de grasa 8 5 13
176
b2 = 8 - 6 = 2,
Ejemplo 1:
Un empresario, fabricante de artculos de cuero ha decidido lanzar un nuevo
producto de bolsas de piel para damas. El distribuidor Alda de lnea de cartera,
bolsas y bolsones acepta comprar todas las bolsas que fabrique la empresa. Las
operaciones necesarias para la fabricacin de las bolsas son las siguientes:
1)
2)
Coser
3)
Terminar
4)
Inspeccionar y embalar
Corte y
teido.
7/10
1
603
TIEMPO DE PRODUCCIN
Costura
Terminado
Insp. y
Emb.
1/2
1
1/10
5/6
2/3
1/4
600
708
135
Utilidad
10
9
-
177
Sujeto a:
7/10 X1 + 10X2 <= 630
1/2X1 + 5/6X2 <= 600
X1 + 2/3X2 <= 708
1/10X1 + 1/4X2 <= 135
X1, X2 >= 0
El desarrollo en Lindo es:
OBJECTIVE FUNCTION VALUE
1)
7662.147000
VARIABLE
X1
X2
ROW
2)
3)
4)
5)
VALUE
538.418091
253.107346
SLACK OR SURPLUS
0.000000
120.711861
0.000000
17.881355
REDUCED COST
0.000000
0.000000
DUAL PRICES
4.331450
0.000000
6.967985
0.000000
N ITERATIONS = 0
RANGES IN WHICH THE BASIC IS UNCHANGED:
178
CURRENT COEF
X1
X2
10.000000
9.000000
ALLOWABLE
INCREASE
3.432836
5.285714
ALLOWABLE
DECREASE
3.700000
2.300000
CURRENT RHS
2)
3)
4)
5)
630.000000
600.000000
708.000000
135.000000
ALLOWABLE
INCREASE
51.885242
INFINITY
192.000000
INFINITY
ALLOWABLE
DECREASE
134.400009
120.711861
126.599998
17.881355
VARIABLE
X1
X2
VALUE
538.418091
253.107346
179
ROW
2)
3)
4)
5)
REDUCED COST
0.000000
0.000000
SLACK OR SURPLUS
0.000000
120.711861
0.000000
17.881355
/ corte y teido
/ costura
/ terminado
/ inspeccin y embalaje
ROW
2)
3)
4)
5)
DUAL PRICES
4.331450
0.000000
6.967985
0.000000
/ corte y teido
/ costura
/ terminado
/ inspeccin y embalaje
180
Se puede afirmar que una hora adicional de tiempo corte t teido mejora (aumenta)
el valor de la funcin objetivo en $4.33 y una hora adicional de tiempo de
Terminado mejora (aumenta) en $6.967.
En consecuencia, aumenta el tiempo de corte y teido de 630 a 631 horas,
manteniendo constante todos los dems coeficientes del problema, aumentando las
utilidades de la compaa de 7662.1 + 4.33 = 7666.43, similarmente en el caso de
terminado aumenta el tiempo de 7662.1 + 6.97 =7669.07. Los precios duales cero
sealan que aumentar las horas disponibles de estos recursos no mejora el valor de
la funcin objetivo.
Si el precio dual es negativo por ejemplo 4.33 significa que al aumentar el tiempo
de corte y teido de 630 a 631 horas las utilidades disminuiran en $4.33
Ejemplo 2:
Considere el PL siguiente y su tablero ptimo:
ST
2X1 + 5X2 + 4X3 <=18
3X1 7X2 + 3X3 <=10
X1 +
ROW(BASIS)
1 ART
2 slk4
3 X2
4 X3
X1
0.907
0.326
-0.140
0.674
X3 <= 9
X2
0.000
0.000
1.000
0.000
X3
0.000
0.000
0.000
1.000
slk2
1.047
-0.163
0.070
0.163
slk3
0.605
-0.116
-0.093
0.116
slk4
0.000
1.000
0.000
0.000
LD
24.884
4.907
0.326
4.093
181
Solucin:
a)
Para x2
0.907+ (-0.140) Cs2 0
Cs2 6.479
Cs2 -14.957
Cs2 6.505
Para x3
0.907+ (0.674) Cs3 0
Cs3 -1.346
Cs3 -6.423
Cs3 -5.216
182
Valor
0
0.326
4.093
Costo reducido
0.907
0
0
Valor
0
0
4.907
Precio Dual
1.047
0.605
0
RANGOS
Var.
X1
X2
X3
Coef.
3
1
6
Aumento
infinito
6.479
infinito
Disminucin
infinito
14.957
1.346
Cs2 42.302
Cs2 3.505
Cs2 -35.284
Z = 26.996
183
PROBLEMAS RESUELTOS
1.
Beerco fabrica cerveza tipo ALE y BEER, a partir de trigo, lpulo y malta.
Actualmente, se disponen de 40 lb. de trigo, 30 lb. de lpulo y 40 lb. de malta.
Un barril de ALE se vende a 40 dlares y requiere 1 lb. de trigo, 1 lb. de lpulo
y 2 lb. de malta. Un barril de BEER se vende a 50 dlares y se necesitan 2 lb.
de trigo, 1 lb. de lpulo y 1 lb. de malta. Beerco puede vender toda la ALE y
toda la BEER que produce. Suponiendo que la meta de Beerco es maximizar el
ingreso total de las ventas, Beerco tendr que resolver el PL siguiente:
Sujeto a:
ALE + 2 BEER 40
ALE + BEER 30
2ALE + BEER 40
ALE, BEER 0
ALE = barriles de ale producidos y BEER = barriles de cerveza producidos.
En la tabla se muestra un cuadro ptimo para este PL.
Z
BEER
S2
ALE
Z
1
0
0
0
ALE
0
0
0
1
BEER
0
1
0
0
S1
20
2/3
-1/3
-1/3
S2
0
0
1
0
S3
10
-1/3
-1/3
2/3
LD
1200
40/3
10/3
40/3
Sujeto a:
Para k = 5, r = 3
Z5 - C5 + Y35CALE 0
10 + (2/3)CALE 0
CALE -15
185
Para k = 5, r = 1
Z5 - C5 + CBEER 0
10 + (1/3)CBEER 0
CBEER 30
b1 -20
10/3 + (-1/3)b1 0
b1 10
186
40/3 + (-1/3)b1 0
b1 40
40/3 + (0)b2 0
b2 -
10/3 + (1)b2 0
b2 -10/3
40/3 + (0)b2 0
b2 -
187
40/3 + (-1/3)b3 0
b3 40
10/3 + (-1/3)b3 0
b3 10
40/3 + (2/3)b3 0
b3 -20
Z
S1
S2
S3
Z
1
0
0
0
ALE
- 40
1
1
2
BEER
- 50
2
1
1
MALT
-50
0.5
3
3
S1
0
1
0
0
S2
0
0
1
0
S3
0
0
0
1
LD
0
40
30
40
188
-1
Zj - Cj = CVBB aj - Cj
20
0.5
3
10
2/3
-1
A 3 = B a3 =
1/ 3
1/ 3
- 50 = 10 + 30 - 50 = -10
0
1
0
0.5
2 / 3
11 / 6
11 / 6
MALT LQ
-10
-2/3
11/6
11/6
S1
20
2/3
-1/3
-1/3
1 / 3
1 / 3
2/3
3 =
Z
CERV.
S2
ALE
Z
1
0
0
0
ALE
0
0
0
1
BEER
0
1
0
0
S2
0
0
1
0
S3
10
-1/3
-1/3
2/3
LD
1200
40/3
10/3
40/3
Dado que el valor del coeficiente de Malt Liquor es < 0 (= -10), esto indica que la
base no sera optima; por lo tanto, podra usar el algoritmo Simplex a fin de no
tener coeficientes negativos en el rengln o, con lo que la variable Malt Liquor
entrara a la base, convirtindose en variable de decisin.
2. Radioco fabrica dos tipos de radios. El nico recurso escaso que se necesita para
producir los radios es la mano de obra. Actualmente la compaa tiene dos
trabajadores. El trabajador 1 est dispuesto a trabajar hasta 40 horas a la semana,
y se le para 5 dlares la hora. El trabajador 2 est dispuesto a trabajar hasta 50
horas a la semana, y se le paga 6 dlares la hora. En la tabla siguiente se dan los
precios, as como los recursos requeridos para fabricar cada tipo de radio.
RADIO 1
RADIO 2
189
Precio
Recursos
(dlares)
Requeridos
25
Trabajador 1: 1 hora
a).
Precio
(dlares)
22
Recursos
Requeridos
Trabajador 1: 2 horas
Trabajador 2: 2 horas
Trabajador 2: 1 hora
Materia prima:
Costo: 5 dlares
Materia prima
costo: 4 dlares
siguiente tabla)
b). Para qu valores del precio de un radio tipo 1, la base actual permanece
ptima?
c). Para qu valores del precio de un radio tipo 1, la base actual permanece
ptima?
d). Si el trabajador 1 estuviera dispuesto a trabajar solamente 30 horas a la
semana, permanecera ptima la base actual?
e). Si el trabajador 2 estuviera dispuesto a trabajar hasta 60 horas a la semana,
f).
Sujeto a:
190
X1 + 2X2 40
2X1 + X2 50
X1, X2 0
En la tabla se muestra un cuadro ptimo para este PL.
X1
0
1
0
Z
X1
X2
X2
0
0
1
S1
1/3
-1/3
2/3
S2
4/3
2/3
-1/3
LD
80
20
10
Solucin:
a) Calculando:
Trabajador 1: X1 + 2X2 40
Trabajador 2: 2X1 + X2 50
Desarrollo del PL con la ayuda del LINDO
MAX 3X1 + 2X2
SUBJECT TO:
X1 + 2X2 <= 40
2X1 + X2 <= 50
LP OPTIMUM FOUND AT STEP
80.000000
VALUE
20.000000
10.000000
REDUCED COST
0.000000
0.000000
191
ROW
2)
3)
SLACK OR SURPLUS
0.000000
0.000000
DUAL PRICES
0.333333
1.333333
N ITERATIONS = 2
RANGES IN WHICH THE BASIC IS UNCHANGED:
VARIABLE
X1
X2
ROW
2)
3)
ALLOWABLE
DECREASE
2.000000
0.500000
ALLOWABLE
DECREASE
15.000000
30.000000
1/3 + (-1/3)C1 0
C1 1
k=5
4/3 + (2/3)C1 0
C1 -2
C1 = 3 + 1 = 4
C1 = 3 - 2 = 1
P1 - 22 = 4
P1 - 22 = 1
P1 = 26
P1 = 23
192
c) Valores del precio de un radio tipo 1,en donde la base actual permanece optima
1/3 + (2/3) C2 0
C2 -1/2
4/3 + (-1/3) C2 0
C2 4
C2 = 2 - 1/2
P2 - 20 = 2 - 1/2
P1 = 21.5
C2 = 2 + 4
P2 - 20 = 2 + 4
P1 = 26
b1 60
40 + 60 = 100
10 + (2/3)b1 0
b1 -15
40 - 15 = 25
193
20 + (2/3)b2 0
b2 -30
50 + 30 = 80
10 + (-1/3)b2 0
b2 30
50 - 30 = 20
De otra manera:
Z = 80 - 2x4/3 = 77.33
194
3 3
3
2
28 1
3 0
33 3
Por lo tanto Radioco no tendra que producir radios tipo 3.
AUTOMVIL
DAS EN LA
DAS EN LA
TONELADAS DE
MQUINA TIPO 1
MQUINA TIPO 2
ACERO
0.8
0.6
195
CAMIN
0.7
REPORTE EN LINDO:
MAX 300X1 + 400X2 - 50M1
SUBJECT TO
2) 0.8X1 + X2 - M1 0
3)
M1 98
4) 0.6X1 + 0.7X2
73
5) 2X1+ 3X2
260
6)
X1
88
7)
X2
26
END
LP OPTIMUM FOUND AT STEP 1
196
32540.000000
VARIABLE
X1
X2
X3
ROW
2)
3)
4)
5)
6)
7)
VALUE
88.000000
27.599998
98.000000
SLACK OR SURPLUS
0.000000
0.000000
0.879999
1.200003
0.000000
1.599999
REDUCED COST
0.000000
0.000000
0.000000
DUAL PRICES
400.000000
350.000000
0.000000
0.000000
-20.000000
0.000000
N ITERATIONS = 1
RANGES IN WHICH THE BASIC IS UNCHANGED:
VARIABLE
X1
X2
X3
ROW
2)
3)
4)
5)
ALLOWABLE
DECREASE
INFINITY
25.000000
350.000000
ALLOWABLE
DECREASE
1.599999
1.599999
0.879999
1.200003
197
6)
7)
88.000000
26.000000
1.999999
1.599999
3.000008
INFINITY
Solucin:
a). Si los automviles contribuyeran con 310 dlares, se estara adicionando $ 10 a
la utilidad.
Observando el reporte en LINDO, vemos que 10 est dentro del rango
admisible para un incremento de X1 (automviles); por lo tanto la F.O. seguira
siendo ptima.
Nuevo valor objetivo
c). Carco no utiliza todo el recurso disponible de acero; por lo tanto, no le interesa
comprar 1 ton extra de acero.
0
d). Si Carco tuviera que producir por lo menos 86 automviles (dos automviles
menos que el planteamiento original). Veremos que una disminucin en 2 est
dentro del rango permisible; por lo tanto:
Nueva utilidad = 32540 + (-2)(-20) = 32580
198
El reporte en LINDO nos arroja una F.O. de 32631 dlares, que es mayor a la
F.O. original, lo que significa que aumentara nuestras utilidades.
Por lo tanto, es recomendable producir jeep.
199
Sujeto a:
e) Supngase que se puede utilizar una libra de materia prima para fabricar
0.8 lb. del producto 3. El producto 3 se vende a 24 dlares/lb. y por
procesar 1 lb. de materia prima de 0.8 lb. del producto 3 requiere 7 horas
de mano de obra. Tendra que producir WIVCO algn producto 3?
MAX
13 P1 + 13.2P2 10RM
SUBJECT TO
- P1P2+RM 0
2P1+3P2 200
RM 90
0.33P2 40
200
END
LP OPTIMUM FOUND AT STEP
90.000000
VARIABLE
P1
P2
RM
VALUE
70.000000
20.000000
90.000000
REDUCED COST
0.000000
0.000000
0.000000
ROW
2)
3)
4)
5)
SLACK OR SURPLUS
0.000000
0.000000
0.000000
33.400002
DUAL PRICES
-12.600000
0.200000
2.600000
0.000000
N ITERATIONS = 3
RANGES IN WHICH THE BASIC IS UNCHANGED:
OBJ COEFFICIENT RANGES
VARIABLE
CURRENT COEF
P1
P2
RM
13.000000
13.200000
-10.000000
ROW
2)
3)
4)
ALLOWABLE
INCREASE
0.200000
1.300000
INFINITY
ALLOWABLE
DECREASE
0.866667
0.200000
2.600000
ALLOWABLE
DECREASE
10.000000
20.000000
23.333334
201
5)
40.000000
INFINITY
33.400002
Solucin:
a.
sta disminucin
SLACK OR SURPLUS
0.000000
DUAL PRICES
-12.600000
Por lo que WIVCO estar dispuesto a pagar 12.60 dlares por rentar dicha
mquina
d. De igual manera, la respuesta la obtendremos directamente observando los
precios duales del reporte LINDO:
ROW
3)
SLACK OR SURPLUS
0.000000
DUAL PRICES
0.200000
Por lo que WIVCO estar dispuesto a pagar 0.20 dlares (20 centavos $) por
una hora adicional de mano de obra.
202
Sujeto a:
203
Sujeto a:
204
1)
4250.000000
VARIABLE
A1
A2
L
ROW
2)
3)
4)
5)
6)
VALUE
25.000000
20.000000
350.000000
SLACK OR SURPLUS
0.000000
0.000000
0.000000
15.000000
40.000000
REDUCED COST
0.000000
0.000000
0.000000
DUAL PRICES
75.000000
12.000000
2.500000
0.000000
0.000000
N ITERATIONS = 3
RANGES IN WHICH THE BASIC IS UNCHANGED:
VARIABLE
A1
A2
L
ROW
2)
3)
4)
5)
6)
Solucin:
205
a) La mxima cantidad que estar dispuesto a pagar el granjero Leary por una
hora adicional de mano de obra ser menos de 2.5 dlares.
ROW
4)
SLACK OR SURPLUS
0.000000
DUAL PRICES
2.500000
...
SLACK OR SURPLUS
0.000000
DUAL PRICES
75.000000
206
6 A1 + 10 A2 + 3 A3 - L 0
L 350
5 A1 140
4 A2 120
A1, A2, L 0
Lo cual nos da una nueva F.O. de $ 4,350.0, que es una utilidad mayor a la
anterior de 4,250.0, por lo que si se puede cultivar cebada.
6. Con rubes y zafiro Zales Jewelers producen dos tipos de anillos. Un anillo tipo
1 requiere 2 rubes, 3 zafiros y 1 h. de trabajo de un joyero. Un anillo tipo 2
requiere 3 rubes, 2 zafiros y 2 h de trabajo de un joyero. Cada anillo tipo 1 se
vende a 400 dlares y cada anillo tipo 2 a 500 dlares. Se pueden vender todos los
anillos producidos por Zales. Actualmente Zales dispone de 100 rubes, 120 zafiros
y 70 horas de trabajo de un joyero. Se puede compras ms rubes a un costo de 100
dlares el rub. La demanda del mercado requiere una produccin por lo menos de
20 anillos tipo 1 y por lo menos 25 anillos tipo. Para maximizar la ganancia Zales
tendr que resolver el PL siguiente:
= Anillos tipo 1 producidos
= Anillos tipo 2 producidos
= Nmero de rubes comprados
Sujeto a:
207
a) Suponga que cada rub cuesta 190 dlares, en lugar de 100 dlares
Todava comprara Zales rubes? Cul sera la nueva solucin ptima
para el problema?
b) Suponga que Zales solamente tuviera que producir 23 anillos tipo 2 Cul
sera la utilidad de Zales ahora?
c) Cul es la mxima cantidad que tendra que estar dispuesto a pagar Zales
por otra hora de trabajo a un joyero?
d) Cul es la mxima cantidad que tendra que estar dispuesto a pagar Zales
por otro zafiro?
e) Zales considera producir anillos tipo 3. Cada anillo tipo 3 puede venderse a
550 dlares y requiere 4 rubes, 2 zafiros y 1 hora de trabajo de un joyero.
Tendra que producir Zales anillos tipo 3?
MAX 400X1 + 500X2 - 100R
SUBJECT TO
2X1 + 3X2 R 100
3X1 + 2X2 120
X1 + 2X2 70
X1 >= 20
X2 >= 25
END
LP OPTIMUM FOUND AT STEP
19000.000000
VALUE
20.000000
25.000000
15.000000
SLACK OR SURPLUS
REDUCED COST
0.000000
0.000000
0.000000
DUAL PRICES
208
2)
3)
4)
5)
6)
0.000000
10.000000
0.000000
0.000000
0.000000
100.000000
0.000000
200.000000
0.000000
-200.000000
N ITERATIONS = 1
RANGES IN WHICH THE BASIC IS UNCHANGED:
VARIABLE
X1
X2
R
ROW
ALLOWABLE
DECREASE
100.000000
INFINITY
100.000000
2)
100.000000
15.000000
INFINITY
3)
120.000000
INFINITY
10.000000
4)
70.000000
3.333333
0.000000
5)
20.000000
0.000000
INFINITY
6)
25.000000
0.000000
2.500000
209
Solucin:
a) Suponiendo que cada rub cuesta 190 dlares en lugar de 100 dlares;
entonces habra un incremento de 90 dlares; si observamos el reporte
LINDO, 90 estar dentro del intervalo permisible del aumento por lo que la
nueva solucin seguir siendo ptima.
Nueva solucin ptima = 400(20) + 500(25) - 190(15) = 17650 dlares
b) En caso de que Zales solamente tuviera que producir 23 anillos tipo 2 (2
anillos menos), la F.O. permanecer ptima y la nueva utilidad sera.
Nueva solucin ptima
c) La mxima cantidad que estara dispuesto a pagar Zales por otra hora de
trabajo de un joyero ser: 200 dlares.
ROW
4)
SLACK OR SURPLUS
0.000000
DUAL PRICES
200.000000
d) La mxima cantidad que Zales estara dispuesto a pagar por otro Zafiro
seria 0 dlares.
ROW
3)
SLACK OR SURPLUS
10.000000
DUAL PRICES
0.000000
210
Sujeto a:
2X1 + 3X2 + 4X3 R 100
3X1 + 2X2 + 2X3 120
X1 + 2X2 + X3 <= 70
X1 20
X2 25
X1, X2 0
La nueva F.O. seguir siendo 19000 dlares y no se reducir el anillo tipo
3.
7. SOFA S.A. produce cuatro tipos de losetas, las cuales sern vendidas en la
prxima edicin de la Feria del Hogar. Estos cuatro tipos son:
- Romana
- Esparta
- Sicilia
- Atenas
En la tabla se dan los recursos requeridos para producir una unidad de cada
producto y los precios de venta de cada tipo de loseta. SOFA S.A. dispone
de 5 toneladas de Barbotina (Barro Lquido) y 4600 horas de trabajo .la
empresa debe abastecer su stand con una produccin exacta de 950 unidades
en total .Por estudios de mercado y de aceptacin de productos, se exige que
se produzcan por lo menos 400 unidades de losetas Atenas.
Materia Prima (Barbotina)en Kg./unid
Horas de trabajo (hr/unid)
Precio de venta($/unid)
Romana
3
2
4
Esparta
4
3
6
Sicilia
5
4
7
Atenas
6
6
8
SOFA S.A. necesita conocer la cantidad de losetas a producir por cada tipo a fin
de maximizar ingresos.
MAX: 4X1+6X2+7X3+8X4
211
SUBJECT TO
2)
3)
4)
5)
X1 + X2 + X3 + X4 = 950
X4 >= 400
2X1 + 3X2 + 4X3 + 6X4<= 4600
3X1 + 4X2 + 5X3 + 6X4<= 5000
END
LP OPTIMUM FOUND AT STEP
6900.000000
VARIABLE
X1
X2
X3
X4
ROW
2)
3)
4)
5)
VALUE
0.000000
150.000000
400.000000
400.000000
SLACK OR SURPLUS
0.000000
0.000000
150.000000
0.000000
REDUCED COST
1.000000
0.000000
0.000000
0.000000
DUAL PRICES
2.000000
0.000000
0.000000
1.000000
N ITERATIONS = 3
RANGES IN WHICH THE BASIC IS UNCHANGED:
VARIABLE
X1
X2
X3
X4
ALLOWABLE
DECREASE
INFINITY
0.500000
0.000000
INFINITY
212
ROW
2)
3)
4)
5)
ALLOWABLE
DECREASE
30.000000
150.000000
150.000000
400.000000
213
Z
S1
S2
S3
Z
1
0
0
0
X1
-4
1
0
2
X2
-6
1
0
3
X3
-7
1
0
4
X4
-8
1
1
6
X5
-9
1
0
5
S1
0
1
0
0
S2
0
0
1
0
S3
0
0
0
1
S4
0
0
0
0
LD
0
950
400
4600
214
S4
5000
1
0
1 X
5
5
0
0
a B 1 a
1
0
1
0
0
959 4
1
1
1
2
1
X
5
0
1
Z
X2
S1
X4
X3
Z
1
0
0
0
0
X1
1
2
0
0
-1
X2
0
1
0
0
0
X3
0
0
0
0
1
X4
0
0
0
1
0
X5
-4
-5
0
1
5
S1
0
0
1
0
0
S2
0
1
1
1
-2
S3
0
-2
-1
-1
3
LD
6900
950
400
4600
5000
Dado que el valor del coeficiente de Losetas Npoles <0 (= - 4), esto indica que la
base no sera ptima; por lo tanto, podramos usar el algoritmo Simplex a fin de no
tener coeficientes negativos en el regln (o con lo que la variable de Losetas
Npoles entrara a la base, convirtindose en variables de decisin).
215
Z
1
0
0
0
0
Z
X2
X5
X4
S3
X1
0.5
1.5
-0.5
0
-0.5
X2
0
1
0
0
0
X3
0.5
0.5
0.5
0
0.5
X4
0
0
0
1
0
X5
0
0
1
0
0
S1
2.5
-0.5
1.5
-1
0.5
S2
1.5
-0.5
0.5
0
-0.5
S3
0
0
0
0
1
LD
7325
275
275
400
125
De estos resultados SOFA S.A. puede notar que se deberan producir 275 unidades
de las losetas Npoles y no producir losetas Sicilia; ya que as aumenta la utilidad a
$7325.00
8. SOFA S.A. est planeando dar un acabado especial a sus losetas , el cual
consiste en un fino acabado de oro. Es necesario entonces utilizar ms horas por
cada tipo de losetas : para el tipo Romana , se necesitan 4 horas; para el tipo
Esparta , 3 horas; para el tipo Sicilia, 5 horas; y para el tipo Atenas, 4 horas. La
disponibilidad mxima de horas de acabado es de 4000 Disminuir las utilidades
de la empresa aadiendo este tipo de acabado a las losetas?
Solucin:
Agregamos los valores de las restricciones de las horas de acabado de oro a la
Tabla ptima, con lo que aumentara una variable S 4 de holgura. Esta variable se
asume en el tablero como variable bsica, y se puede notar que las columnas a ij de
cada una de las variables bsicas restantes quedarn alteradas por el ingreso de esta
nueva restriccin. Este problema se solucionar usando el mtodo de
transformaciones de Gauss-Jordan
Tablero ptimo Alterado:
Z
X2
S1
X4
X3
S4
Z
1
0
0
0
0
0
X1
1
2
0
0
-1
4
X2
0
1
0
0
0
3
X3
0
0
0
0
1
5
X4
0
0
0
1
0
4
S1
0
0
1
0
0
0
S2
0
1
1
1
-2
0
S3
1
-2
-1
-1
3
0
S4
0
0
0
0
0
1
LD
6900
950
400
4600
5000
4000
216
Debemos transformar los valores en negrita, que son los que alteran el tablero
ptimo.
Luego de las transformaciones, obtenemos el siguiente tablero ptimo:
Z
1
0
0
0
0
0
Z
X2
S1
X4
X3
S4
X1
1
2
0
0
-1
3
X2
0
1
0
0
0
0
X3
0
0
0
0
1
0
X4
0
0
0
1
0
0
S1
0
0
1
0
0
0
S2
0
1
1
1
-2
3
S3
1
-2
-1
-1
3
-5
S4
0
0
0
0
0
1
LD
6900
300
150
550
100
400
P2
217
PRODUCTO 2
15 dlares
0.75 horas
1.50 horas
2 unidades
8 dlares
0.50 horas
0.80 horas
1 unidad
Precio de Venta
Trabajo requerido
Tiempo de mquina
requerido
Materia prima requerida
PROGRAMA EN LINDO
MAX 15P1+8P2-6OT-1.5RM-A1-A2
SUBJECT TO
2)
3)
4)
5)
6)
7)
8)
P1 - 10A1
<= 50
P2 - 15A2
<= 60
0.75P1 + 0.5P2 - OT <= 160
2P1 + P2 - RM <= 0
RM
<= 400
A1 + A2
<= 100
1.5P1+ 0.8P2
<= 320
END
LP OPTIMUM FOUND AT STEP
OBJECTIVE FUNCTION VALUE
1)
VARIABLE
P1
P2
OT
RM
A1
6900.000000
VALUE
160.000000
80.000000
0.000000
400.000000
11.000000
REDUCED COST
0.000000
0.000000
2.133333
0.000000
0.000000
218
A2
ROW
2)
3)
4)
5)
6)
7)
8)
1.333333
SLACK OR SURPLUS
0.000000
0.000000
0.000000
0.000000
0.000000
87.666664
16.000000
0.000000
DUAL PRICES
0.100000
0.066667
3.866667
6.000000
4.000000
0.000000
0.000000
N ITERATIONS = 5
VARIABLE
P1
P2
OT
RM
A1
A2
ROW
2)
3)
4)
ALLOWABLE
DECREASE
0.533333
0.483333
INFINITY
4.500000
5.333333
7.250000
ALLOWABLE
DECREASE
876.666626
1314.999878
2.500000
219
5)
6)
7)
8)
0.000000
400.000000
100.000000
320.000000
6.666667
6.666667
INFINITY
INFINITY
55.000000
55.000000
87.666667
16.000000
c) Wivco debera pagar por otra unidad de materia prima 4.5 dlares, este valor lo
obtenemos del precio dual de la fila (6), ya que esta fila est referida a las
unidades de materia prima.}
220
d) Con respecto a las horas mquinas, sabemos que Wivco utiliza 1.5 horas en el
producto 1 y 0.8 en el producto 2 (tabla 1), esto est referido en la fila 8,
entonces para saber si se necesitan ms horas mquinas tenemos que ver la
holgura (slack or surplus) de la fila 8, el valor que encontramos es de 16; esto
implica que estn sobrando 16 horas mquina, por lo que Wivco no pagara nada
por hora mquina.
e) Al aumentar las horas-hombre a 45 (45x4=180), encontramos un aumento de 20
horas(45x4-40x4=20), este aumento est en el rango permisible de aumento ; la
fila 4 se refiere a las horas-hombres , y el precio dual de esta fila es 3.866667 , al
calcular la nueva solucin ptima tendremos que esta ser : 2505 dlares, lo que
se puede obtener tanto del reporte LINDO, como de ejecutar:
Z nueva Z C1 2427.667 ( 20)(3.866667) 2505
Precio de Venta
Trabajo requerido
Tiempo de mquina requerido
Materia prima requerida
PRODUCTO 3
17 dlares
2 horas
2 horas
1 unidad
221
10. ABC, puede fabricar los productos A, B y C, para el cual requiere de los
componentes C1, C2 y C3, la cantidad de componentes por unidad de producto as
como el precio de venta unitario se presenta en la siguiente tabla:
PRODUCTO
A
B
C
CANTIDAD
DISPONIBLE
C1
C2
C3
1
3
4
2
2
5
6
2
1
140
130
80
PRECIO DE
VENTA
60
70
80
(BASIS)
ART
XC
SLK 3
XB
SLK 4
8.000
-0.600
1.400
0.800
XA
6.000
-3.200
8.800
4.600
XB
0.000
0.000
0.000
1.000
XC
0.000
1.000
0.000
0.000
SLK 2
18.000
1.400
-1.600
-0.200
SLK 3
0.000
0.000
1.000
0.000
3160.000
8.000
18.000
36.000
222
3160.000000
VARIABLE
XA
XB
XC
ROW
2)
3)
4)
VALUE
0.000000
36.000000
8.000000
SLACK OR SURPLUS
0.000000
18.000000
0.000000
REDUCED COST
6.000000
0.000000
0.000000
DUAL PRICES
18.000000
0.000000
8.000000
N ITERATIONS = 2
VARIABLE
XA
XB
XC
ROW
2)
3)
4)
ALLOWABLE
DECREASE
INFINITY
1.304348
45.000000
ALLOWABLE
DECREASE
20.000000
18.000000
12.857143
223
a) Encuentre el intervalo de los valores del precio de C para los cuales la base
actual permanece ptima.
Como e
Aumento Admisible
1.875
Disminucin Admisible
45
80 45 C C 80 1.875
35 C 81.875
Aumento Admisible
13.33
Disminucin Admisible
12.86
80 12.86 C 3 C 3 80 13.33
67.14 C 3 93.33
Aumento Admisible
224
11.25
Disminucin Admisible
1.304
MODELO
A
B
C
RECURSO
DISPONIBLE
MATERIALE
S
2
4
1
700
M. DE
OBRA
1.5
2
0.5
450
T DE MAQ.
P. DE VENTA
0.5
1.2
0.2
220
80
130
40
225
13825.000000
VARIABLE
XA
XB
XC
ROW
2)
3)
4)
5)
VALUE
150.000000
0.000000
400.000000
SLACK OR SURPLUS
0.000000
25.000000
65.000000
0.000000
REDUCED COST
0.000000
34.000000
0.000000
DUAL PRICES
20.500000
0.000000
0.000000
-3.500000
N ITERATIONS = 2
RANGES IN WHICH THE BASIC IS UNCHANGED:
VARIABLE
XA
XB
XC
ROW
2)
3)
4)
5)
ALLOWABLE
ALLOWABLE
INCREASE
DECREASE
37.500000
3.500000
INFINITY
48.000000
34.000000
INFINITY
20.500000
INFINITY
1.750000
OBJ COEFFICIENT RANGES
CURRENT COEF
ALLOWABLE
DECREASE
400.000000
5.000000
65.000000
50.000000
226
Se pide responder:
a) Cul deber ser el precio de B para hacer atractiva su fabricacin?
Precio de Vta. + Costo reducido
130 34 164
b) Calcule la holgura del tiempo de mquina
0.5XA+1.2XB+0.2XC + H = 220
0.5*150+1.2*0+0.2*400 + H = 220
H = 65
c) Calcule la utilidad total si solo se produce 110 unidades de A.
Z = 13825 + (110 150)*(-3.5) = 13965
d) Calcule la utilidad total si se deseas aumentar en 10 unidades los
MATERIALES, considerando que se deber pagar 20 soles por cada
unidad adicional
Z = 13825 + (20.5 (20 - 15))*10 = 13980
e) Exprese en el modelo los cambios que se deber hacer para incluir lo
solicitado en (d).
MAX
ST
227
228
Se llaman as a los problemas en la cual todas las variables tienen que ser nmeros
enteros.
Ejemplo:
Min Z = 5 X1 + 3 X2
Sujeto a:
5 X1 + 3 X2 6
X1, X2 0, X1, X2 enteros
Se llama as a los problemas en los cuales todas las variables deben ser iguales a 0
1.
Ejemplo:
Min Z = 5 X1 + 3 X2
Sujeto a:
229
5 X1 + 3 X2 6
X1, X2 = 0 1
RESTRICCIONES O BIEN
Frecuentemente se dar la situacin en que se dan restricciones de la forma:
f ( x1 , x 2 ,..., x n ) 0
g ( x1 , x 2 ,..., x n ) 0
f ( x1 , x 2 ,..., x n ) My
Donde:
g ( x1 , x 2 ,..., x n ) M (1 y )
x + y 3, 2 x + 5 y
Solucin:
Escribiendo las restricciones anteriores de la forma () y () se tiene:
230
x+ y 30
2 x + 5 y 12 0
...(1)
...(2)
...(1.1)
...(2.1)
Donde:
yi : es una variable 0 1 y,
M : es un valor suficientemente grande que se escoge para asegurar que
satisfagan las Ec. (1.1) y (2.1)
se
...(1.2)
2 x + 5 y 12 3 (1-yi)
...(2.2)
RESTRICCIONES SI ENTONCES
231
f ( x1 , x2 ,..., xn ) 0 g ( x1 , x2 ,..., xn ) 0
Para lograr esto, es necesario incluir las restricciones siguientes a la formulacin:
g ( x1 , x 2 ,..., x n ) My
f ( x1 , x 2 ,..., x n ) M (1 y )
Donde:
Y: es una variable 0-1
M: es un nmero positivo suficientemente grande para que se cumplan las dos
restricciones anteriores, para todos los valores de x i que satisfacen las otras
restricciones del problema.
EJEMPLOS
1. El gobierno peruano, dentro de sus planes de apoyo y fomento del sector agrario,
est considerando 4 proyectos de irrigacin. El proyecto 1, irrigar aproximadamente
20000 hectreas de terreno, el proyecto 2, 27500 hectreas, el proyecto 3, 15000
hectreas y el proyecto 4, 10000 hectreas. La ejecucin y puesta en marcha del
proyecto 1 tiene un costo de 10 millones de dlares, el proyecto 2, 14 millones de
dlares, el proyecto 3, 8 millones y el proyecto 4, 6 millones. Si adems se sabe que
el presupuesto para la ejecucin de proyectos de inversin en el sector agrario es de
28 millones de dlares, formule un P.E., cuya solucin ayude al gobierno a
maximizar el nmero de hectreas irrigadas.
Solucin:
Como el gobierno puede tomar solo dos decisiones, respecto a cada proyecto de
inversin, empezamos definiendo una variable 0-1.
232
...(2)
X2 + X3 + X4 1
233
3
2
2
1
1
3
3
1
3
3
3
1
1
3
2
3
1
2
DEFENSA
3
2
2
1
2
3
234
A-D
Solucin:
De la tabla 1 se puede observar que los jugadores juegan en las posiciones
siguientes:
Posicin
D
A
C
Jugadores
3,4,5,6,7
1,3,5,7
2,4,6
Sea:
yi = 1,0 : Inicia o no el juego el jugador i (i = 1,2,7)
Max
Z = 3 y 1 + 2 y 2 + 2 y3 + y4 + 2 y 5 + 3 y6 + y 7
Sujeto a:
y 1 + y2 + y3 + y4 + y5 + y6 + y7 = 5
3 y1 + 2 y2 + 2 y3 + y4 + y5 + 3 y6 + 3 y7 10
3 y1 + y2 + 3 y3 + 3 y4 + 3 y5 + y6 + 2 y7 10
y1 + 3 y2 + 2 y3 + 3 y4 + y5 + 2 y6 + 2 y7 10
y 3 + y6 1
y1 - y4 0
y 1 - y5 0
y 2 + y3 1
y 3 + y4 + y5 + y6 + y7 4
y 1 + y3 + y5 + y7 2
y 2 + y4 + y5 1
3. Debido a la contaminacin excesiva del ro Mommis, el estado de Mommis
construir algunas estaciones para el control de la contaminacin. Se est
considerando 3 lugares (Lugares 1,2 y 3). A Mommis le interesa controlar los
niveles de contaminacin de dos contaminantes (Contaminantes 1 y 2), la
legislacin del estado requiere que se eliminen por lo menos 80000 toneladas del
contaminante 1 y por lo menos 50000 toneladas del contaminante 2 del ro. En la
Tabla 2 se encuentran los datos relevantes para este problema.
235
Lugar 1
Lugar 2
Lugar 3
COSTO DE
CONSTRUC.
DE UNA
ESTACION
(dlares)
100000
60000
40000
COSTO DEL
CANTIDAD REMOVIDA
TRATAM. DE
POR TONELADA DE
1 TON. DE
AGUA
AGUA
(dlares)
Contam. 1 Contam. 2
20
0.40 Ton
0.30 Ton
30
0.25 Ton
0.20 Ton
40
0.20 Ton
0.25 Ton
Solucin:
Sea:
yi = 1,0 Se construye o no una estacin en el lugar i (i = 1,2,3)
Xi = N de toneladas de agua tratadas en la estacin i (i = 1,2,3)
Min Z = 100000 y1 + 60000 y2 + 40000 y3 + 20 X1 + 30 X2 + 40 X3
Sujeto a:
0.40X1 + 0.25 X2 + 0.20 X3 80000
0.30 X1 + 0.20 X2 + 0.25 X3 50000
X 1 M 1 y1
X 2 M 2 y2
X 3 M 3 y3
X1, X2, X3 0
Donde:
M1 = Max (80000/0.40 , 50000/0.30) = 200000
M2 = Max (80000/0.25, 50000/0.20) = 320000
M3 = Max (80000/0.20, 50000/0.25) = 400000
236
Calc.
IO
Estr.
Dat
(1)
X
(2)
X
X
(3)
X
MATERIAS
Matemticas
Inv. Oper.
Computacin
Pre-Requisito Ning Ning
X
(6)
CURSOS (i)
Estad
Sim.
adm
Comp
(4)
X
X
(1)
Intr
Prog
(5)
(6)
X
X
(6)
X
Ning
Pred
Total
Reque
r
(7)
X
2
X
2
2
(4)
Sea:
yi = 1,0 : Se estudia o no el curso i (i = 1,2,,7)
237
Min
Z = y 1 + y2 + y3 + y4 + y5 + y6 + y7
Sujeto a:
y1 + y2 + y3 + y4+ y7 2
y2 + y4 + y5+ y7 2
y3 + y5 + y6 2
- y 1 + y4 0
y3 - y6 0
y5 - y6 0
- y 4 + y7 0
5. Una compaa considera la apertura de almacenes en cuatro ciudades: Nueva
York, Los Angeles, Chicago y Atlanta. Cada almacn puede enviar 100 unidades a
la semana. El costo semana fijo para mantener abierto cada almacn es de 400
dlares en Nueva York, de 500 dlares en Los Angeles, de 300 dlares en Chicago,
y de 150 dlares en Atlanta. La regin 1 del pas requiere semanalmente 80
unidades; la regin 2,70 unidades y la regin 3, 40 unidades. En la tabla 3 se
muestran los costos (incluyendo los costos de produccin y de envo) para enviar 1
unidad de la fbrica a una regin. Se desea satisfacer las demandas semanales a un
costo mnimo, sujetas a la informacin anterior y a las restricciones siguientes:
1. Si se abre el almacn en Nueva York, entonces hay que abrir el almacn
en Los ngeles.
2.
3.
20
48
26
40
15
35
50
26
18
238
Atlanta
24
50
35
Solucin:
La tabla 3 se puede re-escribir de la manera siguiente:
Regin 1
(dlares)
HACIA
Regin 2
(dlares)
Regin 3
(dlares)
Envo
(sem.)
Costo Fijo
(semanal)
20
40
50
100
400
48
15
26
100
500
Chicago (3)
26
35
18
100
300
Atlanta
(4)
24
50
35
100
150
Requer./sem.
80
70
40
Ciudad (i)
Sea:
yi = 1,0 Se abre o no un almacn en la ciudad i (i = 1,2,3,4)
Xij = N de unid. Enviadas semanalmente de la ciudad i a la regin j (i = 1,2,3,4; j =
1,2,3)
Min
Sujeto a:
X11 + X12 + X13 100 y1
239
Archivo 1
1
X
Archivo 2
Archivo 3
Archivo 4
Archivo 5
2
X
4
X
8
X
9
X
10
X
X
X
DISCO
5
6
X
X
X
X
X
240
Solucin:
La Capacidad de los Discos la podemos tabular de la manera siguiente:
Capacidad
(K)
DISCO
5
6
1
10
a) Sea:
yi = 1,0 Se selecciona o no el disco i (i = 1,2,,10)
Min Z = 3 y1 + 5 y2 + y3 + 2 y4 + y5 + 4 y6 + 3 y7 + y8 + 2 y9 + 2 y10
Sujeto a:
y1 + y3 + y4 + y5 + y8 + y9 1
y1 + y3 1
y2 + y5 + y7 + y10 1
y3 + y6+ y8 1
y1 + y2 + y4 + y6 + y7 + y9 + y10 1
b) Para esta parte habr que agregar a la formulacin anterior la siguiente
condicin:
Si y3 + y5 > 0 Entonces y2 1
Esta condicin lgica puede escribirse:
1 - y2 M y
y3 + y5 M (1-y)
Donde el mnimo valor que puede tomar M es: M = 2
Por lo tanto:
1 - y2 2 y
y3 + y5 2 (1-y)
241
COSTO
UTILIDAD
(MILES $)
(MILES $)
(MILES $)
departamento
48
40
casas
46
40
Sea:
yj = 1,0 Se construye o no el proyecto recreacional j (j = 1,2)
Xi = El nmero de viviendas de tipo i en miles (i = 1,2: Casa, Dpto.)
Max Z = 8 X2 + 6 X1 - 28000 y1 - 12000 y2
Sujeto a:
X1+ X2 10
(1)
y1 + y2 = 1
(2)
La condicin:
Si se construye una drsena, el nmero de casas en el proyecto tendr que ser por
lo menos el triple del nmero de departamentos, puede escribirse.
Si y2 > 0 Entonces
242
X1 3 X2
3 X2 - X1 M y
y2 M (1-y)
Donde M puede tomar como valor mnimo: M = 3 (10000)
Por lo tanto:
3 X2 - X1 30000 y
y2 30000 (1-y)
(3)
(4)
Obs:
Las restricciones (3) y (4), pueden escribirse como:
3 X2 - X1 30000 (1-y2)
Adems
Teniendo en cuenta que: y1 + y2 = 1, tambin se pueden escribir como:
3 X2 - X1 30000 y1
243
Tabla 5
FF
7
8
4
5
SE
EN
LO
MO
R
2
4
8
4
PF
6
5
2
7
B
5
3
11
5
Solucin:
La tabla 5 la podemos escribir de la manera siguiente:
REGION
FF
PF
SE
EN
LO
MO
Prom. Cheques/da
C. de mantener una
Cta. /ao
X11
X12
X13
X14
90000
X21
X22
X23
X24
90000
X31
X32
X33
X34
90000
X41
X42
X43
X44
90000
50000
40000
30000
20000
Z = 0.15(2555 X11 + 2920 X12 + 1460 X13 + 1825 X14 + 730 X21 + 1460 X22
+ 2920 X23 + 1460 X24 + 2190 X31 + 1825 X32 + 730 X33 + 2555 X34 + 1825 X41 +
1095 X42 + 4015 X43 + 1825 X4 ) 50000 y1 40000y2 30000 y3 20000 y4
Sujeto a:
X11 + X21 + X31 + X41 = 40000
244
PROBLEMAS RESUELTOS
1. Una compaa produce dos productos A y B. Cada unidad de producto A
requiere una hora de servicios de ingeniera y 5 horas de tiempo mquina. Producir
una unidad de producto B requiere 2 horas de servicios de ingeniera y 8 horas de
maquina disponible. Hay 100 horas de ingeniera y 400 horas de tiempo de
maquina disponible. El costo de produccin es una funcin no lineal de la cantidad
producida tal como se da en la tabla 7
Tabla 7
PRODUCTO A
Produccin
Utilidad
PRODUCTO B
Produccin
Utilidad
245
(unidades)
049
50 100
10
8
(unidades)
039
40 100
7
3
Solucin:
Sea:
Xij
j=1,2)
Yi 1,0 ; (i = 1,2)
MAX Z 10 X A1 8 X A 2 7 X B1 3 X B 2
Sujeto a:
(XA1 + XA2) + 2 (XB1 + XB2) 100
5 (XA1 + XA2) + 8 (XB1 + XB2) 400
XA1 50
XA1 50 y1
XA2 50 y1
XB1 40
XB1 40 y2
XB2 60 y2
XA1 , XA2 , XB1 , XB2 0
2. Un urbanizador de bienes races est estudiando varios proyectos estrechamente
interrelacionados. Algunos proyectos solo se pueden llevar a cabo si se cumplen
ciertas condiciones (Tabla 8). Sea R1 la utilidad total de la inversin i y C1 el
costo de hacerlo. Desea maximizar la utilidad total al invertir hasta M dlares.
Formule el problema como un PE. Defina sus variables de decisin.
Tabla 8
PROYECTO
A
CONDICION
Ninguna
246
B
C
D
E
F
G
No si C y solo si E
No si B
Solo si A
No si F y solo si C
No si E y solo si C
Solo si A y B
Solucin:
Sea:
Xi = 1,0 : Se realiza o no el proyecto i (i = A,B,,G)
Max Z = RA XA + + RG XG - (CA XA +.......+ CG XG)
Sujeto a:
XB + X C 1
XE - XB 0
XC XF 1
XA - X D 0
XE + XF 1
XC - XE 0
2 XG - XA - XB 0
XA + + XG 0
3. Un problema que afronta todos los das un electricista consiste en decidir que
generadores conectar. El electricista en cuestin tiene tres generadores con las
caractersticas que se muestran en la Tabla 9. Hay dos periodos del da. En el
primero se necesitan 2900 megawatts. En el segundo 3900 megawatts. Un
generador que se conecte para el primer periodo puede ser usado en el segundo sin
causar un nuevo gasto de conexin. Todos los generadores principales (como son
A, B y C de la tabla 9) son apagados al trmino del da. Formule este problema
como un PLE.
Tabla 9
Generador Costo fijo de
Capacidad mxima
247
A
B
Costos
conexin
3000
2000
1000
5
4
7
en cada periodo en
(MW)
2100
1800
3000
Solucin:
Sea:
Xij = MW utilizados por el generador i en el periodo j (i = A,B,C; j = 1,2)
yi = 1,0 : Se utiliza o no el generador i (i = A,B,C)
Min
Z =
(XC1 + XC2)
Sujeto a:
XA1 + XB1 + XC1 2900
XA2 + XB2 + XC2 3900
XA1 2100 yA
XA2 2100 yA
XB1 1800 yB
XB2 1800 yB
XC1 3000 yC
XC2 3000 yC
XA1 , XA2 , XB1 , XB2 , XC1 , XC2 0
4. La junta de directores de una empresa manufacturera esta estudiando un
conjunto de inversiones sujetas a las siguientes condiciones:
INVERSION
CONDICION
248
1
2
3
4
5
6
7
Ninguna
Solo s 1
Solo s 2
Se har s 1 y 2
No s 1 2
No s 2 3
Solo s 2 y no 3
Sujeto a:
C.2: X2 - X1 0
C.3: X3 - X2 0
C.4: X4 - X1 0
X4 - X 2 0
C.5: X1 - X2 + X5 1
- X1 + X2 + X5 1
C.6: X2 + X3 + X6 2
C.7: X2 + X3 + X7 2
X2 - X7 0
X1 , X2 , X3 , X4 , X5 , X6 , X7 0
INV.2
X2
1
1
X1
1
0
Solo s 1
X 2 X1 0
249
INV.3
INV. 4
INV. 5
INV. 6
0
0
1
0
X3
1
1
0
0
X2
1
0
1
0
X4
1
X1
1
X2
1
1
1
0
0
0
0
X1
1
1
1
1
0
0
0
0
0
0
0
0
1
1
X2
1
1
0
0
0
0
1
1
1
0
0
1
0
1
X5
1
0
1
0
0
1
0
1
X6
1
1
1
1
0
0
X2
1
1
0
0
0
0
X3
1
0
1
0
0
1
Solo s 2
X3 X2 0
Se har s 1 y 2
No s 1 2
X2 + X3 + X6 2
No s 2 y 3
X2 + X3 + X6 2
250
INV. 7
0
0
1
1
0
1
X7
1
1
1
1
0
0
0
X2
1
1
0
0
0
0
1
X3
1
0
1
0
0
1
0
Solo s 2 y no 3
X7 - X2 0
X3 + X7 1
votantes independientes)
El gobernador Blue es demcrata. Suponga que cada elector siempre vota por su
propio partido. Formule un PE para ayudar al gobernador Blue a maximizar el
nmero de demcratas que ganarn una silla en el congreso.
Tabla 6
Ciudad 1
Ciudad 2
Ciudad 3
Ciudad 4
Ciudad 5
Ciudad 6
REPUBLICANOS
80
60
40
20
40
40
DEMOCRATAS
34
44
44
24
114
64
251
Ciudad 7
Ciudad 8
Ciudad 9
Ciudad 10
70
50
70
70
14
44
54
64
Solucin:
Sea:
yij = 1,0 : Pertenece o no la ciudad i al distrito congresional j (i = 1,2,,10;
j = 1,2,,5)
Como cada distrito tiene entre 150000 y 250000 votantes se tienen:
114 y11 + 104 y21 + 84 y31 +.+ 134 y101 150000
114 y11 + 104 y21 + 84 y31 +.+ 134 y101 250000
114 y12 + 104 y22 + 84 y32 +.+ 134 y102 150000
114 y12 + 104 y22 + 84 y32 +.+ 134 y102 250000
114 y13 + 104 y23 + 84 y33 +.+ 134 y103 150000
114 y13 + 104 y23 + 84 y33 +.+ 134 y103 250000
114 y14 + 104 y24 + 84 y34 +.+ 134 y104 150000
114 y14 + 104 y24 + 84 y34 +.+ 134 y104 250000
114 y15 + 104 y25 + 84 y35 +.+ 134 y105 150000
114 y15 + 104 y25 + 84 y35 +.+ 134 y105 250000
y11 + y12 + y13 + y14 + y15 = 1
y21 + y22 + y23 + y24 + y25 = 1
y31 + y32 + y33 + y34 + y35 = 1
..............................................
..............................................
y91 + y92 + y93 + y94 + y95 = 1
y101 + y102 + y103 + y104 + y105 = 1
252
Ciudad 1
Ciudad 2
Ciudad 3
Ciudad 4
Ciudad 5
Ciudad 6
Ciudad 7
Ciudad 8
Ciudad 9
Ciudad 10
REPUBLICANOS
(R )
80
60
40
20
40
40
70
50
70
70
DEMOCRATAS
(D)
34
44
44
24
114
64
14
44
54
64
(R ) (D)
46
16
-4
-4
-74
-24
56
6
26
6
253
X 5 M y5
Finalmente el PE se puede escribir como:
Sea:
yij = 1,0 : Pertenece o no la ciudad i al distrito congresional j (i = 1,2,
,10; j = 1,2,,5)
Xj = 1,0: Gana o no el partido demcrata en el distrito congresional j
yj = 1,0 : Se cumple o no la condicin j
Max
Z = X1 + X2 + X3 + X4 + X5
Sujeto a:
114 y11 + 104 y21 + 84 y31 +.+ 134 y101 150000
114 y11 + 104 y21 + 84 y31 +.+ 134 y101 250000
114 y12 + 104 y22 + 84 y32 +.+ 134 y102 150000
114 y12 + 104 y22 + 84 y32 +.+ 134 y102 250000
254
X 1 M y1
46 y12 + 16 y22 - 4 y32 + + 6 y102 M (1-y2)
X 2 M y2
........
........
........
........
46 y15 + 16 y25 - 4 y35 ++ 6 y105 M (1-y5)
X 5 M y5
255
Distrito 1
Distrito 2
Distrito 3
Distrito 4
Distrito 5
BLANCOS
80
70
90
50
60
NEGROS
30
5
10
40
30
La distancia, en millas, que debe viajar un alumno de cada distrito para llegar a
cada secundaria se muestra en la siguiente tabla:
Distrito 1
Distrito 2
Distrito 3
Distrito 4
Distrito 5
ESC.SEC 1 ESC.SEC 2
1
2
0.5
1.7
0.8
0.8
1.3
0.4
1.5
0.6
256
BNESC(BN,ESC):;
DISESC(DIS,ESC):;
ALUMNO(BN,DIS,ESC):Y,X,DISTA,NUME;
ENDSETS
DATA:
DISTA=1,2,.5,1.7,.8,.8,1.3,.4,1.5,.6,
1,2,.5,1.7,.8,.8,1.3,.4,1.5,.6;
NUME=80,80,70,70,90,90,50,50,60,60,
30,30,5,5,10,10,40,40,30,30;
ENDDATA
MIN=@SUM(ALUMNO:X*DISTA);
@FOR(ESC(K):@SUM(BNDIS(I,J):X(I,J,K))>=150);
@FOR(ESC(K):@SUM(BNDIS(I,J)|
I#EQ#2:X(I,J,K))*4>=@SUM(BNDIS(I,J)|I#EQ#1:X(I,J,K)));
@FOR(ALUMNO:X=NUME*Y);
@FOR(DISESC(J,K):Y(1,J,K)=Y(2,J,K));
@FOR(BNDIS(I,J):Y(I,J,1)+Y(I,J,2)=1);
@FOR(ALUMNO:@BIN(Y));
END
2. State University debe comprar 1100 computadoras de tres vendedores. El
vendedor 1 cobra $500 por computadora, ms un costo de transporte de $5000.
El vendedor 2 cobra $350 ms $4000, el vendedor 3 cobra $250 ms $6000. El
vendedor 1 vender a lo ms 500 computadoras a la universidad, el vendedor 2
vender a lo ms 900 computadoras, y el vendedor 3 vender a lo ms 400
computadoras. Formule una PE para minimizar el costo de la compra de las
computadoras necesarias.
EN LINGO
257
DATA:
CV=500,350,250;
CF=5000,4000,6000;
LIM=500,900,400;
ENDDATA
MIN=@SUM(F:CV*X+CF*Y);
@FOR(F(I):X(I)<=LIM(I)*Y(I));
@SUM(F:X)=1100;
@FOR(F:@BIN(Y));
END
ALGEBRAICAMENTE:
MIN
5000Y(1)+4000Y(2)+6000Y(3)+500X(1)+350X(2)+250X(3)
SUBJECT TO
2]- 500 Y( 1) + X( 1) <=
0
3]- 900 Y( 2) + X( 2) <=
0
4]- 400 Y( 3) + X( 3) <=
0
5] X( 1) + X( 2) + X( 3) =
END
INTE
3
1100
258
AO 1
AO 2
AO 3
AO 4
AO 5
PLANTA 1
PLANTA 2
PLANTA 3
PLANTA 4
Costo
Operacin
1.5
0.8
1.3
0.6
Costo
Reapertura
1.9
1.5
1.6
1.1
Costo por
cierre
1.7
1.2
1.3
0.8
i=Planta j=ao;
SETS:
PLA/1..4/:CAPL;
ANO/1..5/:CAPAC;
FUN(PLA,ANO):Y,X,Z,COSOP,COSCIE,COSABR;
ENDSETS
DATA:
CAPL=70,50,60,40;
CAPAC=80,100,120,140,160;
COSOP=1.5,1.5,1.5,1.5,1.5,.8,.8,.8,.8,.8,1.3,1.3,1.3,1.
3,1.3,.6,.6,.6,.6,.6;
COSCIE=1.7,1.7,1.7,1.7,1.7,1.2,1.2,1.2,1.2,1.2,1.3,1.3,
1.3,1.3,1.3,.8,.8,.8,.8,.8;
COSABR=1.9,1.9,1.9,1.9,1.9,1.5,1.5,1.5,1.5,1.5,1.6,1.6,
1.6,1.6,1.6,1.1,1.1,1.1,1.1,1.1;
ENDDATA
MIN=@SUM(FUN:X*COSOP+Y*COSCIE+Z*COSABR);
@FOR(ANO(J):@SUM(FUN(I,J):X(I,J)*CAPL(I))>=CAPAC);
@FOR(FUN(I,J)|J#EQ#1:X(I,J)=1);
@FOR(FUN(I,J)|J#EQ#1:Z(I,J)=0);
259
@FOR(FUN(I,J)|J#NE#5:X(I,J)-Y(I,J)=X(I,J+1));
@FOR(FUN(I,J):Y(I,J)<=X(I,J));
@FOR(FUN(I,J)|J#NE#1:X(I,J)-Z(I,J)=X(I,J-1));
@FOR(FUN(I,J):Z(I,J)<=X(I,J));
@FOR(FUN:@BIN(Y));
@FOR(FUN:@BIN(X));
@FOR(FUN:@BIN(Z));
END
3. Se disponen de cuatro camiones para entregar leche a cinco tiendas de
comestibles. En la siguiente tabla se muestran la capacidad y el costo diario de
operacin de cada camin.
Nota: He supuesto que ese costo diario es por CADA operacin, es decir, si el
camin 1 va a dos tiendas, el costo de operacin del camin 1 ser 45+45=90.
Se puede satisfacer la demanda de una tienda de comestibles mediante un solo
camin pero un mismo camin puede entregar leche a ms de una tienda. La
demanda diaria de cada una de las tiendas es la siguiente: tienda 1 demanda 100
galones, tienda 2=200galones, tienda3=300 galones, tienda 4=500 galones, tienda 5
demanda 800 galones. Formule una PE que se puede usar para minimizar el costo
diario para satisfacer las demandas de las cinco tiendas.
Camin 1
Camin 2
Camin 3
Camin 4
CAPACIDA
COSTO
D (galones) OPERACIN ($)
400
45
500
50
600
55
1100
60
! Xij=galones transportados
i=camin
j=tienda;
MODEL:
SETS:
CA/1..4/:CAP,OPE;
TI/1..5/:DEM;
260
MILK(CA,TI):X,Y,COS;
ENDSETS
DATA:
CAP=400,500,600,1100;
COS=45,45,45,45,45,50,50,50,50,50,55,55,55,55,55,60,60,
60,60,60;
DEM=100,200,300,500,800;
ENDDATA
MIN=@SUM(MILK:COS*Y);
@FOR(TI(J):@SUM(CA(I):X(I,J))=DEM(J));
@FOR(MILK(I,J):X(I,J)<=CAP(I)*Y(I,J));
@FOR(CA(I):@SUM(TI(J):X(I,J))<=CAP(I));
@FOR(MILK:@BIN(Y));
END
SECUENCIACIN EN LINGO
Sea la siguiente matriz de tiempos:
MAQUINA
J
TRABAJO I
M1
A
3
B
7
C
5
M2
-3
4
M3
8
-3
3
--M2
261
262
@FOR(GOL(I,J)|J#EQ#1:T(3,1)+5<=T(2,1)+500*(1-Y(3)));
@FOR(GOL(I,J)|J#EQ#2:T(2,2)+3<=T(3,2)+500*Y(4));
@FOR(GOL(I,J)|J#EQ#2:T(3,2)+4<=T(2,2)+500*(1-Y(4)));
@FOR(GOL(I,J)|J#EQ#3:T(1,3)+8<=T(3,3)+500*Y(5));
@FOR(GOL(I,J)|J#EQ#3:T(3,3)+3<=T(1,3)+500*(1-Y(5)));
@FOR(B:@BIN(Y));
END
FORMULACIN ALGEBRAICA
MIN
TT( 1)
SUBJECT TO
2]- T( 1, 1) + T( 1, 3) <= - 8
5] T( 1, 1) - TT( 1) <= - 3
10] T( 1, 1) - TT( 1) <= - 3
14]- T( 2, 1) + T( 2, 2) <= - 3
17] T( 2, 1) - TT( 1) <= - 7
26] T( 3, 1) - T( 3, 3) <= - 5
29]- T( 3, 2) + T( 3, 3) <= - 3
32] T( 3, 2) - TT( 1) <= - 4
41]- 500 Y( 1) + T( 1, 1) - T( 2,
44] 500 Y( 1) - T( 1, 1) + T( 2,
47]- 500 Y( 2) + T( 1, 1) - T( 3,
50] 500 Y( 2) - T( 1, 1) + T( 3,
53]- 500 Y( 3) + T( 2, 1) - T( 3,
56] 500 Y( 3) - T( 2, 1) + T( 3,
59]- 500 Y( 4) + T( 2, 2) - T( 3,
62] 500 Y( 4) - T( 2, 2) + T( 3,
65]- 500 Y( 5) + T( 1, 3) - T( 3,
68] 500 Y( 5) - T( 1, 3) + T( 3,
END
INTE
INTE
INTE
INTE
INTE
Y(
Y(
Y(
Y(
Y(
1)
1)
1)
1)
1)
1)
2)
2)
3)
3)
<=
<=
<=
<=
<=
<=
<=
<=
<=
<=
- 3
493
- 3
495
- 7
495
- 3
496
- 8
497
1)
2)
3)
4)
5)
263
Variable
80
15.00000
2
Value
TT( 1)
15.00000
Y( 1)
1.000000
Y( 2)
1.000000
Y( 3)
1.000000
Y( 4)
0.0000000
Y( 5)
0.0000000
T( 1, 1)
12.00000
T( 1, 2)
0.0000000
T( 1, 3)
0.0000000
T( 2, 1)
5.000000
T( 2, 2)
0.0000000
T( 2, 3)
0.0000000
T( 3, 1)
0.0000000
T( 3, 2)
11.00000
T( 3, 3)
8.000000
264
0.0000000
0.0000000
0.0000000
0.0000000
0.0000000
0.0000000
0.0000000
0.0000000
D( 1, 1)
0.0000000
D( 1, 2)
0.0000000
D( 1, 3)
0.0000000
D( 2, 1)
0.0000000
D( 2, 2)
0.0000000
D( 2, 3)
0.0000000
D( 3, 1)
0.0000000
D( 3, 2)
0.0000000
D( 3, 3)
0.0000000
0.0000000
Como se puede ver, el tiempo mnimo de ejecucin de los trabajos es 15.
265
6.2 ANEXO
LENGUAJE DE MODELADO LINGO
Qu es LINGO?
LINGO es una herramienta simple para utilizar la potencialidad de la optimizacin
lineal y no lineal para formular problemas muy grandes de una manera concisa,
resolverlos y analizar su solucin. La optimizacin le ayuda a encontrar la
respuesta que representa la mejor solucin; obtiene la mayor utilidad, respuesta o
felicidad; o logra el menor costo, desperdicio o disconformidad. A menudo estos
problemas significan hacer el uso ms eficiente de sus recursos- incluyendo
dinero, tiempo, maquinaria, personal, inventario y mucho ms. Los problemas de
optimizacin se clasifican a menudo como lineales y no lineales, dependiendo si
las relaciones entre las variables son o no lineales.
LINGO es un lenguaje de modelado matemtico diseado para formular y
resolver problemas de programacin lineal, programacin entera y
programacin no lineal.
Lenguaje de modelado de LINGO
Sintaxis de LINGO
La sintaxis que se utiliza en este programa es muy sencilla. Para el nombre de las
Variables y otros identificadores se establece que pueden tener 32 caracteres como
Mximo, Deben comenzar con una letra seguido de letras, dgitos o _.
LINGO no distingue entre maysculas y minsculas.
Con respecto a las sentencias:
Todas las sentencias deben terminar en un punto y coma.
Para darle un nombre a la funcin objetivo o a las restricciones, estos se deben
colocar entre corchetes.
Para declarar la funcin objetivo debemos colocar las palabras reservadas MAX o
MIN, (aparecern resaltadas en azul) seguidas del signo = .
Los comentarios deben comenzar con un signo ! , los cuales aparecen
resaltados en verde. Al igual q las sentencias los comentarios finalizan con un
punto y coma.
266
267
J = 2:
J=3
268
269
Objective value:
160.0000
Variable
CAPACIDAD( F1)
CAPACIDAD( F2)
DEMANDA( C1)
DEMANDA( C2)
DEMANDA( C3)
C( F1, C1)
C( F1, C2)
C( F1, C3)
C( F2, C1)
C( F2, C2)
C( F2, C3)
X( F1, C1)
X( F1, C2)
X( F1, C3)
X( F2, C1)
X( F2, C2)
X( F2, C3)
Value
Reduced Cost
30.00000
0.0000000
20.00000
0.0000000
10.00000
0.0000000
25.00000
0.0000000
15.00000
0.0000000
2.000000
0.0000000
4.000000
0.0000000
6.000000
0.0000000
7.000000
0.0000000
10.00000
0.0000000
1.000000
0.0000000
5.000000
0.0000000
25.00000
0.0000000
0.0000000
10.00000
5.000000
0.0000000
0.0000000
1.000000
15.00000
0.0000000
@BIN
@FREE
@BND
270
271
@FOR(ITEMS: @BND(QL, Q, QU)): Fija los lmites de todas las variables del
atributo Q en QL y QU (A QL y QU deben habrsele asignado valores en la
seccin de datos)
OPERADORES LOGICOS
LINGO tiene nueve operadores lgicos:
NOT (no), EQ (igual), NE (no igual), GT (mayor que), GE (mayor igual), LT
(menor igual), LE (menor igual), AND (y) y OR (o) que se utilizan para
comparar valores, la forma de usar es: #operador#.
Polmero A
4
3
6
500
Base
6
9
2
1100
272
Solucin:
Producto
airtex
extendex
resistex
inventario
Polmero
A
4
3
6
500
ganancia
7
7
6
Sea:
Xi: la cantidad de ingredientes del PRODUCTO i(i=airtex,extendex,resistex)que se
puede usar.
3
FUNCION OBJETIVO:
X
i 1
* ganacia(i )
sea del
Producto airtex, extendex,
resistex.
Entonces: MAX Z = 7*X1 + 7*X2 + 6*X3
SUJETO A:
3
X
i 1
COMPROMISO(i )
X
i 1
, donde requerimiento (i , j)
es la cantidad
de producto de cada tipo
de ingrediente.
para J=1------4X 1 + 3X2 + 6X3 <= 500*16;
para J=2------2X 1 + 2X2 + 3X3 <= 425*16;
273
EXTENDEX
DATA:
COMPROMISO=1000,500,400;
GANANCIA=7,7,6;
INVENTARIO =500,425,650,1100;
X=4,2,4,6,
3,2,2,9,
6,3,5,2;
ENDDATA
!FUNCION OBJETIVO MAXIMIZAR LA UTILIDAD;
MAX=@SUM(PRODUCTO:GANANCIA*NIVEL);
!RESTRICCION DEL INVENTARIO;
@FOR(INGREDIENTE(I):@SUM(PRODUCTO(P):X(P,I)*NIVEL(P))<=INVEN
TARIO(I)*16);
!RESTRICCION DEL COMPROMISO;
@FOR(PRODUCTO:NIVEL>=COMPROMISO);
END
274
NUTRIENTE(mg / 100g)
Tiacina
Tiamina
Protena
Hierro
Vitna c
Grasa
Espagueti
5000
1.1
1.4
0.18
0.0
5000
Pavo
Papas
29300
5300
1.8
0.5
5.4
0.9
0.06
0.06
0.0
10
5000
7900
Espinacas
3000
2.2
0.5
0.07
28
300
Pastel
4000
1.2
0.60
0.15
3.0
14300
Hierro
Espagueti
5000
1.1
NUTRIENTE(mg/100g)
Tiacina Tiamin Vitna c
a
1.4
0.18
0
Pavo
Papas
Espinacas
29300
5300
3000
1.8
0.5
2.2
5.4
0.9
0.5
0.06
0.06
0.07
Pastel
Minimo
4000
63000
1.2
10
0.6
15
0.15
1
Grasa
Maximo
5000
300
0
10
28
5000
7900
300
300
200
100
3
50
14300
0
100
275
Sujeto a:
Restriccin de la cantidad de grasa total que debe haber :
5
X 6 X i * CGi 0 ; donde
i 1
X
i 1
potaje ya sea de
Espagueti, pavo, papas, espinacas y pastel
Por lo tanto:
Para i=1--X1<=300/100
Para i=2--X2<=300/100
Para i=3--X3<=200/100
Para i=4--X4<=100/100
Para i=5--X5<=100/100
Restriccin de la cantidad de nutriente por pataje que debe haber como
mnimo:
5
X
i 1
* nutriente j ,i MINIMO j
Donde:
NUTRIENTEj,i . j ,es la cantidad de nutriente ya sea de protena, hierro, tiacina,
tiamina, vitna C ,grasa por cada potaje i ya sea: Espagueti, pavo, papas ,espinacas,
pastel.
y MINIMOj es la cantidad de potaje como mximo que debe existir.
276
Por lo tanto:
Para i=1,2,3,4,5 Y j=1
X1*5000 + X2*29300 + X3*5300 + X4*3000 + X5*4000 >=63000
Para i=1,2,3,4,5 Y j=2
X1*1.1 + X2*1.8 + X3*0.5 + X4*2.2 + X5*1.2 >=10
Para i=1,2,3,4,5 Y j=3
X1*1.4 + X2*5.4 + X3*0.9 + X4*0.5 + X5*0.6 >=15
Para i=1,2,3,4,5 Y j=4
X1*0.18 + X2*0.06+ X3*0.06 + X4*0.07 + X5* 0.15 >=1
Para i=1,2,3,4,5 Y j=5
X1*5000 + X 2*5000 + X3*7900 + X4*300 + X5*14300 >=0
SETS:
POTAJE/ESPAGUETI PAVO PAPAS ESPINACAS PASTEL/:NIVEL,MAXIMO;
NUTRIENTE/PROTEINA
HIERRO
TIACINA
TIAMINA
VITAMC
GRASA/:MINIMO;
PONU(POTAJE,NUTRIENTE):REQ;
ENDSETS
DATA:
MAXIMO=300,300,200,100,100;
MINIMO=63000,10,15,1,50,0;
REQ=5000,1.1,1.4,0.18,0,5000,
29300,1.8,5.4,0.06,0,5000,
5300,0.5,0.9,0.06,10,7900,
3000,2.2,0.5,0.07,28,300,
4000,1.2,0.6,0.15,3,14300;
ENDDATA
MIN=GRASA;
GRASA=@SUM(POTAJE(I):REQ(I,6)*NIVEL(I));
@FOR(POTAJE(I):NIVEL(I)<=MAXIMO(I)/100);
@FOR(NUTRIENTE(J):@SUM(POTAJE(I):REQ(I,J)*NIVEL(I))>=MINIMO(J));
END
HACIENDO CORRER EL PROGRAMA CON LINGO 10.0
277
Objective value:
Total solver iterations:
54800.00
3
Variable
Value
Reduced Cost
NIVEL( ESPAGUETI)
3.000000
0.000000
NIVEL( PAVO)
2.833333
0.000000
NIVEL( PAPAS)
2.000000
0.000000
NIVEL( ESPINACAS)
1.000000
0.000000
NIVEL( PASTEL)
0.6666667
0.000000
A ASI SUCESIVAMENTE.
Jess Mara
15
Callao
17
San Luis
22
Los Olivos
12
ventanilla
30
Villa el salvador
25
chorrillos
21
278
SOLUCION:
El problema nos menciona q debemos determinar la distribucin ptima de
pescado o tambin dicho la cantidad de toneladas pescado q debe ir de cada
proveedor a cada sucursal de la empresa, este dato debe reflejarse en la funcin
objetivo (FO).
Xij: cantidad de pescado distribuido (en toneladas) por el proveedor i hacia destino
j donde: (i va de 1 a 3 y j va de 1 a 4)
En esta expresin i (por convencin representa las filas) representa a los
proveedores y j representa a las sucursales de la empresa CUATRO MARIAS
Ya tenemos el conjunto de variables con que vamos a trabajar pero eso no es
suficiente para determinar la FO, para nuestro caso nos pide determinar la
distribucin ptima para tener el menor costo posible debido a la distribucin del
pescado.
Proveedor
Ventanilla
V. salvador
chorrillos
Jess Maria
6
4
8
Callao
2
9
8
sucursal
San Luis
6
5
1
Los Olivos
7
3
5
Para armar la funcin objetivo necesitamos relacionar los datos costo unitario por
tonelada de pescado por nmero de toneladas pescados trasportados esto nos dara
es costo total de transporte.
En la tabla relacionamos primero al primer proveedor ventanilla (rojo) con la
primera sucursal Jess Mara (verde) este dato nos refleja el costo por tonelada de
trasporte desde ventanilla hacia Jess Mara (azul) lo cual nos da la relacin 6*X 11
este mismo paso es para los dems datos con lo cual tendramos:
3
MIN Z = 6*X11 +2*X12+ 6*X13+ 7*X14 +4*X21+ 9*X22 +5*X23 +3*X24 +8*X31
+8*X32 +1*X 335*X34
Determinando las restricciones:
279
proveedor
ventanilla
V. Salvador
Chorrillos
Demanda
Sucursal
Callao
S. Luis
2
6
9
5
8
1
17
22
Jess M.
6
4
8
15
Los Olivos
7
3
5
12
disponibilidad
30
25
21
A la tabla de costos hemos hecho unos aadidos los cuales son la disponibilidad y
la demanda, la que nos ayudaran a determinar las restricciones del problema
Restriccin de la demanda:
Para la determinacin de las restricciones se debe tener muy en cuenta la relacin
entre los datos, para el caso de la demanda ella est relacionada directamente con
los proveedores ya q refleja la cantidad de toneladas q requiere cada sucursal.
Entonces la restriccin de la demanda va estar basada segn sucursal j.
3
Para j=1-----
X
i 1
i ,1
disponibilidad (1)
15 (Jess Maria)
3
Para j=2-----
X
i 1
i,2
disponibilidad (2)
Para j=3-----
X
i 1
22
i ,3
(San Luis)
Para j=4-----
X
i 1
disponibilidad (3)
i,4
disponibilidad ( 4)
280
(Los Olivos)
30 (ventanilla)
2
i=2
i=3
X31+ X32 +X33 21 (Chorrillos)
Al igual que la demanda la disponibilidad est en funcin a la cantidad de
toneladas trasportadas as q no debe estar relacionada con ninguna dato que refleje
costos.
Problema N 3: (usando 3 subndices):
281
DISTRITO
ESTUDIANTES
DE GRUPOS
DE GRUPOS
NO
MINORITARIOS MINORITARIOS
1
2
3
50
50
100
200
250
150
La corte local a decidido que cada una de las dos escuelas de segunda
enseanza de la ciudad (Cooley y walt whitman) debe tener aproximadamente
(ms o menos 5%) el mismo porcentaje de estudiantes de minoras, que la
ciudad entera. En la tabla B se da las distancias entre los distritos escolares y
las escuelas. Cada escuela debe tener entre 300 y 500 estudiantes. Utilice la
programacin lineal para determinar la asignacin de los estudiantes a cada
escuela para minimizar la distancia total que tienen que viajar los estudiantes
para llegar a ella.
Tabla B:
WALT
COOLE WHITMA
DISTRITO
Y
N
1
1
2
2
2
1
3
1
1
Solucin:
Primero vamos a encontrar la funcin objetivo, la escuela busca minimizar la
distancia total recorrida por sus estudiantes desde su distrito a la escuela y cuantos
estudiantes son mayoras y minoritarios, entonces vamos a llamar a la variable
estudiantes i,j,k, donde i: estudiantes del distrito i (i=1,2,3) que pertenecen al
282
Las Restricciones:
La primera restriccin va ser con respecto a cantidad de alumnos de los dos tipos
minora y mayora en los distritos i, segn la encuesta realizada matemticamente
lo expresaramos as.
For i ( For j
283
!colegios ;
SETS:
DIST/1..3/:;
TIPO/1..2/:;
COLE/1..2/:;
DT(DIST,TIPO):CANT;
DC(DIST,COLE):DISTA;
ALUMNOS(DIST,TIPO,COLE):X;
ENDSETS
DATA:
CANT=50,200,50,250,100,150;
DISTA=1,2,2,1,1,1;
ENDDATA
MIN=@SUM(ALUMNOS:DISTA*X);
!ALUMNOS POR DISTRITO Y POR TIPO;
@FOR(DT(I,J):@SUM(ALUMNOS(I,J,K):X(I,J,K))=CANT(I,J));
@FOR(COLE(K):@SUM(ALUMNOS(I,J,K):X(I,J,K))>300);
@FOR(COLE(K):@SUM(ALUMNOS(I,J,K):X(I,J,K))<500);
@FOR(COLE(K):@SUM(ALUMNOS(I,J,K)|
J#EQ#1:X(I,J,K))>0.2*(@SUM(DT(I,J):X(I,J,K))));
@FOR(COLE(K):@SUM(ALUMNOS(I,J,K)|
J#EQ#1:X(I,J,K))<0.3*(@SUM(DT(I,J):X(I,J,K))));
END
EL MODELO ALGEBRAICO:
284
MODEL:
[_1] MIN= X_1_1_1 + 2 * X_1_1_2 + 2 * X_1_2_1
+ X_2_1_1 + 2 *
X_2_1_2 + 2 * X_2_2_1 + X_2_2_2 + X_3_1_1 + 2
+ 2 * X_3_2_1 +
X_3_2_2 ;
[_2] X_1_1_1 + X_1_1_2 = 50 ;
[_3] X_1_2_1 + X_1_2_2 = 200 ;
[_4] X_2_1_1 + X_2_1_2 = 50 ;
[_5] X_2_2_1 + X_2_2_2 = 250 ;
[_6] X_3_1_1 + X_3_1_2 = 100 ;
[_7] X_3_2_1 + X_3_2_2 = 150 ;
[_8] X_1_1_1 + X_1_2_1 + X_2_1_1 + X_2_2_1 +
X_3_2_1 <= 300;
[_9] X_1_1_2 + X_1_2_2 + X_2_1_2 + X_2_2_2 +
X_3_2_2 <= 500;
[_10] 0.8 * X_1_1_1 - 0.2 * X_1_2_1 + 0.8 *
0.2 * X_2_2_1 +
0.8 * X_3_1_1 - 0.2 * X_3_2_1 >= 0 ;
[_11] 0.8 * X_1_1_2 - 0.2 * X_1_2_2 + 0.8 *
0.2 * X_2_2_2 +
0.8 * X_3_1_2 - 0.2 * X_3_2_2 >= 0 ;
[_12] 0.7 * X_1_1_1 - 0.3 * X_1_2_1 + 0.7 *
0.3 * X_2_2_1 +
0.7 * X_3_1_1 - 0.3 * X_3_2_1 <= 0 ;
[_13] 0.7 * X_1_1_2 - 0.3 * X_1_2_2 + 0.7 *
0.3 * X_2_2_2 +
0.7 * X_3_1_2 - 0.3 * X_3_2_2 <= 0 ;
END
+ X_1_2_2
* X_3_1_2
X_3_1_1 +
X_3_1_2 +
X_2_1_1 X_2_1_2 X_2_1_1 X_2_1_2 -
285
LIBRERIA
DO
FAB
A
FAB
B
FAB
C
CUZCO
SUPERMERCAD
O
LIBRERIA
Prod
Prod
Prod
Prod
Prod
Prod
Prd
Prod
Prod
Prod
Prod
Prod
13
15
17
11
12
15
14
12
13
15
13
12
10
13
14
12
14
16
13
14
15
11
12
13
12
11
13
10
11
13
11
13
14
12
13
14
Capacidad de produccin
FAB
FAB
FAB C
286
PROD1
75
65
70
PROD2
60
70
80
PROD3
65
75
75
Demanda
TACNA
CUZCO
PROD1
73
67
PROD2
58
72
PROD3
67
74
TACNA
CUZCO
SUPER MERCADO
150
140
LIBRERA
130
150
Solucin:
Xi,j,k,l=cantidad de productos fabricados en la fabrica i(i=A,B,C),en la ciudad si
en TACNA, si y
287
MAX
CZ
P4
I A , J TC K SM L P1
X I , J , K , L * PRECIOI , J , K , L
, donde PRECIO es el
P4
K SM , L P1
I , J ,K ,L
Para I=A,J=P1:
XA,P1_SM,TC + XA,P1,SM,CZ + XA,P1,L,TC + XA,P1,L,CZ <= 75 ;
Para I=A,J=P2:
XA,P2,SM,TC + XA,P2,SM,CZ + XA,P2,L,TC + XA,P2,L,CZ <= 60 ;
Para I=A,J=P3:
XA,P3,SM,TC + XA,P3,SM,CZ + XA,P3,L,TC + XA,P3,L,CZ <= 65 ;
Para I=B,J=P1:
XB,P1,SM,TC + XB,P1,SM,CZ + XB,P1,L,TC + XB,P1,L,CZ <= 65 ;
Para I=B,J=P2:
XB,P2,SM,TC + XB,P2,SM,CZ + XB,P2,L,TC + XB,P2,L,CZ <= 70 ;
288
Para I=B,J=P3:
XB,P3,SM,TC + XB,P3,SM,CZ + XB,P3,L,TC + XB,P3,L,CZ <= 75 ;
Para I=C,J=P1:
XC,P1,SM,TC + XC,P1,SM,CZ + XC,P1,L,TC + XC,P1,L,CZ <= 70 ;
Para I=C,J=P2:
XC,P2,SM,TC + XC,P2,SM,CZ + XC,P2,L,TC + XC,P2,LCZ <= 80 ;
Para I=C,J=P3:
XC,P3,SM,TC + XC,P3,SM,CZ + XC,P3,L,TC + XC,P3,L,CZ <= 75
Demanda:
C
I A, K SM
I , J ,K ,L
289
P4
I A , L P1
I ,J ,K ,L
290
SETS:
! FABRICAS DONDE SE VA A PRODUCIR EL PAPEL;
FABRICAS/A B C/: ;
! PRODUCTOS A SER PRODUCIDO POR LAS FBRICAS;
PRODUCTOS/P1 P2 P3/: ;
! CENTRO DE DISTRIBUCION DE LOS PRODUCTOS;
CDIST/SM L/:;
! CIUDADES DONDE VAN A SER DISTRIBUIDOS LOS PRODUCTOS;
CIUDAD/TC CZ/:;
! REQUERIMIENTO DE PRODUCTOS PARA UNA FABRICA ,EN UNA
CIUDAD,EN UN DETERMINADO SUPERMERCADO;
FPCC (FBRICAS, PRODUCTOS, CDIST, CIUDAD): PRECIO, X;
! CAPACIDAD DE PRODUCCION DE UN TERMENINADOM PRODUCTO
POR FBRICA;
FABPRO (FBRICAS, PRODUCTOS): CAPACIDAD;
! DEMANDA DE PRODUCCION;
PROCIU(PRODUCTOS,CIUDAD):DEMANDA;
! CAPACIDAD DE LOS CENTROS DE DISTRIBUCION;
CDCIUDAD (CDIST, CIUDAD): CAPACCD;
ENDSETS
DATA:
CAPACIDAD= 75,60,65,
65,70,75,
70,80,75;
DEMANDA= 73,67,
58,72,
67,74;
CAPACCD= 150,140,
130,150;
PRECIO= 12,15,17,11,12,15,14,12,13,15,13,12,
10,13,14,12,14,16,13,14,15,11,12,13,
12,11,13,10,11,13,11,13,14,12,13,14;
ENDDATA
!FUNCION OBJETIVO ,MAXIMIZANDO LA UTILIDAD;
[OBJETIVO]MAX =@SUM(FPCC:PRECIO*X);
!RESTRICCION DE LA CAPACIDAD DE DISTRIBUCION ;
@FOR(FABPRO(I,J):@SUM(FPCC(I,J,K,L):X(I,J,K,L))<=CAPACIDAD(I,J));
!RESTRICCION DE LA DEMANDA;
@FOR(PROCIU(J,L):@SUM(FPCC(I,J,K,L):X(I,J,K,L))>=DEMANDA(J,L));
291
! RESTRICCION
DE LA CAPACIDAD DE LOS CENTROS DE
DISTRIBUCION;
@FOR(CDCIUDAD(K,L):@SUM(FPCC(I,J,K,L):X(I,J,K,L))<=CAPACCD(K,L))
;
END
EL MODELO ALGEBRAICO:
MODEL:
[OBJETIVO] MAX= 12 * X_A_P1_SM_TC + 15 * X_A_P1_SM_CZ + 17 *
X_A_P1_L_TC
+ 11 * X_A_P1_L_CZ + 12 * X_A_P2_SM_TC + 15 * X_A_P2_SM_CZ + 14 *
X_A_P2_L_TC + 12 * X_A_P2_L_CZ + 13 * X_A_P3_SM_TC + 15 *
X_A_P3_SM_CZ +
13 * X_A_P3_L_TC + 12 * X_A_P3_L_CZ + 10 * X_B_P1_SM_TC + 13 *
X_B_P1_SM_CZ + 14 * X_B_P1_L_TC + 12 * X_B_P1_L_CZ + 14 *
X_B_P2_SM_TC +
16 * X_B_P2_SM_CZ + 13 * X_B_P2_L_TC + 14 * X_B_P2_L_CZ + 15 *
X_B_P3_SM_TC + 11 * X_B_P3_SM_CZ + 12 * X_B_P3_L_TC + 13 *
X_B_P3_L_CZ +
12 * X_C_P1_SM_TC + 11 * X_C_P1_SM_CZ + 13 * X_C_P1_L_TC + 10 *
X_C_P1_L_CZ + 11 * X_C_P2_SM_TC + 13 * X_C_P2_SM_CZ + 11 *
X_C_P2_L_TC +
13 * X_C_P2_L_CZ + 14 * X_C_P3_SM_TC + 12 * X_C_P3_SM_CZ + 13 *
X_C_P3_L_TC + 14 * X_C_P3_L_CZ ;
[_2] X_A_P1_SM_TC + X_A_P1_SM_CZ + X_A_P1_L_TC + X_A_P1_L_CZ
<= 75 ;
[_3] X_A_P2_SM_TC + X_A_P2_SM_CZ + X_A_P2_L_TC + X_A_P2_L_CZ
<= 60 ;
[_4] X_A_P3_SM_TC + X_A_P3_SM_CZ + X_A_P3_L_TC + X_A_P3_L_CZ
<= 65 ;
[_5] X_B_P1_SM_TC + X_B_P1_SM_CZ + X_B_P1_L_TC + X_B_P1_L_CZ
<= 65 ;
[_6] X_B_P2_SM_TC + X_B_P2_SM_CZ + X_B_P2_L_TC + X_B_P2_L_CZ
<= 70 ;
[_7] X_B_P3_SM_TC + X_B_P3_SM_CZ + X_B_P3_L_TC + X_B_P3_L_CZ
<= 75 ;
292
293
X_C_P3_L_CZ <=
150 ;
END
8260.000
20
Variable
Value
Reduced Cost
PRECIO( A, P1, SM, TC)
12.00000
0.000000
PRECIO( A, P1, SM, CZ)
15.00000
0.000000
PRECIO( A, P1, L, TC)
17.00000
0.000000
PRECIO( A, P1, L, CZ)
11.00000
0.000000
PRECIO( A, P2, SM, TC)
12.00000
0.000000
PRECIO( A, P2, SM, CZ)
15.00000
0.000000
PRECIO( A, P2, L, TC)
14.00000
0.000000
PRECIO( A, P2, L, CZ)
12.00000
0.000000
PRECIO( A, P3, SM, TC)
13.00000
0.000000
A ASI SUCESIVAMENTE.
294
Componente
Cuerpo
Base
Armado
Costo de Adquisicin
( en minutos)
10
20
16
Costo de Produccin ( en
minutos)
8
20
16
Componente
Cuerpo
Base
Armado
Tiempo
Preparacin
2
5
4
de Tiempo
Molde
4
2
5
de Tiempo
de
Fabricacin
2
4
5
295
Disponibilidad
49200
49200
49200
Nombre
Costo
Coef1
Coef2
FO
Rest1
Rest2
x
296
DATA:
n=6;
m=3;
ENDDATA
SETS:
v1/1.. n/:c,x;
v2/1..m/:b,e;
m1(v2,v1):a,d;
ENDSETS
DATA:
c,a,d,b,e=@OLE('C:\Samples\coelhos.xls','custo','Rest1'
,'Rest2','Coef1','Coef2');
ENDDATA
MIN=FO;
FO=@SUM(v1(j):c(j)*x(j));
@FOR(v2(i):@sum(v1(j):a(i,j)* x(j))<=b(i));
@FOR(v2(i):@sum(v1(j):d(i,j)* x(j))>=e(i));
DATA:
297
@OLE('C:\Samples\coelhos.xls','x','FO')= x,FO;
ENDDATA
Observacin:
-Intenta poner la Hoja de clculo sin ruta te va salir una venta dicindote OPEN
FILE si no pones la ruta tiene que estar el archivo abierto, para exportar los
resultados.
Si no quieres que el archivo se abra inmediatamente despus de la compilacin
tienes que poner la ruta como en la figura, esta ruta puede ser cualquiera donde tu
decidas guardar tu archivo.
Objective value:
Total
234650.0
iterations:
solver
Reduced Cost
0.000000
0.000000
0.000000
Variable
Value
6.000000
3.000000
FO
234650.0
298
0.000000
1.000000
0.000000
0.000000
0.000000
X( 1)
4675.000
X( 2)
0.000000
X( 3)
6100.000
X( 4)
1425.000
X( 5)
6100.000
X( 6)
0.000000
3.500000
EL PROBLEMA DE LAS BARRAS DE CHOCOLATE
Los requerimientos para la produccin de 3 tipos de barras de chocolate as como
la cantidad de recursos y la utilidad de cada tipo se muestran en el siguiente
cuadro:
MATERIA
PRIMA
B1
B2
B3
CANTIDAD
DISPONIBLE
AZUCAR
50
CHOCOLATE
100
GANANCIA
UNITARIA
SOLUCION.
Xi = Cantidad de barras de tipo i a producir; (i = 1, 2,3)
FO
299
Sujeto a:
X1+ X2 + X3 < 50
2X1+ 3X2 + X3 < 100
SOLUCIN EN LINGO
SETS:
!BARRAS,PRODUCCION Y GANANCIA;
!INGREDIENTES (AZUCAR,CHOCOLATE) Y DISPONIBILIDAD;
!CANTIDAD USADA POR BARRA;
B/B1,B2,B3/:P,G;
IN/A,CH/: D;
CA(IN,B):USO;
ENDSETS
MAX= @SUM(B:P*G);
@FOR(IN(I):@SUM(B(J):USO(I,J)*P(J))<=D(I));
DATA:
G=3,7,5;
D=50,100;
USO= 1,1,1,
2,3,1;
ENDDATA
END
300
1
En EXCEL
Se tiene por ejemplo el conjunto de datos G = 3,7,5 del problema de las
barras de chocolate, se digita en cada casillero, prescindiendo de la coma,
un numero de acuerdo al orden establecido.
Con el ratn se marca las celdas del un conjunto de datos ( 3 7 5)
Con INSERTAR, NOMBRE, DEFINIR del men, se define el nombre del
conjunto (Ejemplo: G)
En LINGO
En la seccin DATA, para importar informacin, se escribe el nombre del
conjunto que se iguala a la funcin @OLE sealando la ruta donde se
ubican los datos y si se desea exportar resultados a un lugar predefinido en
EXCEL se escribe primero la funcin y esta se iguala al nombre de las
celdas definidas en EXCEL
SOLUCION COLOCANDO LOS DATOS EN EXCEL
SETS:
!BARRAS,PRODUCCION Y GANANCIA;
!INGREDIENTES (AZUCAR,CHOCOLATE) Y DISPONIBILIDAD;
!CANTIDAD USADA POR BARRA;
B/B1,B2,B3/:P,G;
IN/A,CH/: D;
CA(IN,B):USO;
ENDSETS
MAX= @SUM(B:P*G);
@FOR(IN(I):@SUM(B(J):USO(I,J)*P(J))<=D(I));
DATA:
G=@OLE('C:\MIS DOCUMENTOS\BARRAS.XLS');
D=@OLE('C:\MIS DOCUMENTOS\BARRAS.XLS');
USO= @OLE('C:\MIS DOCUMENTOS\BARRAS.XLS');
!RESPUESTA DE PRODUCCION DE BARRAS;
@OLE('C:\MIS DOCUMENTOS\BARRAS.XLS')=P;
301
ENDDATA
END
2)
Los pasos que se sigue para crear una base de datos en ACCESS son:
Crear una carpeta en C, por ejemplo: USO DEL LINGO CON ACCESS.
Para crear la primera tabla B se procede como sigue: Crear una tabla
utilizando el asistente/Nuevo/Vista Diseo/Aceptar/en Nombre del
302
Tipo de Datos
Texto
Texto
numrico
IN
II
A
CH
D
50
100
CA
IN
B
A B1
A B2
A B3
CH B1
CH B2
CH B3
USO
1
1
1
2
3
1
303
304