Introduccion A La Modelacion Matematica y Optimizacion 14ENE2011 PDF
Introduccion A La Modelacion Matematica y Optimizacion 14ENE2011 PDF
Introduccion A La Modelacion Matematica y Optimizacion 14ENE2011 PDF
MODELACIÓN MATEMÁTICA
Y OPTIMIZACIÓN
DEDICATORIA
Contenido
SINOPSIS
La presente publicación trata de familiarizar al estudiante con los principales aspectos
de la Investigación de Operaciones (IO), incluyendo su origen histórico, su naturaleza y sus
diversas técnicas. El énfasis principal se da a los modelos de programación lineal, incluyendo
la programación lineal entera y el problema del transporte. Posteriormente, se presenta una
introducción a los modelos de redes. Así, estas notas pueden utilizarse en cursos de
Investigación de Operaciones I, Modelación Matemática e Introducción a la Optimización,
tanto para pregrado como para postgrado.
Se agradece sinceramente a todas las personas que con su esfuerzo como monitores
hicieron posible la recopilación de estas notas de clase: A Héctor Toro, Ana María Pérez, Juan
Carlos Lozano, Jhonatan Arias y a todas las personas que de una u otra forma han colaborado
con la producción de estas notas.
empresa. En los Estados Unidos, el interés de la industria por las técnicas de la IO no fue tan
inmediato y surgió con el advenimiento de las Segunda Revolución Industrial, caracterizada
por la automatización y el advenimiento de los computadores electrónicos digitales.
Como puede observarse, el punto de vista de la IO es muy amplio, pues apunta hacia
cualquier tipo de organización, tratando de resolver los conflictos entre los componentes de la
misma, para obtener el máximo provecho de los recursos disponibles y lograr los objetivos
globales de la empresa. Ese “máximo aprovechamiento de los recursos” se ve reflejado en la
mayoría de los casos que estudia la investigación operacional en la búsqueda de soluciones
óptimas de los problemas.
Por otra parte, el grupo que realiza un estudio específico en alguna organización debe
contener personas expertas en las áreas antes mencionadas y además poseer un conocimiento
10 Introducción a la Modelación Matemática y Optimización
global y profundo de la organización misma, de tal forma que la definición del problema, el
diseño del modelo y la utilización de sus soluciones estén enmarcados dentro del contexto del
sistema real y específico bajo estudio. De esto último puede deducirse que los resultados que
brinda la IO se obtienen en general de mediano a largo plazo, pues el conocimiento inicial del
sistema bajo estudio por parte del grupo puede consumir un tiempo considerable.
Como puede observarse, las posibilidades de aplicación son inmensas y cada día más
especializadas. Las aplicaciones reales de la IO están creciendo cada día, especialmente en
nuestro medio cuando apenas empieza a reconocerse su importancia y utilidad para el diseño y
optimización de sistemas complejos.
Entre las técnicas de la IO que más se han utilizado se resaltan el análisis estadístico
cuyo origen tuvo lugar con el desarrollo de la IO, para posteriormente volverse una técnica
independiente con sus propios desarrollos. Igualmente, la programación lineal y la simulación
han tenido y tienen actualmente un gran desarrollo en diversos campos, como por ejemplo los
sistemas integrados de producción y distribución, los sistemas de manejo de materiales,
algunas aplicaciones en agricultura, los problemas de planeación de las líneas aéreas y muchos
otros más.
Prácticamente cada una de estas técnicas podría dar lugar a cursos completos o
seminarios especializados, con amplia exposición de aplicaciones reales. Esta publicación se
centra en lo que podría denominarse un curso de Introducción a la Optimización o
Investigación de Operaciones I, concentrándose principalmente en los modelos de
programación lineal y brindando una introducción a la teoría de redes y a la programación no-
lineal.
12 Introducción a la Modelación Matemática y Optimización
2. EL PROCESO DE DISEÑO EN
INGENIERÍA
Como ya se ha visto, la IO aplica el método científico para analizar, modelar y
optimizar sistemas complejos reales. Desde el punto de vista de la ingeniería industrial y de
sistemas, esta metodología podría definirse como el “proceso de diseño en ingeniería.” Este
proceso consiste en una serie de pasos bien definidos para encontrar la solución de problemas
reales. La Figura 2.1 muestra una síntesis de este proceso.
Generación de alternativas de
solución
Formulación y
Evaluación y selección de la(s)
validación de un
mejor(es) alternativa(s)
modelo
Especificaciones de la(s)
solución(es)
Implementación, evaluación y
control del diseño utilizado
Los dos pasos siguientes son comunes al proceso bien sea que exista o no un modelo
de decisión. Se trata de especificar claramente la solución con el objeto de que sea
implementada de una forma correcta y de que permita su continua evaluación y control. Estos
pasos pueden identificar desviaciones de la realidad o cambios con respecto de la concepción
original del problema, de tal forma que se puede hacer necesario un proceso continuo de
análisis y refinamiento del modelo, incluso existiendo la posibilidad de reformular el problema
e iniciar el proceso de nuevo.
2.1. MODELOS
Uno de los pasos fundamentales del proceso descrito anteriormente puede ser la
formulación y validación de un modelo de decisión. En general, puede decirse que un modelo
es la abstracción de la realidad, el cual debe poseer dos características básicas:
En todas las disciplinas, de una u otra forma, se utilizan modelos. A continuación se da una
clasificación típica de modelos.
Físicos, los cuales normalmente lucen semejantes a la realidad física, pero a menor
escala, como por ejemplo un avión en miniatura siendo probado en un túnel de
viento.
Otro tipo de clasificación define los modelos como descriptivos, o sea aquéllos que
sólo describen el sistema bajo estudio, o normativos u optimizables, los cuales se formulan
para mejorar u optimizar el sistema bajo estudio. De manera semejante, los modelos se pueden
clasificar en determinísticos o estocásticos (probabilísticos), de acuerdo con su capacidad de
involucrar o no variables aleatorias en el análisis.
Existen otras clasificaciones más detalladas, pero con las dadas aquí es suficiente para
ilustrar los tipos de modelos objeto de esta publicación. En general, los modelos matemáticos
tratados en esta publicación, es decir los de programación lineal, modelos de redes y
programación no-lineal, pueden ser clasificados como modelos matemáticos determinístico-
optimizables.
2.1.2. Ventajas y Desventajas del Uso de Modelos
Por el contrario, existen también algunas desventajas que deben ser evaluadas antes de
emprender la tarea de formulación de un modelo para la solución de un problema. Ellas son:
Los modelos deben siempre considerarse como un soporte o ayuda para la toma de
decisiones y nunca como el reemplazo de las personas que toman las decisiones. Si se formula
un modelo pensando que sus resultados van a resolver todos los problemas de decisión y van a
ser implementados “a ojos cerrados,” se cae en un grave error que puede traer consecuencias
impredecibles. Existen en la realidad muchos factores que no pueden ser involucrados en
modelos matemáticos y que sólo los grupos de personas que toman las decisiones pueden
evaluar. Por lo tanto, desde este punto de vista, los modelos son una poderosa herramienta
para este proceso, pero no el fin mismo de la decisión.
Las tres comunas han llegado a un acuerdo para que cada comuna siembre la misma
proporción de su tierra irrigable disponible. Sin embargo, cualquier combinación de cultivos
puede sembrarse en cualquiera de las comunas. Deben, entonces, planearse cuántas hectáreas
destinar para cada tipo de cultivo en cada comuna. Formule un modelo de programación lineal
para este efecto, de manera que se maximice el retorno neto esperado de todas las comunas.
Una compañía multinacional posee n fábricas situadas en diferentes lugares, desde las
cuales surte a sus m principales distribuidores localizados también en diferentes lugares.
Continuamente las fábricas deben cubrir las demandas de sus distribuidores, teniendo en
cuenta su capacidad de producción. Si ai es la capacidad de producción de la planta i (i = 1, 2,
..., n), bj es la demanda del distribuidor j (j = 1, 2, ..., m), y cij es el costo unitario de transporte
desde la fábrica i hacia el distribuidor j (i = 1, 2, ..., n; j = 1, 2, ..., m), formule un modelo
matemático de programación lineal que permita encontrar los flujos entre cada fábrica y cada
distribuidor, de tal forma que se satisfaga la demanda al costo total mínimo de transporte.
Un salvavidas se encuentra situado en un punto S sobre una playa plana y recta, situado
a 30 m de la orilla, como muestra la Figura 2.2. Dentro del agua, en un punto P, situado 80 m a
la izquierda de S y 100 m dentro del agua, se encuentra un bañista en peligro de ahogarse. El
salvavidas puede correr a una velocidad lineal constante de 9 m/seg sobre la tierra y nadar a
una velocidad lineal constante de 2.5 m/seg. Formule un modelo de programación matemática
para determinar la ruta que debe seguir el salvavidas de tal forma que llegue lo más pronto
posible donde se encuentra la persona en peligro. Discuta detalladamente el mayor número
posible de supuestos que tuvo que considerar para obtener el modelo formulado.
P
agua
100 m
80 m 30 m
tierra
S
Figura 2.22.2.
Figura El El
problema del
problema del salvavidas
salvavidas
18 Introducción a la Modelación Matemática y Optimización
Una empresa posee dos plantas, P1 y P2, que van a alimentar a una nueva bodega, W, la
cual a su vez proveerá a tres centros de consumo, C1, C2 y C3. El producto A se produce en la
planta P1 y el producto B se produce en la planta P2. Las plantas enviarán sus productos
respectivos hacia la bodega, para ser distribuidos hacia los centros de consumo. El volumen
total, el costo unitario de transporte y la localización de cada instalación se muestran en la
Tabla 2.3. El costo de transporte se considera proporcional al volumen transportado y a la
distancia recorrida, medida como la distancia euclidiana entre los dos puntos, o sea que la
distancia entre los puntos (x1, y1) y (x2, y2) viene dada por:
Un hotel turístico muy demandado recibe una serie de ofertas para ocupación de uno de
sus mejores conjuntos de habitaciones entre Mayo 1 y Agosto 31. Las ofertas varían en cuanto
al precio de las habitaciones, ya que existen diferentes tipos de clientes, y en cuanto al día de
llegada y al día de salida de los turistas. Por ejemplo, una oferta puede ser por $800,000 entre
Junio 15 y Junio 18, mientras que otra oferta puede ser por $1,000,000 entre Junio 16 y Junio
20. Como el ejemplo muestra, no se permite superposición de ofertas, ya que se supone que
cada una de ellas ocuparía completamente las habitaciones disponibles. Por lo tanto, de las dos
ofertas mostradas en el ejemplo, sólo una puede escogerse. También, pueden existir períodos
de tiempo donde no exista ninguna oferta. Formule un modelo de redes que permita encontrar
la mejor forma de seleccionar las ofertas entre Mayo 1 y Agosto 31, de tal forma que se
maximice el ingreso total del hotel para este período.
2.3. ALGORITMOS PARA LA SOLUCIÓN DE
MODELOS
Una vez formulado un modelo matemático, éste debe resolverse para cumplir con los
pasos de evaluación de las alternativas y selección de la(s) alternativa(s) óptima(s). El proceso
que se utiliza para lograr este objetivo se denomina algoritmo. Un algoritmo es, por lo tanto,
un conjunto finito de pasos, especialmente diseñado para la solución de un problema.
(2.1)
Es importante comentar, finalmente, que existen problemas que en general son NP-
completos, pero que para algunos de sus casos particulares dejan de serlo. Un ejemplo típico
es la programación lineal entera. En general, un problema de programación lineal entera es
NP-completo. Sin embargo, cierta clase de problemas de transporte y el problema de
asignación, a pesar de ser problemas de programación lineal entera, no son NP-completos, ya
que su estructura particular permite diseñar algoritmos eficientes para encontrar su solución
óptima. A medida que se vayan desarrollando los diversos temas de esta publicación, se
aclarará el tipo de problema tratado y las características de sus diversos algoritmos de
solución.
2.3.3. Tipos de Algoritmos
Finalmente, existen sistemas expertos que pueden utilizarse para resolver cierto tipo de
problemas. En general, estos sistemas resuelven problemas operacionales, pero no son
adecuados para resolver problemas tácticos o estratégicos de mayor nivel.
Ejercicios 2.1
2) Una compañía petrolera está planeando construir un oleoducto para llevar petróleo crudo
desde un pozo hasta un punto donde se embarcará en tanques y será transportado a la
refinería. La figura siguiente muestra la disposición geográfica del pozo y del punto de
embarque.
Punto de embarque
Río 25 Km
Pozo de petróleo
320 Km
Por la ribera del río donde está ubicado el pozo de petróleo: $72/Km
Por la ribera del río donde está ubicado el punto de embarque: $90/Km
Atravesando el río (por cualquier parte): $150/Km
3) Un barco de carga tiene tres bodegas: en la proa, en el centro y en la popa. Los limites de
capacidad de cada bodega son:
Usted como propietario del barco puede aceptar el total o una parte de los artículos que se
ofrecen para el transporte. Las características de estos artículos son las siguientes:
1
WinQSB (“Quantitative Systems for Business para Windows”), versión 1.00, por Yih-Long Chang.
ARTICULO CANTIDAD VOLUMEN UTILIDAD
EN TON. PIES CUB/TN $/TN
A 6000 60 6000
B 1000 50 8000
C 2000 25 5000
Para preservar el equilibrio del barco el tonelaje transportado en cada bodega debe ser
proporcional a su capacidad en toneladas. Formule un modelo matemático que le permita
determinar la mejor forma de cargar el barco para obtener la utilidad máxima.
4) Su compañía posee tres plantas, en las cuales elabora un componente pequeño para un
producto industrial. La compañía comercializa el producto a través de cinco distribuidores
en el país. Los pronósticos de ventas indican que los requerimientos mensuales por
distribuidor son los siguientes:
Distribuidor 1 2 3 4 5
Demanda mensual 2700 2700 9000 4500 3600
Planta 1 2 3
Capacidad 4500 9000 11250
Costo Unitario
de produccion $ 60 30 54
Los costos de envío a distribuidores desde las plantas se muestran en el cuadro siguiente,
en $/unidad:
Distribuidor 1 2 3 4 5
Planta 1 1.5 2.1 3.3 4.5 4.8
Planta 2 2.4 1.8 3.0 3.6 4.5
Planta 3 3.0 2.7 2.7 3.0 4.8
Formule un modelo matemático que le permita determinar dónde producir los
componentes y la forma de despacharlos hacia los distribuidores.
Se trata de obtener una mezcla de una tonelada (1000 Kg), cuyo contenido del factor f sea
por lo menos del 18% y con la condición que las materias primas B y C no constituyan
más del 20% de la mezcla. Formule un modelo matemático que le permita determinar
cuánto de cada materia prima debe utilizar para la mezcla con el mínimo costo posible.
6) Una industria de papel produce pulpa la cual puede vender al mercado local o utilizar para
fabricar papel blanco o cartón. Una tonelada de cartón requiere 0.7 ton de pulpa, mientras
que una tonelada de papel blanco consume 0.9 ton de pulpa (el cartón y el papel blanco
requieren de otras materias primas que no se consideran en este problema). La pulpa se
produce a partir de bagazo de caña de azúcar, con un rendimiento del 40%. Se dispone de
260,000 ton/año de bagazo. Las instalaciones para producir pulpa tienen capacidad para
250 ton/día. La máquina de cartón trabaja a una velocidad efectiva de 200 ton/día y la de
papel blanco a 150 ton/día. La producción se hace durante 335 días al año, ya que el resto
de tiempo se dedica a mantenimiento.
Cada tonelada de papel blanco producida arroja al río 10 unidades de contaminación; una
tonelada de cartón arroja 6 unidades y cada tonelada de pulpa arroja 20 unidades. Se
permite un máximo total de 1,000,000 unidades de contaminación/año arrojada al río.
Las utilidades netas por tonelada de pulpa, cartón y papel blanco son $50, $60 y $80,
respectivamente.
Formule un modelo de PL que permita estimar el mejor plan de producción anual. Asuma
que todo lo que se produce puede venderse.
3. FORMULACIÓN DE MODELOS DE
PROGRAMACIÓN LINEAL
3.1. ORÍGENES Y ANTECEDENTES HISTÓRICOS
La Programación Lineal (PL) es una de las técnicas de la investigación de Operaciones.
Muchos autores consideran que ha sido uno de los más importantes avances científicos del
presente siglo y, de hecho, su gran aplicación y la magnitud de todos los problemas que ha
resuelto, así lo confirman. Mediante la utilización de la PL se han logrado ahorros millonarios
en las organizaciones que la han aplicado.
Una vez definidas las actividades y la función objetivo, se pasa a estudiar las
restricciones del sistema, las cuales son el reflejo de los recursos limitados de que se dispone.
El desarrollo de cada actividad consumirá parte de los recursos y la cuantificación total de la
contribución de todas las actividades conformará las restricciones del sistema. Estas
restricciones se expresan como funciones matemáticas lineales de igualdad o desigualdad.
2
Para mayores detalles sobre esta discusión, consultar Chvátal (1983), Capítulo 4, pág. 45–52.
acotamiento requiere de la solución de un gran número de subproblemas de programación
lineal continúa.
La obtención del modelo se basa en la observación cuidadosa del sistema bajo estudio.
Este paso es fundamental para que los resultados eventuales del modelo tengan validez.
3.3.1. Proporcionalidad
Este no es el caso más común en la realidad. Sin embargo, dentro de ciertos rangos
normales de operación de los sistemas, existe comúnmente la proporcionalidad. Fuera de estos
rangos es probable que haya que utilizar otras técnicas, tales como la programación no-lineal.
3.3.2. Aditividad
El hecho de que exista la proporcionalidad para todas las actividades no garantiza que
todas las funciones sean lineales, ya que puede existir cierto grado de dependencia entre las
actividades, lo que ocasionaría la aparición de no-linealidades. Lo que plantea la aditividad,
por lo tanto, es que las contribuciones de cada actividad a la función objetivo y a las
restricciones sean independientes de otras actividades.
A manera de ilustración, supóngase que una empresa produce los artículos A y B, los
cuales producen utilidades unitarias de $30/unidad y $20/unidad, respectivamente, e
independientemente de las cantidades que se produzcan. Así, la función objetivo se escribiría
como:
(3.1)
Supóngase ahora que el hecho de producir una mayor cantidad del artículo B, hace que
los costos de producción del artículo A se incrementen. Así, por ejemplo, la utilidad unitaria
del artículo A se vería disminuida mediante cierto factor y podría expresarse como $(30 –
0.0001XB)/unidad. La función objetivo sería por lo tanto:
– (3.2)
Cabe la pregunta de por qué la PL se utiliza tan ampliamente si estos supuestos casi
nunca se cumplen en la práctica?
Algunos autores hablan de otros dos supuestos, la certeza y la divisibilidad (Ver Hillier
y Lieberman (1997)). La divisibilidad se aplica fundamentalmente a problemas de números
enteros. Existen muchos casos prácticos en los cuales las variables de decisión o actividades
sólo tienen significado físico si sus valores son enteros. Por ejemplo, si una compañía de
aviación está analizando cuántos aviones comprar de cierto tipo, es obvio que una respuesta de
2.33 aviones no es satisfactoria y deberá darse una respuesta entera en su lugar, la cual no
necesariamente corresponde a alguno de los enteros vecinos. Para problemas enteros se han
desarrollado técnicas de solución especiales, las cuales se siguen investigando actualmente.
Una introducción a este tópico se presentará posteriormente.
La certeza, más que un supuesto, es una característica propia de los modelos de PL. Se
dice entonces que todos los valores de los parámetros que alimentan el modelo tienen valores
determinísticos, o sea que no se consideran como variables aleatorias. Una demanda de un
producto, por ejemplo, se obtiene como un valor promedio estimado, pero su distribución
probabilística propia no se considera en el modelo. En la mayoría de los casos prácticos,
ningún parámetro o variable puede considerarse determinístico. En la mayoría de los casos, los
parámetros y variables tienen un rango de validez entre los cuales pueden oscilar sus valores y,
por lo tanto, se constituyen en variables aleatorias cuya distribución probabilística puede o no
ser conocida. Actualmente están en investigación y desarrollo técnicas de programación lineal
estocástica y programación lineal con conjuntos difusos, donde se consideran algunos
parámetros de los modelos como variables aleatorias3.
3
Ver, por ejemplo, Taha (1998), pág. 807.
Dado que el desarrollo de la programación lineal estocástica es muy incipiente, el
problema de la incertidumbre de los parámetros de los modelos de PL aún se resuelve de una
forma relativamente sencilla, haciendo uso de los análisis de sensibilidad. Normalmente se
obtiene la solución óptima utilizando los mejores estimados de los parámetros y,
posteriormente, se analizan los posibles efectos de su variación sobre la solución óptima
original. Si estos efectos no son significativos, no se justifica invertir mayor esfuerzo, tiempo
y dinero en su estimación más precisa. Si por el contrario, se observa que el efecto de cambios
pequeños en los valores de ciertos parámetros es significativo, debe prestársele especial
atención a su estimación más precisa y a su posibilidad de cambio real en el sistema bajo
estudio. Igualmente, existen otras tendencias actuales tales como la combinación de modelos
matemáticos con modelos de simulación4, los cuales representan de una manera más fiel los
sistemas y sus variabilidades, y otros como el análisis de diversos escenarios aleatorios con la
ayuda de modelos de optimización5.
Una compañía productora de elementos eléctricos tiene durante este mes un sobrante
en su capacidad total de producción, el cual quiere utilizar para la manufactura de dos artículos
de rápida venta, los transformadores de 40 VA y los transformadores de 75 VA. Por su
experiencia, se han reunido los siguientes datos:
a) Variables de decisión
4
Ver, por ejemplo, Hicks (1999).
5
Ver, por ejemplo, Escudero et al. (1999).
30 Introducción a la Modelación Matemática y Optimización
Las actividades en las que está interesada esta empresa en el momento son: Producir
transformadores de 40 VA y/ó de 75 VA. Por lo tanto, las variables de decisión son:
b) Función Objetivo
Obsérvese que es muy importante definir las unidades de la función objetivo. En este
caso la función “U” debe ser MAXIMIZADA.
c) Restricciones
Las restricciones surgen por la limitación de los recursos. En este caso sólo se dispone
de cierto número de hr.hombre y de hr.máquina para distribuir entre las dos actividades
definidas. Cada actividad “gasta” cierta fracción de esos recursos. Con base en lo anterior, las
restricciones son:
Función Objetivo
Sujeto a:
Restricciones
Restricciones Obvias
Debido a la difícil situación económica por la que atraviesa el mundo, suponga que
usted decidió alimentarse solo con huevos y arroz. Asuma que el cuerpo humano como
mínimo debe disponer diariamente de 2700 Kcalorías, 56 gr de proteínas y 1.4 mg de vitamina
B1. Usted sabe que el contenido que brindan los alimentos seleccionados y su costo
aproximado son:
HORAS- HORAS-
UTILIDAD HORAS-
TRANSFOR MÁQUINA 1 MÁQUINA 2
NETA HOMBRE POR
MADOR POR POR
UNITARIA [$] UNIDAD
UNIDAD UNIDAD
40 VA 400 1 1,0 1,0
75 VA 700 2 1/3 1,4 1,0
32 Introducción a la Modelación Matemática y Optimización
Nota: Una porción se refiere a una cantidad estándar normal para el consumo humano.
Por ejemplo, una porción de huevo = 1 huevo y una porción de arroz = 100 gr de arroz.
a) Variables de decisión
Sean:
c) Restricciones:
Note que en este caso las restricciones son de mayor o igual, pues se trata de
requerimientos mínimos. Además no aparece inicialmente la restricción por enteros (variables
continuas).
La industria tiene en sus bodegas rollos semejantes, pero de 114 pulg. de ancho, y en
cantidad suficiente y decide utilizarlos para el pedido, cortándolos en los diferentes anchos
solicitados. ¿Cuál es la mejor forma de cortar los rollos de 114 pulg. de ancho para satisfacer
el pedido y minimizar el desperdicio de papel?
Xi = Número de rollos de 114 pulg. de ancho a cortar según el patrón i (i = 1, 2, .., 12).
b) Función Objetivo
c) Restricciones
34 Introducción a la Modelación Matemática y Optimización
Las restricciones en este caso surgen simplemente de la satisfacción del pedido, así
Note que lo que interesa es simplemente satisfacer el pedido y por ello las restricciones
son de igualdad. Sin embargo, la solución del modelo varía notablemente con restricciones de
(ver pág. 26- A). El modelo matemático será la reunión de la función objetivo y las
respectivas restricciones.
2. Con restricciones de
X 3 100
X 4 50
X 7 150
X 9 300
Dmín = 1800 pul. Aquí sobran rollos, pero se cumple con el pedido. Comente sobre las
implicaciones de estas dos soluciones y su relación con casos reales.
3.4.4. Modelo Nº 4: Carga de un Avión
El propietario del avión tiene posibilidad de llevar parte de la carga o toda la que se le
ofrece (si tiene capacidad). Esta carga y sus características son las siguientes:
VOLUMEN EN UTILIDAD
CANTIDAD
CLASE DE METROS OBTENIDA POR
TIPO OFRECIDA PARA
CARGA CUBICOS POR SU TRANSPORTE
LLEVAR (EN TON.)
TON. EN $/TON.
1 Herram. 20 1,0 250.000
2 Libros 15 2,0 280.000
3 Flores 8 10,0 500.000
4 Artesanías 10 6,0 360.000
a) Variables de decisión
En este tipo de problema se empieza a ver la necesidad de definir las variables con dos
subíndices, pues la decisión debe hacerse en dos etapas: qué tipo de mercancía y a dónde
llevarla. Así:
Es decir, que existirán 4 × 3 = 12 variables de decisión, a saber: X11, X12, X13, X21, X22,
X23, X31, X32, X33, X41, X42 y X43.
Si se considera más fácil, las variables podrían reenumerarse, por ejemplo, Y1, Y2,.....,
Y12 ó incluso X1, X2,....., X10, X11, X12. Sin embargo, la utilización de los subíndices siempre
deja mayor claridad en el significado físico de las variables.
36 Introducción a la Modelación Matemática y Optimización
b) Función Objetivo
c) Restricciones
En este caso existen cuatro grupos diferentes de restricciones, así: Por la capacidad en
peso de cada compartimiento del avión, por la capacidad volumétrica de los mismos, por la
cantidad de tipo de carga ofrecida para su transporte y por la condición de equilibrio del avión.
Estos grupos de restricciones son:
Restricciones obvias:
d) Modelo de Programación Lineal
Utilidad máxima: 1.605 x 107 [$]. El avión debería cargarse de la siguiente forma (en
toneladas):
X11 = 5 X32 = 8
X12 = 12 X41 = 4.8
X21 = 6.2 X43 = 5.2
X23 = 8.8
Holguras:
Un fabricante debe cumplir un contrato a cuatro meses durante los cuales varían los
costos de producción. El costo de almacenamiento de unidades producidas en un mes
determinado y no vendidas en ese mes es de $10 por unidad y por mes. Se dispone de la
siguiente información:
Este modelo ilustra el hecho de que pueden existir formas diferentes de definir las
actividades. Algunas de estas formas pueden ser mejores y más manejables que las otras y la
eficiencia del modelo obtenido depende en gran parte del “arte” de definirlas correctamente.
38 Introducción a la Modelación Matemática y Optimización
a) Variables de decisión
b) Función Objetivo
La función objetivo tiene dos componentes: Los costos de producción y los costos de
almacenamiento.
Para encontrar la expresión para estos costos, es necesario ilustrar el “balance de las
unidades” a través del tiempo, así:
MES
1 2 3 4
DESCRIPCIÓN
Inventario inicial 0 X1-20 X1+X2-50 X1+X2+X3-100
Producción X1 X2 X3 X4
Ventas 20 30 50 40
Inventario Final X1-20 X1+X2-50 X1+X2+X3-100 X1+X2+X3+X4-140
c) Restricciones:
Obvias:
d) Modelo matemático:
Sujeto a:
Nótese que la condición j i indica el significado lógico de que una unidad sólo puede
ser vendida después de haber sido producida; así, las variables de decisión son:
X11, X12, X13, X14, X22, X23, X24, X33, X34 y X44,
40 Introducción a la Modelación Matemática y Optimización
b) Función Objetivo
En este caso no hay necesidad de construir el flujo de las unidades para encontrar la
función objetivo, pues ya lo tienen intrínseco. Así:
En la expresión anterior, el primero término representa el costo de las unidades que han
estado almacenadas durante un mes, el segundo término las que han estado almacenadas
durante dos meses, y el último durante tres meses. Así la función objetivo será:
c) Restricciones:
Obvias:
41
d) Modelo matemático:
Sujeto a:
X11 = 20 X11 = 20
X13 = 20 X12 = 20
X22 = 30 X22 = 10
X33 = 30 X23 = 20
X44 = 40 X33 = 30
X44 = 40
Obsérvese la completa correspondencia entre las dos soluciones. Nótese que este
modelo aunque tenga igual número de restricciones que el de la primera forma de solución,
tiene 10 variables, mientras que el primero tiene 4. Esto podría ser una desventaja de manejo
del modelo. Sin embargo, este último permite no sólo tener las unidades que se deben producir
en un mes determinado (tal como el primer modelo), sino también cuando deben venderse. O
sea que, indudablemente el segundo modelo es superior al primero.
42 Introducción a la Modelación Matemática y Optimización
Una compañía multinacional posee “n” fábricas en diferentes partes del mundo, desde
las cuales surte a sus principales distribuidores, “m” en total, localizados en diferentes países
(en algunos posee fábricas y distribuidores simultáneamente). El costo de transporte desde las
fábricas hacia los distribuidores es diferente para cada uno de ellos. Continuamente las
fábricas deben cubrir las demandas de sus distribuidores, teniendo en cuenta su capacidad de
producción.
Si: a1, a2, a3, ..., an son las capacidades de producción de las plantas 1, 2, ..., n,
respectivamente; b1, b2, b3, ..., bm son las demandas de los distribuidores 1, 2, ..., m,
respectivamente, y Cij (i =1, 2, 3,...,n; j =1, 2, 3,..., m) son los costos unitarios de transporte
desde la fábrica i hacia el distribuidor j, formúlese un modelo de Programación Lineal para
satisfacer la demanda de todos los distribuidores sin violar la capacidad de producción de las
plantas, con el costo total mínimo de transporte.
a. Variables de decisión
b. Función Objetivo
c. Restricciones
Obvias:
a. Lograr un Valor Presente Neto mínimo de mil millones de pesos (Utilidad a largo
plazo).
b. Mantener el recurso laboral actual de 100 empleados (Nivel de empleo).
c. Sostener la inversión de capital en el nuevo equipo de 400 millones de pesos
(Inversión inicial).
CONTRIBUCIÓN UNITARIA
UNIDAD DE
PRODUCTO PRODUCTO PRODUCTO
OBJETIVO MEDIDA
1 2 3
Decenas de
(a) 15 8 12
millones de $
Cientos de
(b) 3 2 1
empleados
Decenas de
(c) 4 5 7
millones de $
¿Cuáles deben ser las tasas de producción de cada producto para que los objetivos se
cumplan de la mejor forma posible?
44 Introducción a la Modelación Matemática y Optimización
Si se definen las actividades como Xi = Tasa de producción del producto i (i=1,2, 3),
las metas a cumplir serían las siguientes (en su orden: Objetivos (a), (b) y (c)):
Obsérvese que las variables Yi (i =1, 2, 3,) pueden ser positivas o negativas. Por
ejemplo, si Y1> 0, indica que la utilidad ha sobrepasado los 1000 millones de pesos, pero si
Y1< 0, entonces la utilidad ha sido inferior a esa cifra y el objetivo no se habría cumplido.
Dado que manejar variables libres (las que pueden tomar valores positivos) no es
conveniente y debido a la naturaleza del problema, se reemplaza cada variable Yi (i = 1, 2, 3)
por la diferencia de dos variables positivas, así:
Sujeto a:
45
Obsérvese que hay variables que aparecen en las restricciones mas no en la función
objetivo, y más aún, las variables de decisión no aparecen directamente en la función objetivo.
Esto puede parecer sorprendente, pero en realidad las variables de decisión están
estrechamente relacionadas con las variables auxiliares de la función objetivo (recordar la
definición de estas últimas). La solución óptima de este modelo, aplicando el método simplex,
es:
O sea que el producto 2, de acuerdo con la solución óptima, no debería producirse. Las
metas de utilidad a largo plazo y de inversión inicial se cumplen a cabalidad, produciéndose
100 decenas de millones de utilidad es invirtiéndose inicialmente 40 decenas de millones de
pesos. La meta de nivel de empleo no puede ser cumplida y, dado que
Y2 Y2 Y2 , Y2 5.087 0.0000 5.0877 , el nivel de empleados debe ser
aumentado en 5.0877 cientos, para poder cumplir con las otras dos metas. Obsérvese que el
puntaje óptimo P* = 15.2632 se obtiene de penalizar con 3 puntos por cada 10 empleados de
más en la función objetivo.
KILÓMETROS NÚMEROS
TIPO DE BOMBARDERO
POR GALÓN DISPONIBLE
Pesado 3,0 48
Mediano 4,5 35
PROBABILIDAD DE SER
DISTANCIA
DESTRUÍDA POR:
PLANTA DESDE LA
BOMB. BOMB.
BASE EN Km.
PESADO MEDIANO
1 675 0,10 0,08
2 720 0,20 0,16
3 810 0,15 0,12
4 900 0,25 0,20
a. Variables de decisión:
Sean Xij = Número de bombarderos tipo “i” a ser enviados a la planta “j”;
b. Función objetivo:
En este caso hay dos formas de plantear la función objetivo: Maximizar la probabilidad
de éxito ó Minimizar la probabilidad de fracaso.
Dado que un solo bombardero puede destruir la planta, la probabilidad de éxito sería
difícil de expresar, pues sería la suma de las probabilidades de todas las posibles
combinaciones. Es más sencillo minimizar la probabilidad de fracaso, pues el fracaso total de
47
la misión significa que todos los bombarderos fallen en todas las plantas y recuérdese que con
sólo destruir una planta se para la producción de tanques del enemigo. La probabilidad de
fracaso sería la intersección (producto) de todas las probabilidades de fracaso de cada
bombardero a cada planta (se considera independiente la acción de cualquier bombardero con
respecto a la de cualquier otro). Así, la expresión para la función objetivo sería:
Como puede observarse, así definida la función objetivo no sería lineal. Pero se puede
“linearizar” fácilmente por medio de la función log P. Dada la naturaleza de la función
logarítmica (Función biyectiva), es equivalente minimizar P que minimizar log P; así, la
función objetivo se convierte en:
La base del logaritmo sacado es indiferente, pues cualquier base logarítmica se puede
cambiar a otra multiplicando por una constante. Así, la función objetivo obtenida es lineal.
c. Restricciones:
Un avión cualquiera debe tener combustible para ir a la planta, volver y tener una
reserva de 150 galones; por ejemplo, los aviones tipo pesado que se envíen a la planta 1
utilizarán la siguiente cantidad de combustible:
Simplificando, se obtiene:
Restricciones Obvias:
Sujeto a:
La solución óptima de este modelo de PL entera es: X14 = 43 y X24 = 34 con Zmáx =
19.9572, ó, equivalentemente, una probabilidad de falla de P = e 19.9572 2.15 109 . O sea
que si se envían 43 aviones pesados y 34 aviones medianos, todos a la planta 4 del enemigo, la
probabilidad de falla de la misión es casi cero.
Ejercicios 3.1.
3. Una compañía de alquiler de camiones dispone de dos tipos de vehículos el tipo A que
posee 20 pies cúbicos de espacio refrigerado y 40 pies cúbicos de espacio no refrigerado.
El tipo B que posee 30 pies cúbicos de espacio refrigerado y la misma cantidad de espacio
no refrigerado. Una fábrica de alimentos debe transportar 900 pies cúbicos de producto
refrigerado y 1200 píes cúbicos de producto no refrigerado. ¿Cuántos camiones de cada
tipo se deben alquilar, si el camión A se alquila a $30 por milla y el camión B a $40 por
milla, de tal forma que se minimice el costo total del transporte por milla recorrida?
(Prado, Hernando, Notas de Clase, 1978.)
4. Una factoría posee dos minas. La mina A produce diariamente 1 tonelada de material de
alta calidad, 3 toneladas de calidad intermedia y 5 toneladas de baja calidad; la mina B
produce diariamente 2 toneladas de cada una de las tres calidades. La compañía necesita,
para su posterior procesamiento, al menos 100 toneladas de material de alta calidad, 150
de mediana calidad y 180 toneladas de baja calidad. ¿Cuántos días debe operarse sobre
cada mina para satisfacer las necesidades de la compañía si el costo diario de explotación
es de $200.000 en cualquier mina? (Prado, Hernando, Notas de Clase, 1978.)
5. Se dispone de n alimentos. Cada unidad del alimento j contiene Aij del nutriente i en sus
unidades características y tiene un costo de Pj pesos por unidad. Formule un modelo de
PL cuya solución proporcione los niveles Bi mínimos requeridos diariamente del nutriente
i de tal manera que el costo sea mínimo (problema de la dieta). Trate de obtener una dieta
real típica en la ciudad de Cali, a costo mínimo, investigando los requerimientos mínimos
de energía, proteínas, vitaminas, minerales y otros nutrientes (basado en una dieta de 2000
Kcal). Incluya alimentos de los cuatro grupos principales, a saber: Lácteos, carnes,
50 Introducción a la Modelación Matemática y Optimización
vegetales y frutas y harinas y cereales. Investigue igualmente los precios actuales de estos
alimentos y determine el valor promedio diario de una dieta balanceada. Resuelva el
problema, primero, sin restricciones de raciones mínimas de los alimentos y luego
considerando cantidades mínimas y máximas de dichos alimentos.
8. Formule el problema anterior si además del gran número de rollos de 60” de ancho, se
dispone también de 50 rollos de 80” de ancho.
9. Se hace un pedido a una papelería de 800 rollos de papel corrugado de 30” de ancho, 500
rollos de 45” de ancho y 1000 rollos de 56‟ de ancho. Si la papelería recibe de su
proveedor rollos de 108” de ancho únicamente y en cantidad suficiente, ¿cómo deben
cortarse estos para cumplir con el pedido en forma exacta y para minimizar el papel
desperdiciado? (Prado, Hernando, Notas de Clase, 1978.)
10. Formule el problema anterior si se le permite cortar rollos en exceso de cualquier ancho
(adicionales a los del pedido), los cuales generan un costo de inventario despreciable y
podrían ser consumidos en un futuro. Compare los resultados en cuanto al desperdicio
mínimo se refiere. Desprecie los costos de inventario que se pueden generar.
11. Usted recibe una orden de producción en su taller que requiere del corte de láminas de
acero de cierto calibre. Se dispone de 250 láminas rectangulares de 120 240 cm y de 70
láminas rectangulares de 100 150 cm que habían sobrado de una orden anterior. La
orden consiste en 75 láminas rectangulares de 50 80 cm y 120 láminas cuadradas de 100
51
100 cm. Formule un modelo de PL que le ayude a decidir la mejor forma de cortar las
láminas disponibles para cumplir con la orden. Asuma como desperdicio cualquier retal
de lámina del cual no pueda obtenerse ninguna de las láminas requeridas en la orden.
Asuma igualmente que sólo pueden realizarse cortes paralelos a las caras de las láminas.
12. Una industria produce dos artículos distintos A y B. La elaboración de una unidad del
articulo A cuesta $20 por concepto de mano de obra y de una unidad del articulo B, $10.
Cada unidad de A utiliza $10 de materia prima y cada unidad de B $30. El desgaste del
equipo se considera proporcional a la producción. Por cada unidad producida de A, el
equipo se desgasta $5 y por cada unidad producida de B, el equipo se desgasta $1. Se
cuenta con un presupuesto de $100.000 para salarios, de $180.000 para materia prima y
no conviene que el desgaste de los equipos exceda de $40.000. Determinar la cantidad que
debe producirse de cada artículo para obtener la máxima utilidad si el beneficio por cada
artículo A es de $8 y por cada artículo B es de $5.
13. Un ebanista dispone de dos tipos diferentes de madera, tiene 1500 pies de tabla tipo A y
1000 pies de tabla tipo B. También dispone de 800 horas hombre para efectuar el trabajo.
La demanda estimada es la siguiente: cuando menos 40 mesas, 130 sillas exactamente, 30
escritorios exactamente y no más de 10 estantes. Las cantidades de madera A y B y las
horas hombre que se requieren para la elaboración de cada uno de los artículos anteriores
vienen dadas en la tabla siguiente:
Determine cuantas unidades de cada uno de los artículos debe producir el ebanista con el
fin de cumplir su demanda y de obtener la máxima utilidad. (Prado, Hernando, Notas de
Clase, 1978.)
14. Un fabricante de aparatos de televisión tiene facilidades para ensamblar dos tipos de
televisor: el tipo A a color y el tipo B a blanco y negro. El televisor de tipo A se vende a
$49.500 la unidad y ocasiona un costo de producción de $26.800. El de tipo 5 se vende en
$23.400 y ocasiona un costo de producción de $11.190. La planta tiene capacidad diaria
para fabricar hasta 50 pantallas para televisión en colores. No se pueden comprar pantallas
a otros proveedores. Cada TV A requiere 18 horas hombre para el ensamblaje del chasis y
el tipo B requiere 8 hr.hombre. La planta emplea 225 hombres con un turno diario de 8
horas en el departamento de ensamblaje. Cada aparato tipo A requiere 1.6 horas hombre
para su armado completo, mientras que el tipo B requiere 1 hora hombre. La planta
emplea 30 hombres con turnos de 8 horas diarias para esta labor. Cada TV de color
requiere 2 horas hombre para inspección final mientras que cada TV en blanco y negro
requiere 0.5 horas hombre. La planta emplea 20 inspectores de tiempo completo y uno de
medio tiempo. ¿Cuántos televisores de cada tipo deben producirse para maximizar la
utilidad neta total? (Prado, Hernando, Notas de Clase, 1978.)
52 Introducción a la Modelación Matemática y Optimización
Se trata de obtener una mezcla de una tonelada, cuyo contenido del factor f sea por lo
menos del 18% y con la condición que las materias primas B y C no constituyan más del
20% de la mezcla, con el mínimo costo posible.
16. Un barco tiene tres bodegas: en la proa, en el centro y en la popa. Los limites de
capacidad de cada bodega son:
Los dueños del barco pueden aceptar el total o una parte de los artículos que se ofrecen
para el transporte. Las características de estos artículos son las siguientes:
Para preservar el equilibrio del barco el tonelaje transportado en cada bodega debe ser
proporcional a su capacidad en toneladas. Determine la distribución de carga para obtener
la máxima utilidad por el transporte. Formule este mismo problema desde el punto de
vista del propietario de la carga y compare resultados.
17. Se está diseñando una nave espacial que lleve y traiga astronautas a Marte. Esta nave
tendrá tres compartimientos, cada uno con su propio sistema que permite vivir en él. El
elemento clave de cada uno de estos sistemas es una pequeña unidad oxidante que
provoca un proceso químico para producir oxígeno. Estas unidades no se pueden probar
con anticipación y solo algunas de ellas tienen éxito en el proceso químico. Por esto, es
importante tener unidades de repuesto para cada sistema. Como los requerimientos son
diferentes para cada compartimiento, las características de las unidades para cada uno
varían. Se debe tomar una decisión sobre el número de unidades que se deben incluir en
cada compartimiento, tomando en cuenta limitaciones de diseño y la cantidad total de
espacio, peso y costo que puede asignarse a estas unidades para toda la nave. La siguiente
53
tabla resume estas limitaciones, al igual que las características individuales para cada
compartimiento:
Si todas las unidades fracasan en uno ó dos de los compartimientos, los astronautas
podrán ocupar el ó los restantes y continuar su viaje espacial, pero con algunas pérdidas
en la cantidad de información científica que puedan obtener. No obstante, si todas las
unidades fracasan, todavía tienen la manera de regresar en la nave a salvo, pero el viaje
completo sería un fracaso total a un gran costo. El objetivo es entonces minimizar la
probabilidad de que todas las unidades fallen, sujeto a las limitaciones anteriores y a la
restricción adicional de que cada compartimiento deberá tener una probabilidad no mayor
que 0.05 de que todas sus unidades fallen. Formule un modelo de PL para este problema.
(Hillier y Lieberman, segunda edición, 1989, pág. 273)
18. Cierto hacendado dispone de los siguientes recursos para emplearlos en la próxima
cosecha:
Esta tierra es propia para sembrar maíz, millo o fríjol. Se supone que tiene a su
disposición hombres suficientes. Los costos de producción son los siguientes: tractor e
implementos $50/hora, mano de obra $5/hora, alquiler del terreno por la cosecha
$200/hectárea. El hacendado ha acordado retornar el 150% del dinero no invertido
exclusivamente en los cultivos a la Caja Agraria, donde efectuó el préstamo. Los datos de
los cultivos posibles son los siguientes y vienen dados por hectárea:
19. Una refinería de petróleo tiene un oleoducto que la alimenta, el cual tiene una capacidad
de 50000 barriles de petróleo diarios. La refinería debe suministrar por lo menos 27000
barriles por día a las bombas vendedoras de la Compañía (barriles de gasolina). La
refinería utiliza dos procesos de refinación con los siguientes rendimientos por barril de
petróleo procesado:
54 Introducción a la Modelación Matemática y Optimización
20. Una fábrica de automóviles y camiones consta de los departamentos que a continuación se
enuncian:
El Departamento A puede estampar por mes las planchas necesarias para 25000
automóviles ó 35000 camiones ó las correspondientes combinaciones de automóviles y
camiones. El Departamento B puede armar por mes 33333 motores de automóvil ó 16667
motores de camión ó una combinación de ambos. El Departamento O puede montar y
terminar 22500 automóviles/mes. El Departamento D puede montar y terminar 25000
camiones/mes. Cada automóvil deja una utilidad neta de $40 u.p. y cada camión de $55
u.p. Planifique la producción para obtener la máxima utilidad.
21. La firma “PAVIMENTOS S.A.” está licitando por un contrato para la construcción de la
calzada de una carretera. Las especificaciones dadas indican que debe tener un mínimo de
12 cm. de espesor y un máximo de 48 cm. Debe, a su vez, construirse de concreto, asfalto,
gravilla o cualquier combinación de los tres, siempre y cuando la resistencia total sea al
menos equivalente a la que tendría una calzada de 9 cm. de concreto. La firma ha
establecido que 3 cm. de asfalto son tan fuertes como 1 cm. de concreto y que 6 cm. de
gravilla son tan resistentes como 1 cm. de concreto. Su costo estimado para un metro
cuadrado y un centímetro de espesor para el concreto es de $2000, para el asfalto es $700
y para la gravilla es de $300. Formule el problema de Programación Lineal que le permita
a la firma saber cual es la mejor combinación para la calzada y conocer su costo.
22. Una compañía de aviación tiene el proyecto de comprar aviones de vuelos cortos,
medianos y largos. Los aviones de vuelos cortos tienen un precio de US$3.500.000 cada
uno; los de vuelos medianos US$5.000.000 cada uno; y los de vuelos largos
US$6.700.000 cada uno. Se tiene un presupuesto máximo de US$150.000.000 para el
55
23. Una compañía de aviación que opera con base en una terminal central, tiene 8 aviones
tipo 1, 15 aviones tipo 2 y 12 aviones tipo 3, disponibles para vuelos diarios. La capacidad
en toneladas es de: 45 para los de tipo 1, 7 para los de tipo 2 y 4 para los de tipo 3. La
compañía despacha sus aviones a las ciudades A y B. Las demandas en toneladas en cada
ciudad son: 200 para la ciudad A y 180 para la ciudad B. El exceso de toneladas
suministrado a una ciudad no tiene costo adicional. Cada avión puede hacer solamente un
vuelo al día. El costo del envío de un avión hacia cada ciudad está dado por la siguiente
tabla:
24. Una compañía de aviación debe decidir a cuantas nuevas azafatas emplear y entrenar
durante los 6 meses siguientes. Las necesidades expresadas como el número de azafatas ·
hora · vuelo son 8000 en enero, 9000 en febrero, 8000 en marzo, 10000 en abril, 9000 en
mayo y 12000 en junio. Para que una azafata pueda ser puesta en vuelo regular toma un
mes su entrenamiento. Una azafata aprendiz requiere 100 horas de supervisión por una
azafata experimentada durante el mes de entrenamiento; por lo tanto, durante dicho mes
se disponen de 100 horas menos de azafatas para el servicio regular de vuelos. Cada
azafata experimentada puede trabajar hasta 150 horas en el mes y la compañía dispone de
60 azafatas a comienzos de enero. Si el máximo tiempo disponible de azafatas
experimentadas excede a la necesidad de un mes de vuelos y de requisitos de
entrenamiento, las azafatas regulares trabajan menos de 150 horas, pero ninguna es
despedida. Cada mes aproximadamente, el 10% de las azafatas experimentadas se retiran
de la compañía, para casarse. Una azafata experimentada cuesta a la compañía $800
mensuales y una aprendiz $400. Formule un modelo de Programación Lineal para la
situación de empleo y entrenamiento.
56 Introducción a la Modelación Matemática y Optimización
Note usted que el período 1 sigue inmediatamente al periodo 6. Cada policía trabaja 8
horas consecutivas. El departamento de policía busca un programa de trabajo diario que
emplee el menor número de policías en el departamento, teniendo presente cada uno de
los requerimientos anotados. (Prado, Hernando, Notas de Clase, 1978.)
26. Una cervecería posee tres plantas localizadas en tres ciudades diferentes, con las cuales
surte el consumo del país dividido en cuatro zonas: zona norte, zona central, zona
suroriental y zona del pacífico, a las cuales se denominarán respectivamente A, B, C y D.
Mensualmente la planta I produce 1.300.000 litros, la planta II 700.000 litros y la planta
III 300.000 litros. La demanda mensual en las diferentes zonas de consumo es como
sigue: zona A 600.000 litros, zona B 400.000 litros, zona C 400.000 litros y zona D
500.000 litros. Los costos combinados de producción y el transporte de las diferentes
plantas a los centros de consumo, se muestran en la tabla siguiente:
27. Un individuo cuyo negocio es mezclar whisky importa tres grados A, B y C. Los combina
de acuerdo con recetas que especifican los porcentajes máximo y mínimo de los grados A
y C en la mezcla. Estos porcentajes se dan en la tabla siguiente:
La provisión de los tres whiskies básicos, junto con sus costos, se presenta en la tabla que
sigue:
57
Proyecte una política de producción que maximice las ganancias. (Prado, Hernando,
Notas de Clase, 1978.)
28. Una empresa siderúrgica produce tres aleaciones diferentes. El diagrama de proceso como
sigue:
Aleación 1
CAJA DE
RECOCIDO
Aleación 2
RECOCIDO MOLINOS
Aleación 3 CONTINUO CONTINUO
Se deben determinar las cantidades de cada aleación que deben producirse dentro de las
limitaciones del volumen de ventas y las capacidades de las máquinas con el fin de
maximizar las ganancias. Los datos sobre las capacidades y utilidades se presentan en las
tablas siguientes:
Los rollos de cada aleación son de 400 pies de longitud y pesan 4 toneladas. Formule un
modelo de programación lineal del cual pueda obtenerse una política de producción para
la siderúrgica (Asuma que 1 mes = 4 semanas).
29. Una planta de productos químicos fabrica dos productos A y B, los cuales tienen que
pasar por cuatro centros de proceso: 1, 2, 3 y 4, según se muestra con las líneas continuas
en la figura:
A
CENTRO 2
CENTRO 1 CENTRO 4
B
CENTRO 3
Cada centro puede manejar solamente el paso de un producto a la vez. Sí hay capacidad
disponible en el centro 3, es posible enviar el producto A a través de 3, en lugar de hacerlo
dos veces a través del centro 2, pero esto es más costoso. Con la información de la tabla,
¿cómo deberá programarse la producción para hacer máximas las ganancias?. Se entiende
por Programa de Producción la especificación de las siguientes cantidades:
(1) La cantidad de materia prima utilizada para A diariamente, por el curso regular.
(2) La cantidad diaria de materia prima usada para A diariamente, por el curso opcional y
(3) La cantidad diaria de materia prima usada para B.
1 500 90 300
B 3 480 85 250
4 400 80 240
59
Los centros 1 y 4 trabajan hasta 16 horas al día; los centros 2 y 3 trabajan hasta 12 horas
al día. Una restricción final la proporcionan las facilidades de envío que limitan la
producción diaria de A y B a un total de 2500 galones. Formule el modelo de
programación lineal correspondiente.
30. Una refinería produce dos tipos de gasolina, la regular y la extra, y otros productos
derivados del petróleo, cuyos precios de venta por galón son $2.500 para la gasolina
regular, $3.200 para la gasolina extra y se puede considerar un promedio de $1.100 para
otros productos. Ambos tipos de gasolina y los otros productos se fabrican a partir de dos
clases de petróleo crudo, el nacional y el importado. Las tablas siguientes muestran las
especificaciones principales que debe cumplir cada tipo de gasolina y las características
de cada tipo de petróleo crudo.
La refinería puede emplear tres procesos de refinación con las siguientes características de
eficiencia y costo:
Asumiendo que en cada proceso los componentes de una mezcla de petróleos crudos
contribuyen al octanaje general y a la presión de vapor proporcionalmente a su cantidad
en galones, formule un modelo de PL que le permita determinar la mejor programación de
producción de la refinería.
31. Una compañía manufactura secadoras y lavadoras automáticas de ropa para consumo
nacional. La gerencia de producción ha calculado que por cada secadora requiere 1.5
horas de mano de obra y por cada lavadora 2 horas. La fuerza laboral de la empresa en el
último cuatrimestre de este año se estima en 5000 horas. La gerencia no desea modificar
la fuerza laboral en más de un 10% por periodo. Las metas de venta en los próximos
cuatro periodos del próximo año se ha fijado así:
Los costos asociados con cada esquema de producción deberán incluir costos de
producción (excluyendo mano de obra), costo de almacenamiento y costo de mano de
obra. La sección de costos de la compañía ha desarrollado la siguiente tabla que describe
los costos unitarios para cada periodo:
Establezca esquemas de producción que minimicen los costos anuales para satisfacer las
demandas estimadas para ambos productos.
32. Un inversionista puede elegir entre las actividades A ó B ó ambas, disponibles al principio
de cada uno de los próximos 5 años. Cualquier cantidad invertida y recuperada en el
futuro puede ser reinvertida en cualquier alternativa disponible. Cada peso que invierte en
A al comienzo de cada año le produce $1.40 dos años más tarde. Cada peso invertido en B
al comienzo de un año le produce 1.70 tres años después. Además, las actividades C y D
están disponibles una sola vez en el futuro, C al comienzo del segundo año y D al
comienzo del quinto año. Cada peso invertido en C genera $1.60 en 2 años. Cada peso
invertido en D le produce $1.30 un año después. El inversionista dispone hoy de
$100.000. Formule un modelo de programación lineal que le permita determinar la mejor
forma de inversión a lo largo de los cinco años para maximizar el capital dispone al final
61
del quinto año (comienzo del sexto año). (Adaptado de Hillier y Lieberman, 1989, pág.
50.)
33. La gerencia de una compañía aérea debe decidir con respecto a la asignación de aviones a
rutas aéreas. El cuadro siguiente indica la capacidad máxima en número de pasajeros, el
número de aviones de cada tipo actualmente disponible en la empresa, el máximo número
de vuelos que realiza cada tipo de avión por ruta específica y el número esperado de
pasajeros en cada ruta.
Los costos de operación por vuelos en las cuatro rutas y el costo de oportunidad (pérdida
– ganancia) originado al no servir a un pasajero se muestran así:
Formule un modelo de programación lineal para asignar los aviones a las rutas a costo
mínimo. (Adaptado de Taha, 1981, pág. 37–38; Taha, 1998, pág. 60.)
34. Usted ha sido contratado como asesor del departamento de planificación nacional grupo
de desarrollo energético. El grupo está considerando la construcción de una planta nuclear
y de una planta desalinizadora de agua para un complejo agroindustrial. La planta nuclear,
un reactor adaptado a una planta de vapor, es la única fuente de electricidad. La inversión
en la planta es del orden de $9.000 millones de pesos por cada 1000 megavatios de
capacidad y los costos de operación anuales para esta capacidad son de $1350 millones. El
agua para las plantas industriales y la agricultura será provista por una planta
desalinizadora, con costo de inversión de $1500 millones por cada 1000 millones de
metros cúbicos de agua al año y costo de operación anual de $450 millones. Se estima que
para esta producción de agua se necesitan 300 megavatios de potencia por año. En la
región hay 120.000 hectáreas aptas para la agricultura. Existen fundamentalmente dos
cultivos que pueden ser desarrollados a escala industrial: Cultivo A y cultivo B. Los
fertilizantes se obtendrán básicamente de lo producido por la industria de la región. La
siguiente tabla muestra la inversión, costo de operación, agua, fertilizantes y un estimado
de los ingresos generados por una extensión determinada dedicada a cada uno de los
cultivos.
62 Introducción a la Modelación Matemática y Optimización
Suponga que para este tipo de proyecto se dispone de financiación por prestamos
internacionales de hasta $30.000 millones de pesos a ser pagaderos en 10 años a una tasa
anual de 15%. Analice el problema y si es posible plantee un modelo que permita la
evaluación de posibles alternativas.
35. Suponga que en cierta población se pretende hacer inversiones cuantiosas en el cultivo de
aguacate, lima reina, mango y zapote prieto. Se persiguen dos objetivos: reducir el
desempleo rural y aumentar las exportaciones que vendrán a equilibrar la balanza de
pagos de la nación. Se sabe que la producción promedio de cada árbol está dada por la
siguiente tabla:
36. La Fábrica Nacional de Cerveza tiene plantas ubicadas en Bogotá, Cali, Barranquilla,
Medellín y Palmira. La Fábrica Nacional de Latas, una subsidiaria, tiene plantas ubicadas
en Pereira, Bucaramanga y Cartagena. La demanda mensual de latas de cerveza se
pronostica así:
Los fletes son una función de las distancias que existen entre las plantas productoras de
cerveza y las plantas productoras de latas. Estos fletes son:
37. Se supone que la Secretaria de Educación Pública se ha trazado un plan a 4 años, para
entrenar maestros urbanos y rurales de enseñanza primaria. El entrenamiento de ambas
clases de maestros tarda 1 año y los costos de entrenamiento están dados por:
Bajo tales circunstancias, ¿Cuántos maestros de cada tipo deberían prepararse por año,
cuantos deberían cambiar de calidad por año y cuantos maestros Jubilados deberían
contratarse por año, tal que el costo total se minimice y se atenga a la demanda
pronosticada?
38. Una siderúrgica debe decidir cuántas libras de acero puro y cuántas de chatarra utilizar en
la preparación de una aleación para un cliente. El costo por libra de acero puro es de $3 y
de chatarra $6 (por las impurezas), la demanda del cliente es de por lo menos 5 libras y
aceptaría más si se le ofrece. La disponibilidad de acero puro es 4 libras y la de chatarra
de 7 libras. La relación entre chatarra y el acero no puede exceder de 7/8. La fábrica tiene
18 horas disponibles para el proceso; una libra de acero puro requiere 3 horas mientras
que la chatarra solo requiere 2 horas. Formule un modelo de programación lineal que
determine la carga óptima del horno.
39. Una empresa fabrica dos modelos de productos: Z-1200 y Z-1500. Los requerimientos de
producción y las disponibilidades están mostradas a continuación:
Los beneficios unitarios logrados a la venta de los modelos Z1200 y Z-1500 son de $50 y
$40, respectivamente. Encuentre el número óptimo de cada producto a producir. Si la
corporación está produciendo actualmente 30 unidades del modelo Z-1200 y 12O
unidades del modelo Z-1500, ¿cuánto beneficio está perdiendo o está dejando de ganar?
(Adaptado de Thierauf y Grosse, 1972, pág. 273–274.)
40. El dueño de una editorial está imprimiendo un nuevo libro y tiene las alternativas de
empastarlo con cartón fino y/o cartulina corriente. La venta de un libro con pasta de
cartón fino le genera un beneficio de $30, mientras que la de un libro con cartulina es de
$7.50. El tiempo requerido para empastar un libro con cartón fino es de 3 minutos y con
cartulina ordinaria 2 minutos. El tiempo total disponible para empastado es de 800
minutos. Se estima que las ventas serán de no más de 6000 copias para el libro empastado
con cartón fino y no más de 10000 copias para el empastado con cartulina corriente.
Formule un modelo de programación lineal para este enunciado.
42. Una compañía de paqueteo y envíos desea contratar operadores extra para las festividades
navideñas, pero debido a las limitaciones de espacio el número de empleados no puede
exceder de 10. De experiencias anteriores se sabe que un operador hombre maneja 3000
cartas diarias ó 1000 paquetes diarios, mientras que una mujer puede manejar 4000 cartas
ó 500 paquetes. Sabiendo que al menos llegarán 36000 cartas diarias y 10000 paquetes
por día y que un hombre y una mujer ganan diariamente $250 y $220 respectivamente,
¿cuántos hombres y mujeres deberán ser empleados?. ¿Cómo se afecta la respuesta si el
gobierno ordena el pago de $250 diarios a cada mujer? Suponga que el número de
operadores en condiciones normales es de 2 hombres y 2 mujeres quienes reciben salarios
de $200 y $180 semanal. ¿Qué problemas encuentra usted en la formulación de este
problema? ¿Cómo resolverlos para obtener un modelo satisfactorio?
Actualmente se dispone de 100 unidades MT1 y 200 unidades de MT2. Se desea conocer
en que forma se puede utilizar la información disponible para decidir el esquema de
producción de AS1 y BC1 para tener el máximo número de productos Z terminados.
44. Una compañía posee tres plantas, en las cuales elabora un componente pequeño para un
producto industrial. La compañía comercializa el producto a través de cinco distribuidores
en el país y el precio al distribuidor es de $75 la unidad incluyendo el costo de
distribución. Los pronósticos de ventas indican que los requerimientos mensuales por
distribuidor son los siguientes:
Distribuidor 1 2 3 4 5
Demanda mensual 2700 2700 9000 4500 3600
Planta 1 2 3
Capacidad 4500 9000 11250
Costo Unitario
de producción $ 60 30 54
67
Distribuidor 1 2 3 4 5
Planta 1 1.5 2.1 3.3 4.5 4.8
Planta 2 2.4 1.8 3.0 3.6 4.5
Planta 3 3.0 2.7 2.7 3.0 4.8
45. Una empresa papelera posee un centro de recolección de desperdicios los cuales somete a
diferentes tratamientos de tal manera que pueda producir materia prima para la venta. De
acuerdo a la mezcla de los materiales utilizados es posible producir tres tipos de calidades
diferentes del producto. Para la mezcla existe cierta flexibilidad y se han especificado
estándares de calidad que indican los niveles máximo y mínimo en porcentaje (por peso)
de los materiales que se permiten en cada tipo de producto. Las especificaciones se dan en
la tabla siguiente junto con el costo de amalgamado por libra y el precio de venta por
libra.
La compañía tiene que determinar cuanto debe producir de cada tipo de producto y la
mezcla exacta de materiales que debe utilizar para cada tipo de tal manera que se
68 Introducción a la Modelación Matemática y Optimización
maximice el beneficio total por semana (ventas totales menos costos totales de
amalgamamiento y tratamiento).
PROYECTO
CONTRATISTA
Pl P2 P3
C1 28 32 36
C2 36 28 30
C3 38 34 40
Formule un modelo de programación lineal que permita determinar como asignar los
proyectos para minimizar los costos totales de todos ellos. Se asume que a cada contratista
se le asignará un solo proyecto.
48. Cierta empresa es un contratista grande que realiza trabajos en techos. Puesto que el
precio de las tejas varía con las estaciones del año, la compañía trata de acumular
existencias cuando los precios están bajos y almacenarlas para su uso posterior. La
compañía cobra el precio corriente en el mercado por las tejas que ínstala, sin importar
cuando las haya adquirido. La tabla mostrada abajo refleja lo que la compañía ha
proyectado como costo, precio de venta y demanda para las tejas durante las próximas
cuatro temporadas. Cuando las tejas se compran, se incurre en un costo de manejo de 6
u.p. por cada 1000 piezas, así como también en un costo de almacenamiento de 12 u.p.
por cada 1000 piezas por cada temporada en la que se almacenan (suponga que se incurre
en el costo de almacenamiento solo para las tejas que se almacenan para períodos su venta
en períodos posteriores).
Lo máximo que se puede guardar en el almacén son 220.000 piezas; Esto incluye el
material que se compra para utilizarlo en el mismo periodo. La compañía ha fijado una
política que señala que no se conservan materiales más de cuatro temporadas. Plantee un
modelo de programación lineal para que la compañía pueda maximizar sus ganancias
netas para las cuatro temporadas.
49. Una compañía compra y vende maíz en efectivo. Posee una bodega con capacidad de
1.000 toneladas. El inventario inicial a Enero 01 es de 200 toneladas y $7 millones en
caja. El precio estimado del maíz por tonelada para el primer trimestre es el siguiente:
70 Introducción a la Modelación Matemática y Optimización
PRECIO DE PRECIO DE
MES COMPRA VENTA
($/ton) ($/ton)
Enero 8.550 9.300
Febrero 9.150 9.750
Marzo 8.700 8.850
El maíz es entregado en el mes de compra y no puede ser vendido sino hasta el mes
siguiente o los meses que siguen. La compra y venta se hace estrictamente al contado
contra-entrega. La compañía desea tener un inventario final de 400 toneladas al terminar
el trimestre. Formule un modelo de PL que le permita determinar las políticas de compra
y venta en cada mes que maximicen la ganancia neta total en el trimestre. Ignore los
posibles costos de inventario. (Adaptado de Moskowitz y Wright, 1982, pág. 289–290.)
50. Una empresa pública debe llevar a cabo cinco proyectos, pero solo han cotizado tres
contratistas para realizar los cinco proyectos. Las cotizaciones se muestran en la tabla
siguiente:
PROYECTO
CONTRATISTA
Pl P2 P3 P4 P5
C1 65 37 42 29 29
C2 59 39 50 29 31
C3 62 46 33 24 31
Formule un modelo de Programación lineal para cada uno de los casos planteados a
continuación, suponiendo que siempre se busca minimizar el costo total de realización de
los proyectos respectivos.
a) Cada contratista puede llevar a cabo máximo dos proyectos y todos los proyectos
deben ser ejecutados.
b) Cada contratista solo puede ejecutar un solo proyecto y, por lo tanto, solo podrán
llevarse a cabo tres de los cinco proyectos.
c) No hay límite en el número de proyectos que cada contratista puede llevar a cabo y
deben ejecutarse todos los proyectos.
51. Una gran empresa constructora posee 2000 hectáreas de tierra de primera clase, pero no
urbanizada. En el pasado no existía mucha regulación a nuevas urbanizaciones en el sitio
donde se encuentra la tierra, pero hoy en día las cosas han cambiado. Debido a la falta de
desagües por alcantarillado, se utilizan muchos tanques sépticos, la mayoría instalados en
forma inadecuada. Como las tierras rodean un lago, con el paso de los años, la filtración
de los tanques sépticos ha provocado un severo problema de contaminación de agua.
Para controlar la degradación más profunda en la calidad del agua, los funcionarios del
municipio presentaron y aprobaron algunos reglamentos estrictos aplicables a todas las
urbanizaciones del futuro:
71
a) Solo se pueden construir casas para una, dos y tres familias, donde las unifamiliares
constituyen cuando menos el 50% del total.
b) Para limitar el número de tanques Sépticos, se requieren tamaños de lote mínimos de
5, 7.5 y 10 hectáreas para casas de una, dos y tres familias.
c) Se deben establecer áreas de recreo de 2.5 hectáreas cada una, a razón de un área para
cada 200 familias.
d) Para preservar la ecología del lago, no se puede extraer agua del subsuelo para uso en
la casa o el jardín.
53. Usted dispone de cinco proyectos en los cuales puede invertir. Los proyectos se realizarán
totalmente o no se realizarán, pero no se podrán emprender parcialmente. Los datos de
flujo de dinero se muestran en la tabla siguiente, donde un valor negativo representa una
salida de dinero o inversión en el proyecto, y un valor positivo representa un retorno o una
entrada de dinero producto del proyecto:
AÑO Proyectos:
1 2 3 4 5
0 $1.000 $2.500 $3.000 $3.300 $4.000
1 $500 $1.000 $2.000 $800 $1.000
2 $800 $1.500 $1.000 $1.500 $2.000
3 $1.700 $1.800 $4.000 $1.500 $5.000
4 $1.700 $2.000 $2.000 $1.000 $0
5 $0 $1.500 $5.000 $500 $0
54. La puerta de un vehículo consta de tres partes de metal cuya pintura se puede hacer
mediante tres diferentes procesos piscina de pintura, pistola o cuarto de baño de pintura.
Se dispone de la siguiente información:
55. Una firma está planeando la construcción de nuevos edificios en cuatro lugares de una
ciudad (1, 2, 3 y 4). En cada uno de los lugares puede haber tres posibilidades para el
diseño de los edificios (A, B y C). También existe la opción de no construir ningún
edificio en cualquiera de los lugares 1, 2, 3 y 4. La tabla siguiente muestra cada una de las
opciones, su inversión y su ingreso anual:
Opción A1 B1 C1 A2 B2 C2 A3 C3 A4 B4 C4
Inversión 12 20 24 30 69 45 10 28 42 50 55
Ingreso 2 4 6 5 7 11 5 12 16 19 22
Anual
56. Una cadena de tiendas pequeñas está planeando entrar a un nuevo mercado y debe
determinar donde localizar nuevas tiendas en cierta zona. La figura muestra un mapa de
las principales calles en la zona de interés.
A B C D E
1 1 2
3 4
2
5 6
3
7 8 9
4
10
5
1 milla
1 milla
Las calles adyacentes están separadas 1 milla. Las localizaciones, marcadas de 1 a 10, son
los lugares potenciales donde se pueden ubicar nuevas tiendas. La distancia entre un punto
y otro se mide con base en la red de calles. Por ejemplo, la distancia entre el punto 4 y el
punto 6 es de 3 millas; entre el punto 1 y el punto 8 es de 4 millas. El costo de
construcción de una tienda depende del lugar donde se haga, de acuerdo a la tabla
siguiente:
Localización 1 2 3 4 5 6 7 8 9 10
Costo 100 80 90 50 80 90 100 70 90 80
a) No se pueden abrir dos tiendas sobre la misma calle. Por ejemplo, no se pueden abrir
las tiendas 5 y 6 simultáneamente, ni por ejemplo las tiendas 4 y 8 simultáneamente.
b) Las tiendas deben estar localizadas de tal forma que estén separadas entre sí por 3
millas o más. Por ejemplo, las tiendas 4 y 5 no podrían abrirse simultáneamente por
este motivo.
c) Todo punto en el que se intercepten dos calles (o sea A1, B1, ... hasta E5.) no debe
estar alejado más de 3 millas de alguna de las nuevas tiendas (3 millas es aceptable,
pero más no).
57. Una compañía manufactura tres productos, A, B y C. Cada unidad de A requiere 1 hora de
servicio de ingeniería, 8 horas de mano de obra directa y 4 libras de material. Para
producir una unidad de producto B se necesitan 3 horas de ingeniería, 3 horas de mano de
obra directa y 3 libras de material. Cada unidad de producto C requiere de 2 horas de
ingeniería, 4 horas de mano de obra directa y 2 libras de material. Se dispone de 80 horas
de ingeniería, 800 horas de mano de obra directa y 300 libras de material cada mes. Las
utilidades unitarias de cada producto se han determinado como muestra la tabla siguiente:
58. Una fábrica necesita cumplir un contrato de suministro de 1200 puertas dentro de un plazo
de cinco semanas, así:
Al final de la semana 1 2 3 4 5
Núm. de puertas a entregar 100 200 300 400 200
Formule un modelo de programación lineal que permita planear la producción a costo mínimo.
59. Un ingenio produce azúcar morena, azúcar refinada, azúcar pulverizada y mieles con el
jarabe de la caña de azúcar. La compañía puede procesar hasta 4000 toneladas de jarabe a
la semana y tiene un contrato para entregar un mínimo de 25 toneladas semanales de cada
tipo de azúcar y mieles. El proceso de producción se inicia fabricando azúcar morena y
mieles con el jarabe. Una tonelada de jarabe produce 0.3 toneladas de azúcar morena y 0.1
toneladas de mieles. Las mieles así obtenidas están listas para la venta. Después, el azúcar
76 Introducción a la Modelación Matemática y Optimización
Las utilidades netas de azúcar morena, azúcar refinada, azúcar pulverizada y mieles son
de 150, 200, 230 y 35 $/tonelada, respectivamente. Formule un modelo de programación
lineal que le permita determinar el programa de producción semanal de cada tipo de
azúcar y de las mieles, para obtener la máxima utilidad neta. (Taha, 1998, pág. 55)
60. La tabla siguiente muestra la demanda requerida de cierto producto por los próximos doce
meses.
MES 1 2 3 4 5 6 7 8 9 10 11 12 TOTAL
Demanda 10 62 12 13 15 129 88 52 12 160 238 41 1200
0 4 4
61. A su empresa ha llegado un pedido de 3500 unidades de cierto producto, el cual puede
usted manufacturar en cualquiera de cuatro máquinas, con las siguientes características:
62. Una planta manufacturera debe cumplir con la siguiente demanda en unidades en los
próximos 12 meses:
Los cambios en la fuerza laboral, bien sea en su aumento o disminución, no pueden ser
mayores que 40 trabajadores por mes, debido al convenio realizado con la empresa de
empleos temporales. Adicionalmente, estos cambios cuestan $300 por cada empleado
contratado y $420 por cada empleado que se despida, así sea en forma temporal.
El trabajo en horas extras es también limitado. Cada trabajador puede producir hasta 20
unidades/mes en tiempo normal de trabajo, pero no puede producir más de 6
unidades/mes en tiempo extra. Además, los costos por unidad producida en hora extra
exceden en $20/unidad al costo de producción en tiempo normal.
63. Una corporación ha decidido producir tres productos nuevos. En este momento cinco de
sus plantas tienen exceso de capacidad de producción. El costo unitario de fabricación del
primer producto será de 31, 29, 32, 28 y 29 $/unidad en las plantas 1, 2, 3, 4 y 5,
respectivamente. El costo de fabricación del segundo producto será de 45, 41, 46, 42 y 43
$/unidad en las plantas respectivas 1, 2, 3, 4 y 5. El costo unitario de fabricación del tercer
producto será de 38, 35 y 40 $/unidad en las plantas 1, 2 y 3, respectivamente, mientras
que las plantas 4 y 5 no pueden fabricar este producto. Los pronósticos de ventas indican
que la producción diaria de cada uno de los tres productos debe ser de 6000, 10000 y 8000
unidades. Las plantas 1, 2, 3, 4 y 5 tienen capacidades para producir 4000, 6000, 4000,
10000 y 3500 unidades diarias, independientemente del producto o combinación de
productos que se quiera. Supóngase que cualquier planta que tiene capacidad y puede
fabricarlos, podrá producir cualquier combinación de productos en cualquier cantidad.
Formule este problema como un modelo de transporte para determinar cómo asignar los
78 Introducción a la Modelación Matemática y Optimización
nuevos productos a las plantas con el fin de minimizar el costo total de fabricación.
(Adaptado de Hillier y Lieberman, segunda edición, 1989, pág. 227.)
64. Supóngase que Inglaterra, Francia y España producen todo el trigo, cebada y avena del
mundo. La demanda mundial de trigo requiere que se dediquen 125 millones de acres a la
producción de este cereal. De igual manera, se necesitan 60 millones de acres para cebada
y 70 millones de acres para avena. La cantidad total de tierra disponible en los tres países
es de 70 millones de acres en Inglaterra, 110 millones de acres en Francia y 80 millones
de acres en España. El número de horas de mano de obra necesarias en Inglaterra, Francia
y España para producir un acre de trigo es 18, 13 y 16 horas, respectivamente. La
producción de un acre de cebada requiere 15 horas de mano de obra en Inglaterra, 12 en
Francia y 12 en España. El número de horas de mano de obra necesarias para producir un
acre de avena es de 12 en Inglaterra, 10 en Francia y 16 en España. El costo de mano de
obra por hora en los tres países es $3.00, $2.40 y $3.30 para la producción de trigo, $2.70,
$3.00 y $2.80 para la producción de cebada y $2.30, $2.50 y $2.10 para la producción de
avena. El problema es asignar el uso de la tierra en cada país de manera que se cumpla
con los requerimientos de alimentación en el mundo y se minimice el costo total de mano
de obra. Formule este problema como un problema de transporte construyendo la tabla de
costos y requerimientos apropiada. (Adaptado de Hillier y Lieberman, segunda edición,
1989, pág. 228.)
66. Considere la tabla de la página siguiente, la cual muestra la información básica para
construir un problema real de dieta humana, y resuelva el caso que se presenta a
continuación.
79
Vamos a asumir que se trata de una dieta promedio con los siguientes requerimientos
mínimos diarios:
6
Desarrollado con base en la información mostrada en Chvátal, Vasek, Linear Programming, W. H. Freeman
and Company, New York, 1983, pág. 3–5, 182–187.
80 Introducción a la Modelación Matemática y Optimización
Igualmente, definimos los tamaños de las porciones de cada alimento de acuerdo con la
tabla siguiente:
Manzanas 1 mediana
Brócoli ½ taza
Zanahorias 1 mediana
Naranjas 1 mediana
Jugo de naranja 1 taza
Papas 1 mediana
Espinaca ½ taza
Tomates 1 mediano
Pan 2 tajadas
Cereal (hojuelas) 100 gr
Avena 120 gr
Pasta 140 gr
Arroz 100 gr
Mantequilla 1 cucharada
Azúcar ½ taza
Manteq. de maní 1 cucharada
a) Adicionar a la lista de alimentos presentada cinco (5) nuevos alimentos a los cuales
pueda conseguirle la información nutricional y su correspondiente costo, por porción.
(Se supone que cada grupo de trabajo del curso encuentre información diferente)
b) Investigar los costos de los alimentos mostrados en la tabla anterior en nuestro medio,
por porción (Valores aproximados). Recuerde que para formular un problema real,
81
este costo debe calcularse por porción de cada alimento. Por ejemplo, si una libra
(500 gr) de carne para freír cuesta aproximadamente $8000, entonces el costo de cada
porción de 100 gr sería de 8000/5 = $1600.
c) Formular un modelo de programación lineal que le permita encontrar una dieta diaria
de costo total mínimo, que satisfaga todos los requerimientos mínimos diarios
mostrados anteriormente. Este modelo puede hacerlo en forma simbólica y compacta,
definiendo nombres para los datos anteriores. En otras palabras, para formular el
modelo no se necesita escribir explícitamente cada término de la función objetivo ni
cada restricción, sino su representación compacta en forma de sumatorias.
e) Analizar los problemas de sabor y gusto de las personas para replantear el modelo
formulado de tal forma que raciones de alimentos seleccionados siempre aparezcan
en la dieta. Resolverlo nuevamente y analizar la solución.
82 Introducción a la Modelación Matemática y Optimización
En la mayoría de los casos, el método que más se aplica es el simplex. Sin embargo,
para lograr la comprensión del método simplex es necesario estudiar inicialmente dos formas
de solución: El método gráfico y el método de enumeración de soluciones básicas. Existen
otros métodos denominados algoritmos de punto interior que no se estudian aquí.
Este método se ilustrará encontrando la solución del modelo Nº 1 anterior. Este modelo
servirá de base para el estudio y la comparación de todos los métodos.
Lo primero que hay que observar es que el modelo al cual se llegó, tiene sólo dos
variables de decisión, (X1 y X2). Esto hace pensar que es posible representar las restricciones
en un eje de coordenadas (X1, X2). El modelo encontrado anteriormente es el siguiente:
Sujeto a:
X2), las cuales las satisfacen simultáneamente. Este conjunto de parejas son las soluciones
factibles y forman la denominada región factible (Figura 3.1)
negatividad R4 y R5, las cuales limitan la región al primer cuadrante. La intersección de todos
los semiplanos definidos por las desigualdades constituye la región factible.
Los infinitos puntos en los límites y dentro de la región factible son las soluciones
factibles del modelo. Por lo general, sólo una de estas infinitas soluciones es la solución
óptima, o sea aquella que produce el valor óptimo de la función objetivo (máximo o mínimo).
Por ejemplo, se observa claramente que X1 = 400 y X2 = 200 es una solución factible,
pues está dentro de la región factible. Esta solución produciría un valor de la función objetivo
U = $300.000. Otra solución factible es X1 = 500 y X2 = 300, la cual produciría U = $410.000.
O sea que esta última es mejor que la primera, pues produce un mayor valor de la función
objetivo “U”, ya que ésta se está maximizando.
Se puede demostrar que la región factible definida por restricciones lineales es una
región convexa. En una región convexa, si se escogen cualquier par de puntos dentro o en los
límites de la región y se unen por una línea recta, entonces toda la línea recta queda
comprendida dentro de la región. La propiedad de convexidad es un aspecto fundamental para
la caracterización de los problemas de optimización. En general, puede afirmarse que un
problema convexo (programación convexa) es más fácil de resolverse que otro que no lo es.
Los vértices de la región factible son los interceptos de los segmentos de reata límites
de la región. En este caso los vértices son:
V1 (0,0)
V2 (0,600)
V3 (350,450)
V4 (700,200)
V5 (900,0)
Todos los puntos de esta recta que queden dentro de la región factible son soluciones
factibles que producen obviamente un valor de U = 280.000 (por ejemplo X1 = 350; X2 = 200).
Esta recta es obviamente paralela a la primera, pues tiene la misma pendiente igual a
4/7 y, de nuevo, todos los puntos de ella que están dentro de la región factible son soluciones
factibles que producen un valor de U = 420.000.
O sea que la función objetivo “U” aumenta de valor a lo largo de una familia de rectas
paralelas con pendiente –4/7. Es aquí donde se puede concluir la forma de encontrar la
85
X1 = 350
X2 = 450
U* = 455.000 (U* = U óptima)
Análisis de la solución:
involucran parejas de ellas; por ejemplo, el siguiente modelo puede ser resuelto utilizando el
método gráfico:
Sujeto a:
Y se utiliza el método gráfico independientemente para cada Zi, con sus respectivas
restricciones.
En el caso de tres variables de decisión, aún puede ser utilizado el método gráfico en
algunas ocasiones donde resulte fácil apreciar la región factible y los “planos objetivos”. En
este caso la región factible es un poliedro convexo, formado por el corte de los diferente
planos que constituyen las restricciones; el crecimiento de la función objetivo se puede
observar mediante el desplazamiento de planos paralelos. De nuevo, el último punto del
poliedro que toque el plano al abandonar la región factible será la solución óptima.
Sujeto a:
Sin embargo, el método gráfico resulta inadecuado en la mayoría de los casos de tres
variables de decisión, pues se necesita gran habilidad en el dibujo y una gran imaginación
espacial.
87
De los dos ejemplos resueltos por el método gráfico, puede sacarse una conclusión
muy importante, la cual es válida para cualquier modelo de Programación Lineal: La solución
óptima se encuentra en uno de los vértices de la región factible. Esta afirmación, que puede
ser planteada como un teorema, es de fundamental importancia, pues de infinitas soluciones
posibles (definidas por la región factible), sólo basta con mirar aquellas definidas por los
vértices de la región factible.
Este caso ocurre cuando la recta de la función objetivo es paralela a uno de los lados de
la región factibles. Obviamente, aquí la recta de la función objetivo toca infinitos puntos al
abandonar la región factible. Estos infinitos puntos están en el segmento de la recta que une
dos vértices determinados (sin embargo, también en este caso la solución óptima esta en los
88 Introducción a la Modelación Matemática y Optimización
vértices de la región factible). Todos ellos son solución óptima y producen el mismo valor de
la función objetivo.
En este caso una ó varias variables pueden crecer (ó decrecer) sin restricción,
produciendo valores cada vez mayores (ó menores) en la función objetivo. Es el caso en el que
la región factible es “abierta”, probablemente debido a que faltó incluir alguna restricción en el
modelo.
Sujeto a:
Ejercicio: Dibuje la región factible del modelo anterior y compruebe que es “abierta”.
Incluya una restricción adicional que permita tener solución óptima única.
En este caso no existe en realidad una región factible, pues la “intersección de las
restricciones” es el conjunto vacío. Esta situación se presenta cuando hay inconsistencia en las
restricciones y el modelo debe ser replanteado. Por ejemplo el modelo:
89
Sujeto a:
Ejercicio: Trate de dibujar la región factible del modelo y compruebe que ésta no existe
y que por lo tanto el modelo no tiene solución óptima.
Función Objetivo:
90 Introducción a la Modelación Matemática y Optimización
Sujeto a:
Igualdades Restrictivas
En el modelo anterior los bi deben ser mayores ó iguales a cero. Por lo general, se
cumple que m < n y, por lo tanto, el sistema de igualdades tiene infinitas soluciones (las
infinitas soluciones de la región factible).
Se transforma a:
Se puede transformar a:
Se transforma en:
91
Las variables de holgura, obviamente, deben cumplir con las restricciones de nó-
negatividad.
c. Si alguno de los bi es menor que cero, entonces se multiplica por (-1) a ambos
lados de la desigualdad (o igualdad) y posteriormente se le suma (ó se le resta) la variable de
holgura (ó de exceso), dependiendo del signo resultante de la desigualdad. Si se trata desde un
comienzo de una igualdad, lógicamente no es necesario adicionar variable de holgura ó
exceso.
Xj es variable libre (Xj puede ser mayor, igual ó menor que cero)
Sujeto a:
f. Donde aparezca la variable X4, debe cambiarse por (X4+ X4), X4+ y X4 mayores
ó iguales que cero.
g. La función objetivo se cambia a maximización.
Sujeto a:
Con base en lo anterior, se pueden hacer las siguientes definiciones sobre cualquier
modelo de programación lineal en forma estándar:
Solución Factible:
Una solución factible es cualquier solución del sistema n × m que satisfaga las
condiciones de no-negatividad (o sea aquella solución en la cual todas las variables son
positivas ó cero).
Solución Básica:
se les denomina variables básicas y se dice que constituyen una base del sistema de
ecuaciones.
Es una solución básica que cumple con las condiciones de no-negatividad (en este caso
las variables básicas son positivas. (Si existe por lo menos una variable básica igual a cero, se
dice que la solución es degenerada)
Solución Óptima:
Con base en las definiciones anteriores, se pueden enunciar los tres principales
teoremas de optimalidad (sin demostración), a saber:
a. Las soluciones factibles conforman un conjunto convexo, cuyos vértices son las
soluciones básicas factibles.
b. Si el sistema tiene una solución factible, entonces existe por lo menos una solución
básica factible.
c. Si la función objetivo tiene un óptimo (máximo ó mínimo) finito, entonces existe
por lo menos una solución óptima la cual es una solución básica factible.
O sea, que de las infinitas soluciones factibles, “sólo” basta con inspeccionar las
básicas factibles y una de éstas será una solución óptima del problema de PL. Igualmente, de
la afirmación anterior queda claro por qué se decía en el desarrollo del método gráfico que la
solución óptima estaba en uno de los vértices de la región factible.
Dado que el método de enumeración de soluciones básicas debe inspeccionar todas las
soluciones básicas, identificar la factible y escoger la mejor, debe primero determinarse el
número posible de soluciones básicas. Como en el sistema de “m” igualdades con “n”
variables, “sobran” (n – m) variables a las cuales se les asigna el valor de cero, el número
94 Introducción a la Modelación Matemática y Optimización
posible de soluciones básicas resulta ser igual a todas las posibles formas de escoger las (n –
m) variables de las n totales, o sea:
Aquí ya pueden observarse las limitaciones de este método, pues supóngase que se
tiene un 30 modelo
de 30 variables con 15 restricciones (el cual es un modelo modesto y
relativamente
15 pequeño), entonces tendrían que ser resueltos = 155,117,520 sistemas de
ecuaciones lineales de 15 × 15 para determinar la solución óptima. Por lo tanto, la importancia
de estudiar este método radica simplemente en la formación de los fundamentos para la
comprensión del método simplex.
Sujeto a:
Sujeto a:
Obsérvese que como las restricciones son todas , sólo fue necesario sumar a cada
restricción su correspondiente variable de holgura.
Solución
X1 X2 X3 X4 X5 Factible? Valor de U
Básica
1 0 0 1400 980 900 x 0
2 0 600 0 140 300 x 420.000
3 0 700 233 1/3 0 200
4 0 900 -700 -280 0
5 1400 0 0 -420 -500
6 980 0 420 0 -80
7 900 0 500 80 0 x 360.000
/Solución
8 350 450 0 0 100 x 455.000
óptima
9 525 375 0 -70 0
10 700 200 233 1/3 0 0 x 420.000
La forma de asignar el valor de cero a dos variables en cada solución debe ser lógica y
ordenada, con el objeto de no repetir ó dejar de tener en cuenta algunas de las soluciones
básicas.
Obsérvese que las soluciones básicas No. 3, 4, 5, 6, y 9 no son factibles, pues por lo
menos una de las variables no cumple con las restricciones de no-negatividad. Además, es
importante notar que las soluciones 1, 2, 7, 8 y 10 son factibles y corresponden precisamente a
los vértices V1, V2, V5, V3 y V4, respectivamente, de la región factible mostrada en la figura
3.1 anterior. La solución óptima es la No. 8, pues es la que produce el valor máximo de la
función objetivo U. La solución óptima completa viene dada por:
X1 = 350
X2 = 450
X3 = 0
X4 = 0
X5 = 100
U* = 455.000
Sujeto a:
Sujeto a:
Solución Valor de
X1 X2 X3 X4 X5 Factible?
Básica U
1 0 0 0 10 -18
2 0 0 2 0 ? Sist. Inconsistente
3 0 0 18 10 0 x 36
4 0 10 0 0 12 x 10
5 0 6 0 4 0 x 6
6 0 10 -12 0 0
7 10 0 0 0 2 x 30
8 9 0 0 1 0 x 27
9 10 0 -2 0 0
10 12 -2 0 0 0
La tabla anterior sugiere que la solución óptima del modelo es la solución básica No. 3,
pero en realidad, como se demostró anteriormente, la función objetivo es no acotada y, por lo
tanto, el modelo no tiene solución óptima. Por esta razón, este método puede fallar en estos
casos.
En la página siguiente se presenta un taller para la preparación del primer examen del
curso. Se sugiere su solución completa.
97
Ejercicio 3.2.
1. Una refinería de petróleo puede comprar dos tipos de petróleo crudo. El tipo 1 cuesta
$57.500/barril y el tipo 2 cuesta $50.600/barril. Cada barril de petróleo crudo de cualquier
tipo ya refinado de produce tres tipos de productos terminados: gasolina para avión,
gasolina normal y kerosén. Un barril de petróleo crudo tipo 1 produce 0.45 barriles de
gasolina para avión, 0.18 barriles de gasolina normal, 0.30 barriles de kerosén y el resto
son desechos cuya eliminación cuesta $2.300/barril. Por otra parte, un barril de petróleo
crudo tipo 2 produce 0.35 barriles de gasolina para avión, 0.36 barriles de gasolina
normal, 0.20 barriles de kerosén y el resto son desechos cuya eliminación cuesta
$3.450/barril. La refinería ha firmado un contrato para proveer por lo menos 1.260.000
barriles de gasolina para avión, 900.000 barriles de gasolina normal y 300.000 barriles de
kerosén. Formule un modelo de programación lineal que le permita determinar cuántos
barriles de petróleo crudo de cada tipo debe comprar la refinería para cumplir con la
demanda establecida a costo total mínimo.
2.
a) Resuelva el problema anterior mediante el método gráfico e interprete la solución.
b) Replantee la función objetivo del modelo asumiendo que los precios de venta de cada
producto son los siguientes: Gasolina para avión $134.600/barril, gasolina normal
$75.000/barril y kerosén $40.500/barril. Asuma que ahora las demandas a satisfacer
son máximas. Resuelva de nuevo el problema mediante el método gráfico e interprete
la solución obtenida.
3. Una gran empresa distribuidora tiene que suministrar 800 toneladas/mes de productos a
sus clientes y está pensando en renovar completamente su flota de transporte. La empresa
tiene un presupuesto disponible de $1.150 millones (incluyendo los ingresos por venta del
equipo actual) para invertir en tres tipos de camiones. La tabla siguiente muestra la
capacidad, el costo de compra, el costo operativo y el número máximo de viajes/mes de
cada tipo de camión.
4. Una empresa productora estima con muy buena precisión que la demanda para las
próximas cuatro semanas es de 600.000 tornillos pequeños y 400.000 tornillos grandes.
Estos tornillos pueden producirse en dos máquinas distintas. La máquina 1 está disponible
56 horas/semana y la máquina 2 está disponible 84 horas/semana. Los requerimientos de
costos y tiempo para producir cada tamaño de tornillo en cada máquina y el precio de
venta de cada tamaño de tornillo se muestran a continuación:
Por motivos de sabor, la cantidad mínima de cada ingrediente en el dulce debe ser de 0.50
onzas.
7. Se está analizando cómo combinar dos tipos de fertilizantes, de tal forma que la mezcla
contenga al menos 200 unidades de un cierto compuesto. El fertilizante A contiene 8
unidades del compuesto/Kg de fertilizante, mientras que el fertilizante B contiene 25
unidades del compuesto/Kg de fertilizante. El fertilizante A cuesta $1.000/Kg, mientras
que el fertilizante B cuesta $3.125/Kg. Ambos fertilizantes se producen a partir de dos
materias primas M1 y M2. Cada Kg de fertilizante A consume 6 litros de M1 y 3 litros de
M2. Cada Kg de fertilizante B consume 12 litros de M1 y 2 litros de M2. Se dispone de
216 litros de M1 y 60 litros de M2. Finalmente, por razones de calidad, la mezcla debe
contener a lo más 16 Kg de fertilizante A.
4. EL MÉTODO SIMPLEX
Los autores presentan el método simplex en diversas formas: tabular, algebraica,
matricial, etc. De todas estas maneras de presentación, tal vez la más adecuada y la que más
ventajas ofrece para sentar en el estudiante una mejor base teórica, es la forma matricial. Sin
embargo, posteriormente se hará énfasis en las aplicaciones resueltas por computador, ya que
es en últimas el instrumento que se utiliza para resolver cualquier modelo de Programación
Lineal.
M aximizar Z 5X1 3X 2
sujeto a :
3X1 5X 2 S1 15 (3.2)
5X1 2X 2 S 2 10
(X1 , X 2 , S1 , S 2 ) 0
Luego, el modelo en su forma estándar presenta n = número de variables = 4 y m =
número de restricciones = 2.
n 4
El número máximo de soluciones básicas = 6
m 2
La tabla siguiente muestra las soluciones básicas y la solución óptima:
Solución
Básica No X1 X2 S1 S2 Factible Valor de Z
1 0 0 15 10 Si 0
2 0 3 0 4 Si 9
3 0 5 -10 0 No –
4 5 0 0 -15 No –
5 2 0 9 0 Si 10
235/19
6 20/19 45/19 0 0 Si
12.4
Por lo tanto, la solución óptima es la No. 6. La tabla anterior debe tenerse presente para
el desarrollo que sigue.
102 Introducción a la Modelación Matemática y Optimización
M AX (M IN) Z CX
sujeto a :
(3.3)
AX b
X0
donde,
X1
X
2
.
C C 1 C 2 C 3 ·······
C n 1n X (incluye variables de holgura y/ó exceso Sk)
.
.
X n n1
C 5 0; X X 1 X 2 S1 S2 ;
T
3 0
3 5 1 0 15
A ; b
5 2 0 1 10
103
A las variables que quedan formando un sistema de ecuaciones lineales m×m, se les
denomina VARIABLES BÁSICAS, ya que sus columnas forman la matriz base, así:
1 0 S1 15
0 1 S 10
2
1
S1 1 0 15 15
S 2 0 1 10 10
La cual corresponde a la solución básica No 1 (en este caso se trata de una solución
básica factible).
3 0 X 1 15
5 1 S 10
2
1
X 1 3 0 15 1 / 3 0 15 5
S 5 1 10 5 / 3 1 10 15
2
XB
A B R X C CB C R
X R
M AX (M IN) Z C B X B C R X R
sujeto a :
(3.4)
BX B RX R b
X0
En palabras, lo que hace el método SIMPLEX es partir de una solución básica inicial
factible, con una base inicial asociada, la cual, en general, es la matriz idéntica de orden m.
Posteriormente, se pasa a otra solución básica factible, cambiando una sola columna de la
base, de tal forma que se logre hacer crecer lo más posible a la función objetivo Z, si se está
maximizando, o decrecer, si se está minimizando.
105
El paso a otra solución básica factible se logra cambiando de base, de tal forma que
una de las variables básicas sale de la base, pasando a ser no básica, y una de las variables no
básicas entra a la base, volviéndose básica y entrando a ocupar el lugar de la variable que salió
de la base. En otras palabras, la diferencia entre dos bases sucesivas está en una sola columna.
El proceso de cambio de base continúa hasta que se detecta que no se puede mejorar
más la función objetivo (minimizando o maximizando) y es cuando se ha logrado encontrar la
solución óptima del problema.
Z C B X B C R X R (3.5)
BX B RX R b (3.6)
X B YX R X B (3.7)
CB X B CB YX R CB X B
Z CB X B CB YX R (3.8)
Z Z CR X R CB YX R
106 Introducción a la Modelación Matemática y Optimización
O mejor:
Z Z (CB Y CR )X R (3.9)
Así, la ecuación (3.9) muestra el cambio que tendría la función objetivo en una
iteración, al cambiar de base. Por lo tanto se puede tratar de controlar este cambio para que sea
lo más grande posible. Se puede entonces modificar esta ecuación observando que se trata de
una ecuación escalar. La notación se modifica así:
Los términos C B Y j son los productos escalares entre el vector C B y las columnas de la
matriz Y, y se le denomina Z j . Así, la ecuación (3.9) puede transformarse a:
Z Z ( Z j C j )X j (3.10)
j
Como sólo una de las variables no básicas X j es la que va a entrar a la base y tomará
seguramente un valor positivo, entonces la sumatoria se reduce a un sólo término, ya que el
resto de variables no básicas seguirá siendo igual a cero. Por lo tanto, le ecuación (3.10) se
transforma a:
Z Z ( Z k C k )X k (3.11)
7
En realidad puede entrar a la base cualquier variable cuyo (Z jCj) sea negativo para maximización, pero si se
escoge el menor es más probable llegar al óptimo (máximo en este caso) más rápidamente. Un razonamiento
semejante se presenta si se está MINIMIZANDO.
107
Como debe garantizarse que la próxima solución básica sea factible, entonces:
x1 y1k xk 0
x y x 0
2 2k k
xm -ymk xk 0
Por lo tanto:
x
xk mín s ; ysk 0
ysk
y sk debe ser mayor que cero, puesto que xk debe tomar también un valor no negativo.
Por otra parte, como la variable que va a salir de la base pasará a tomar valor cero,
entonces:
xs xs ysk xk 0 , donde xs es la variable a salir de la base. Así:
x
xk s ; y sk 0
y sk
108 Introducción a la Modelación Matemática y Optimización
Por lo tanto, la variable que entra a la base toma el valor del cociente mínimo de la
variable que sale de la base.
Este cociente se efectúa entre los valores actuales de la solución y los términos de la
columna k de la matriz Y que sean positivos solamente (recuérdese que la columna k es la
columna asociada a la variable xk , previamente escogida para entrar a la base). El valor que
toma la variable que va a entrar a la base es precisamente el del cociente mínimo, o sea
xk mín . Es importante notar que este criterio es válido tanto para problemas de
maximización como para problemas de minimización, pues su origen es el de obtener una
solución básica que sea factible.
Si La solución
¿Se cumple el CRITERIO DE básica factible
PARADA? actual es una
SOLUCIÓN
ÓPTIMA
No
1 0 15 3 5 S X
B ; C B 0 0; C R 5 3; b ; R ; X B 1 ; X R 1
0 1 10 5 2 S2 X2
S1 S2 X1 X2
Entonces:
1
11 0 15 15
XB B b
0 1 10 10
Z CB X B 0
8
Los ( Z j C j ) de las variables básicas siempre son iguales a cero, como se verá posteriormente. Si existe
algún ( Z j C j ) de una variable no básica igual a cero, se concluye que existen múltiples soluciones óptimas,
ya que se podrá cambiar de base sin modificar el valor de la función objetivo.
110 Introducción a la Modelación Matemática y Optimización
Z 2 C B Y2 0 05 2 0
T
Z1 C1 0 5 5
ENTRA A LA BASE LA VARIABLE x1 9
Z 2 C2 0 3 3
Criterio de salida:
Recuérdese que X B X B YX R , en este caso:
s 15 3 5 x1
XB 1
s2 10 5 2 x2
Columna asociada a X1 , variable a entrar a la base.
15 10
Luego: mín , mín 5 , 2 2
3 5
Por lo tanto sale de la base la variable S2.
1 3 15 0 5 S S
B ; C B 0 5; C R 0 3; b ; R ; X B 1 ; X R 2
0 5 10 1 2 X1 X2
s1 x 1 s2 x 2
1 3 15 1 35 15 9 10
1
1
XB B b
0 5 10 0 15 10 2
Z CB X B 10
9
Obsérvese que X1 es la variable asociada a la columna Y1.
10
En este caso la inversa de la nueva base no se ha obtenido por operaciones fila elementales, como se hace
normalmente en la forma tabular, que se estudiará más adelante.
111
Ahora:
1 3 / 5 0 5 3 / 5 19 / 5
Y B 1R 2 / 5
(columnas de Y)
0 1 / 5 1 2 1 / 5
S2 X2
Z1 C B Y1 0 5 3 / 5 1 / 5 1
T
Z 2 C B Y2 0 519 / 5 2 / 5 2
T
Z1 C1 1 0 1
Z 2 C2 2 3 1 , ENTRA A LA BASE LA VARIABLE X2 (la
única cuyo Zj – Cj es negativo)
criterio de salida:
s 9 3 / 5 19 / 5 s2
XB 1
x1 2 1 / 5 2 / 5 x2
columna asociada a X2, la variable a entrar a la base.
45 45
Luego, mín , 5
19 19
Por lo tanto, sale de la base la variable S1.
5 3 15 0 1 X S
B ; C B 3 5; C R 0 0; b ; R ; X B 2 ; X R 2
2 5 10 1 0 X1 S1
x2 x1
Ahora:
5 / 19 3 / 19 0 1 3 / 19 5 / 19
Y B 1R (columnas de Y)
2 / 19 5 / 19 1 0 5 / 19 2 / 19
Y1 Y2
Z1 C B Y1 3 5 3 / 19 5 / 19 16 / 19
T
Z 2 C B Y2 3 55 / 19 2 / 19 5 / 19
T
112 Introducción a la Modelación Matemática y Optimización
Z1 C1 16 / 19 0 16 / 19
0 , SE CUMPLE EL CRITERIO DE
Z 2 C2 5 / 19 0 5 / 19
PARADA.
Luego la solución básica factible actual es una solución óptima. Como los únicos (Zj –
Cj) iguales a cero son los de las variables básicas, entonces esta es la única solución óptima del
modelo (3.1).
Cj 5 3 0 0
VARIABLES
CB XB X1 X2 S1 S2 COCIENTE
BÁSICAS
S1 0 15 3 5 1 0 15/3 = 5
S2 0 10 5 2 0 1 10/5 = 2
Zj 0 0 0 0 0 0
Zj C j -5 -3 0 0
S1 0 9 0 19/5 1 -3/5 45/19 2.4
X1 5 2 1 2/5 0 1/5 2/(2/5) = 5
Zj 10 5 2 0 1
Zj C j 0 -1 0 1
X2 3 45/19 0 1 5/19 -3/19
X1 5 20/19 1 0 -2/19 5/19
Zj 5 3 5/19 16/19
Zmáx 235/19
Zj C j 0 0 5/19 16/19
X1 20/19
X 2 45/19
S1 0
S2 0
Z máx 235/19
MAX ( MIN ) Z CX
sujeto a :
(3.12)
AX (, , )b
X0
Si todas las restricciones son de , el método SIMPLEX inicia con la base igual a la
matriz idéntica, formada por las columnas de las variables de holgura. Pero si existe por lo
menos una restricción de = ó de , la matriz idéntica no aparece en forma automática y por lo
tanto debe crearse mediante la adición de variables artificiales, salvo algunas excepciones
(Ver Ejemplos 3.2 y 3.3)
114 Introducción a la Modelación Matemática y Optimización
M aximizar Z 2X1 3X 2
Sujeto a : X1 2X 2 4
X1 X 2 3
X1, X 2 0
La forma estándar de éste modelo es la siguiente:
M aximizar Z 2X1 3X 2
Sujeto a : X1 2X 2 S1 4
X1 X 2 3
X1 , X 2 , S1 0
Como se observa, no aparece la submatriz idéntica de orden 2 en la matriz A, ya que
solo fue necesario adicionar la variable de holgura S1 para obtener la forma estándar del
modelo. Por lo tanto, la submatriz idéntica debe crearse artificialmente así:
M aximizar Z 2X1 3X 2
Sujeto a : X1 2X 2 S1 4
X1 X 2 A 3
X1 , X 2 , S1 0; A variable artificial
Dado que antes de adicionar la variable artificial A ya la igualdad estaba conformada,
la única forma de que el problema original tenga solución es que la variable artificial sea igual
a cero en la solución óptima, lo cual se logra la mayoría de las veces sacándola de la base.
Así, el método debe obligar a que la variable A salga de la base y esto puede conseguirse
penalizándola en la función objetivo con un valor muy grande positivo (si se está
minimizando) o muy pequeño negativo (si se está maximizando). Así, surge el método de la
“Gran M”, el cual transformaría el modelo anterior en:
M aximizar Z 2 X 13 X 2 MA
Sujeto a : X 1 2 X 2 S 4 11
X1 X 2 A3
(X 1 ,X 2 , S) 0; A variable artificial
11
Se ha cambiado S1 por S por existir sólo una variable de holgura en este modelo.
115
Cj 2 3 0 -M
VARIABLES
CB XB X1 X2 S A
BÁSICAS
S 0 4 1 2 1 0 4/2 = 2
A -M 3 1 1 0 1 3/1 = 3
Zj -3M -M -M 0 -M
Z j -C j -M-2 -M-3 0 0
X2 3 2 ½ 1 ½ 0 2/ ½ = 4
A -M 1 ½ 0 -½ 1 1/ ½ = 2
Zj 6-M
3 M
2 2 3
3 M
2 2 -M
Z j -C j - 1 M 0 3 M 0
2 2 2 2
X2 3 1 0 1 1 -1
X1 2 2 1 0 -1 2
Zj 7 2 3 1 1
Z j -C j 0 0 1 1+M
X1 2
X2 1
S 0
A 0 obviamente
Z max. 7
Obsérvese que la variable artificial sale de la base en el segundo tablero. Una vez
salga, no podrá volver a entrar a la base, y puede incluso ignorarse esta columna para
encontrar la solución óptima.12
Debe además recordarse que M siempre es un valor muy grande positivo (M )
y, por lo tanto, su magnitud domina sobre cualquier otro número.
12
Sin embargo, es útil mantener las columnas de las variables artificiales para determinar la inversa de la base
óptima, ya que ellas formaron la matriz idéntica inicial.
116 Introducción a la Modelación Matemática y Optimización
Minimizar Z 2 X 1 3X 2 X 3
Sujeto a : X 1 4 X 2 2 X 3 8
3X1 2 X 2 6
X1, X 2 , X 3 0
Minimizar Z 2 X 1 3 X 2 X 3 MA
Sujeto a : 1 X 1 2 X 2 X 3 S1 4
2
3X1 2 X 2 S2 A 6
X 1 , X 2 , X 3 , S1 , S 2 0; A var iable artificial
Los tableros correspondientes serían los siguientes:
Cj 2 3 1 0 0 M
VAR.
CB XB X1 X2 X3 S1 S2 A
BÁSICAS
X3 1 4 ½ 2 1 -1 0 0 8
A M 6 3 2 0 0 -1 1 2
Zj 4+6M ½ +3M 2+2M 1 -1 -M M
- 3 M 3M -1+2M 0 -1 -M 0
Z j -C j 2
puede entrar a la base sin modificar el valor mínimo de la función objetivo y generar una
solución óptima diferente. Estas soluciones son:
4
X1 2 5
9
X 2 0 5
X 3
3 0 ; 0 1; 0 1 1; ,
S1 0 0
S 0
2 0
Por ejemplo, si = 0.5 y = 0.5, entonces se obtiene la solución óptima (nó básica):
Este hecho hace que sea difícil que un programa de computador detecte la presencia de
soluciones óptimas múltiples, ya que, debido a errores de redondeo, rara vez se tendría un
Z j -C j exactamente igual a cero. Sin embargo, esto se puede lograr estableciendo rangos de
tolerancia para los números muy pequeños considerados iguales a cero.
118 Introducción a la Modelación Matemática y Optimización
En el último tablero, S2 trata de entrar a la base, pero la regla del cociente falla,
entonces como no hay ninguna variable a salir de la base se concluye que hay FUNCIÓN
OBJETIVO NO ACOTADA.13
13
El problema se presenta en la 2ª restricción. La variable X 3 está acotada por debajo pero no por encima y así su
valor puede hacerse tan grande como se desee, haciendo crecer indefinidamente la función objetivo. La variable
de exceso S2 trata de entrar a controlar esto, pero obviamente no puede.
119
Maximizar Z 10 X 1 15 X 2 12 X 3
Sujeto a : 5X 1 3X 2 X 3 144
5 X 1 6 X 2 15 X 3 240
2X1 X 2 X 3 80
X1, X 2 , X 3 0
Forma estándar con variables artificiales:
Maximizar Z 10 X 1 15 X 2 12 X 3 MA
Sujeto a : 5 X 1 3 X 2 X 3 S1 144
5 X 1 6 X 2 15 X 3 S2 240
2X1 X 2 X 3 S 3 A 80
X 1 , X 2 , X 3 , S1 , S 2 , S 3 0; A var iable artificial
En este caso los tableros Simples serían los siguientes:
C j 10 15 12 0 0 0 M
Var.
CB XB X1 X2 X3 S1 S2 S3 A
Básicas
S1 0 144 5 3 1 1 0 0 0 28.8
S2 0 240 -5 6 15 0 1 0 0 NO
A -M 80 2 1 1 0 0 -1 1 40
Zj -80M -2M -M -M 0 0 M -M
Z j Cj -2M-10 -M-15 -M-12 0 0 M 0
X1 10 144/5 1 3/5 1/5 1/5 0 0 0 144
S2 0 384 0 9 16 1 1 0 0 24
A -M 112/5 0 -1/5 3/5 -2/5 0 -1 1 37.3
M 3M 2M
Zj 112M 6 2 2
288 10 5 5 0 M -M
5 5
Zj Cj 0 M
9
3M
10 2M
2
0 M 0
5 5
5
X1 10 24 1 39/80 0 3/16 -1/80 0 0
X3 12 24 0 9/16 1 1/16 1/16 0 0
A -M 8 0 -43/80 0 -7/16 -3/80 -1 1
7 M 21 3M 5
Zj 43M 93
- 8M 528 10 80 8 12 16 8 80 8 M -M
Zj Cj 43M 27
0 7 M 21
3M 5
M 0
80 8 16 8 80 8
120 Introducción a la Modelación Matemática y Optimización
Este método elimina el problema de trabajar con la “Gran M” y los errores de redondeo
asociados. El método comprende las siguientes fases:
FASE I: Trata de encontrar una solución básica factible inicial. Aquí siempre se
minimiza la suma de las variables artificiales, independientemente de si el problema es de
MAXIMIZACIÓN o de MINIMIZACIÓN, sujeto a las restricciones del problema original.
Hay dos posibilidades:
Minimizar Z 4X1 X 2
Sujeto a : 3X1 X 2 3
4 X1 X 2 6
X1 2 X 2 4
( X1, X 2 ) 0
El modelo se transforma como sigue para iniciar el método de las dos fases:
14
Ver Taha (1997), página 92
121
M inimizar Z ' A1 A2
sujeto a : 3 X 1 X 2 A1 3
4 X 1 3X 2 S1 A2 6
X1 2X 2 S2 4
( X 1 , X 2 , S1 , S 2 ) 0; A1 , A2 Var. Artificiales
Cj 0 0 0 0 1 1
Var.
CB XB X1 X2 S1 S2 A1 A2
Básicas
3/3 = 1
A1 1 3 3 1 0 0 1 0
6/4 =
A2 1 6 4 3 -1 0 0 1
1.5
S2 0 4 1 2 0 1 0 0
4/1 = 4
Zj 9 7 4 -1 0 1 1
Z j Cj 7 4 -1 0 0 0
X1 0 1 1 1/3 0 0 1/3 0 3
A2 1 2 0 5/3 -1 0 -4/3 1 1.2
S2 0 3 0 5/3 0 1 -1/3 0 1.8
Zj 2 0 5/3 -1 0 -4/3 1
Z j Cj 0 5/3 -1 0 -7/3 0
Fin
Cj 4 1 0 0 0 0
Fase I
X1 4 3/5 1 0 1/5 0 3/5 -1/5
X2 1 6/5 0 1 -3/5 0 -4/5 3/5 3
S2 0 1 0 0 1 1 1 -1 NO
1
Zj 18/5 4 1 1/5 0
Z j Cj 0 0 1/5 0
X1 4 2/5 1 0 0 -1/5
X2 1 9/5 0 1 0 3/5
S1 0 1 0 0 1 1
Zj Fin
17/5 4 1 0 -1/5
Fase
Z j Cj 0 0 0 -1/5
II
5. TEORÍA DE DUALIDAD
El problema dual puede definirse en forma matemática. Sin embargo, es adecuado
comenzar por una interpretación económica.
Supóngase que en una empresa se producen dos artículos: 1 y 2. Las utilidades netas
unitarias son respectivamente, 35 y 80 u.p/unidad. Los requerimientos unitarios y la
disponibilidad de las máquinas son:
hr.máq.1 hr.máq.2
ARTICULO
unidad unidad
1 2 5
2 3 4
Disponibilidad 300 500
MAXIMIZAR Z 35 X 1 80 X 2
Sujeto a :
2 X 1 3 X 2 300 (Recurso 1) (Problema Primal)
5X 1 4X 2 500 (Recurso 2)
(X1 , X 2 ) 0
Ahora, los recursos, en este caso las máquinas, pueden ser dedicados a otra actividad
diferente a la de producir los artículos 1 y 2. Para ello, sean:
y el producto:
unid .R2 u. p u. p
5W2
unid . prod.1 unid .R2 unid . prod.1
123
y, si se efectúa la suma:
y, análogamente:
3W1 4W2 80 [u.p/unid. prod. 2]
Para hallar la función objetivo del problema que se está construyendo a partir del
problema primal, puede pensarse así: se desea minimizar el costo total de los recursos
involucrados en las otras actividades. Así, el problema dual es:
Por lo tanto, la solución del problema dual da una valiosa información: el precio
unitario adicional máximo que puede pagarse por cada unidad adicional de recurso escaso.
Esto se concluye de la propia definición de la función objetivo del problema dual, ya que,
como se verá posteriormente, el valor óptimo de la función objetivo del problema primal es
igual al valor óptimo de la función objetivo del problema dual ( o sea Zmáx = Z'mín).
Matricialmente estos dos problemas pueden definirse como:
15
Más adelante se verá que cualquier problema de PL tiene su dual asociado. Sin embargo, para hallarlo debe
estar en la forma mostrada, o sea con función objetivo de maximización y restricciones menor ó igual ().
124 Introducción a la Modelación Matemática y Optimización
CX WT b q.e.d
TEOREMA 4.2: Sean X̂ y Ŵ dos soluciones factibles para los problemas primal y
dual respectivamente. Si se cumple que:
CXˆ Wˆ Tb
entonces X̂ y Ŵ son las soluciones óptimas de los problemas primal y dual
respectivamente.
Demostración:
ˆ , W) . Del Teorema 4.1 se concluye que:
a) Tómese la pareja de soluciones factibles ( X
CXˆ WT b
pero como CX ˆ W ˆ T b , por hipótesis, entonces, W ˆ T b WT b , donde W T es
cualquier solución factible del problema dual. Por lo tanto, W ˆ T b será el mínimo y es entonces
la solución óptima del problema dual.
Var.
CB XB X1 X2 S1 S2 S3
Básicas
X2 700 450 0 1 15/14 -15/14 0
X1 400 350 1 0 -3/2 5/2 0
S3 0 100 0 0 3/7 -10/7 1
Zj 455000 400 700 150 250 0
Z j Cj 0 0 150 250 0
Solución óptima del
problema dual asociado
Z j C j 0, j
O, en forma compacta:
C B B 1 A C
1
Si W C B B , entonces W A C
T T
Con lo que se conseguiría cumplir con las restricciones del dual. Obsérvese que las
condiciones de optimalidad del problema primal son equivalentes a las condiciones de
factibilidad del problema dual.
Falta por demostrar que la solución del dual propuesta satisface las condiciones de no
negatividad y el valor óptimo de la función objetivo. En particular para las variables de
holgura, las cuales forman la matriz idéntica, se cumple que:
C B B 1I 0
1
pues sus coeficientes en la función objetivo son iguales a cero. Luego, si W C B B ,
T
se cumple que:
WT 0 ,
o sea que se trata de una solución factible del problema dual. Pero, será la óptima??
Z1 W T b C B B 1b C B X B Z
Luego, de acuerdo con el Teorema 4.2 dicha solución debe ser óptima. (q.e.d)
MIN Z 3X 1 5X 2
Sujeto a : X1 4
2 X 2 12
3 X 12 X 2 18
(X1, X 2 ) 0
Para obtener el problema dual, el modelo debe llevarse a la forma dada en el modelo
generalizado (4.1). La función objetivo se multiplica por (-1) y queda entonces convertida a
una función de Maximización. Las restricciones de basta con multiplicarlas por (-1) y las
restricciones de igualdad se reemplazan por dos restricciones de desigualdad: una de y otro
de . Nótese que aquí no interesa que queden valores negativos en los términos del lado
derecho (Vector b)
Maximizar U 3 X 1 5 X 2
Sujeto a :
X1 4
2 X 2 12
2 X 2 12
3 X 1 2 X 2 18
X1, X 2 0
Obsérvese que las variables duales W2 y W3 aparecen siempre en la forma (W2 W3) y
por lo tanto podrían ser reemplazadas por una sola variable W= W2 W3 no restringida en
signo.
MAX Z 4 X 1 3X 2
Sujeto a : X1 6
X2 8
X1 X2 7
3 X 1 X 2 15
X2 1
(X1, X 2 ) 0
Obsérvese que el problema tiene dos variables de decisión y cinco restricciones. Así,
las bases serían de orden m = 5. Pero, si se plantea el problema dual asociado, tendría cinco
variables duales y dos restricciones, por lo tanto se trabajaría con bases de orden 2, lo cual es
más manejable manual y computacionalmente. Obviamente, esta aplicación adquiere mucho
más sentido en grandes problemas que se presentan en la vida real.
5 1 l
W1 W2 W5 0 W3 W4 Z mín 25
2 2
X 1 4 X 2 3 Zmáx 25
Z j Cj
; y kj 0 y se escoge el valor MINIMO(para problemas de MINIMIZACION)
y kj
Z j Cj
; y kj 0 y se escoge el valor MINIMO(para problemas de MAXIMIZACION)
y kj
donde k es el subíndice asociado a la variable que va a salir de la base.
c) La reducción de Gauss-Jordan es semejante a la del algoritmo SIMPLEX normal.
Ejemplo 5.1: aplicar el algoritmo SIMPLEX dual para resolver el siguiente modelo de
programación lineal:
MIN Z 16 X 1 11X 2 15 X 3
Sujeto a : 2X1 X 2 X 3 3
X 1 2 X 2 3X 3 5
(X1, X 2 , X 3 ) 0
Obsérvese que si se fuera a resolver mediante el algoritmo SIMPLEX normal, se
necesitarían dos variables artificiales para crear la base inicial. El procedimiento mediante el
método SIMPLEX dual es como sigue:
MIN Z 16 X 1 11X 2 15 X 3
Sujeto a : 2 X 1 X 2 X 3 S1 3
X 1 2 X 2 3X 3 - S2 5
( X 1 , X 2 , X 3 , S1 , S 2 ) 0
MIN Z 16 X 1 11X 2 15 X 3
Sujeto a : - 2 X 1 - X 2 - X 3 S1 -3
- X 1 - 2 X 2 - 3X 3 S 2 -5
( X 1 , X 2 , X 3 , S1 , S 2 ) 0
Cj
16 11 15 0 0
Variables
CB XB X1 X2 X3 S1 S2
Básicas
S1 0 -3 -2 -1 -1 1 0
S2 0 -5 -1 -2 -3 0 1
Zj 0 0 0 0 0 0 Se satisfacen
Z j Cj -16 -11 -15 0 0 condiciones de
16 5.5 5
optimalidad
S1 0 -4/3 -5/3 -1/3 0 1 -1/3
X3 15 5/3 1/3 2/3 1 0 -1/3
Zj 25 5 10 15 0 -5
Z j Cj -11 -1 0 0 -5 Se mantiene
6.6 3 15 optimalidad
X2 11 4 5 1 0 -3 1
X3 15 -1 -3 0 1 2 -1
Zj 29 10 11 15 -3 -4
Z j Cj -6 0 0 -3 -4 Se mantiene
2 NO 4
optimalidad
X2 11 7/3 0 1 5/3 1/3 -2/3
Se logro la
X1 16 1/3 1 0 -1/3 -2/3 1/3
factibilidad:
Zj 31 16 11 13 -7 -2 SOLUCIÓN
Z j Cj 0 0 -2 -7 -2 ÓPTIMA
6. ANÁLISIS DE SENSIBILIDAD
Dado un problema de PL de la forma:
MAX Z CX
Sujeto a : AX b
X0
Así, los cambios en el vector b pueden afectar la factibilidad y los cambios en los
coeficientes de la función objetivo y de algunas de las columnas de la matriz A pueden afectar
la optimalidad.
Todos estos cambios presuponen que la matriz base óptima B va a permanecer sin
cambios, ya que si se afecta la base puede dañarse la factibilidad y la optimalidad
simultáneamente, y puede ser preferible volver a resolver el problema desde un comienzo.
133
En este caso los Z j no cambian pues el CB permanece constante; sólo cambian los C j .
Si C j va a cambiar a C j , entonces debe cumplirse que Z j C j 0 (Maximización)
* *
para que la solución siga siendo óptima. Si se rompe el criterio de optimalidad, entonces se
sigue iterando con el algoritmo SIMPLEX normal a partir de ese punto.
En este caso, dado que cambia el vector C B , cambian los C j . Por lo tanto, éstos deben
recalcularse para investigar las condiciones de optimalidad.
criterio de optimalidad se rompe, debe seguirse iterando con el algoritmo SIMPLEX normal.
(interpretación económica del problema dual), es válido en el rango para el cual la base
óptima sigue siendo la misma.
b * b b
B 1 b * B 1 b B 1 b
Y j B 1a *
*
j
Y se recalcula:
Z j CB Yj
* *
dicha condición, la solución actual sigue siendo óptima; si no se cumple, entonces se continúa
iterando con el algoritmo SIMPLEX normal, hasta obtener las condiciones de optimalidad.
Entonces se calcularía:
Yn1 B 1a
n 1
Z n1 C B Yn 1
Obsérvese que este caso puede tratarse como si la columna original an 1 hubiera sido
de ceros y estuviera cambiando a los valores dados.
MAX Z 3 X 1 5 X 2
Sujeto a : X 1 4
2 X 2 24
3 X 1 2 X 2 18
(X1, X 2 ) 0
MAX Z 3 X 1 5 X 2
Sujeto a : X 1 S1 4
2X 2 S2 24
3X 1 2 X 2 S 3 18
( X 1 , X 2 , S1 , S 2 , S3 ) 0
Cj 3 5 0 0 0
Variables
CB XB X1 X2 S1 S2 S3
Básicas
S1 0 4 1 0 1 0 0
X2 5 9 3/2 1 0 0 ½
S2 0 6 -3 0 0 1 -1
Zj 45 15/2 5 0 0 5/2
Z j Cj 9/2 0 0 0 5/2
Cj 3 5 0 0 0 0
Variables
CB XB X1 X2 S1 S2 S3 S4 Entra como
Básicas variable básica
S1 0 4 1 0 1 0 0 0 Obsérvese que
X2 5 9 3/2 1 0 0 ½ 0 X2 es básica,
S2 0 6 -3 0 0 1 -1 0 pero su
columna es [0 1
S4 0 24 2 3 0 0 0 1 0 3], luego,
S1 0 4 1 0 1 0 0 0 debe corregirse
X2 5 9 3/2 1 0 0 ½ 0 mediante
Gauss-Jordan
S2 0 6 -3 0 0 1 -1 0
S4 0 -3 -5/2 0 0 0 -3/2 1
Zj 45 15/2 5 0 0 5/2 0
Z j Cj 9/2 0 0 0 5/2 0
1.8 1.66
Aquí se toma el
S1 0 4 1 0 1 0 0 0 en valor absoluto
X2 5 8 2/3 1 0 0 0 1/3 (Maximización)
S2 0 8 -4/3 0 0 1 0 -2/3 NUEVA
S3 0 2 5/3 0 0 0 1 -2/3 SOLUCIÓN
ÓPTIMA
Zj 40 10/3 5 0 0 0 5/3
Z j Cj 1/3 0 0 0 0 5/3
Obsérvese que el nuevo valor óptimo de la función objetivo es Z*máx = 40 < 45 (Valor
anterior), como es de esperarse al añadir una nueva restricción que vuelva al problema no-
factible.
1. Incrementar la capacidad del taller T3, introduciendo una nueva máquina que
añade 6000 hr/año de capacidad con un costo de $1‟020.000 anuales.
2. Una inversión en publicidad que vale $700.000 anuales y que podría cambiar las
demandas en 20, 14 y 10 respectivamente (o sea que las nuevas demandas serían
30, 22 y 20, respectivamente).
3. Una investigación que cuesta $1‟200.000 al año y que permitiría disminuir el
tiempo requerido en el taller T3 para el producto A3, de 1000 a 500 hr/unidad.
4. Un esfuerzo técnico por reducir los costos, aumentando el beneficio unitario del
producto A3 en $60.000/unidad. El costo de amortización durante la vida
remanente de fabricación del producto sería de $400.000/año.
Nota: las alternativas NO son mutuamente excluyentes, o sea que podría estudiarse
combinación de ellas.
El lector debe comprobar que la solución óptima única viene expresada en el siguiente
tablero final óptimo:
Cj 3 4 4.8 0 0 0 0 0 0
Var.
CB XB X1 X2 X3 S1 S2 S3 S4 S5 S6
Básicas
S1 0 7000 0 0 0 1 0 -0.4 140 -500 0
Solución S2 0 5000 0 0 0 0 1 -0.8 -20 -1000 0
óptima
degenerada X3 4.8 0 0 0 1 0 0 0.001 -0.6 0 0
(X3=0 X1 3 10 1 0 0 0 0 0 1 0 0
y básica) X2 4 8 0 1 0 0 0 0 0 1 0
S6 0 10 0 0 0 0 0 -0.001 0.6 0 1
Zj 62 3 4 4.8 0 0 .0048 0.12 4 0
Z j Cj 0 0 0 0 0 .0048 0.12 4 0
17
Se ha dividido la función objetivo Z entre 50.000 ( Z Z / 50000 ) para trabajar más fácilmente.
140 Introducción a la Modelación Matemática y Optimización
w1=0
w2=0
w3=(0.0048)(50000)=240 18
w4=(0.12)(50000)=6.000
w5=(4)(50000)=200.000
w6=0
Es importante notar que algunos solucionadotes como WinQSB pueden dar una
solución diferente a la anterior, ya que en el penúltimo tablero simplex se presenta un empate
en la regla de salida.
18
Obsérvese que la solución del problema dual debe multiplicarse también por 50.000 , al igual que el valor real
óptimo de la función objetivo.
19
Tómese B 1 del tablero óptimo de la sección anterior y compruébese este resultado.
141
ˆ
X B [4600 200 6 10 8 4]T
Obsérvese que la nueva solución sigue siendo factible y, por lo tanto, es también
óptima. Se mantienen los niveles de producción de A1 y A2 y ahora si debe producirse A3
(X3=6 unidades/año de A3).
Por lo tanto, esta alternativa por sí sola es adecuada (el valor óptimo anterior era de
$3‟100.000/año).
Nota: Del análisis dual se concluyó que por cada hora adicional de T3 la función
objetivo se incrementaría en $240/hr·T3. En este caso esto se cumple ya que la base óptima no
cambia, y:
$240
600hr··T 3 1'440.000 $ / Año
hr··T 3
3'100.000 $ / Año
4'540.000 $ / Año
que era lo que se había obtenido.
Este análisis del problema dual será válido si la composición básica de la solución
óptima NO cambia, o sea si la solución permanece factible. Además, obsérvese que se está
pagando:
1'020.000 $ / Año
$170 / hr $240 / hr
6000 hr / Año
que es el máximo adicional que estaba dispuesto a pagarse por cada hora adicional en
T3, lo cual confirma aún más lo dicho anteriormente.
Una pregunta que podría surgir es: cuánto es lo máximo que se pueden aumentar las
hr/Año disponibles en T3 para que se siga manteniendo la factibilidad y por lo tanto la función
objetivo crezca a razón de $240/hr·T3?
7000 0.4 K
5000 0.8 K
ˆ 0.001K
X B X B X B
10
8
10 0.001K
Pero, como para que se conserve la factibilidad se requiere que X B 0 , entonces:
En este caso:
Luego se rompe la factibilidad y habría que seguir iterando con el algoritmo SIMPLEX
DUAL; el lector debe comprobar que la nueva solución óptima es:
X 1 10 X 2 13 X 3 0 S1 4500 S 2 S 3 0 S 4 20 S 5 9 S 6 20
*
Z máx (82)(50.000) 4'100.000 $ / Año
Costo Publicidad 700.000 $ / Año
$3'400.000/Año
143
Obsérvese que aquí NO se cumple lo previsto por la solución del dual, ya que cambia
base óptima al romperse la factibilidad. De acuerdo al dual la función objetivo debía
incrementarse en:
Cj 3 4 4.8 0 0 0 0 0 0
Var.
CB XB X1 X2 X3 S1 S2 S3 S4 S5 S6
Básicas
S1 0 7000 0 0 700/3 1 0 -1/6 0 -500 0
S2 0 5000 0 0 -100/3 0 1 -5/6 0 -1000 0
X1 3 10 1 0 5/3 0 0 1/600 0 0 0
S4 0 0 0 0 -5/3 0 0 -1/600 1 0 0
X2 4 8 0 1 0 0 0 0 0 1 0
S6 0 10 0 0 1 0 0 0 0 0 1
Zj 62 3 4 5 0 0 1/200 0 4 0
Z j Cj 0 0 0.2 0 0 1/200 0 4 0
Y luego si se itera con la nueva columna. El lector debe comprobar que la nueva
solución óptima es:
X 1 1 23 X 2 8 X 3 10 S1 3833 13 S 2 1166 23 S 3 0 S 4 8 13 S5 S6 0
*
Z máx (85)(50.000) 4'250.000 $ / Año
Costo alternativa 1'200.000 $ / Año
$3'050.000/Año
X 1 0 X 2 8 X 3 6 S1 5600 S 2 5200 S 3 0 S 4 10 S 5 0 S 6 4
*
Z máx (68)(50.000) 3'400.000 $ / Año
Costo alternativa 400.000 $ / Año
$3'000.000/Año
Luego, por sí sola esta alternativa no sería adecuada. Trate el lector de verificar
combinaciones de alternativas(primero parejas, luego ternas y luego las cuatro) por
computador y llegar a la conclusión de que la mejor alternativa combinada que puede hacerse
es:
X 1 3 13 X 2 8 X 3 10
*
Z máx 5'100.000 $ / Año
Costo alternativa 1'420.000 $ / Año
$3'680.000/Año
La forma estándar:
M AX Z CX
sujeto a :
AX b
X0
pasa a:
Z CX 0
AX b
X0
1 C Z 0
0 A X b
El problema es determinar una solución básica tal que Z sea máximo; por lo tanto Z
siempre será una variable básica. Así, una base cualquiera para este caso sería:
1 C 1 C 2 ··· ··· Ck
0
0
: Se le denomina B1
B1 pues se aumentó su
: B
dimensión en 1.
1 CB :
B1
0
0 B 0
147
1 C B B 1 - C j - C j C B B 1 a j
1
Y1j B1 a (1)
j
0 B 1 a j B 1 a j
1
B1 b (1)
(1)
XB
Pero como los coeficientes de las variables de holgura en la función objetivo son
iguales a cero, entonces:
1 0
1
B1
0 I m
Luego:
0
Ib b
(1) (1) (1)
XB
b
En cualquier iteración:
1 1 C B B 1 0 C B B 1b Z
XB 1 1
0 B b B b X B
Luego, se procede a detectar la variable que sale, recordando que Z siempre es básica.
Para lograr esto, se calcula la nueva columna Yk1 , si k es la variable a entrar. Recuérdese que
este cálculo se realiza así:
Z CK
YK1 B11 AK1 K
YK
149
Y así, se puede aplicar la regla del cociente, saliendo de la base aquella variable cuyo
sea mínimo, donde:
X X
Br B i ; Yik 0
Yrk Yik
En este caso se hará por Gauss-Jordan. Si el lector está interesado en conocer el otro
método puede consultar a Taha (1991)20.
20
TAHA, Hamdy. Investigación de operaciones. Alfa y Omega. 2ª edición, México, 1991. Páginas 273-282.
150 Introducción a la Modelación Matemática y Optimización
MAX Z
Sujeto a :
Z - 400 X 1 - 700 X 2 0
7
X1 X 2 S1 1400
3
X 1 1.4 X 2 S 2 980
X1 X 2 S 3 900
( X 1 , X 2 , S1 , S 2 , S3 ) 0
B1 b1 b 2 b3 b4
Inicialmente, y dado que las variables básicas iniciales son Z, S1, S2, S3, coincide con la
matriz idéntica de orden 4. Por notación:
B1 1 2 3 4
1
1 siempre será el vector unitario ya que Z siempre será variable básica. Luego, el
tablero inicial simplificado es:
Inversa de la base
(1)
Base XB Y21
1 2 3 4
Z 1 0 0 0 0 -700
S1 0 1 0 0 1400 7/3 600
S2 0 0 1 0 980 1.4 700
S3 0 0 0 1 900 1 900
400
1
Var X 1 : Z 1 C1 1 0 0 0 400
1
1
700
7/3
Var X 2 : Z 2 C 2 1 0 0 0 700
1.4
1
151
Luego la variable que entra es X2. Por lo tanto, debe calcularse Y21 , para poder
determinar la variable que sale (por esa razón aparece en la tabla mostrada arriba).
Dado que Y2(1) debe convertirse en el vector unitario, entonces se realizan las
operaciones fila elementales necesarias. El resultado aparece en el segundo tablero:
Inversa de la base
(1)
Base XB Y21
1 2 3 4
Z 1 300 0 0 420000 -100
X2 0 3/7 0 0 600 3/7 1400
S2 0 -3/5 1 0 140 2/5 350
S3 0 -3/7 0 1 300 4/7 525
0
1
var . S1 : Z 2 C 2 1 300 0 0 300
0
0
Entonces:
1 300 0 0 100
400
0 3 7 0 0 1 37
Y11 B11 A11 3
0 1 0 1 2
5 5
0 3 0 1 1 4
7 7
0
1
var . S1 : Z 1 C1 1 150 250 0 150
0
0
0
0
var . S 2 : Z 2 C 2 1 150 250 0 250
1
0
Aunque este método no parece tener ventajas en la forma manual mostrada, sus
principales beneficios radican en los grandes ahorros de cálculo computacional para
problemas de gran tamaño que se presentan en la práctica.
153
Ejercicios 7.1
1. Encuentre por inspección una solución óptima del siguiente modelo de PL:
Maximizar z 3 x1 4 x2 6 x3
sujeto a :
2 x1 3 x2 4 x3 16
( x1 , x2 , x3 ) 0
2. Cuáles de las siguientes afirmaciones sobre el método Simplex son Falsas ó Verdaderas
(Justifique su respuesta):
c. El criterio de salida del método Simplex garantiza que la solución básica siguiente sea
factible.
d. En un tablero óptimo, todos los (zj – cj) de las variables básicas deben ser iguales a
cero.
e. La solución óptima del problema dual asociado a un problema primal aparece siempre
en la fila de los (zj – cj) debajo de la submatriz idéntica inicial.
cj 7 3 5 0 0
Var. CB XB x1 x2 x3 S1 S2
Bás.
x1 100 - 2 1/2
7/2
x3 250 1 14/ 1/6
5
zj
z j cj
c. Si usted quisiera encontrar otra solución básica factible óptima diferente a la dada en el
tablero, ¿qué variable entraría a la base y cuál saldría de ella?
d. ¿Si la inversa de la base inicial se formó con las columnas de las variables S1 y S2, cuál
es la solución óptima del problema dual asociado?
a. Cuando en cualquier tablero Simplex existe una variable candidata a entrar a la base,
pero ninguna variable puede salir de ella porque el criterio de salida falla, ¿qué se
puede concluir acerca del problema original?
b. Si en un tablero Simplex óptimo existe una variable básica artificial con valor positivo,
¿qué se puede decir del problema original? ¿Qué puede decirse si la variable básica
artificial es igual a cero? (Solución degenerada).
5. Considere los siguientes resultados de un programa que resuelve modelos de PL, por
ejemplo el WinQSB. El problema original es un problema de producción semejante al de
los transformadores, con 4 productos candidatos a ser producidos (representados por las
variables x1, x2, x3 y x4), y 3 recursos utilizados para su producción (hr.hombre, materia
prima y hr.máq):
Restricción Valor Dir. Valor del lado Holgura ó Precio Mínimo Máximo
Lado izq. Derecho Exceso Sombra RHS RHS
(RHS)
Hr.hombre 2150 2500 350 0 2150 M
Mat_Prima 800 800 0 1.5 700 1050
Hr.máq. 1250 1250 0 2.0 1000 1375
Responda las siguientes preguntas en forma breve y precisa y justifique sus respuestas:
e. Usted ha estimado que el costo normal por cada hora en la máquina ($/hr.máq.) es de
$18/hr.máq. Dado que se trata de un recurso escaso, tal como lo muestra la solución
arriba, usted está considerando la posibilidad de alquilar ciertas horas máquina a un
costo de $21/hr.máq. ¿Debería usted aceptar esta oferta de alquiler? Recuerde que
usted está alquilando para usted para tener más recursos.
Var. CB XB x1 x2 x3 S1 S2 S3
Bás.
x3 0 0.4 0.2
x2 0.5 0.1 0.3
S3 0.5 0.7 0.1
zj
zj c j
a. Complete totalmente el tablero anterior y compruebe que se trata del tablero óptimo del
problema original.
b. Suponga que le ofrecen 100 unidades adicionales de cada uno de los tres recursos a un
costo adicional de $1.0/unidad. ¿Cuál de los tres recursos debería comprarse en
primera instancia? ¿Qué información faltaría para tomar una decisión definitiva acerca
de la compra de estos recursos?
c. Suponga que se aumenta la utilidad unitaria del producto 1 (x1) de 2 a $3/unidad. ¿Cree
usted que entonces ahora sería rentable producirlo? ¿Por qué?
156 Introducción a la Modelación Matemática y Optimización
e. Suponga que existe una restricción adicional por transporte, en la cual las cantidades
sumadas de los tres productos a producir no puede ser mayor que 280 unidades. ¿Cuál
sería entonces la nueva solución óptima?
sujeto a:
X1 + X2 + X3 25
2X1 - X2 + 3X3 10
(X1, X2, X3) 0
sujeto a:
0.4X1 + 0.35X2 + 0.5X3 + 0.7X4 = 0.55
0.6X1 + 0.65X2 + 0.5X3 + 0.3X4 = 0.45
(X1, X2, X3, X4) 0
sujeto a:
3X1 - 3X2 + 2X3 + 8X4 50
5X1 + 6X2 - 4X3 - 4X4 40
-4X1 + 2X2 - X3 - 3X4 -20
(X1, X2, X3, X4) 0
Minimizar Z = X1 - X2 + 3X3
sujeto a:
X1 - 2X2 + X3 5
2X1 - 3X3 1
-4X2 + 5X3 8
(X1, X2, X3) 0
sujeto a:
X1 + 2X2 + X3 430
3X1 + 2X3 460
X1 + 4X2 420
(X1, X2, X3) 0
Nota: Trate independientemente cada una de las alternativas anteriores [(b) a (f)].
158 Introducción a la Modelación Matemática y Optimización
Ejercicio 7.2
1. Usted posee 50 acres de tierra en la que puede plantar cualquier cantidad de maíz, soya y
lechuga. La siguiente tabla muestra la información relevante con respecto de la
producción, la ganancia neta y los requerimientos de agua de cada cultivo.
Dado que se tiene disponible 100.000 litros de agua y que la compañía se ha comprometido
a vender al menos 5.120 Kg de maíz, se desarrolló el siguiente modelo de PL, donde Xi =
Acres a sembrar de cada tipo de cultivo i = 1, 2, 3:
a) Con la misma cantidad de dinero, qué sería mejor, ¿adquirir 10 acres de terreno adicionales
ó comprar 13.000 litros adicionales de agua? (Sólo puede escoger una de las dos
alternativas).
b) Manteniendo la misma cantidad de tierra y condición de demanda del maíz, usted desea
obtener una ganancia de $21.000. ¿Cuántos litros adicionales de agua necesitaría para
lograrlo?
i) Reconstruya el tablero final Simplex. Observe que debió haberse adicionado una
variable artificial A en la restricción de demanda de maíz (además de la variable de
exceso S3), y, por lo tanto, la tercera columna de la matriz inversa anterior
corresponde a la columna de dicha variable artificial y las dos primeras a las
variables de holgura S1 y S2 adicionadas en las dos primeras restricciones de tierra y
agua.
iii) Compruebe el rango de sensibilidad dado por el programa del recurso tierra en la
función objetivo.
f) Formule el problema dual del modelo primal formulado en este punto y encuentre su
solución óptima sin resolverlo desde el comienzo.
sujeto a:
a) Sería de $6.625
161
b) Sería de $9.125
c) Sería de $5.125
d) Permanecería igual a $5.875
e) Sólo podría determinarse resolviendo el problema desde el tablero inicial Simplex
8. OPTIMIZACIÓN NO LINEAL
Los problemas de optimización no-lineal son en general difíciles de resolver, es decir,
son problemas NP-Hard. Más aún, no existen algoritmos generales para resolver cualquier
problema de programación no-lineal. Dependiendo de la estructura del problema en particular,
se han desarrollado algunos algoritmos eficientes, como por ejemplo Programación
Cuadrática, Programación Convexa, Programación Geométrica, etc.
Uno de los problemas de la optimización no-lineal surge del hecho de que la región
factible resultante no necesariamente es convexa. Además, el óptimo no necesariamente ocurre
en uno de los vértices de la región factible.
Sujeto a:
10
9
8
7
6
5
4
3
2
1
1 2 3 4 5 6 7 8 9 10
163
Ejercicios 8.1
Minimizar
R/. Mínimo
Se van a tratar funciones de varias variables en general. El caso para una variable es un
caso especial, para el cual las ecuaciones y criterios siguen siendo válidos.
Condiciones necesarias. Una condición necesaria para que sea un punto extremo ó
estacionario (máximo, mínimo ó punto de silla) de es que:
ó, equivalentemente
Cómo los signos del menor principal. De orden k es , se puede concluir que IH
( ) es negativa definida y por lo tanto es un punto máximo de .
Ejercicio 8.2
a)
b)
c)
d)
e)
Donde:
=(
g=
Donde
IP =
|Δ|=0
Deben ser
a) Negativas si es un punto máximo.
b) Positivas si es un punto mínimo.
167
Minimizar
Sujeto a:
= (2, 2, 1, 1)
= (2, -2, 1, 1)
= (2.8, 0, 1.4, 1.4)
En este caso m=1 y n=3. Se evaluarán los últimos (n-m) = 2 menores principales de
, iniciando con el menor principal de orden 2m+1=3. Para cada punto estacionario, esto
daría lo siguiente:
a) = (2, 2, 1, 1)
b) Para = (2, -2, 1, 1) se obtienen los mismos resultados que para el caso
(a), y por lo tanto el punto es un mínimo (compruébelo)
c) Para = (2.8, 0, 1.4, 1.4), es:
d)
Cuyas raíces son µ=2 y µ=8/9, ambas mayores que cero, reforzando el hecho de que
ambos son puntos mínimos.
Cuyas raíces son µ=2 y µ=-0.8, con lo que se concluye que este punto no es un punto
extremo.
Ejercicios 8.3
1) Resuelva:
Minimizar
Sujeto a:
Minimizar =
Sujeto a:
Donde C > 0.
3) Resuelva:
Minimizar =
Sujeto a:
Considérese el problema:
S= y =
Dadas las restricciones g( ) 0, una condición necesaria para optimizar es que sea
no negativa (no positiva) para problemas de maximización (minimización). Ademas, tomando
las derivadas parciales de L( , se obtiene:
Por todo lo anterior, las condiciones necesarias de Kuhn Tucker para que y sea
un punto estacionario del problema de maximización son:
( )
Maximizar (ó Minimizar) Z =
Sujeto a
Cóncava
Lineal
Minimización Convexa Convexa
Cóncava
Lineal
Minimizar =
Sujeto a:
=
=
=
=
=
(
173
Ejercicio 8.4
a) Maximizar
Sujeto a:
b) Minimizar
Sujeto a:
174 Introducción a la Modelación Matemática y Optimización
Nótese que la ventaja de este método es que la función a optimizar no requiere ser
unimodal. Su desventaja radica en el probable número muy grande de números aleatorios que
deben generarse y en que al final del proceso no se conoce que tan cerca se puede estar del
óptimo real, o sea que no se conoce el intervalo de incertidumbre final.
Tomando los intervalos anterior y posterior al punto donde se obtuvo el mejor valor. Si
el mejor valor de la función objetivo se obtiene en un extremo, entonces sólo se toma el sub-
175
intervalo anterior (ó posterior) al extremo dado. Se repite así el proceso con el nuevo intervalo
de incertidumbre Δ, hasta donde sea necesario.
x1 x2
XL = a XR = b
a+b
2
Dado que la función f(x) es unimodal sobre el intervalo [a, b], se pueden presentar tres
casos:
En el intervalo [0,5]
Utilice Δ = 0.001
ITER
No
0 0 5 2.4995 2.5005R 21.35938 21.39063
1 0 2.5005 1.24975L 1.25075 5.35961 5.358676
R
2 1.24975 2.5005 1.874625 1.875625 8.3335272 8.345171
R
3 1.24975 1.875625 1.5621875 1.5631875 5.86457007 5.8690465
4 1.24975 1.5631875 1.40596875 1.40696875R 5.40090202 5.40245366
5 1.24975 1.40696875 1.327859375 1.328859375R 5.3317479 5.3320018
6 1.288804688 1.327859375 1.288804688L 1.289804688 5.33408798 5.3337343
7 1.308332032 1.327859375 1.308332032L 1.309332032 5.3299532 5.32989992
8 1.308332032 1.327859375 1.318095704 1.319095704R 5.33010099 5.33020048
9 1.308332032 1.319095704 1.313213868 1.314213868R 5.32984076 5.32986364
10 1.308332032 1.314213868 1.31077295L 1.31177295 5.329850525 5.329835273
11 1.31077295 1.314213868 1.311993409 1.312993409R 5.32983401 5.329837814
12 1.31077295 1.311993409 1.31138318L 1.31238318 5.329839632 5.329833635
13 1.31138318 1.312993409
Obsérvese que en la iteración No. 13 se puede escoger el punto medio del intervalo de
incertidumbre [1.31138318, 1.312993409], el cual es x* 1.312188.
Dado que Δ = 0.001, no se puede esperar una precisión mayor que Δ. El valor más
exacto de este mínimo en [0, 5] es x*= 1.31224996
entonces es el valor de r que maximiza a h(r). como h(r) es una función de una sola
variable, las técnicas de búsqueda vistas anteriormente son muy útiles aquí, teniendo en
cuenta que h(r) debe ser unimodal.
Ejemplo:
Max.
Se tiene que:
Primera iteración:
h(r) = (1-2r, 1) = -2
Segunda iteración:
Sea = + r(0 , 1) =
178 Introducción a la Modelación Matemática y Optimización
h(r) = -2 , obteniéndose y
El lector debe comprobar que las siguientes iteraciones producen los resultados
mostrados en la siguiente tabla:
Iteración (K)
0 -1,1 (-2,0) -
1 (0, 1)
Óptimo
3
2
1
2
1 1 3 2
2 2
179
Ejercicio 9.1
Max , donde:
A=
Considere .
180 Introducción a la Modelación Matemática y Optimización
Max(min) Z = CX + DX
Sujeto a:
AX
AX
Donde:
X=
C= (
b=
A es una matriz m x n =
D es una matriz cuadrada n x n
Maximizar Z = CX + DX
G(X)= X- 0
Sean
=
U= (
Nótese que:
Sean:
Nótese que excepto para las condiciones y =0, las demás ecuaciones son
lineales en X, , u y S. O sea que el problema se reduce a resolver un conjunto de ecuaciones
lineales satisfaciendo las condiciones adicionales y =0 .
Como las condiciones de Kuhn Tucker son también suficientes en este caso, la
solución factible que satisface todas estas condiciones tiene que ser la solución optima única.
182 Introducción a la Modelación Matemática y Optimización
Para resolver este sistema se utiliza la fase I del método de las dos fases, con la
condición adicional =0 y . Esto se controla evitando que y sean
positivas (básicas) simultáneamente, y que sean positivas (básicas) simultáneamente.
Ejemplo:
Maximizar Z =
Sujeto a:
Max Z = +
Sujeto a:
(1 , 2) 2
(
Para resolver estas ecuaciones se agregan dos variables artificiales y de tal forma
que la submatriz idéntica de orden 3 aparezca en las condiciones de Kuhn Tucker:
183
Min r = +
0 0 0 0 0 1 1 0
Var. Bás. Ѳ
← 1 4 4 2 1 -1 0 1 0 0 4÷4 =1
1 6 2 4 2 0 -1 0 1 0 6÷2= 3
0 2 1 2 0 0 0 0 0 1 2÷1 =2
14 10 6 6 3 -1 -1 1 1 0
6 6 3 -1 -1 0 0 0
Nótese que puede entrar a la base pues = 0 (nó-básica). Luego el segundo tablero
es:
0 0 0 0 0 1 1 0
Ѳ
Var. Bás.
0 1 1 0 0 0 1÷ =2
Como se encontró una solución factible sin variables artificiales en la base, la solución
óptima del problema original es:
*= , *= *= 4
= = +
CTU(y)= +
+ = 0 y* =
= =
Para 1, 2, … , n, sean:
= Demanda
= Costo de preparación
= Costo de manejo / unidad por unidad de tiempo
= Cantidad de pedido
= Requerimiento de área (o volumen) de almacenamiento por cada unidad en
inventario
A = Area (ó volumen) máximo de almacenamiento disponible para todos los n
artículos.
Si se supone que no hay faltantes (“backorders”), el modelo matemático sería el
siguiente:
Sujeto a:
*= ; i= 1,2,…, n
186 Introducción a la Modelación Matemática y Optimización
L( =
Recuérdese de las condiciones de Kuhn Tucker que en este caso para minimización,
<0.
*=
*= = = 11.55 Unidades
*= = = 20.00 Unidades
*= = = 24.49 Unidades
Ejercicio 10.1
2) Serie de problemas 20.3c (Pág 770-771): No. 2. Resuelva además el siguiente problema:
Sujeto a:
6) Serie de problemas 11.3c (Pág 451-452): No. 3. Nota: Utilice los resultados del problema
No. 4 para aproximar a .
189
Proceso de Esmaltado:
Las ventas anuales de aluminio terminado son de 1,000 ton en la planta D y de 1,200 ton en la
planta E.
190 Introducción a la Modelación Matemática y Optimización
DE HACIA B C D E
j = 1 2 3 4
A 1 400 2,010 510 1,920
i = B 2 10 630 220 1,510
C 3 1,630 10 620 940
DE HACIA D E
k = 1 2
B 1 220 1,510
j = C 2 620 940
D 3 0 1,615
E 4 1,465 0
Wj
Xij B
Yjk
A
A T
C D l e
u r
B m m
i i
n n
D E i a
o d
C o
E
Variables de Decisión
Función Objetivo
Mina A : 420( X AB X AC X AD X AE )
Mina B : 360( X BB X BC X BD X BE )
Mina C : 540( X CB X CC X CD X CE )
Costo anual de transporte de bauxita desde las minas de bauxita hacia las plantas de alúmina
($/año):
Costo anual de transporte de alúmina desde las plantas de alúmina hacia las plantas de
esmaltado ($/año):
Planta B : 3,000,000WB
Planta C : 2,500,000WC
Planta D : 4,800,000WD
Planta E : 6,000,000WE
Restricciones
Mina A : X AB X AC X AD X AE 36,000
Mina B : X BB X BC X BD X BE 52,000
Mina C : X CB X CC X CD X CE 28,000
Estas restricciones expresan que todo el flujo anual de bauxita que sale de cada mina no puede
exceder su capacidad anual de explotación. Note que aquí no hay decisiones de cierre de
minas involucradas, es decir, que se asume que las minas siempre van a estar abiertas.
193
Planta B : X AB X BB X CB 40,000WB
Planta C : X AC X BC X CC 20,000WC
Planta D : X AD X BD X CD 30,000WD
Planta E : X AE X BE X CE 80,000WE
Estas restricciones aseguran que no se puede recibir flujo alguno de bauxita desde las
minas, si la planta de alúmina destino no ha sido abierta. Las restricciones han sido
formuladas de acuerdo con la capacidad de cada planta de alúmina de procesar la bauxita
entrante. Otra forma en la que estas restricciones podrían escribirse es de acuerdo con la
capacidad de cada planta de alúmina de producir alúmina saliente. Sin embargo, la
información dada en el problema no permite formularlas de esta forma.
Este ejemplo es útil para ilustrar que las restricciones de capacidad en una planta
productora pueden escribirse con base en las materias primas entrantes o también con base en
los productos terminados salientes. En el caso de los ingenios azucareros, para ilustrar con
otro ejemplo, podríamos escribir las restricciones de capacidad de procesamiento de la caña de
azúcar que llega al ingenio o podríamos definir las restricciones con base en la capacidad de
producción de azúcar. Obsérvese finalmente que las anteriores restricciones, en combinación
con el conjunto de restricciones (5) descrito más adelante, aseguran que no se puede despachar
cantidad alguna de alúmina hacia las plantas de esmaltado, si la planta origen no ha sido
abierta, completando así la lógica del modelo.
Estas restricciones expresan que todo el flujo anual de alúmina que llega a cada planta de
esmaltado no puede exceder su capacidad anual de procesamiento de alúmina. De nuevo, son
restricciones de capacidad dadas por la materia prima entrante a la planta.
4) Por ventas anuales de aluminio terminado en cada planta de esmaltado (Ton de aluminio
terminado/año):
Estas restricciones aseguran que la demanda anual proyectada en cada planta de esmaltado se
va a satisfacer en forma total. Nótese que se aplica el factor de rendimiento de alúmina en
aluminio terminado (0.4 toneladas de aluminio terminado por cada tonelada de alúmina). Esta
relación sencilla permite en este modelo no formular variables adicionales para el aluminio
194 Introducción a la Modelación Matemática y Optimización
terminado. En extensiones de este modelo, tal como la planteada en el Problema No. 1 de los
Ejercicios 12.1, es probable que se requieran formular variables de decisión adicionales.
Estas restricciones aseguran que la cantidad anual de bauxita que entra a cada planta de
alúmina, afectada por su correspondiente rendimiento para producir alúmina, es igual a la
cantidad anual de alúmina que sale de dicha planta hacia las plantas de esmaltado. Al
formular estas restricciones se asume, por una parte, la linealidad del proceso de
transformación de bauxita en alúmina y, por la otra, que no hay diferenciación de calidad de la
alúmina saliente. En la práctica, es posible que las cosas no sean así de sencillas y que cada
materia prima tenga calidad diferente y produzca diferentes calidades de un mismo producto
terminado.
6) Por límites en los valores de cada una de las variables (restricciones obvias):
X ij 0 i, j
Y jk 0 j , k
W j 0, 1 j
Sujeto a:
X AB X AC X AD X AE 36,000
X BB X BC X BD X BE 52,000
X CB X CC X CD X CE 28,000
195
X AB X BB X CB 40,000WB 0
X AC X BC X CC 20,000WC 0
X AD X BD X CD 30,000WD 0
X AE X BE X CE 80,000WE 0
X ij 0 i, j
Y jk 0 j , k
W j 0, 1 j
Este modelo contiene 24 variables de decisión, de las cuales 4 son binarias y el resto son
continuas, y 15 restricciones (sin incluir las restricciones obvias). Los modelos matemáticos
de cadenas de abastecimiento internacionales que se formulan en la práctica pueden contener
decenas o cientos de miles de variables con decenas o cientos de miles de restricciones,
llegando incluso a millones de ellas. La dificultad para su solución eficiente depende
principalmente del número de variables enteras (binarias) involucradas en el modelo y de su
estructura.
La solución de este modelo puede encontrarse con la ayuda de software comercial especial
para ello, mediante el solver de Excel™ o incluso mediante un software más especializado
como AMPL, a pesar de tratarse de un modelo pequeño. La hoja electrónica que resuelve este
modelo se muestra en la Figura 2. En general, cualquier modelo de programación lineal
entera-mixta relativamente pequeño (por ejemplo, con un máximo de 240 celdas variables) se
puede resolver en formatos parecidos a los que se presentan en la figura mencionada.
196 Introducción a la Modelación Matemática y Optimización
La celda objetivo corresponde a la celda del costo total $C$35. Las celdas cambiantes son
todas las variables de decisión, incluyendo todos los flujos (celdas B10 a la B29, es decir, las
variables Xij y Yjk) y las variables binarias de apertura de las plantas de alúmina (celdas B30 a
la B33, es decir, las variables Wj).
Las restricciones corresponden primero a aquéllas que son de desigualdad de
correspondientes a las capacidades de las minas de bauxita, las plantas de alúmina y las
plantas de esmaltado (columnas D a la L). Obsérvese que estas restricciones se pueden
escribir en forma compacta mediante la expresión $D$37:$L$37 $D$38:$L$38, lo que
198 Introducción a la Modelación Matemática y Optimización
significa que está tomando el valor de cada una de estas nueve restricciones y lo está haciendo
menor al recurso disponible en cada caso. El segundo grupo de restricciones hay que
escribirlo aparte porque son restricciones de igualdad y corresponden a las dos de demanda de
aluminio terminado, la cual se va a satisfacer exactamente, y a las cuatro de balance que por su
naturaleza siempre son igualdades (éstas corresponden a la expresión $M$37:$R$37
$M$38:$R$38 de las restricciones del solver). Al solver se le debe decir cuáles de las
variables o celdas cambiantes son binarias. Para este efecto se tiene la expresión Apertura =
binario en solver. Nótese que aquí no aparecen las referencias de las celdas, ya que se ha
utilizado, para ilustrar, la definición de un „nombre‟ en Excel y se han bautizado las celdas
B30 a B33 correspondientes a las variables binarias de apertura de plantas de alúmina con el
nombre de “Apertura”. Esto no es necesario hacerlo, pero muchas veces es conveniente para
efectos de identificación y fácil recordación de las variables.
El resto de variables o celdas cambiantes, por defecto, las tomará reales mayores ó iguales
que cero (aunque debe tenerse cuidado de que esto en realidad esté definido así en la ventana
que se abre al pulsar las „Opciones…‟ de la ventana del solver). Finalmente, cuando se le da el
comando „Resolver‟ al solver, se obtiene entonces la solución óptima mostrada en la Figura 2
anterior.
B
40,000
El valor de la función objetivo indica que el costo total de logística mínimo es de 87,455,600
$/año. Este costo incluye los costos de explotación y producción, los costos de transporte y
los costos fijos. Esto se logra abriendo las plantas de alúmina ubicadas en B, C y D (variables
binarias WB = WC = WD = 1) y cerrando la planta ubicada en E (variable binaria WE = 0), la
199
cual al parecer es una planta muy ineficiente, tiene un costo fijo muy alto y/o tiene problemas
de localización. Los flujos óptimos son los que se muestran en las celdas B10 a B29 de la
Figura 2, ilustrados en forma gráfica en la Figura 4. Se sugiere al lector comprobar que los
balances de flujo tanto en las plantas de alúmina como en las plantas de esmaltado se cumplen.
En cuanto a los recursos sobrantes es importante notar que la mina A apenas está siendo
explotada en un 4.63% de su capacidad (1,667 de 36,000 ton/año disponibles). Muy
probablemente análisis adicionales pueden revelar que debería considerarse el cierre de esta
mina o la redistribución de su materia prima. Además, su bauxita es la de menor calidad de
acuerdo con su rendimiento de alúmina. En cuanto a las plantas de alúmina abiertas, la B y la
C trabajan al 100% de su capacidad y la D tiene una capacidad sobrante del 54.44% (sobran
16,333 ton/año de capacidad de las 30,000 disponibles). Las plantas de esmaltado también
están trabajando holgadamente: La planta D trabaja al 62.5 % de su capacidad y la E al
42.86%.
Se sugiere al lector comprobar que los siguientes análisis de sensibilidad son correctos:
# CONJUNTOS PRINCIPALES
# PARÁMETROS
# VARIABLES DE DECISIÓN
# FUNCIÓN OBJETIVO
# RESTRICCIONES
# subject to pl_abierta_b:
# w["B"] = 1;
203
# subject to pl_abierta_c:
# w["C"] = 1;
# subject to pl_abierta_d:
# w["D"] = 1;
# subject to pl_abierta_e:
# w["E"] = 1;
# CONJUNTOS PRINCIPALES
set MINAS:= A B C;
set PLALU:= B C D E;
set PLESM:= D E;
# PARÁMETROS
param ctran_al: D E :=
B 220 1510
C 620 940
D 0 1615
E 1465 0;
param ctran_b: B C D E :=
option show_stats 1;
option solution_precision 0;
option omit_zero_rows 1;
option omit_zero_cols 1;
option display_precision 6;
option display_round 1;
option display_width 50;
# COMANDO DE SOLUCIÓN:
solve;
printf "\n\n*************************************\n";
printf "RESULTADOS DEL PROBLEMA DE LA BAUXITA\n";
printf "*************************************\n\n";
Los resultados que se obtienen con el anterior archivo de comandos son los siguientes, los
cuales coinciden totalmente con lo que ya habíamos encontrado anteriormente:
*************************************
RESULTADOS DEL PROBLEMA DE LA BAUXITA
*************************************
w [*] :=
B 1.0
C 1.0
D 1.0
x :=
A D 1666.7
B B 40000.0
B D 12000.0
C C 20000.0
207
y :=
B D 1440.0
B E 1760.0
C E 1240.0
D D 1060.0
cap_exp.slack [*] :=
A 34333.3
C 8000.0
cap_exp [*] :=
B -660.0
cap_prodal.slack [*] :=
D 16333.3
cap_prodal [*] :=
B -196.4
C -436.4
cap_proc_alu.slack [*] :=
D 1500.0
E 4000.0
5. Si se obliga a la planta de alúmina E a estar abierta: Se abren todas las plantas, pero no se
produce nada en la planta de alúmina E. Simplemente se paga su costo fijo y el costo total
anual asciende a 93,455,600 $/año.
7. Si sólo se abren las plantas B, C y E, se obtiene una solución óptima con un costo mínimo
de 103,457,000 $/año.
Una red consta de un conjunto de nodos unidos por arcos. La notación empleada es
(N,A), donde:
N= Conjunto de nodos
A= Conjunto de Arcos
1 3
2 4
Normalmente, un arco tiene una capacidad asociada, la cual limita el posible flujo a
través de el. Por ejemplo, se puede hablar de flujo de petróleo a través de un oleoducto, de
flujo de productos a través de una cadena de abastecimiento, o de flujo de vehículos a lo largo
de una red de careteras.
A consecuencia de lo anterior, una red puede ser dirigida u orientada si permite flujo
positivo en un sentido de sus arcos y flujo cero en el sentido contrario. Una red puede también
ser cíclica cuando presenta al menos un ciclo, o sea cuando existe una ruta desde un nodo
hacia sí mismo.
1 2 4
3 5
Una red es conectada cuando cada dos nodos distintos están unidos por lo menos por
una ruta.
Un árbol es una red conectada acíclica, incluyendo un subconjunto de los nodos de la red. Un
árbol de expansión es un árbol que incluye todos los nodos de la red. Por ejemplo, para la red
mostrada en la página anterior, se muestra un árbol y un árbol de expansión:
1 3 1 3
2 2 4
Árbol de expansión mínima
En este algoritmo, se escoge un nodo cualquiera para iniciar y se definen dos conjuntos
de nodos en cada iteración, los que ya han sido conectados y los que aún no. Cada iteración
determina el arco de longitud minima para unir uno de los nodos que ya han sido conectados
con los que no lo han sido. El algoritmo termina cuando todos los nodos han sido conectados,
resultando así el árbol de expansión mínima. Formalmente, el algoritmo es como sigue:
21
“Minimun Spanning Tree”
211
Sean:
Paso 0: Hacer = , N
Paso 1: Escoger cualquier nodo ϵ N. Hacer = { } y, por lo tanto, N–{ .
Hacer K=2.
Paso K: Seleccionar un nodo j* ϵ que produce el arco más corto hacia un nodo en
el conjunto conectado . Unir permanentemente j* con y eliminarlo de , o sea:
,
Si ⟹ Parar. De lo contrario, hacer k=k+1 y repetir el paso k.
Ejemplo: Se muestran 6 poblaciones, las cuales deben ser unidas por una red de
carreteras con la longitud mínima de carreteras pavimentadas. Los arcos mostrados
representan todas las conexiones posibles con sus respectivas longitudes en km. Utilizando el
algoritmo del árbol de expansión mínima, determine cómo deben unirse las poblaciones.
50
2 4
13
15
0
48
25
17
32
1 6
50 52
12
75
0
100
3 5
Por conveniencia, se escoge el nodo 1 como el punto de partida del algoritmo. Por
conveniencia, los nodos que ya han sido conectados se muestran con doble círculo.
2 4
15
25
32
1 6
50
12
0
3 5
Iteración 1
212 Introducción a la Modelación Matemática y Optimización
50
2 4
15
25
48
32
1 6
50
12
0
3 5
Iteración 2
Las líneas delgadas muestran los arcos candidatos en cada iteración para unir el
conjunto de nodos que ya han sido conectados con el suplemento.
Las líneas punteadas representan el arco de longitud mínima y las líneas gruesas
representan los arcos permanentes que unen nodos ya conectados.
Nótese que en la iteración 1, por ejemplo, j*=2, y en la iteración 2, j*=4 y que daría
unido al nodo 1 , mediante el arco (1 , 4). En la página siguiente se muestra el resto de las
iteraciones.
50
2 4
13
0
15
25
48
32
1 6
50
17
12
0
3 5
Iteración 3
50
2 4
13
0
15
25
48
32
1 6
50
17
52
100
3 5
Iteración 4
213
2 4
15
25
48
32
1 6
50
17
75
100
3 5
Iteración 5
2 4
15
25
48
32
1 6
17
3 5
Iteración 6
(Árbol de expansión mínima)
O sea que la longitud mínima de carreteras pavimentadas que uniría a las poblaciones es igual
a 15+25+31+17+48 km = 137 km. Aunque visualmente el árbol de la iteración 6 pareciese
cíclico (lo cual sería una contradicción), la figura siguiente aclara este aspecto:
1
32
15
25
2 6
4
48
17
5 Árbol de expansión
mínima
3
214 Introducción a la Modelación Matemática y Optimización
Este problema consiste en determinar la ruta mas corta entre un punto de origen y un
punto de destino de la red (ver algunas aplicaciones en Taha, pág. 222-225).
Dos algoritmos para resolver este problema en redes cíclicas o acíclicas son:
El algoritmo de Dijkstra.
Algoritmo de Floyd.
El primero determina la ruta más corta entre el nodo origen y cada uno de los otros
nodos de la red. L algoritmo de Floyd permite determinar la ruta más corta entre dos nodos
cualesquiera de la red.
ALGORITMO DE DIJKSTRA:
Sea:
[ ]=[ ], .
Una clasificación temporal se puede reemplazar con otra cuando se encuentra una ruta
más corta. Cuando es evidente que no se puede encontrar una ruta más corta, la clasificación
temporal se cambia a permanente. El algoritmo completo es como sigue:
Paso 0: Clasificar el nodo de origen (nodo 1) como permanente [0,-]. Haga i=1.
Paso i:
rompiendo con los empates en forma arbitraria. Haga i=r, clasifique el nodo r como
permanente y repita el paso i.
Ejemplo: Determine la ruta más corta desde el nodo 1 hacia todos los nodos de la red
2 12
7
8
5 9
(Nodo 8
Origen) 1 3 9
5 7
7
6
13
4 6
siguiente:
La forma más fácil de resolver los problemas de la ruta mas corta es hacerlo
directamente en la red, copiando las clasificaciones sucesivas. Cuando una clasificación se
vuelve permanente, por notación, se encierra en un rectángulo. La red de la página siguiente
vá mostrando las clasificaciones sucesivas.
[7 , 1]
2 12 [12 , 4]
7
8
5 9
8
1 3
9
[0 , -] 7 [21 , 5]
5
[8 , 1]
7 6
13 6
4
[17 , 3]
[5 , 1]
Las rutas más cortas desde el nodo de origen (nodo 1) y cualquier nodo se encuentran
partiendo del nodo destino y devolviéndose con la información dada por las clasificaciones
permanentes. La tabla siguiente resume los resultados:
216 Introducción a la Modelación Matemática y Optimización
13. BIBLIOGRAFÍA
AHUJA, Ravindra K., Thomas L. Magnanti y James B. Orlin, Network Flows: Theory,
Algorithms, and Applications, Prentice Hall, Englewood Cliffs, New Jersey, 1993. (Excelente
texto para profundización en Teoría de Redes)
DUQUE M., Ramón, Modelos Lineales: Formulación, Universidad del Valle, Cali,
1987. (Colección de problemas formulados y resueltos para consulta)
VIDAL, Carlos J. y M. Goetschalckx, "A Global Supply Chain Model with Transfer
Pricing and Transportation Cost Allocation," European Journal of Operational Research, Vol.
129, pp. 134 – 158, 2001. (Artículo de referencia)