Cap 4 Dualidad y Sensibilidad
Cap 4 Dualidad y Sensibilidad
Cap 4 Dualidad y Sensibilidad
CAPÍTULO 4
DUALIDAD Y SENSIBILIDAD EN PL
1. INTRODUCCIÓN
En el capítulo anterior se describieron las características de los modelos
de programación lineal, así como los diferentes caminos a partir de los
cuales encontrar la solución: resolución gráfica y algoritmo Simplex.
En este capítulo, se describen dos técnicas relacionadas con la
programación lineal: la dualidad y el análisis de sensibilidad.
En la primera parte se desarrolla la teoría asociada a la dualidad: cómo
se obtiene el dual de un programa lineal, la interpretación del concepto
de precio sombra y una serie de teoremas y resultados útiles para la
interpretación de un modelo lineal.
En la segunda parte se muestran las posibilidades del análisis de
sensibilidad en la programación lineal. Se tratará de analizar cómo varía
la solución del modelo (tanto el valor de la función objetivo como el
valor de las variables de decisión) en función de dos conjuntos de
parámetros del modelo: los coeficientes de la función objetivo y los
términos independientes de las restricciones.
5 y1 + 4 y2 + 6 y3 8
Observemos que, con los mismos datos del problema original, hemos
formulado otro PL que brindará información sobre el “valor” que los
recursos tienen para la empresa, lo que en economía se conoce como
precio sombra o valor marginal del recurso.
Si ahora resolvemos ambos problemas, observamos que los valores
óptimos son iguales; esto era de esperarse ya que la empresa no
aceptaría, por la venta de sus insumos, menos dinero del que podría
obtener si los utiliza en la fabricación de sus productos.
EL PROBLEMA DUAL
Decimos entonces que, para cada problema de programación lineal
existe siempre asociado al mismo, otro problema lineal. A este nuevo
programa, se lo puede emplear para obtener la solución del problema
original y además sus variables proporcionan información útil acerca de
la solución óptima del problema lineal original.
DUALIDAD Y SENSIBILIDAD 111
Cabe aclarar que Bazaraa et. al. (1981) demuestran que el dual del dual
es el primal, por lo cual las definiciones dadas se pueden aplicar al revés
y los términos “primal” y “dual” son relativos al marco de referencia que
se seleccione.
No
Canónica ≤ → ≥0
Negativa
Restricción No Variable
≥ → ≤0 No Positiva
Canónica
No
Igualdad = → n/r
Restringida
No
≥0 → ≥ Canónica
Negativa
No
Variable No Positiva ≤0 → ≤ Restricción
Canónica
No
n/r → = Igualdad
Restringida
DUALIDAD Y SENSIBILIDAD 113
Veamos un ejemplo:
Z * (B´Y*)
= = yi *
bi bi
Por ejemplo, si la i-ésima restricción representa la disponibilidad de b i
unidades de insumo para elaborar un producto y Z representa la
contribución total a las utilidades, entonces yi (variable dual) representa
el incremento en las utilidades por adicionar una unidad del i-ésimo
insumo.
Si en cambio, la i-ésima restricción representa la demanda de al menos
bi unidades producidas y Z representa el costo total de producción,
entonces yi es el costo incremental de producir una unidad más del i-
ésimo producto.
Económicamente puede interpretarse al vector de variables duales Y*
como un vector de precios sombra para el vector del lado derecho, es
decir que es el precio justo o valoración interna de los recursos.
3. ANÁLISIS DE SENSIBILIDAD
Uno de los supuestos sobre los que está basada la PL es el de
certidumbre. Es decir que el modelo supone que todos los parámetros
que en él intervienen se conocen con exactitud.
Nosotros sabemos que en los problemas que se nos presentan
diariamente existe un grado de incertidumbre o aleatoriedad en los
datos que poseemos. Por ejemplo, podemos estimar que la
disponibilidad de horas de mano de obra mensuales es en promedio de
500, pero cada mes en particular la cantidad real de horas disponibles
pueden no ser exactamente 500, aunque sí un valor muy aproximado.
En general, al modelizar, se utilizan estimaciones estadísticas de los
parámetros del modelo y luego se los trabaja como valores ciertos.
Debido a esto se hace necesario realizar un análisis de la sensibilidad de
la solución del problema. Esto es, estudiar los efectos que tienen en la
solución óptima del problema, variaciones que puedan producirse en los
valores de los parámetros.
DUALIDAD Y SENSIBILIDAD 115
➢ un análisis gráfico
➢ estudio de los intervalos de sensibilidad
➢ cálculo de los intervalos.
La ecuación explícita
Supongamos que tenemos el siguiente problema: de la recta de
isoutilidad es:
Max Z = 12 x1 +20 x2 x2=(Z/c2) – (c1/c2) x1
Sa : 10 x1 300 de donde, cualquier
cambio en los
12 x2 360
coeficientes cj,
15 x1 + 10 x2 600 modificarán la
pendiente de la recta.
x1 y x2 ≥ 0
96
90
84
116 78
CAPÍTULO 4
72
66
54
48
42
36
A
30
24
18
x2 Payoff: 12.0 x1 + 20.0 x2 = 840.0
120
12
114
6
108
0
0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170
102
:: 15.0
: 10.00.0
x1 x1
+x1+0.0
+ 12.0
10.0x2
x2 =x2==360.0
300.0
600.0
Optimal Decisions(x1,x2): (20.0, 30.0)
: 10.0x1 + 0.0x2 <= 300.0
96
Figura 1
: 0.0x1 + 12.0x2 <= 360.0
: 15.0x1 + 10.0x2 <= 600.0
90
60
54
48
42
36
A
30
24 x2
120
18
B Payoff: 20.0 x1 + 20.0 x2 = 1000.0
114
12
108
6
102
0
0 96 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170 1
Figura 2
84 + 10.0x2 <= 600.0
: 15.0x1
78
54
48
42
36
A
30
Payoff: 35.0 x1 + 20.0 x2 = 1350.0
24
18
B
12
0
0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 17
Figura 3
Veamos, con el mismo ejemplo, que sucede cuando varían los valores
120
vértice actual deja de ser óptimo en este caso pasa del vértice A al C.
96
las variables básicas eran x1, x2 y S2, siendo las no básicas S1y S3.
78
54
48
42
x2
36 120
A C
30 114
Payoff: 12.0 x1 + 20.0 x2 = 920.0
24 108
18 102
12 96
6 90
0 84
0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 17
78
Optimal Decisions(x1,x2): (26.7, 30.0)
: 10.0x1 + 0.0x2 <= 300.0 Figura 4
72: 0.0x1 + 12.0x2 <= 360.0
: 15.0x1 + 10.0x2 <= 700.0
diferentes. 48
42
36
D A
30
Payoff: 12.0 x1 + 20.0 x2 = 760.0
24
18
12
0
0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160
Figura 5
Optimal Decisions(x1,x2): (13.3, 30.0)
: 10.0x1 + 0.0x2 <= 300.0
: 0.0x1 + 12.0x2 <= 360.0
90
84
118 CAPÍTULO 4
78
72
36
A E
30
24
18
12
6
x2
120
0
114
0 10 20 30 40 50 60 70 80 90 100 110 120 130 140
: 10.0x1 +
Figura 6
108Optimal Decisions(x1,x2): (30.0, 30.0)
0.0x2 <= 300.0
102: 0.0x1 + 12.0x2 <= 360.0
60
42
36
: 0.0 x1 + 12.0 x2 = 360.0
A
30
x2
120
24
114
18
108
12
102
6
: 15.0 x1 + 10.0 x2 = 600.0
96
0
0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 17
90
Figura 7
Optimal Decisions(x1,x2): (20.0, 30.0)
: 10.0x1 + 0.0x2 <= 400.0
: 84
0.0x1 + 12.0x2 <= 360.0
: 15.0x1 + 10.0x2 <= 600.0
36
30
A
24 : 10.0 x1 + 0.0 x2 = 200.0
18
: 0.0 x1 + 12.0 x2 = 360.0
12
0
0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 15
96
90
78
72
36
A
30
18
: 0.0 x1 + 12.0 x2 = 360.0
12
0
0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160
Figura 9
: 0.0x1 + 12.0x2 <= 360.0
: 15.0x1 + 10.0x2 <= 600.0
REGLA DE 100%
La validez de los cambios informados por el análisis de sensibilidad son
ceteris paribus, es decir, se deben analizar uno por vez. No obstante,
existe una regla práctica, conocida como regla del 100%, la cual
sostiene que "para considerar cambios simultáneos se deben sumar los
porcentajes de cambio tanto de los incrementos como de las
disminuciones permisibles; si la suma de los cambios porcentuales
no excede el 100%, la solución óptima no se modificará". Esto es
válido tanto para cambios en el vector de términos independientes de
las restricciones como en los coeficientes que preceden a las variables
en la FO.
Para que la solución actual siga siendo óptima, ningún valor cj - zj debe
hacerse positivo.
(c j ) (
- zj = c j - z j + λkjΔck )
(c j ) ( )
- zj = c j - z j - λkjΔck para j = 1, 2, ...,n
(c j ) ( )
- zj = c j - zj - λkjΔck 0 para j = 1, 2, ...,n
1. Restricciones no limitantes
Si la restricción es del tipo el intervalo de variación de bi será [Si , ],
donde Si es el valor de la variable de holgura correspondiente a la
i-ésima restricción.
Si la restricción es del tipo el intervalo de variación de bi será [- , Si],
donde Si es el valor de la variable de excedente correspondiente a la
i-ésima restricción.
En estos casos no cambiará el valor de la función objetivo, pero sí el
valor de la holgura correspondiente.
2. Restricciones limitantes
Cuando hicimos el análisis de los elementos de la tabla simplex
definimos a las tasas de sustitución (ij) como el sacrificio que se debía
hacer de la variable básica xi para poder incrementar en una unidad la
variable no básica xj.
Cuando la variable no básica analizada es una variable de holgura
debemos interpretar ese incremento como dejar libre o dejar de usar
una unidad del recurso (bi ) de la restricción a la cual corresponde dicha
holgura. Es decir que, estas tasas nos mostrarán el sacrificio de cada
variable básica si disminuimos en una unidad un valor del lado derecho.
Por el contrario, si aumentamos en una unidad un bi , la interpretación
de las ij será exactamente lo opuesto. De esta manera podemos usar a
las tasas de sustitución de las variables de holgura correspondiente al
lado derecho analizado, para determinar en cuánto puede aumentar o
disminuir su valor sin que la solución actual deje de ser factible
Para determinar el intervalo de sensibilidad usamos la condición de
factibilidad de la solución y las tasas de sustitución de las variables de
holgura. El procedimiento para determinar este intervalo en caso de
restricciones del tipo es el siguiente:
124 CAPÍTULO 4
1 1j 0
0
2 2j
3 + bi 3j 0 (1)
m mj 0
Solución actual
1 1j 0
2j 0
2
3 − bi 3j 0
(2)
m m j 0
ck − zk
Si ck − zk es positivo, convendrá introducir este nuevo producto o
actividad y si es negativo nos indicará la magnitud del cambio que
debería realizarse para poder introducirlo.
PRODUCTO
Alfombra Alfombra Alfombra Alfombra
Disponibilidad
I II III IV
Materia Prima 3 4 8 6 22000
(Kg/unid)
Hrs. Sección 8 2 4 2 28000
Teñido
Hrs. Sección 4 6 2 4 8000
Tejidos
Contribución 40 60 30 10
($/unid) Tabla 1
Definición de Variables:
x1, x2, x3 y x4 indican las unidades de las alfombras I, II, III y IV a
fabricar, respectivamente.
x1, x2, x3 y x4 0
DUALIDAD Y SENSIBILIDAD 127
1) 105000.0
Disminución que se produce en z
VARIABLE VALUE REDUCED COST por cada unidad en que se
incrementa la variable
X1 0.000000 0.500000
X2 500.000000 0.000000
X3 2500.000000 0.000000
X4 0.000000 35.000000
3) 17000.000000 0.000000
4) 0.000000 9.000000
Objetivo 0 105000
Celdas de variables
Nombre Valor original Valor final
Variables Alf. I 0 0
Variables Alf. IV 0 0
Restricciones
Nombre Valor de la celda Estado Holgura
Análisis de Sensibilidad
donde bi representa el incremento del valor del lado derecho que en
este caso es negativo.
cj 40 60 30 10 0 0 0
cj Base VLD x1 x2 x3 x4 S1 S2 S3
30 x3 2500 0,05 0 1 0,5 0,15 0 -0,10
0 S2 17000 6,5 0 0 -1 -0,5 1 0
60 x2 500 0,65 1 0 0,5 -0,05 0 0,20
zj 105000 40,5 60 30 45 1,5 0 9
cj - zj -0,5 0 0 -35 -1,5 0 -9
Tabla 2
x4= 10
S1= 0
S3= 0
Z = 105000 - (35)10 = $104650
Se trata de una solución factible no básica.
Para calcular los intervalos de sensibilidad, tanto de los coeficientes de
la función objetivo como de los valores del lado derecho, necesitamos la
tabla óptima de simplex:
cj 40 60 30 10 0 0 0
cj Base VLD x1 x2 x3 x4 S1 S2 S3
30 x3 2500 0,05 0 1 0,5 0,15 0 -0,10
0 S2 17000 6,5 0 0 -1 -0,5 1 0
60 x2 500 0,65 1 0 0,5 -0,05 0 0,20
zj 105000 40,5 60 30 45 1,5 0 9
cj - zj -0,5 0 0 -35 -1,5 0 -9
Tabla 3
i) Intervalo de sensibilidad para los coeficientes de las variables:
Como la variable x1 es no básica, el intervalo de sensibilidad de su
coeficiente se determina como:
− , c - z − , 0, 5
j j o sea
(c j ) ( )
- z j = c j - z j - λkjΔck
Entonces,
-0,5 - 0, 65Δc2 0
−35 − 0,5Δc2 0
−1,5 − (−0, 05)Δc2 0
−9 − 0,20Δc2 0
Δc2 −0,7692
Δc2 −70
Δc2 30
Δc2 −45
DUALIDAD Y SENSIBILIDAD 133
1 1j 0
2 2j 0
3 + bi 3j 0
0
m mj
2500 0,15
17000 + Δb -0,5 0
1
500 -0, 05
Max (z) = 50 x1 + 10 x2
Sa
20 x1 + 10 x2 = 40
10 x1 + 10 x2 ≤ 60
20 x1 + 10 x2 ≥ 50
x1 , x2 ≥0
Cj 50 10 0 -M -M 0
CB Base VLD x1 x2 S2 A1 A3 S3
50 x1 1 1 0 0 2/30 -1/30 1/30
0 S2 30 0 0 1 -1/3 -1/3 1/3
10 x2 2 0 1 0 -1/30 2/30 -2/30
zj 70 50 10 0 3 -1 1
cj - zj 0 0 0 -M -M -1
Tabla 5
1 2 / 30 0
30 + Δb -1 / 3 0
1
2 -1 / 30 0
2
1 + Δb1 0
30
-1
30 + Δb1 0
3
-1
2 + Δb1 0
30
Despejando nos queda que b1 ≥ -15 y b1 ≤ 60, es decir que b1 puede
disminuir hasta en 15 e incrementarse hasta en 60 o lo que es lo mismo
el valor de b1 puede estar entre [25, 100].
1 -1 / 30 0
30 + Δb -1 / 3 0
3
2 2 / 30 0
1 1 / 30 0
30 - Δb 1 / 3 0
3
2 -2 / 30 0
1
1 - Δb3 0
30
1
30 - Δb3 0
3
2
2 + Δb3 0
30
En ambos casos, despejando nos queda que b3 ≥ -30 y b3 ≤ 30,
es decir que b3 puede disminuir hasta en 30 e incrementarse hasta
en 30 o lo que es lo mismo el valor de b 3 puede estar entre [20, 80].
1 2 / 30 5 / 3
30 + 10 -1 / 3 = 80 / 3
2 -1 / 30 5 / 3
1 2 / 30 1 / 3
30 -10 -1 / 3 = 100 / 3
2 -1 / 30 7 / 3
ACTIVIDADES DE AUTOEXAMEN
ACTIVIDAD 1
RESPONDA LAS SIGUIENTES PREGUNTAS
ACTIVIDAD 2
EXPLIQUE SI LAS SIGUIENTES AFIRMACIONES SON VERDADERAS O FALSAS:
ACTIVIDAD 3
Explique:
a) ¿Qué representa el costo reducido (reduced cost) que aparece en los
informes de computadora sobre la solución de un PL?
b) ¿Cuál es la utilidad de esta información?
c) ¿En qué parte de la tabla simplex se encuentra?
ACTIVIDAD 4
En base al problema de la SuperMovil SA y su tabla óptima de simplex:
ACTIVIDAD 5
En la empresa Amarras SA Juan es gerente de producción y está
tratando de decidir cuantos ganchos para trailer debe hacer para usar
un metal de desperdicio. Tiene tres tipos de metal y puede hacer
cualquiera de tres tipos de ganchos. En la tabla siguiente se
proporcionan los datos necesarios:
max 13 G1 + 16 G2 + 14 G3
sa
4 G1 + 5 G2 + 6 G3 <= 950 Hierro acanalado
6 G1 + 3 G2 + 5 G3 <= 800 Hierro plano
4 G1 + 8 G2 + 6 G3 <= 1150 Hierro redondo
G3 >= 40 Producción mínima
Celdas de variables
Nombre Valor original Valor final
Variables G1 0 57,5
Variables G2 0 85
Variables G3 0 40
Restricciones
Nombre Valor de la celda Estado Holgura
H.Acanaldo 895 No vinculante 55
H. Plano 800 Vinculante 0
H. Redondo 1150 Vinculante 0
Prod. Mínima 40 Vinculante 0
DUALIDAD Y SENSIBILIDAD 141
Celdas de variables
Final Reducido Objetivo Permisible Permisible
Nombre Valor Coste Coeficiente Aumentar Reducir
Variables G1 57,5 0 13 19 1,7273
Variables G2 85 0 16 10 2,375
Variables G3 40 0 14 1,0556 1E+30
Restricciones
Final Sombra Restricción Permisible Permisible
Nombre Valor Precio Lado derecho Aumentar Reducir
H.Acanaldo 895 0 950 1E+30 55
H. Plano 800 1,1111 800 165 258,75
H. Redondo 1150 1,5833 1150 110 510
Prod. Mínima 40 -1,0556 40 41,25 40
ACTIVIDAD 6
La siguiente tabla corresponde a un PL de maximización canónico:
Cj 40 60 50 0 0 0
CB Base VLD x1 x2 x3 S1 S2 S3
600 1/2 0 1 0 0 1/2
200 9/4 1 0 1/2 0 -1/4
200 -5/4 0 0 -1/2 1 1/4
Zj
Cj - Z j
ACTIVIDAD 7
Caso: Fábrica de bolsos y carteras “Sureñas”
Rodrigo es un pequeño empresario que se dedica a la fabricación de
carteras y bolsos femeninos.
En este momento está analizando el lanzamiento de dos nuevos
modelos de bolsos. En su confección utiliza cuero ecológico, herrajes,
cierres, hilo de seda reforzado y una tela especial con diseños originales.
Cj 0 0 0 0
Cj Básicas VLD X1 X2 X3 S1 S2 S3 S4
24 1 0 0 0 0 0 -1
428 0 0 0 0 1 -0.625 -0.125
23 0 0 1 -1 0 0.075 -0.125
160 0 1 0 0 0 0.125 0.625
Zj
Cj - Zj