UIA2
UIA2
UIA2
MODELOS INDUSTRIALES
UNIDAD
2
PROGRAMACIÓN LINEAL Y DUAL
Se prohíbe la reproducción total o parcial de esta obra, por cualquier medio,
electrónico o mecánico sin la autorización de la Dirección de Educación a
Distancia (UDB Virtual).
En la programación lineal, el adjetivo lineal significa que todas las funciones matemáticas
del modelo deben ser funciones lineales. La palabra programación no se refiere a términos
computacionales; en esencia es sinónimo de planeación (Hillier & Gerald J. Liberman,
2010).
1
La programación lineal es un caso especial de la programación matemática, en donde todas
las funciones que contempla el modelo son líneas: siempre tendremos una función objetivo
lineal a optimizar (maximizar o minimizar), sujeta a restricciones lineales individuales. Las
variables del modelo, que son continuas, únicamente pueden introducirse valores no
negativos. Si bien puede parecer que estos supuestos quitan realismo al problema porque
el analista este limitado al uso de ecuaciones que quizás no son frecuentes en el mundo
real, las técnicas de programación lineal se utilizan en un amplio espectro de campos como
se describió anteriormente.
2
5.1.1¿En qué consiste la programación lineal?
Según Salazar et al.(2014) la programación lineal consiste en encontrar los valores de unas
variables que maximizan o minimizan un único objetivo sujeto a una serie de restricciones.
Las principales características de la PL (Programación Lineal) son:
Con este video demostrativo lograras ampliar las características de los modelos a
resolver con la programación lineal.
Goal Project. (8 de agosto de 2017). Programación Lineal Introducción 01. [Video].
https://www.youtube.com/watch?v=fxVHxXuhBLY
A continuación, se analizarán con más detalles las características y lo que ocurre si una o
varias de ellas no se cumplen en los modelos que involucran la PL.
En primer lugar, cabe destacar que en la PL todas las funciones utilizadas tanto en el
objetivo como en las restricciones son lineales. Es decir, las restricciones consisten en la
suma de variables multiplicadas por sus respectivos parámetros, siendo esta función
menor, igual o mayor que un determinado recurso. El objetivo también es lineal, si bien
desconocemos a priori su valor. En caso de que tanto su valor como una o más restricciones
no sean lineales, sería necesario introducir métodos de programación no lineal, son muchos
más complejos de resolver y cuya optimalizad no siempre está garantizada.
3
En segundo lugar, la PL considera que las variables de decisión son continuas. Desde el
punto de vista matemático de obtención de soluciones, esta característica no ofrece
problemas. En muchas situaciones la interpretación económica de la solución de un
problema de PL no tiene sentido si obtenemos fracciones en las variables. Por ejemplo, si
estamos asignando trabajadores a tareas, no tiene sentido un resultado que en un momento
determinado asigne 3.4 trabajadores a una determinada tarea. Sin embargo, se opta por
redondear al entero más próximo se puede cometer un grave error.
En tercer lugar, los modelos de programación lineal consideran que hay un único objetivo
a maximizar o minimizar. Muchas veces podemos tener que resolver problemas que tienen
más de un objetivo. Por ejemplo, en un enfoque queremos maximizar la cobertura de un
determinado servicio, mientras que por otro enfoque queremos reducir los costos generales.
Ambos objetivos son conflictivos, en el sentido de que aumentar la cobertura significaría un
aumento en la necesidad de recursos con el consecuente incremento de costos del sistema.
Esta conflictividad se resuelve utilizando métodos de programación multicriterio o
multiobjetiva.
Sin embargo, existen situaciones en las que uno o más parámetros tienen un componente
estocástico, o en palabras menos técnicas, tienen una variabilidad (en algunos casos puede
ser representada por una distribución estadística). Si esto acontece, la PL ya no es un buen
instrumento para la obtención de soluciones.
4
5.1.2 Formulación de modelos
A continuación, se presentan algunos ejemplos de los problemas con los cuales se puede
encontrar una organización y como la programación lineal puede expresarlos
matemáticamente.
Una clínica de salud ha decidido ampliar su servicio de urgencias (abierto las 24 horas) con
la necesidad de nuevo personal de enfermería. La gerencia de la clínica ha estimado las
necesidades mínimas de personal por tramos de horarios para cubrir las urgencias que se
presenten. Se definieron 6 tramos de 4 horas. La necesidad mínima de personal en cada
tramo se idéntica en la Tabla 1. El departamento de recursos humanos ha informado a
gerencia que los contratos laborales han de ser de ocho horas seguidas. El problema se
centra en encontrar el número mínimo de personal necesario para cubrir la demanda.
Tabla 1.
Necesidades de personal por tramos de horarios.
Tramos de horarios
𝒊 1 2 3 4 5 6
0:00-4:00 4:00-8:00 8:00-12:00 12:00-16:00 16:00-20:00 20:00-24:00
Personal 𝑵𝒊 8 4 2 6 4 5
Fuente: Nota. Elaboración propia
Inicialmente se tienen que definir las variables del modelo que queremos desarrollar. Como
debemos establecer el número de personal en cada turno, definimos 𝑁𝑖 como la cantidad
de personal por turno 𝑖, en donde 𝑖 = 1, … ,6. Es decir, existe una variable correspondiente
a cada turno.
5
Las restricciones del modelo tienen que reflejar la necesidad de que la cantidad de personal
en el periodo 𝑖 más el número de personas que entraron a trabajar en el turno 𝑖 − 1 sea
suficiente para cubrir las necesidades del turno 𝑖(𝑁𝑖 ). Esta situación queda reflejada en la
Tabla 2. En esta tabla, un trabajador que entra a trabajar, por ejemplo, a las 4:00, trabajará
en los turnos 2 y 3, por lo tanto, contribuirá a cubrir las necesidades de estos dos turnos.
En otras palabras, el turno 𝑖 estará siendo atendido por 𝑋𝑖−1 y 𝑋𝑖 . Por lo tanto, tendremos
que 𝑋𝑖−1 + 𝑋𝑖 (el personal que trabaja durante el turno 𝑖) tiene que ser, como mínimo, igual
a 𝑁𝑖 , que es el número mínimo de personal de enfermería necesario para este turno. En
términos matemáticos la restricción es la siguiente:
𝑋𝑖−1 + 𝑋𝑖 > 𝑁𝑖
En términos generales existirá una restricción por cada horario, es importante que se
considera la analogía como horario de entrada.
𝑀𝑖𝑛 𝑍 = 𝑋1 + 𝑋2 + 𝑋3 + 𝑋4 + 𝑋5 +𝑋6
C.S.R
𝑋6 +𝑋1 > 8
𝑋1 +𝑋2 > 4
𝑋2 +𝑋3 > 2
𝑋3 +𝑋4 > 6
𝑋4 +𝑋5 > 4
𝑋5 +𝑋6 > 5
𝑋1 , 𝑋2 , 𝑋3 , 𝑋4 , 𝑋5 , 𝑋6 > 0
6
Tabla 2
Necesidades de personal
Tramos de horarios
X 1 2 3 4 5 6
0:00-4:00 4:00-8:00 8:00-12:00 12:00-16:00 16:00-20:00 20:00-24:00
Personal Xi 8 4 2 6 4 5
00:00 𝑋1 𝑋1
04:00 𝑋2 𝑋2
08:00 𝑋3 𝑋3
12:00 𝑋4 𝑋4
16:00 𝑋5 𝑋5
20:00 𝑋6 𝑋6
Personal 𝑵𝒊 8 4 2 6 4 5
En la tabla 2 se observa una matriz de distribución del personal en cada turno, las variables
contempladas en cada columna de tramos de horarios y la disponibilidad de personal
(Personal 𝑁𝑖 ) forman las restricciones anteriormente planteadas.
7
El equipo médico ha estimado el tiempo medio que necesita cada paciente de tipo A y de
tipo B en cada uno de los servicios pre‐quirúrgico (PQ), quirúrgico (QI) y postoperatorio
(PO). La experiencia en un hospital similar muestra que por cada tres pacientes de tipo A
que llegan al hospital como mínimo llega uno de tipo B. Por otra parte, se ha estimado el
costo de cada paciente en los diferentes servicios. En la tabla 3 se muestra los datos del
problema, se considera que la capacidad ociosa es en horas mensuales y el costo por
paciente en dólares.
Tabla 3
Estimación de horas de las cirugías A y B
Horas necesarias de cirugía Capacidad
Salas
A B Ociosa
Sala PQ 1 3 144
Sala QI 3 2 162
Sala PO 4 2 -
Costo $20 $25
Fuente: Nota. Elaboración propia
Como el servicio postoperatorio (PO) aún no existe, el gerente argumenta que para justificar
su creación tiene que utilizarse durante un mínimo de 135 horas al mes. Por otra parte, el
presupuesto mensual asignado a las nuevas cirugías es de $1,500. El gerente quiere saber
el número máximo de pacientes que podrán ser operados al mes.
Primero definimos las variables del modelo. Sean 𝑋1 y 𝑋2 el número total de pacientes por
mes que pueden ser tratados con la cirugía A y B respectivamente. A continuación, se
presentan las restricciones.
8
Se ha establecido que en la sala PQ se disponen de 144 horas. En otras palabras, la
utilización de esta sala no puede sobrepasar las 144 horas. Como cada uno de los pacientes
de tipo A y de tipo B consumen 1 hora y 3 horas en esta sala respectivamente, el número
total de horas mensuales consumidas en PQ para los dos tipos será igual a 𝑋1 + 3𝑋2. Este
número tiene que ser inferior o igual a las 144 horas. La restricción será la siguiente:
𝑋1 + 3𝑋2 ≤ 144
El mismo razonamiento puede ser utilizado para determinar el número límite de horas en la
sala QI. Como el total de horas consumidas será igual a 3𝑋1 + 2𝑋2 , y hay un máximo de 162
horas disponibles, la restricción sobre QI será:
3𝑋1 + 2𝑋2 ≤ 162
El gerente ha determinado que, para viabilizar los nuevos tratamientos, se tiene que ocupar
la nueva sala PO durante un mínimo de 135 horas al mes. Como el número de horas
mensuales que se utilizará en PO es igual a 4𝑋1 + 2𝑋2, tendremos que:
La experiencia en otros hospitales muestra que, por cada 3 pacientes de tipo A, viene
como mínimo un paciente de tipo B. Matemáticamente, esto se expresa como:
X1
≤ 𝑋2
3
Es equivalente a:
𝑋1 − 3𝑋2 ≤ 0
9
Finalmente, el gasto mensual realizado en las dos cirugías no puede exceder $1,500. Como
cada paciente de tipo A y de tipo B cuesta $20 y $25 respectivamente, el gasto total mensual
será de 20𝑋1 + 25𝑋2 , cantidad que no puede exceder $1,500 tendremos que:
Por último, se necesita formular el objetivo. El gerente quiere saber el número máximo de
enfermos de tipo A y de tipo B que se puede atender cada mes. Simplemente, tendremos
que, si Z es este número, el objetivo se expresará como:
𝑀𝑎𝑥 𝑍 = 𝑋1 + 𝑋2
𝑀𝑎𝑥 𝑍 = 𝑋1 + 𝑋2
C.S.R
𝑋1 + 3𝑋2 ≤ 144
3𝑋1 + 2𝑋2 ≤ 162
4𝑋1 + 2𝑋2 ≥ 135
𝑋1 − 3𝑋2 ≤ 0
20𝑋1 + 25𝑋2 ≤ 1,500
𝑋1 , 𝑋2 > 0
10
Problema de transporte
Se desea enviar productos a dos clientes ubicados en Honduras, los destinos son San
Pedro Sula y Tegucigalpa, desde tres bodegas diferentes ubicadas en El Salvador, San
Miguel, San Salvador y Santa Ana. Los costos de transporte unitarios se muestran en la
tabla 4, así como las cantidades de producto con los que cuenta cada almacén y las
unidades que necesita cada cliente. Determinar el modelo de transporte que representa
cada situación.
Tabla 4
Costos de transporte
Bodegas San Pedro Sula Tegucigalpa Oferta
San Miguel $0.42 $0.51 20,000
San Salvador $0.45 $0.48 10,000
Santa Ana $0.47 $0.45 25,000
Demanda 30,000 30,000
Fuente: Nota. Elaboración propia
En este caso, el supuesto no se cumple, debido a que Demandas > Ofertas; por tanto, se
requiere un punto de oferta ficticio con cinco mil unidades, que en el contexto del problema
será producto que los almacenes no podrán cumplir. En la tabla 5 se observa el punto de
oferta ficticio de 5,000 unidades para cumplir el aspecto de 𝑂𝑓𝑒𝑟𝑡𝑎𝑠 = 𝐷𝑒𝑚𝑎𝑛𝑑𝑎𝑠.
11
Tabla 5
Costos de transporte con oferta artificial
Bodegas San Pedro Sula Tegucigalpa Oferta
San Miguel $0.42 $0.51 20,000
San Salvador $0.45 $0.48 10,000
Santa Ana $0.47 $0.45 25,000
Artificial $0.00 $0.00 5,000
Demanda 30,000 30,000
Fuente: Nota. Elaboración propia
𝑥11 = 𝑐𝑎𝑛𝑡𝑖𝑑𝑎𝑑 𝑑𝑒 𝑝𝑟𝑜𝑑𝑢𝑐𝑡𝑜𝑠 𝑎 𝑒𝑛𝑣𝑖𝑎𝑟 𝑑𝑒 𝑏𝑜𝑑𝑒𝑔𝑎 𝑒𝑛 𝑆𝑎𝑛 𝑀𝑖𝑔𝑢𝑒𝑙 𝑎 𝑆𝑎𝑛 𝑃𝑒𝑑𝑟𝑜 𝑆𝑢𝑙𝑎
𝑥12 = 𝑐𝑎𝑛𝑡𝑖𝑑𝑎𝑑 𝑑𝑒 𝑝𝑟𝑜𝑑𝑢𝑐𝑡𝑜𝑠 𝑎 𝑒𝑛𝑣𝑖𝑎𝑟 𝑑𝑒 𝑏𝑜𝑑𝑒𝑔𝑎 𝑒𝑛 𝑆𝑎𝑛 𝑀𝑖𝑔𝑢𝑒𝑙 𝑎 𝑇𝑒𝑔𝑢𝑐𝑖𝑔𝑎𝑙𝑝𝑎
𝑥21 = 𝑐𝑎𝑛𝑡𝑖𝑑𝑎𝑑 𝑑𝑒 𝑝𝑟𝑜𝑑𝑢𝑐𝑡𝑜𝑠 𝑎 𝑒𝑛𝑣𝑖𝑎𝑟 𝑑𝑒 𝑏𝑜𝑑𝑒𝑔𝑎 𝑒𝑛 𝑆𝑎𝑛 𝑆𝑎𝑙𝑣𝑎𝑑𝑜𝑟 𝑎 𝑆𝑎𝑛 𝑃𝑒𝑑𝑟𝑜 𝑆𝑢𝑙𝑎
𝑥22 = 𝑐𝑎𝑛𝑡𝑖𝑑𝑎𝑑 𝑑𝑒 𝑝𝑟𝑜𝑑𝑢𝑐𝑡𝑜𝑠 𝑎 𝑒𝑛𝑣𝑖𝑎𝑟 𝑑𝑒 𝑏𝑜𝑑𝑒𝑔𝑎 𝑒𝑛 𝑆𝑎𝑛 𝑆𝑎𝑙𝑣𝑎𝑑𝑜𝑟 𝑎 𝑇𝑒𝑔𝑢𝑐𝑖𝑔𝑎𝑙𝑝𝑎
𝑥31 = 𝑐𝑎𝑛𝑡𝑖𝑑𝑎𝑑 𝑑𝑒 𝑝𝑟𝑜𝑑𝑢𝑐𝑡𝑜𝑠 𝑎 𝑒𝑛𝑣𝑖𝑎𝑟 𝑑𝑒 𝑏𝑜𝑑𝑒𝑔𝑎 𝑒𝑛 𝑆𝑎𝑛𝑡𝑎 𝐴𝑛𝑎 𝑎 𝑆𝑎𝑛 𝑃𝑒𝑑𝑟𝑜 𝑆𝑢𝑙𝑎
𝑥32 = 𝑐𝑎𝑛𝑡𝑖𝑑𝑎𝑑 𝑑𝑒 𝑝𝑟𝑜𝑑𝑢𝑐𝑡𝑜𝑠 𝑎 𝑒𝑛𝑣𝑖𝑎𝑟 𝑑𝑒 𝑏𝑜𝑑𝑒𝑔𝑎 𝑒𝑛 𝑆𝑎𝑛𝑡𝑎 𝐴𝑛𝑎 𝑎 𝑇𝑒𝑔𝑢𝑐𝑖𝑔𝑎𝑙𝑝𝑎
𝑥41 = 𝑐𝑎𝑛𝑡𝑖𝑑𝑎𝑑 𝑑𝑒 𝑝𝑟𝑜𝑑𝑢𝑐𝑡𝑜𝑠 𝑎 𝑒𝑛𝑣𝑖𝑎𝑟 𝑑𝑒 𝑏𝑜𝑑𝑒𝑔𝑎 𝑓𝑖𝑐𝑡𝑖𝑐𝑖𝑜 𝑎 𝑆𝑎𝑛 𝑃𝑒𝑑𝑟𝑜 𝑆𝑢𝑙𝑎
𝑥42 = 𝑐𝑎𝑛𝑡𝑖𝑑𝑎𝑑 𝑑𝑒 𝑝𝑟𝑜𝑑𝑢𝑐𝑡𝑜𝑠 𝑎 𝑒𝑛𝑣𝑖𝑎𝑟 𝑑𝑒 𝑏𝑜𝑑𝑒𝑔𝑎 𝑓𝑖𝑐𝑡𝑖𝑐𝑖𝑜 𝑎 𝑡𝑒𝑔𝑢𝑐𝑖𝑔𝑎𝑙𝑝𝑎
En este problema se busca minimizar el costo de transporte, dado a que se conoce el costo
de transportar una unidad de cada bodega a cada cliente, la función objetico se determina
de la siguiente forma:
12
Las restricciones del problema se plantean de la siguiente forma:
Modelo matemático:
𝑀𝑖𝑛𝑖𝑚𝑖𝑧𝑎𝑟 𝑧 = 0.42𝑥11 + 0.51𝑥12 + 0.45𝑥21 + 0.48𝑥22 + 0.47𝑥31 + 0.45𝑥32
C.S.R
𝐵𝑜𝑑𝑒𝑔𝑎 𝑑𝑒 𝑆𝑎𝑛 𝑀𝑖𝑔𝑢𝑒𝑙: 𝑥11 + 𝑥12 = 20,000
𝐵𝑜𝑑𝑒𝑔𝑎 𝑑𝑒 𝑆𝑎𝑛 𝑆𝑎𝑙𝑣𝑎𝑑𝑜𝑟: 𝑥21 + 𝑥22 = 10,000
𝐵𝑜𝑑𝑒𝑔𝑎 𝑑𝑒 𝑆𝑎𝑛𝑡𝑎 𝐴𝑛𝑎: 𝑥31 + 𝑥32 = 25,000
𝐵𝑜𝑑𝑒𝑔𝑎 𝑓𝑖𝑐𝑡𝑖𝑐𝑖𝑎: 𝑥41 + 𝑥42 = 5,000
𝐶𝑙𝑖𝑒𝑛𝑡𝑒 𝑒𝑛 𝑆𝑎𝑛 𝑃𝑒𝑑𝑟𝑜 𝑆𝑢𝑙𝑎: 𝑥11 + 𝑥21 + 𝑥31 + 𝑥41 = 30,000
𝐶𝑙𝑖𝑒𝑛𝑡𝑒 𝑒𝑛 𝑡𝑒𝑔𝑢𝑐𝑖𝑔𝑎𝑙𝑝𝑎: 𝑥12 + 𝑥22 + 𝑥32 + 𝑥42 = 30,000
𝑥𝑖𝑗 ≥ 0 𝑝𝑎𝑟𝑎 𝑖 = 1, … ,4; 𝑗 = 1, … ,2
13
En este video observaras la formulación de un modelo de programación lineal del tipo
maximización.
Universidad Continental – Modalidad a Distancia. (19 de septiembre de 2017).
Programación lineal. [Video]. https://www.youtube.com/watch?v=AB20K1IQLVo
14
Clase 6| Método gráfico
6.Método gráfico
6.1 Método Gráfico
El método grafico se considera el primer método de solución usado para resolver modelos
de programación lineal. Este método se limita a que conste de dos variables de decisión y
un numero finito de restricciones lineales (Salazar et al., 2014).
Para la solución grafica de problemas de programación lineal con dos variables, lo que se
debe hacer es trazar un eje de coordenadas cartesianas, para graficar las desigualdades
dadas por el problema, después encontrar el área de soluciones factibles y proceder a
graficar la función objetivo para conocer el valor optimo (maximizar o minimizar) que será
la solución del problema.
15
6.1.1 Problema de maximización
Quimes, Inc. Es una empresa pequeña que fabrica una variedad de productos químicos.
En un proceso de producción particular se utilizan tres materias primas para elaborar dos
productos: un aditivo para combustible y una base para solventes.
El aditivo se vende a las compañías petroleras y se utiliza en la producción de gasolina y
otros combustibles. La base para solventes se vende a una variedad de compañías de
productos químicos y se usa en artículos de limpieza para el hogar y la industria. Las tres
materias primas se mezclan para formar el aditivo para el combustible y la base para
solventes como en la tabla 6, en la que se muestra que una tonelada de aditivo para
combustible es una mezcla de 0.4 ton de material 1 y 0.6 ton de material 3, mientras que
una tonelada de base para solventes es una mezcla de 0.5 ton de material 1, 0.2 ton de
material 2 y 0.3 ton de material 3.
Tabla 6
Requerimientos de material por tonelada
Producto
Aditivo para combustible Base para solventes
Material 1 0.4 0.5
Material 2 0.2
Material 3 0.6 0.3
Se utilizan 0.6 ton de material 3 en cada tonelada de aditivo para
combustible
Fuente: Nota. Elaboración propia
La producción de Quimes está restringida por una disponibilidad limitada de las tres
materias primas. Para el periodo de producción actual cuenta con las siguientes cantidades
de cada materia prima:
16
Tabla 7
Cantidad disponible de materia prima
Cantidad disponible
Material
para la producción
Material 1 20 ton
Material 2 5 ton
Material 3 21 ton
Fuente: Nota. Elaboración propia
El departamento de contabilidad analizó las cifras de producción, asignó todos los costos
relevantes y se llegó a precios para ambos productos que generarían una contribución a
las utilidades de $40 por cada tonelada de aditivo para combustible producido y $30 por
cada tonelada producida de base para solvente.
17
Para la definición de la función objetivo se debe considerar que la contribución a las
utilidades de la organización proviene de la producción de F toneladas de aditivo para
combustible y S toneladas de base para solvente. Como la organización gana $40 por cada
tonelada de aditivo para combustible producida y $30 para cada tonelada de base para
solvente producida, Quimes ganará $40F de la producción de aditivo para combustible y
$30S de la producción de base para solvente. Por tanto:
Cada tonelada de aditivo para combustible que se produce utiliza 0.4 ton de material 1. Por
tanto, se utilizan 0.4F ton de material 1 para producir F ton de aditivo para combustible.
Cada tonelada de base para solvente que Quimes produce utiliza 0.5 ton de material 1, así
que se emplean 0.5S ton de material 1 para producir S ton de base para solvente. Por
consiguiente, el número de toneladas de material 1 utilizado para producir F ton de aditivo
para combustible y S ton de base para solvente es:
0.4𝐹 + 0.5𝑆 ≤ 20
Restricción 2:
El aditivo para combustible no utiliza material 2, pero cada tonelada de base para solvente
que Quimes produce utiliza 0.2 ton de material 2, así que se utilizan 0.2S ton de material 2
para producir S ton de base para solvente. Por consiguiente, el número de toneladas de
material 2 empleadas para producir F ton de aditivo para combustible y S ton de base para
solvente es:
0.2𝑆 ≤ 5
18
Restricción 3:
Cada tonelada de aditivo para combustible que Quimes produce utiliza 0.6 ton de material
3. Por lo anterior, se utilizan 0.6F ton de material 1 para producir F ton de aditivo para
combustible.
De manera similar, cada tonelada de base para solventes que la empresa produce utiliza
0.3 ton de material 3, así que se emplean 0.3S ton de material 1 para producir S ton de
base para solvente. Por consiguiente, el número de toneladas del material 3 empleadas
para producir F ton de aditivo para combustible y S ton de base para solventes es:
0.6𝐹 + 0.3𝑆 ≤ 21
𝐹≥0 𝑦 𝑆≥0
Al escribirse forma abreviada:
𝐹, 𝑆 ≥ 0
Modelo matemático:
La grafica de la figura 1 tiene los valores de F en el eje horizontal y los de S en el eje vertical.
Cualquier punto de la gráfica puede identificarse por medio de los valores de F y S, los
cuales indican la posición del punto a lo largo de los ejes horizontal y vertical,
respectivamente. Esto nos dice que cada punto en la gráfica corresponde a una solución
posible. La solución 𝐹 = 0 𝑦 𝑆 = 0 se conoce como el origen. Debido a que tanto F como S
deben ser no negativos, la gráfica de la figura 1 solo muestra las soluciones en que 𝐹 ≥
0 𝑦 𝑆 ≥ 0.
0.4𝐹 + 0.5𝑆 ≤ 20
Para mostrar todas las soluciones que satisfacen esta relación, empezamos trazando la
gráfica de la recta que corresponde a la ecuación
0.4𝐹 + 0.5𝑆 ≤ 20
20
Figura 1
Gráfica con dos soluciones para el problema de 2 variables
Trazamos la gráfica de esta ecuación al identificar dos puntos que satisfagan esta ecuación
y luego trazar una recta que pase por los puntos.
0.4𝐹 + 0.5𝑆 ≤ 20
21
¿Puede identificar todas las soluciones que satisfacen esta restricción? Primero observe
que cualquier punto de la recta 0.4𝐹 + 0.5𝑆 = 20 debe cumplir con las restricciones. Pero,
¿Dónde están las soluciones que satisfacen 0.4𝐹 + 0.5 < 20? Considere dos soluciones
(𝐹 = 10, 𝑆 = 10) y (𝐹 = 40, 𝑆 = 30).
La figura 2 muestra que la primera solución está por debajo de la recta de restricción y la
segunda solución está por encima de la misma. ¿Cuál de estas soluciones satisface la
restricción del material 1? Para (𝐹 = 10, 𝑆 = 10) tenemos
Las 31 toneladas son mayores que las 20 toneladas disponibles, así que la solución 𝐹 =
10, 𝑆 = 10 no satisfacen la restricción.
Figura 2
Recta de restricción del material 1
Por lo anterior, solo se necesita evaluar una solución para determinar cuál lado de una recta
de restricción proporciona soluciones que satisfacen la restricción. A continuación, se
identifican todas las restricciones que satisfacen la restricción del material 2:
0.2𝑆 ≤ 5
23
Figura 3
Soluciones que satisfacen la restricción del material 2
Asimismo, podemos determinar las soluciones que satisfacen la restricción del material 3.
La figura 4 muestra el resultado.
24
Figura 4
Soluciones que satisfacen la restricción del material 3.
Ahora existen 3 gráficas separadas que indican las soluciones que satisfacen cada una de
las tres restricciones. En un problema de programación lineal, debemos identificar las
soluciones que satisfacen todas las restricciones simultáneamente. Para encontrar estas
soluciones, podemos trazar las tres restricciones en una gráfica y observar el área que
contiene los puntos que sí satisfacen todas las restricciones de forma simultánea. Las
gráficas de las figuras 2, 3 y 4 pueden superponerse para obtener una gráfica con las tres
restricciones.
La figura 5 muestra esta gráfica de restricciones combinadas; la región sombreada incluye
cada punto de solución que satisface todas las restricciones de forma simultánea. Como
las soluciones cumplen con todas las restricciones de forma simultánea se llama soluciones
factibles, la región sombreada se llama región de soluciones factible.
Cualquier punto en el límite de la región factible o dentro de ésta es un punto de solución
factible para el problema de programación lineal.
25
Figura 5
Región factible del problema
Una vez identificada la región factible, seguimos con el método de solución gráfica y obtener
la solución gráfica y obtener la solución óptima para el problema de Quimes. Es importante
mencionar que la solución óptima para un problema de programación lineal es la solución
factible que proporciona el mejor valor posible de la función objetivo.
Un método para encontrar la solución óptima sería evaluar la función objetivo para cada
solución factible; la solución óptima por ende es la que produce el valor más grande. La
dificultad con este método es que el número finito de soluciones factibles vuelve imposible
la evaluación de todas las soluciones factibles.
Por lo anterior, este procedimiento de prueba y error no se puede utilizar para identificar la
solución óptima.
26
En vez de tratar de calcular la contribución a las utilidades para cada solución factible,
seleccionamos un valor arbitrario para la contribución a las útil para la contribución a las
utilidades e identificamos todas las soluciones factibles que producen el valor seleccionado.
Por ejemplo, ¿cuáles soluciones factibles proporcionan una contribución a las utilidades de
$240? Estas soluciones están determinadas por los valores de F y S en la región factible
que proporcionara la función objetivo
Esta expresión es simplemente la ecuación de una recta. Por tanto, todas las soluciones
factibles (F, S) que producen una contribución a las utilidades de $240 deben estar en la
recta. Aprendimos antes en esta sección cómo trazar la gráfica de una recta de restricción.
Al trazar la recta que pasa por estos dos puntos identificamos todas las soluciones que
tienen una contribución a las utilidades de $240. En la figura 6 se presenta una gráfica de
esta recta de utilidades, la cual muestra que un número infinito de combinaciones de
producción factible proporcionara una contribución a las utilidades de $240.
El objetivo es encontrar la solución factible que produzca la más alta contribución a las
utilidades, así que proseguimos con la selección de contribuciones a las mayores utilidades
y encontramos las soluciones que producen los valores establecidos.
27
Figura 6
Recta de utilidades de $240 para el problema
El objetivo es encontrar la solución factible que produzca la más alta contribución a las
utilidades, así que proseguimos con la selección de contribuciones a las mayores utilidades
y encontramos las soluciones que producen los valores establecidos. Por ejemplo, ¿cuáles
soluciones proporcionan una contribución a las utilidades de $720? ¿Cuáles soluciones
proporcionan una contribución a las utilidades de $1,200? Para responder estas preguntas
debemos determinar los valores de F y S que están en las rectas de utilidades:
28
¿Podemos encontrar una solución factible que produzca una contribución a las utilidades
de incluso mayor? Observe con atención la figura 7 y haga algunas observaciones
generales sobre las rectas de utilidades. Deberá poder identificar las propiedades
siguientes: 1) las rectas de utilidades son paralelas entre sí, y 2) Las rectas de utilidades
con contribuciones a las utilidades mayores están más alejadas del origen.
Figura 7
Rectas de utilidades paralelas
Como las rectas de utilidades son paralelas y aquellas que son mayores están más alejadas
del origen, podemos obtener soluciones que producen valores cada vez mayores para la
función objetivo al continuar alejando la recta de utilidades del origen, pero manteniéndola
paralela a las demás rectas de utilidades. Sin embargo, se llegará a un punto en que el
alejamiento coloque a la recta de utilidades completamente fura de la región factible.
29
Dado que los puntos fuera de la región factible son inaceptables, el punto de la región
factible que se encuentra en la recta de utilidades mayor es una solución óptima. El punto
de la solución óptima, se muestra en la figura 8. Los valores óptimos de las variables de
decisión son los valores de F y S en este punto.
Figura 8
Solución óptima para el problema
Remítase a la figura 5 y observe que el punto de solución óptima del problema está en la
intersección de las rectas de restricción del material 1 y del material 3. Es decir, la solución
óptima está tanto en la recta de restricción del material 1,
0.4𝐹 + 0.5𝑆 = 20
0.6𝐹 + 0.3𝑆 = 21
30
Por tanto, los valores de las variables de decisión F y S deben satisfacer ambas ecuaciones
de manera simultánea. Al utilizar la ecuación 0.4𝐹 + 0.5𝑆 = 20 y calcular F obtenemos
𝐹 = 50 − 1.25𝑆
𝐹 = 50 − 1.25(20)
= 50 − 25 = 25
La ubicación exacta del punto de solución óptima es 𝐹 = 25 y 𝑆 = 20. Este punto de solución
proporciona las cantidades de producción óptima para Quimes a 25 toneladas de aditivo
para combustible y 20 toneladas de base para solvente y produce una contribución a las
utilidades de 40(25) + 30(20) = $1600. Para un problema de programación lineal con 2
variables de decisión, usted puede determinar los valores exactos de estas variables para
la solución óptima al utilizar primero el procedimiento grafico para identificar el punto de
solución óptima y luego resolver las dos ecuaciones simultaneas asociadas a ese punto.
31
Visualiza el siguiente video, te será de mucha ayuda para reforzar el concepto de resolución
de problemas de minimización.
Abel Esteban Ortega Luna. (22 de marzo de 2014). Programación Lineal Minimización.
[Video]. https://www.youtube.com/watch?v=MK29boPJDcg
32
Clase 7| Método simplex
7. Método simplex
7.1 Método simplex
El método gráfico presentado anteriormente demuestra que la solución óptima de un
modelo de programación lineal está siempre asociada a un punto extremo o esquina del
espacio de soluciones como se muestra en la figura 9. Esta idea conduce precisamente a
la creación del método simplex. En términos generales lo que hace el método simplex es
trasladar la definición geométrica del punto extremo a una definición algebraica.
Figura 9
Conjunto de intersecciones
33
La búsqueda se hace siempre a través de los lados del polígono (o de las aristas del
poliedro, si el número de variables es mayor). Cómo el número de vértices (y de aristas) es
finito, siempre se podrá encontrar la solución.
El método del simplex fue creado en 1947 por el matemático George Dantzig, para
resolver problemas de programación lineal en los que intervienen más de dos
variables. El álgebra matricial y el proceso de eliminación de Gauss-Jordán para
resolver un sistema de ecuaciones lineales constituyen la base del método simplex
(Taha, 2004, p. 71).
3. Iguale la función objetivo a cero (tenga en cuenta la tabla 8 para completar los
coeficientes objetivos).
34
Condición de Factibilidad: tanto en los problemas de maximización como de
minimización, la variable saliente es la variable básica actual, con la menor razón
(con denominador positivo distinto de cero) que resulta al dividir los valores del lado
derecho entre el valor respectivo de la columna de entrada.
Tabla 8
Criterios para convertir las desigualdades en igualdades.
Tipo de Restricción Problema de Maximizar Problema de Minimizar
Restricción Menor o igual ≤: Variable de holgura toma Variable de holgura toma
Requiere una variable de coeficiente cero en la ecuación coeficiente cero en la
holgura positiva. objetivo. ecuación objetivo.
Variable de holgura toma Variable de holgura toma
Restricción Mayor o Igual ≥:
coeficiente cero en la ecuación coeficiente cero en la
Requiere una variable de
objetivo. La variable artificial ecuación objetivo. La variable
holgura negativa y la suma de
toma coeficiente -M en la artificial tomo coeficiente +M
una variable artificial.
ecuación objetivo. en la ecuación objetivo.
Restricción de igualdad: La Variable artificial tomo La Variable artificial tomo
Requiere que se le sume una coeficiente -M en la ecuación coeficiente +M en la ecuación
variable artificial. objetivo. objetivo.
Fuente: Nota. Elaboración propia
35
El siguiente video muestra de manera introductoria la resolución de un
problema maximización por medio del método simplex.
Universidad Continental – Modalidad a Distancia. (19 de septiembre de 2017).
Resolución de problemas de PL con el Método Simplex. [Video].
https://youtu.be/8BieHhoc6-M
Se introduce una variable de holgura por cada una de las restricciones, para convertirlas en
igualdades, observe que todas las restricciones son del tipo ≤ , esto implica de acuerdo a
la tabla 8 que se debe agregar una variable de holgura positiva.
2𝑥1 + 𝑥2 + 𝑆1 = 18
2𝑥1 + 3𝑥2 + 𝑆2 = 42
3𝑥1 + 𝑥2 + 𝑆3 = 24
36
3. Igualar la función objetivo a cero
𝑍 − 3𝑥1 − 2𝑥2 = 0
En las columnas aparecerán todas las variables del problema y, en las filas, los coeficientes
de las igualdades obtenidas, una fila para cada restricción y la primera fila con los
coeficientes de la función objetivo:
Tabla 9
Iteración n1 del problema
Base Variable de decisión Variable de holgura Valores solución
𝑥1 𝑥2 𝑆1 𝑆2 𝑆3
Z -3 -2 0 0 0 0
𝑆1 2 1 1 0 0 18
𝑆2 2 3 0 1 0 42
𝑆3 3 1 0 0 1 24
Fuente: Nota. Elaboración propia
a. Para escoger la variable de decisión que entra en la base, nos tenemos que
ubicar en la primera fila, la de los coeficientes de la función objetivo y
escogemos la variable con el coeficiente negativo mayor (el valor absoluto). En
nuestro caso, la variable 𝑥 de coeficiente - 3.
37
Por tanto, lo que va a determinar el final del proceso de aplicación del método
simplex, es que en la primera fila no existan elementos negativos debido a que
estamos desarrollando un problema de maximización. La columna de la variable
que entra en la base se llama columna pivote (columna sombreada en la tabla
9).
b. Para encontrar la variable de holgura que tiene que salir de la base, se divide
cada término de la última columna (valores solución) por el término
correspondiente de la columna pivote, siempre que estos últimos sean mayores
que cero. En nuestro caso:
18 42 24
= 9, = 21 𝑦 =8
2 2 3
Si hubiese algún elemento menor o igual que cero no se realiza dicho cociente.
En el caso de que todos los elementos fuesen menores o iguales a cero,
entonces tendríamos una solución no acotada y no se puede seguir con el
proceso de iteración. El término de la columna pivote que en la división anterior
dé lugar al menor cociente positivo, en nuestro caso el 3, ya que 8 es el menor
24
( = 8), indica la fila de la variable de holgura que sale de la base,
3
𝑆3 . Esta fila se llama fila pivote (la última fila sombreada en la tabla 1). Si al
calcular los cocientes, dos o más son iguales, indica que cualquiera de las
variables correspondientes puede salir de la base.
38
También se puede hacer utilizando el siguiente esquema:
Los nuevos coeficientes de 𝑥1 (recordemos que es nuestra variable que entra), se obtienen
dividiendo todos los coeficientes de la fila 𝑆3 en la tabla 9 por el pivote operacional, 3, que
es el que hay que convertir en 1 (ver calculo efectuado en tabla 10, fila 𝑥1 ).
𝑁𝑢𝑒𝑣𝑎 𝑓𝑖𝑙𝑎 = 𝑣𝑖𝑒𝑗𝑎 𝑓𝑖𝑙𝑎 − (𝑐𝑜𝑒𝑓𝑖𝑐𝑖𝑒𝑛𝑡𝑒 𝑑𝑒 𝑙𝑎 𝑣𝑖𝑒𝑗𝑎 𝑓𝑖𝑙𝑎 𝑒𝑛 𝑙𝑎 𝑐𝑜𝑙𝑢𝑚𝑛𝑎 𝑑𝑒 𝑙𝑎 𝑣𝑎𝑟𝑖𝑎𝑏𝑙𝑒 𝑒𝑛𝑡𝑟𝑎𝑛𝑡𝑒)
∗ (𝑁𝑢𝑒𝑣𝑎 𝑓𝑖𝑙𝑎 𝑑𝑒𝑙 𝑝𝑖𝑣𝑜𝑡𝑒)
Una vez calculada la fila pivote proveniente de la tabla 9 en la tabla 10, Puede observar
dicho calculo en la fila 𝑥1 . Se realizará el cálculo de la nueva fila correspondiente a la
variable 𝑆2 .
Vieja fila de s 2 3 0 1 0 42
- - - - - -
Coeficiente 2 2 2 2 2 2
X X X X X X
Nueva fila pivote 1 1/3 0 0 1/3 8
= = = = = =
Nueva fila de s 0 7/3 0 1 -2/3 26
39
Tabla 10
Iteración n2 del problema
Variable de decisión Variable de holgura
Base Valores solución
𝑥1 𝑥2 𝑆1 𝑆2 𝑆3
Z 0 -1 0 0 1 24
𝑆1 0 1/3 1 0 -2/3 2
𝑆2 0 7/3 0 1 -2/3 26
𝑥1 1 1/3 0 0 1/3 8
Fuente: Nota. Elaboración propia
Como en los elementos de la primera fila hay una variable negativa -1, significa que no
hemos llegado todavía a la solución óptima. Se debe repetir el proceso:
b. Para calcular la variable que sale, dividimos los términos de la última columna
entre los términos correspondientes de la nueva columna pivote; como
el menor cociente positivo es 6, tenemos que la variable de holgura quesale es
𝑆1 .
2 26 78 8
= 6, = 𝑦 =8
1⁄3 7⁄3 7 1⁄3
40
Operando de forma análoga a la anterior obtenemos la nueva tabla:
Tabla 11
Iteración n3 del problema
Variable de decisión Variable de holgura
Base Valores solución
𝑥1 𝑥2 𝑆1 𝑆2 𝑆3
Z 0 0 3 0 -1 30
𝑥2 0 1 3 0 -2 6
𝑆2 0 0 -7 0 4 12
𝑥1 1 0 -1 0 1 6
Fuente: Nota. Elaboración propia
En los elementos de la primera fila existen valores negativos, -1, significa que no hemos
llegado a la solución óptima. Se debe repetir el proceso:
b. Para calcular la variable que sale, dividimos los términos de la última columna entre
los términos correspondientes de la nueva columna pivote, y como el menor
cociente positivo es 3, tenemos que la variable de holgura quesale es 𝑆2 .
6 12 6
= −3, =3𝑦 =6
−2 4 1
41
Obtenemos la tabla:
Tabla 12
Finalización de las iteraciones
Variable de decisión Variable de holgura
Base Valores solución
𝑥1 𝑥2 𝑆1 𝑆2 𝑆3
Z 0 0 5/4 0 0 33
𝑥2 0 1 -1/2 0 0 12
𝑆3 0 0 -7/4 0 1 3
𝑥1 1 0 -3/4 0 0 3
Fuente: Nota. Elaboración propia
Como todos los coeficientes de la fila de la función objetivo son positivos, hemos llegado a
la solución óptima. La solución óptima viene dada por el valor de Z en la columna de los
valores solución en nuestro caso: 33. En la misma columna se puede observar el vértice
donde se alcanza, observando las filas correspondientes a las variables de decisión que
han entrado en la base: D(3,12).
Para ampliar el tema del método simplex, en este video podrás visualizar la resolución de
un ejercicio sencillo de minimización por el método simplex.
Pedro Orazzi. (30 de agosto de 2020). 9 método simplex minimizar ejemplo. [Video].
https://www.youtube.com/watch?v=_H8F07Mw6Ls
42
7.2 El método de la “M”
En este apartado se abordarán los problemas que introducen las otras formas de
restricciones funcionales (≥ o =), con ellas será necesario identificar una solución inicial
básica factible que no es necesariamente el origen. Anteriormente esta solución inicial
se encontraba en forma muy conveniente al hacer que las variables de holgura fueran
las variables básicas iniciales, donde cada una era igual a la constante no negativa del
lado derecho de la ecuación correspondiente.
𝑀𝑖𝑛𝑖𝑚𝑖𝑧𝑎𝑟 𝑍 = 4𝑥1 + 𝑥2
C.S.R
3𝑥1 + 𝑥2 = 3
4𝑥1 + 3𝑥2 ≥ 6
𝑥1 + 2𝑥2 ≤ 18
𝑥1 , 𝑥2 ≥ 0
43
1. Se aplica la técnica de las variables artificiales, introduciendo una variable no negativa
(en este caso será 𝑅1 ) en la ecuación 3𝑥1 + 𝑥2 = 3, como si fuese una variable de
holgura:
3𝑥1 + 𝑥2 + 𝑅1 = 3
2. En la restricción 4𝑥1 + 3𝑥2 ≥ 6, se aplica de igual forma la técnica, se resta una variable
de holgura, para este caso será 𝑆1 y se suma una variable artificial, 𝑅2 ya que posee
una desigualdad del tipo ≥.
4𝑥1 + 3𝑥2 − 𝑆1 + 𝑅2 = 6
3. Se asigna una penalización enorme a la función objetivo por el hecho de agregar las
variables 𝑅1 , 𝑅2 ≥ 0, en las ecuaciones anteriores.
Donde M simbólicamente representa un número positivo muy grande. Este método que
fuerza a R1 𝑦 R 2 hasta el nivel de R1 , R 2 = 0 en la solución óptima se llama método de la M.
Para el caso de minimización, penalizamos a la variable artificial, colocándola en la función
objetivo con un coeficiente +𝑀 (en el lado derecho el coeficiente es positivo).
Con los artificios realizados anteriormente se debe encontrar la solución óptima para el
problema real aplicando el método simplex al problema artificial.
44
Con el modelo en su forma canónica pasamos a colocar los coeficientes a la tabla simplex:
Tabla 13
Tabla simplex inicial
Variables de
Variables no básicas Valores
Base decisión
solución
𝑥1 𝑥2 𝑠1 𝑠2 𝑅1 𝑅2
z -4 -1 0 0 -M -M 0
𝑅1 3 1 0 0 1 0 3
𝑅2 4 3 -1 0 0 1 6
𝑠2 1 2 0 1 0 0 4
Fuente: Nota. Elaboración propia
Recordemos que la condición que debe cumplir toda tabla del método simplex nos dice que
toda variable básica debe tener un 1 en la intersección de su renglón y columna
correspondiente y cero en los demás renglones incluido el renglón de Z en otras palabras,
que toda variable que sea básica solamente debe aparecer en el renglón de la restricción
que representa. Para hacer cero los coeficientes -M, utilizamos el renglón de
𝑅1 𝑦 𝑅2 respectivamente como renglón pivote multiplicándolo por M y sumando el resultado
al renglón de Z en cada operación.
(M) (3 1 0 0 1 0 3)
= 3M M 0 0 M 0 3M
45
Para la variable 𝑅2 de igual manera se multiplica el renglón respectivo a esta variable por
M:
(M) (4 3 -1 0 0 1 6)
= 4M 3M -M 0 0 M 6M
𝑅1 (𝑀) 3M M 0 0 M 0 3M
𝑅2 (𝑀) 4M 3M -M 0 0 M 6M
Z -4 -1 0 0 -M -M 0
7M-4 4M-1 -M 0 0 0 9M
Tabla 14
Tabla simplex inicial para realizar iteración 1
Variables de
Variables no básicas Valores
Base decisión
solución
𝑥1 𝑥2 𝑠1 𝑠2 𝑅1 𝑅2
z 7M-4 4M-1 -M 0 0 0 9M
𝑅1 3 1 0 0 1 0 3
𝑅2 4 3 -1 0 0 1 6
𝑠2 1 2 0 1 0 0 4
Fuente: Nota. Elaboración propia
Tabla 15
Iteración 2
Variables de
Variables no básicas Valores
Base decisión
solución
𝑥1 𝑥2 𝑠1 𝑠2 𝑅1 𝑅2
1 + 5𝑀 4 − 7𝑀
z 0 -M 0 0 4+2M
3 3
𝑥1 1 1/3 0 0 1/3 0 1
𝑅2 0 5/3 -1 0 -4/3 1 2
𝑠2 0 5/3 0 1 -1/3 0 3
Fuente: Nota. Elaboración propia
46
Tabla 16
Iteración 3
Variables de
Variables no básicas Valores
Base decisión
solución
𝑥1 𝑥2 𝑠1 𝑠2 𝑅1 𝑅2
4 − 5𝑀 −1 − 5𝑀
z 0 0 1/5 0 5 18/5
5
𝑥1 1 0 1/5 0 3/5 -1/5 3/5
𝑥2 0 1 -3/5 0 -4/5 3/5 6/5
𝑠2 0 0 1 1 1 -1 1
Fuente: Nota. Elaboración propia
Tabla 17
Tabla final del problema
Variables de
Variables no básicas Valores
Base decisión
solución
𝑥1 𝑥2 𝑠1 𝑠2 𝑅1 𝑅2
7 − 5𝑀
z 0 0 0 -1/5 -M 17/5
5
𝑥1 1 0 0 -1/5 2/5 0 2/5
𝑥2 0 1 0 3/5 -1/5 0 9/5
𝑠2 0 0 1 1 1 -1 1
Fuente: Nota. Elaboración propia
En la tabla 17 se observa que en la fila de Z todos los valores son negativos, esto nos dice
2 9
que ya se ha optimizado el problema, donde 𝑥1 = 5 𝑦 𝑥2 = 5, generan un valor de
47
Clase 8| Método dual simplex
8.Metodo dual simplex
8.1 Definición del problema
El problema de programación lineal dual que se define a partir de un problema original que
denominaremos como primal, comparte con el los mismos coeficientes tanto de la función
objetivo de las restricciones, pero en diferente posición como más adelante se especificara.
Se debe considerar otros aspectos importantes del análisis de dualidad que es necesario
tener presentes, por ejemplo, para el modelo primal el objetivo es minimizar y las
restricciones son del tipo ≥, mientras que para el modelo dual el objetivo es maximizar y las
restricciones son del tipo ≤, lo anterior se ilustra en la tabla 1 donde se aprecian estas
diferencias en las representaciones matriciales de estos modelos.
48
Tabla 18
Problema primal y dual asociado a un mismo problema real
Modelo Primal Modelo Dual
𝑍𝑚𝑖𝑛 = 𝐶𝑋 𝑍𝑚𝑎𝑥 = 𝐵𝑌
Sujeto a: Sujeto a:
𝐴𝑋 ≥ 𝐵 𝐴𝑇 𝑌 ≤ 𝐶
𝑋≥0 𝑌≥0
Fuente: Nota. Elaboración propia
Problema primal
𝑍𝑚𝑖𝑛 = 𝑐1 𝑥1 + 𝑐2 𝑥2 + ⋯ + 𝑐𝑛 𝑥𝑛
𝑆𝑢𝑗𝑒𝑡𝑜 𝑎:
𝑎11 𝑥1 + 𝑎12 𝑥2 + ⋯ + 𝑎1𝑛 𝑥𝑛 ≥ 𝑏1
𝑎21 𝑥1 + 𝑎22 𝑥2 + ⋯ + 𝑎2𝑛 𝑥𝑛 ≥ 𝑏2
…
𝑎𝑚1 𝑥1 + 𝑎𝑚2 𝑥2 + ⋯ + 𝑎𝑚 𝑥𝑛 ≥ 𝑏𝑚
𝑥1 , 𝑥2 , … , 𝑥𝑛 ≥ 0
49
Problema dual
𝑍𝑚𝑎𝑥 = 𝑏1 𝑦1 + 𝑏2 𝑦2 + ⋯ + 𝑏𝑚 𝑦𝑚
𝑆𝑢𝑗𝑒𝑡𝑜 𝑎:
𝑎11 𝑦1 + 𝑎21 𝑦2 + ⋯ + 𝑎𝑚1 𝑦𝑚 ≥ 𝑐1
𝑎12 𝑦1 + 𝑎22 𝑦2 + ⋯ + 𝑎𝑚2 𝑦𝑚 ≥ 𝑐2
𝑎13 𝑦1 + 𝑎23 𝑦2 + ⋯ + 𝑎𝑚3 𝑦𝑚 ≥ 𝑐3
…
𝑎1𝑛 𝑦1 + 𝑎2𝑛 𝑦2 + ⋯ + 𝑎𝑚𝑛 𝑦𝑚 ≥ 𝑐𝑛
𝑦1 , 𝑦2 , … , 𝑦𝑚 ≥ 0
Cabe resaltar que si el modelo primal tiene 𝑛 variables y 𝑚 restricciones, el modelo dual
tendrá 𝑚 variables y 𝑛 restricciones, como se puede observar en los problemas primal y
duales desarrollados. La solución óptima de un problema corresponde a la solución del otro,
como ya se ha mencionado.
Cada solución óptima en el problema primal tiene una solución básica optima
complementaria en el problema dual, en donde los valores respectivos de las funciones
objetivos son iguales, 𝑍𝑚𝑎𝑥 = 𝑍𝑚𝑎𝑥 y el valor de las restricciones del problema primal serán
los coeficientes de las variables artificiales.
Esto quiere decir que los valores de las restricciones del modelo primal son los valores que
se encuentran en el reglón 𝑅0 en las columnas de las variables artificiales de la tabla simplex
del modelo dual.
A continuación, se presenta el algoritmo para generar el modelo dual a partir del problema
primal, para lo cual utilizaremos la tabla primal dual:
50
8.2 Método dual simplex
Partiendo de un modelo con el objetivo de minimizar, se utiliza la llamada tabla primal dual
para trasladar el problema a maximización y volver al uso del mismo algoritmo. Si tenemos
el modelo primal de minimización:
Problema primal
𝑍𝑚𝑖𝑛 = 𝑐1 𝑥1 + 𝑐2 𝑥2 + ⋯ + 𝑐𝑛 𝑥𝑛
Sujeto a:
3. En la primera columna a partir del segundo renglón se escriben los nombres de las
variables del problema dual (m variables artificiales).
51
4. En el primer renglón se escriben los nombres de las variables del problema primal
(n variables).
Tabla 19
Tabla para relacionar los modelos dual primal.
Primal
𝒙𝟏 𝒙𝟐 … 𝒙𝒏 ≥
Dual
𝑦1
𝑦2
…
𝑦𝑚
≤
Fuente: Nota. Elaboración propia
6. Se escribe cada uno de los coeficientes de las restricciones del problema primal en
forma horizontal, ocupando los renglones de la tabla.
8. De esta tabla podemos obtener el modelo dual, lo único que debemos hacer es leer
el modelo de manera vertical y los coeficientes de la función objetivo se
obtienen de la última columna.
52
A continuación, se desarrolla un modelo atendiendo los puntos anteriores:
𝑥1 + 4𝑥2 + 6𝑥3 ≥ 9
7𝑥1 + 3𝑥2 + 2𝑥3 ≥ 7
𝑥1 , 𝑥2 , … , 𝑥𝑛 ≥ 0
El problema primal tiene 3 variables y 2 restricciones, por lo tanto, el problema dual tendrá
2 variables y 3 restricciones. El tamaño de la tabla es de 2 + 2 = 4 renglones y 3 + 2 = 5
columnas para este caso.
3. En la primera columna a partir del segundo renglón se escriben los nombres de las
variables del problema dual (2 variables artificiales).
4. En el primer renglón se escriben los nombres de las variables del problema primal.
53
Tabla 20
Tabla inicial de los modelos dual-primal
Primal
𝒙𝟏 𝒙𝟐 𝒙𝟑 ≥
Dual
𝑦1
𝑦2
≤
Fuente: Nota. Elaboración propia
Tabla 21
Coeficientes del modelo primal en el ultimo reglón
Primal
𝒙𝟏 𝒙𝟐 𝒙𝟑 ≥
Dual
𝑦1
𝑦2
≤ 5 8 1
Fuente: Nota. Elaboración propia
6. Se escribe cada uno de los coeficientes de las restricciones del problema primal en
forma horizontal, ocupando los renglones de la tabla.
Tabla 22
Coeficientes de las restricciones del problema primal
Primal
𝒙𝟏 𝒙𝟐 𝒙𝟑 ≥
Dual
𝑦1 1 4 6
𝑦2 7 3 2
≤ 5 8 1
Fuente: Nota. Elaboración propia
54
7. En la última columna se escriben las cantidades limitantes de las restricciones del
modelo primal.
Tabla 23
Cantidades limitantes de las restricciones en la última columna
Primal
𝒙𝟏 𝒙𝟐 𝒙𝟑 ≥
Dual
𝑦1 1 4 6 9
𝑦2 7 3 2 7
≤ 5 8 1
Fuente: Nota. Elaboración propia
8. De esta tabla podemos obtener el modelo dual, lo único que debemos hacer es
leerlo y de manera vertical, y los coeficientes de la función objetivo se obtienen de
la última columna.
Tabla 24
Obtención del modelo dual a partir de la tabla.
Primal
𝒙𝟏 𝒙𝟐 𝒙𝟑 ≥
Dual
𝑦1 1 4 6 9
𝑦2 7 3 2 7
≤ 5 8 1
𝑍𝑚𝑎𝑥 = 9𝑦1 + 7𝑦2
Sujeto a:
𝑦1 + 7𝑦2 ≤ 5
4𝑦1 + 3𝑦2 ≤ 8
6𝑦1 + 2𝑦2 ≤ 1
𝑦1 , 𝑦2 , 𝑦3 ≥ 0
El resultado obtenido es el modelo dual con el objetivo de maximizar y, así como las
restricciones del tipo ≤.
𝑍𝑚𝑎𝑥 = 9𝑦1 + 7𝑦2
55
Sujeto a:
𝑦1 + 7𝑦2 ≤ 5
4𝑦1 + 3𝑦2 ≤ 8
6𝑦1 + 2𝑦2 ≤ 1
𝑦1 , 𝑦2 , 𝑦3 ≥ 0
4𝑥1 + 2𝑥2 ≥ 80
3𝑥1 + 4𝑥2 ≥ 90
𝑥1 , 𝑥2 ≥ 0
56
La tabla primala dual está dada por:
Tabla 25
Tabla de relación dual
Primal
𝒙𝟏 𝒙𝟐 ≥
Dual
𝑦1 4 2 80
𝑦2 3 4 90
≤ 12 7
Fuente: Nota. Elaboración propia
2𝑦1 + 4𝑦2 ≤ 7
𝑦1 , 𝑦2 ≥ 0
Ahora se utiliza el método dual simplex para resolver el modelo dual. La tabla inicial de
este modelo es:
Tabla 26
Tabla inicial de problema
Variables
Z 𝒚𝟏 𝒚𝟐 𝒔𝟏 𝒔𝟐 Solución
básicas
Z 1 -80 -90 0 0 0
𝒔𝟏 0 4 3 1 0 12
𝒔𝟐 0 2 4 0 1 7
Fuente: Nota. Elaboración propia
57
El primer elemento pivote de la tabla inicial es:
Tabla 27
Iteración 1 del modelo
Variables
Z 𝒚𝟏 𝒚𝟐 𝒔𝟏 𝒔𝟐 Solución
básicas
Z 1 -80 -90 0 0 0
𝒔𝟏 0 4 3 1 0 12
𝒔𝟐 0 2 4 0 1 7
Fuente: Nota. Elaboración propia
Tabla 28
Iteración 2 del modelo
Variables
Z 𝒚𝟏 𝒚𝟐 𝒔𝟏 𝒔𝟐 Solución
básicas
Z 1 -35 0 0 22.5 157.5
𝒔𝟏 0 2.5 0 1 -0.75 6.75
𝒚𝟐 0 0.5 1 0 0.25 1.75
Fuente: Nota. Elaboración propia
Tabla 29
Iteración 3 del modelo
Variables
Z 𝒚𝟏 𝒚𝟐 𝒔𝟏 𝒔𝟐 Solución
básicas
Z 1 -35 0 0 22.5 157.5
𝒔𝟏 0 2.5 0 1 -0.75 6.75
𝒚𝟐 0 0.5 1 0 0.25 1.75
Fuente: Nota. Elaboración propia
58
De manera similar con el primer pivote, con operaciones básicas entre renglones se
resuelve la tabla simplex:
Tabla 30
Finalización de iteraciones
Variables
Z 𝒚𝟏 𝒚𝟐 𝒔𝟏 𝒔𝟐 Solución
básicas
Z 1 0 0 14 12 252
𝒚𝟏 0 1 0 0.4 -0.3 2.7
𝒚𝟐 0 0 1 -0.2 0.4 0.4
Fuente: Nota. Elaboración propia
Tabla 31
Resultados de las variables
Variables
Z 𝒚𝟏 𝒚𝟐 𝒔𝟏 𝒔𝟐 Solución
básicas
Z 1 0 0 14 12 252
𝒚𝟏 0 1 0 0.4 -0.3 2.7
𝒚𝟐 0 0 1 -0.2 0.4 0.4
Por lo que al transferir la solución de la tabla simplex a las variables originales del
problema primal se tiene:
𝒔𝟏 = 𝒙𝟏 = 14
𝒔𝟐 = 𝒙𝟐 = 12
𝒁𝒎𝒂𝒙 = 𝒁𝒎𝒊𝒏 = 252
59
En el próximo video se desarrollará un problema de minimización con el método dual
simplex, te ayudará a reforzar la temática.
Red Tutores Programación Lineal. (3 de noviembre de 2018). Minimización por el
método simplex dual. [Video]. https://www.youtube.com/watch?v=R6qj-A6VXvQ
60
Referencias citadas en UNIDAD 2
61
Glosario de los términos citados en la
UNIDAD 2
62