PLC4
PLC4
ANÁLISIS DE SENSIBILIDAD
4.1 INTRODUCCIÓN
Toda solución a un problema de toma de decisiones se basa en determinados parámetros que se
presumen como fijos. En un ambiente real los datos utilizados en el modelo a menudo tienen que
estimarse y por lo tanto están sujetos a inexactitudes, o bien pueden cambiar con el tiempo
debido a la dinámica del negocio.
Una solución óptima sólo es óptima respecto al modelo específico que se está usando para
representar al problema real y, tal solución, se transforma en una guía confiable para la acción
únicamente después de que se ha verificado que se comporta con igual propiedad para otras
representaciones razonables del problema.
El análisis de sensibilidad o análisis post-óptimo es un conjunto de actividades que sirven
para estudiar y determinar qué tan sensible es la solución a los cambios en las hipótesis. Tales
cambios comprenden:
Cambios en los coeficientes de la función objetivo.
Cambios en los términos independientes de las restricciones.
Cambios en los coeficientes tecnológicos.
Agregado de una nueva variable.
Agregado de una nueva restricción.
Es importante destacar que el análisis de sensibilidad supone que los coeficientes varían
sólo uno a la vez manteniendo los restantes parámetros tal cual se plantearon en la forma
original.
73
Capítulo 4
La resolución gráfica de la Figura 4-1 nos muestra que mientras la pendiente del funcional
oscile entre la pendiente de la restricción debida a las hora s de mano de obra y la correspondiente
a la demanda máxima de cántaros, el plan diar io de producción se mantendrá óptimo.
Ma
te r i
n
aP
r im
a
Luego,
pendiente del funcional c1 /c 2 4/9
1 c
pendiente de la restricción de mano de obra 1 / 2 1 0 0 c1 22,5
2 45
pendiente de la restricción de demanda 0
por lo tanto, mientras la contribución marginal de las vasijas varíe entre $0 y $22,5 el plan de
producción actual seguirá siendo óptimo.Lógicamente el valor de z* dependerá del valor que
asuma c1.
Un análisis similar para la contribución marginal de los cántaros nos permitirá concluir que
ésta puede oscilar entre $40 e infinito, sin alterar el punto de óptimo.
Veamos ahora cómo generalizar este estudio cuando se dispone de la solución óptima de un
modelo de programación lineal resuelto por el método Simplex. En el Ejemplo 1-1 dicha
solución está dada por la tabla siguiente.
Tabla de Óptimo
20 45 0 0 0 xi/yij
ci Ai A1 A2 A3 A4 A5 xi (yij>0)
20 A1 1 0 1 0 –2 10
0 A4 0 0 –3 1 9/2 45/2
45 A2 0 1 0 0 1 15
zj 20 45 20 0 5
z = 875
cj – zj 0 0 –20 0 –5
74 Norma Torrent
Análisis de Sensibilidad
c j c i y ij c k y kj 0 j I N c k y kj c j z jj I N
iI
B
si y kj 0, c k (c j z j ) / y kj si y kj 0, c k (c j z j ) / y kj
por lo tanto,
cj zj cj zj
max y kj 0 c k min y kj 0
y kj y kj
20 5
c1 0 c1 22,5
1 2
coincidente con el análisis gráfico efectuado. Análogamente
cj zj cj zj
max y 2 j 0 c 2 min y 2 j 0
y2j y2j
5
c 2 40 c 2
1
Notemos que en el ejemplo presentado las variables no básicas en la solución óptima son
variables de holgura, no teniendo por tanto mayor sentido estudiar las variaciones en sus
coeficientes económicos.
75
Capítulo 4
óptima o no básica. Excepto en algunos casos especiales los citados problemas hacen que sea
dificultoso determinar sistemáticamente el efecto de tales cambios sobre el óptimo, s iendo
recomendable en estas circunstancias volver a resolver el problema.
Nos limitaremos entonces al estudio de variaciones en los coeficientes tecnológicos no
básicos.
Para una solución óptima c j z j 0 Aj (caso de maximización). Frente un cambio en los
coeficientes tecnológicos de la variable no básica xk, bastará con analizar el correspondiente
c k z k . Sea Ak' el nuevo vector de coeficientes, tendremos
77
Capítulo 4
Por supuesto que la tabla anterior ya no es una tabla Simplex dado que A1 , A4 y A2 son
vectores básicos y sus correspondientes coeficientes en la restricción adicional deben s er nulos.
Realizando operaciones elementales sobre el cuerpo central de la tabla obtenemos :
20 45 0 0 0 0 xi/yij
ci Ai A1 A2 A3 A4 A5 A6 xi (yij>0)
20 A1 1 0 1 0 –2 0 10
0 A4 0 0 –3 1 9/2 0 45/2
45 A2 0 1 0 0 1 0 15
0 A6 0 0 –1 0 2 1 –2
zj 20 45 20 0 5 0
cj – zj 0 0 –20 0 –5 0
Como puede apreciarse se verifica la condición de optimización pero la solución actual es
no factible. En estos casos para obtener el nuevo óptimo se emplea la tabla precedente como
tabla de inicio para la aplicación del método Dual Simplex que estudiaremos en el próximo
capítulo.
x 1 2 x 2 40
x 1+
2x 2
2 x 1 x 2 50 20 ÓPTIMO
=z 0
x1 , x 2 0
x j : unidades semanales del modelo j; 0 25 40 x1
j = 1 (Oxford), 2 (clásico)
78 Norma Torrent
El gráfico indica que la solución se encuentra en la intersección de las rectas
x1 2x 2 40 y 2x1 x 2 50, resolviendo el sistema tenemos x1* 20 y x 2* 10, por lo tanto
z * 80. Llevando el programa a su forma estándar resulta inmediato que x3* 0 y x4* 0.
Para obtener la correspondiente tabla de óptimo calculamos la inversa de la base.
1
1 1 2 1 / 3 2 / 3
A1 , A 2
2 1 2 / 3 1 / 3
Luego,
3 2 0 0
ci Ai A1 A2 A3 A4 xi
3 A1 1 0 –1/3 2/3 20
2 A2 0 1 2/3 –1/3 10
zj 3 2 1/3 4/3
z = 80
cj – zj 0 0 –1/3 –4/3
b) Variación en c1
c1 corresponde a una variable básica, entonces
cj zj cj zj
max y1 j 0 c1 min y1 j 0
y1 j y1 j
4/3 1/ 3
c1 2 c1 1 1 c1 4
2/3 1/ 3
Observemos que hubiésemos obtenido idéntico resultado basándonos en el análisis gráfico.
En efecto, mientras la pendiente del funcional oscile entre la pendiente de la restricción debida a
las horas disponibles de Thelma y la correspondiente a las horas disponibles de Luisa, el plan
semanal de producción se mantendrá óptimo. Es decir, 2 c1 / 2 1 / 2 1 c1 4 . El
precio de venta de los jeans Oxford podrá oscilar entonces entre $23 y $26 y el valor de z*
dependerá del valor que asuma c1.
c) Variación en b1
x *i actual x *i actual
max
ri1 0 Δb 1 min rik 0
ri1 ri1
–1
donde ri1 el i-ésimo elemento en la columna 1de B .
10 20
b1 15 b1 60 25 b1 100
2/3 1/ 3
Para 25 b1 100 la base actual se mantiene óptima. Si b1 = 60 la nueva solución es
1 / 3 2 / 3 40 20 x 1 (1 / 3)20 40 / 3
x *B B 1b Nueva
2 / 3 1 / 3 50 x 2 (2 / 3)20 70 / 3
x1* 13,33, x 2* 23,33, x3* 0, x4* 0 y z * 86,66.
A igual conclusión hubiéramos arribado analizando el gráfico. Si b1 cambia, la nueva
solución óptima seguirá teniendo como variables básicas a x1 y x2 sólo si se encuentra en la
intersección de las rectas correspondientes a la carga horaria de ambas costureras (siempre que
esta intersección pertenezca al conjunto de soluciones posibles). Así, b1 puede incrementarse
hasta que esta línea de restricción pase por el punto (0; 50), y puede disminuir hasta que la
misma toque el punto (25; 0). De esta forma, para 25 b1 100 la base (A1, A2) se mantendrá
óptima y solución óptima dependerá del valor asumido por b1 (en este caso b1 = 60).
Ahora bien, dado que el nuevo plan de producción debe estar en términos valores enteros,
* * * *
redondeando obtenemos x1 14, x 2 22, x3 2(horas ociosas de Luisa), x4 0 y z * 86.
79
Capítulo 4
* * * *
Otras posibilidades hubiesen sido (x1 14, x 2 23), (x1 13, x 2 24) sin embargo tales
puntos no satisfacen las restricciones (no pertenecen a la región factible). 18
d) Agregado de una nueva variable
Para que el nuevo modelo tenga posibilidad de integrar la base óptima (c Nuevo z Nuevo )
debe ser no negativo. Luego
1 / 3 2 / 3 2 2 / 3 2 / 3
y iNuevo B 1A Nuevo z Nuevo c i y iNuevo 3, 2 10 / 3
2 / 3 1 / 3 2 2 / 3 2 / 3
c Nuevo z Nuevo 0 c Nuevo 10 / 3
c Nuevo Pr ecio Venta Costo Materia Prima Costo Costureras Pr ecio Venta 25
en consecuencia el precio de venta del nuevo modelo deberá ser mayor o igual a $28,33.
El resultado precedente nos dice que los valores implícitos de cada restricción son los
valores zj de la tabla de óptimo correspondientes a las columnas de la base inicial identidad.19
Para entender la importancia y características de este conce pto trabajaremos sobre algunos
ejemplos.
Analicemos nuevamente el Ejemplo 1-1.
Max z 20x 1 45x 2 (maximizar contribución marginal)
s. a
x 1 2 x 2 40 (restricción debida a la disponibilidad de horas de mano de obra)
3x 1 1,5x 2 75 (restricción debida a la disponibilidad de materia prima)
x 2 15 (restricción debida a la demanda máxima de cántaros)
x1 ; x 2 0
Tabla de Óptimo
20 45 0 0 0 xi/yij
ci Ai A1 A2 A3 A4 A5 xi (yij>0)
20 A1 1 0 1 0 –2 10
0 A4 0 0 –3 1 9/2 45/2
45 A2 0 1 0 0 1 15
zj 20 45 20 0 5
z = 875
cj – zj 0 0 –20 0 –5
18
Véase apartado 1.9 en Capítulo 1.
19
Véase Cómo identificar B-1 en una tabla Simplex, apartado 3.3 en Capítulo 3.
80 Norma Torrent
Análisis de Sensibilidad
incremente en $20.20
Observemos además que en el óptimo el costo reducido de x3 es
z *
(c 3 z 3 )* 0 z *3 z *3 20 por lo que si incrementamos x3 en una unidad z*
x3
disminuirá en $20. Asimismo, pretender incrementar x3 equivale a restarle capacidad al recurso
b1 ya que, en la situación actual, el hecho que x3 asuma un valor positivo significa reducir a
40 x3 las horas diarias de mano de obra destinadas a las actividades x1 y x 2 . De esta forma
cualquier incremento en x3 equivale a un decremento en b1 . Luego
z * z *
u1* z *3 20 .
b1 x3
Cuando el término independiente de una restricción de representa la cantidad disponible
de un determinado recurso, el valor implícito se denomina también precio sombra, costo
marginal o costo de oportunidad del recurso asociado a la restricción (ya que existe la
oportunidad de mejorar el valor de z* al aumentar la disponibilidad del recurso).
Los costos de oportunidad se utilizan frecuentemente para saber cuál es la cantidad máxima
que estaríamos dispuestos a pagar por una unidad adicional de un determinado recurso.
Concretamente, en nuestro ejemplo, ¿cuánto estaríamos dispuestos a pagar por una hora
adicional de mano de obra?
Dado que el costo marginal de dicho recurso es $20, podríamos pagar hasta $20, por encima
de su valor actual. Lógicamente en el límite el aumento de z será nulo.
Pero cuidado, si deseamos incrementar la cantidad de horas de mano de obra lo primero que
debemos hacer es realizar el análisis de sensibilidad para saber si dicho incremento pertenece al
intervalo de variación permitida. Si esto ocurre el precio sombra mantendrá su valor, en caso
contrario puede cambiar.
Mediante un razonamiento análogo podremos concluir que un incremento en b2 no variará el
valor actual de z (la materia prima es un recurso ocioso, x*4 22,5 ).
En cuanto a la tercera restricción, u*3 5 nos dice que un incremento de una unidad en la
demanda máxima de cántaros producirá una mejora de $5 en el funcional. En este caso, a
diferencia de los anteriores, la restricción no se debe a los recursos disponibles sino a un
requisito o requerimiento a satisfacer. El valor implícito, en vez de costo marginal (o de
oportunidad) asociado a la restricción, se denomina entonces valor marginal asociado al
requerimiento.
Vayamos ahora el Ejemplo 3-2.
20
En este caso el valor implícito puede interpretarse también como el beneficio que deja de percibir la organización
por no disponer de una unidad más del recurso correspondiente.
81
Capítulo 4
Para este problema u1* z 3 0 revela que al incrementar en una unidad el término
independiente de la primera restricción, el valor del funcional no experimentará variación alguna
mientras que, si adicionamos una unidad a b2 el valor de z disminuirá en 20 unidades
( u *2 z 4 20 ).
21
Véase Efecto espejo, apartado 3-5 en Capítulo 3.
82 Norma Torrent
Análisis de Sensibilidad
83
Capítulo 4
Problema I: Max z = x1 + x2
Problema II: Max z = 2x1 + x2
Problema III: Max z = x1 + 3x2
Los tres programas tienen la misma región factible dada por:
x1 x 2 4
x1 1
x1 , x 2 0
Calcularemos el óptimo de cada uno de ellos y examinaremos el efecto en z* al incrementar
en una unidad el término independiente de la segunda restricción.
z=2
z=2
z=
z=
x 1+
x 1+
x
x1
1+
+x
x 2=
x 2=
x2
2=
=z
z0
z0
z0
Finalmente cabe acotar que si una solución óptima es degenerada una variación en bi puede
causar que la base actual no se factible ( z * / bi puede no ser derivable en algunos puntos) en
cuyo caso los valores implícitos dejan de ser válidos.
84 Norma Torrent
Análisis de Sensibilidad
Tales reglas nos permitirán saber qué sucede si más de un coeficiente económico, o un
término independiente, varían dentro de sus respectivos intervalos de sensibilidad.22
puesto que 140% > 100% no podemos asegurar que x* siga siendo el óptimo.
Observemos que si cambiara un solo coeficiente de la función objetivo x* seguirá siendo el
óptimo si rk 1 (o equivalentemente rk 100% ). La regla 100% es una generalización de esta
idea.
22
El lector interesado en la demostración de las regles 100% puede consultar Bradley, S., A. Hax y T. Magnanti,
Applied Mathematical Programming, MA: Addison-Wesley Publishing Company, 1977.
85
Capítulo 4
Desde el punto de vista geométrico, el conjunto convexo (en este caso no acotado) de la
figura es sólo un subconjunto de una región de sensibilidad más amplia para los cambios
posibles en ambos coeficientes económicos (cambios fuera de los intervalos de sensibilidad
calculados en 4-2 y 4-3). Por lo tanto, permanecer dentro de la región sombreada en la figura es
sólo una condición suficiente (no necesaria) para que x* siga siendo el óptimo.
Pueden obtenerse resultados similares para cambios simultáneos de los términos
independientes.
86 Norma Torrent
Análisis de Sensibilidad
necesaria), ya que si se viola, los precios sombra pueden cambiar o seguir iguales.
Veamos el efecto que provoca en el Ejemplo 1-1 disponer, simultáneamente, de 46 horas de
mano de obra y 85kg de arcilla, por día.
Los resultados de la sección 4.3 nos muestran que para 10 Δb1 7,5 (manteniendo fijo
b2 ) y para 22,5 Δb2 (manteniendo fijo b1 ), la base actual seguirá siendo óptima. De
esta forma,
r1 6/7,5 y r2 10 / r1 r2 0,8 0 0,8
80% < 100%, en consecuencia los cambios mantienen la base óptima. Dado que el costo
marginal de b1 es igual $20 y el de b2 igual a $0, el incremento en z estará dado por 620 +
100 =120 es decir el nuevo z * será igual a 875 + 120 = 925. Los valores de las variables
1
básicas pueden calcularse mediante B* bnuevo .
Por último es conveniente aclarar que si bien las reglas del 100% pueden extenderse a
programas lineales de cualquier magnitud, a medida que aumenta la dimensión del problema la
región de sensibilidad se reduce y los resultados obtenidos no son tan útiles para la gestión. En
tales casos el análisis paramétrico que veremos a continuación es la herramienta adecuada.
87
Capítulo 4
Expresando a c 2 en función de λ ,
20 45(1+) 0 0 0 xi/yij
ci Ai A1 A2 A3 A4 A5 xi (yij>0)
20 A1 1 0 1 0 –2 10
0 A4 0 0 –3 1 9/2 45/2
45(1+) A2 0 1 0 0 1 15
zj 20 45(1+) 20 0 5+45
z = 875+675
cj – zj 0 0 –20 0 –5–45
23
Limitaremos nuestro análisis al estudio de variaciones lineales. Si las variaciones son no lineales seguimos
estando frente a un problema de programación lineal pero el análisis se torna bastante complicado.
88 Norma Torrent
Análisis de Sensibilidad
89
Capítulo 4
si m j 0, h j / m j si m j 0, h j / m j
- hj
Para λ max m j 0 la nueva solución se obtendrá introduciendo el vector para el cual
mj
- hj - hj
resulte c j z j max m j 0 . En forma similar para λ min m j 0 ingresará a la base el vector
mj mj
- hj
cuyo c j z j min m j 0 .
mj
90 Norma Torrent
Análisis de Sensibilidad
91
Capítulo 4
que mientras 30 b1 47,5 el costo marginal de cada hora de mano de obra es de $20.
Lógicamente este costo marginal refleja cómo cambia z * en función de b1 .
Para entre –40 y –10 ( 0 b1 30 ), x*B (A2 , A4 , A5 ) 1 (40 , 75, 15)T . En este
caso.
1/2 0 0
A 2 , A 4 , A 5 1
3/4 1 0
1 / 2 0 1
luego, x*2 20 0,5 , x*4 45 0,75 , x*5 5 0,5 y z * 900 22,5 .
Por ejemplo, para 20 , x* (0; 10; 0; 60; 5)T y z * 450. El costo marginal de
b1 es ahora $22,5.
Para ( b1 47,5 ) un análisis similar a los anteriores nos indicará que
x1* 17,5 , x*2 15 y z * 1.025 , en tanto que x*3 7,5 . Notemos que al ser x*3
una variable básica el costo marginal de b1 es nulo.
El siguiente cuadro resume los resultados obtenidos.
92 Norma Torrent
Análisis de Sensibilidad
-10 7,5
b1 30 47,5 z
x1 0 10+ 17,5 1.025
x2 20+0,5 15 15
x3 0 0 - 675
x4 45-0,75 22,5-3 0
x5 -5-0,5 0 0
z 900+22,5 875+20 1.025 0 30 47,5 b1
93