Presentacion Clases V 2
Presentacion Clases V 2
Presentacion Clases V 2
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
I.
Introduccin
I.1. Introduccin
Antecedentes:
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
Durante la Segunda Guerra Mundial, el mando britnico consult a cientficos y tcnicos sobre distintas cuestiones militares: Despliegue de radares. Direccin de operaciones antisubmarinas, de minas, bombardeos y traslado de tropas. El resultado se llam Investigacin de Operaciones Militares, y ms tarde Investigacin Operativa (IO) El MIT contribuy a su puesta en marcha El profesor Morse (MIT) fue pionero en los EE.UU. Fund el Centro OR del MIT y colabor en la fundacin de ORSA.
3
I.
Introduccin
MAESTRIA PYGE
Qu son las ciencias de la gestin (investigacin operativa)? Hoy en da: las ciencias de la gestin y la investigacin operativa suponen el empleo de modelos matemticos para proporcionar pautas que permitan a los gestores tomar decisiones efectivas partiendo de la informacin de la que disponen, o para hallar el modo de ampliar sta, en caso de que sea insuficiente para llegar a la decisin adecuada.
I.
Introduccin
La investigacin operativa en la historia 1947 Proyecto Scoop (Scientific Computation of Optimum Programs), en el que George Dantzig y otros cientficos desarrollan el mtodo simplex de programacin lineal. Dcada de 1950 Aos muy interesantes: progresos matemticos, teora de colas, programacin matemtica. Comparar con: Inteligencia artificial (I.A.) en los 60. Dcada de 1960 Crece el inters: ms progreso, grandes proyectos. Comparar con: Inteligencia artificial en los 80.
MAESTRIA PYGE
I.
Introduccin
Dcada de 1970 poca de desilusin y estancamiento. NPcompleto. Expectativas ms realistas. Dcada de 1980 Gran expansin del uso de computadores personales. Acceso cada vez ms fcil a datos. Se extiende la disposicin de los directivos al empleo de modelos. Dcada de 1990 Uso creciente de sistemas de I.O. Nuevos avances de la tecnologa de I.O; p.ej: ampliaciones de optimizacin y simulacin a hojas de clculo, lenguajes de modelacin, optimizacin a gran escala. Mayor interconexin entre la I.A. y la I.O
6
MAESTRIA PYGE
I.
Introduccin
La investigacin operativa en el ao 2000 CIENTOS de oportunidades para el campo de la I.O. Datos, datos y ms datos Datos de e-business (click stream, compras, otros datos de transacciones, correo electrnico, etc.) El proyecto del genoma humano y su desarrollo Mayor automatizacin en la toma de decisiones Necesidad de mayor coordinacin para la utilizacin de recursos (gestin de la cadena de suministro)
MAESTRIA PYGE
I.
Introduccin
Definicin
MAESTRIA PYGE
Aplicacin del mtodo cientfico por un grupo multidisciplinario personas a la resolucin de un problema. Objetivo El principal objetivo de esta rea de conocimientos consiste en formular y resolver diversos problemas orientados a la toma de decisiones, mediante mtodos cientficos, que optimizan el funcionamiento del proceso analizado, generalmente bajo condiciones que implican la utilizacin de recursos escasos.
I.
Introduccin
Naturaleza
La naturaleza de los problemas abordados pueden ser:
MAESTRIA PYGE
Mtodos
determinsticos: Programacin lineal, programacin entera, transporte, teora de la localizacin o redes, programacin multicriterio, teora de inventarios, etc. (Modelos de Prog. Matemtica)
hbridos: probabilsticos.
Conjugan
mtodos
determinsticos
I.
Introduccin
MAESTRIA PYGE
Uso
la
formulacin
10
I.
Introduccin
MAESTRIA PYGE
Definir el Problema Especificacin de Objetivos Reunir Datos Estimar Valores de Parametros Representacin Idealizada del Problema Estimar el Grado de Acercamiento del Modelo a la Realidad Escoger el modelo que mejor se adapta a los objetivos
2.
Observar el Sistema
3. 4. 5. 6. 7.
Presentar los Resultados y Conclusiones del Estudio a la Organizacin Implantar y Evaluar las Recomendaciones
11
I.
Introduccin
Formular el Problema
MAESTRIA PYGE
Observar el Sistema
Ing. Rodrigo Semprtegui lvarez
I.
Introduccin
MAESTRIA PYGE
Supongamos que se dispone de determinadas piezas para la elaboracin de dos productos finales. Se dispone de 8 piezas pequeas y 6 piezas grandes, que son utilizadas para elaborar sillas (usando 2 piezas pequeas y 1 pieza grande) y mesas (usando 2 piezas de cada tipo). Interesa decidir cuntas sillas y mesas fabricar de modo de obtener la mxima utilidad, dado un beneficio neto de U$ 15 por cada silla y de U$20 por cada mesa fabricada.
13
I.
Introduccin
MAESTRIA PYGE
Posibles soluciones factibles a considerar, esto es soluciones que respetan las restricciones del nmero de piezas disponibles, son por ejemplo, fabricar:
4 sillas, que reportan una utilidad de U$60 1 sillas y 2 mesas , utilidad de U$55
14
I.
Introduccin
Componentes de un modelo matemtico: i) Las variables de decisin, que consiste en definir cules son las decisiones que se debe tomar. En el ejemplo, x: nmero de sillas elaboradas. y: nmero de mesas elaboradas.
MAESTRIA PYGE
Qu puedes decidir?
Ej: cuanto producir; cuanto invertir, y en qu, son variables de decisin
ii) La funcin objetivo del problema, que permita tener el mejor criterio para decidir entre todas las soluciones factibles. En el ejemplo, maximizar la utilidad dada por: z = f(x,y) = 15x + 20y
15
I.
Introduccin
MAESTRIA PYGE
iii) Restricciones del problema, que consiste en definir un conjunto de ecuaciones e inecuaciones que restringen los valores de las variables de decisin a aquellos considerados como factibles. En el ejemplo, respetar la disponibilidad de piezas para la fabricacin de sillas y mesas: Piezas pequeas: 2x + 2y 8 Piezas grandes : x + 2y 6 Tambin se impone restricciones de no negatividad: x,y 0
16
I.
Introduccin
En resumen:
MAESTRIA PYGE
Max sa:
El ejemplo corresponde a un modelo de Programacin Lineal. Si adems restringimos los valores de x e y a nmeros enteros, tendramos un modelo de Programacin Entera. Por otra parte, si hubiese retornos crecientes a escala, deberamos emplear una funcin objetivo no lineal como f(x,y) = cxa + dyb con a,b >1, y tendramos un modelo de Programacin No Lineal.
17
I.
Introduccin
MAESTRIA PYGE
1. Introduccin a la Investigacin de Operaciones, F.S. Hillier y G.J. Lieberman, McGraw Hill, Sexta Edicin, 1997. 2. Investigacin de Operaciones, una introduccin, H.A. Taha, Prentice Hall, Mxico, Sexta Edicin, 1998. 3. Introduction to Management Science, F. Hillier, M. Hillier and G.J. Lieberman. Irwin McGraw-Hill, 1999. 4. Model Operations Research: A practical Introduction. M.W. Carter and C.C.Price. CRC Press, 2000. 5. Practical Management Science: Spreadsheet Modeling and Applications, Winston, W.L., Albright S.C. y Broadie M., International Thomson Publishing Company, 1997.
18
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
Temario:
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
II.1. Introduccin II.2. Definiciones II.3. Suposiciones de la PL II.4. Ejemplos de modelamiento. II.5. Resolucin grfica de problemas. II.6. Anlisis de Sensibilidad. II.7. El Mtodo Simplex. II.8. Dualidad en Programacin Lineal. II.9. Anlisis de Sensibilidad o Post-Optimal
20
MAESTRIA PYGE
matemtica, en donde todas las funciones que participan en el modelo son lineales.
Utilizacin:
Planificacin Gestin de recursos humanos y materiales Transporte Planificacin financiera Organizacin de la produccin. Una extensa gama de problemas que aparecen en las reas de tipo industrial, econmico, administrativo, militar, etc.
21
MAESTRIA PYGE
lineal de x1, x2, x3,. xn, si y solo si para algn conjunto de constantes c1, c2, c3,. cn, existe una funcin f(x1, x2, x3,. xn) = c1 x1+ c2 x2, c3 x3,. cn xn
f(x1, x2, x3,. xn) b es lineal, si y solo si la funcin f(x1, x2, x3,. xn) es lineal y b es cualquier nmero real.
optimizacin, para el cual:
Se busca maximizar o minimizar una funcin lineal de variables de
22
MAESTRIA PYGE
Solucin ptima.- Para un problema de maximizacin en PL la solucin ptima es el punto o conjunto de puntos de la regin factible con mayor valor de la funcin objetivo. Para un problema de minimizacin en PL la solucin ptima es el punto o conjunto de puntos de la regin factible con menor valor de la funcin objetivo.
23
MAESTRIA PYGE
que ser una funcin lineal de las variables de decisin, tiene dos implicaciones:
La contribucin de cada variable de decisin a la funcin objetivo es
proporcional al valor de la variable de decisin. La contribucin a la funcin objetivo por parte de cualquier variable de decisin es independiente de los valores de las otras variables de decisin.
que ser una desigualdad o igualdad lineal de las variables de decisin, tiene tambin dos implicaciones:
La contribucin de cada variable al lado izquierdo de cada
restriccin es proporcional al valor de la variable de decisin. La contribucin de cada variable al lado izquierdo de cada restriccin independiente de los valores de las otras variables de decisin.
24
MAESTRIA PYGE
Suposicin de Certidumbre. La suposicin de certidumbre significa que tiene que conocerse con certidumbre cada parmetro (coeficiente de la funcin objetivo, coeficientes de las variables de las restricciones, lado izquierdo de las restricciones).
25
i) Problema de Transporte (ENERGIA). El problema consiste en decidir cuntas unidades trasladar desde ciertos puntos de origen (plantas, ciudades, etc.) a ciertos puntos de destino (centros de distribucin, ciudades, etc..) de modo de minimizar los costos de transporte, dada la oferta y demanda en dichos puntos. Se suponen conocidos los costos unitarios de transporte, los requerimientos de demanda y la oferta disponible.
26
MAESTRIA PYGE
Por ejemplo, suponga que una empresa posee dos plantas que elaboran un determinado producto en cantidades de 250 y 450 unidades diarias, respectivamente. Dichas unidades deben ser trasladadas a tres centros de distribucin con demandas diarias de 200, 200 y 250 unidades, respectivamente. Los costos de transporte (en $/unidad) son:
27
C.D.1
Orgenes
Destinos
28
MAESTRIA PYGE
Funcin Objetivo: Minimizar el costo total de transporte dado por la funcin: 21x11+25x12+15x13+28x21+13x22+19x23 Restricciones del problema: 1) No Negatividad: xij 0
2) Demanda: CD1 : x11 +x21 CD2 : x12 +x22 CD3 : x13 + x23
MAESTRIA PYGE
3) Oferta : P1 : x11 + x12 + x13 250 P2 : x21 + x22 + x23 450 Las variables de decisin deben aceptar soluciones como nmeros reales para tener un modelo de P.L.
30
MAESTRIA PYGE
ii) Problema de la dieta: este consiste en determinar una dieta de manera eficiente, a partir de un conjunto dado de alimentos, de modo de satisfacer ciertos requerimientos nutricionales. Supongamos que se tiene la siguiente informacin:
Leche Legumbre Naranjas Requerimientos (galon) (1 porcin) (unidad) Nutricionales Niacina Tianina Vitamina C Costo 3,2 1,12 32 2 4,9 1,3 0 0,2 0,8 0,19 93 0,25 13 15 45
31
MAESTRIA PYGE
Funcin Objetivo:
Minimizar el costo total de la dieta, dado por: 2 x1 + 0.2 x2 + 0.25 x3 Restricciones del problema:
MAESTRIA PYGE
Este consiste en hallar una poltica ptima de produccin para satisfacer demandas fluctuantes en el tiempo, de modo de minimizar costos de produccin e inventario, considerando la disponibilidad de diversos recursos escasos. Supongamos que una fabrica puede elaborar hasta 150 unidades en cada uno de los 4 perodos en que se ha subdividido el horizonte de planificacin y se tiene adicionalmente la siguiente informacin:
33
MAESTRIA PYGE
1 2 3 4
6 4 8 9
2 1 2.5 3
Supuestos adicionales: 1) Existe un inventario inicial de 15 unidades. 2) No se acepta demanda pendiente o faltante (es decir, se debe satisfacer toda la demanda del periodo).
34
Variables de decisin:
MAESTRIA PYGE
It : nmero de unidades de inventario al final del periodo t. Funcin objetivo: Consiste en minimizar los costos de produccin y el costo de mantenimiento de inventario. 6x1+ 4x2 + 8x3 + 9x4 + 2I1 + I2 + 2.5I3 + 3I4
Notar que en el ptimo I4 ser 0, as que incluso podramos no incluirla, pero de todos modos la consideramos.
35
MAESTRIA PYGE
1) Restricciones de cotas, que reflejan la capacidad de produccin. xt 150 t Perodo 2) Restricciones de no negatividad xt 0 t Perodo
3) Restricciones de demanda
x1 + I0 I1 = 130 x2 + I1 I2 = 80 Periodo 1 Periodo 2 I0=15
x3 + I2 I3 = 125
x4 + I3 I4 = 195
Periodo 3
Periodo 4
36
MAESTRIA PYGE
Supongamos que un banco dispone de $250 millones para destinar a 4 tipo de crditos ofrecidos, los cuales tienen las siguientes, tasas de crdito:
Primer crdito corriente (PCC) Segundo crdito corriente (SCC) Crdito para el hogar Crdito personal
37
La asignacin de estos crditos, debe satisfacer la siguiente poltica utilizada por la institucin:
MAESTRIA PYGE
El monto asignado a los PCC, debe ser al menos, el 55% del monto asignado a los crditos corrientes, y al menos un 25% del total del dinero prestado.
El SCC, no puede exceder el 30% del total del dinero prestado, por polticas tributarias el inters recibido por el banco no debe exceder a un retorno del 14% sobre el capital prestado. Cunto asignar a cada tipo de crdito, de la manera ms eficiente, respetando la poltica del banco?
38
MAESTRIA PYGE
x4 : Monto asignado al crdito personal. Funcin Objetivo: Se propone maximizar los retornos recibidos en la asignacin, dados por: 0.12 x1 + 0.16 x2 + 0.16 x3 + 0.10 x4 Restricciones del problema: x1 0.55 ( x1 + x2 ) x1 0.25 ( x1 + x2 +x3 + x4 )
MAESTRIA PYGE
v) Problema de mezcla de productos: en este problema una refinera produce 4 tipos de gasolina (gas 1, gas 2, gas 3 y gas 4). Dos caractersticas importantes de cada gasolina son su nmero de performance (NP) y su presin de vapor (RVP), que estn dados por:
NP gas 1 107 RVP 5 Barriles diarios 3814
gas 2
gas 3 gas 4
93
87 108
8
4 21
2666
4016 1300
40
MAESTRIA PYGE
Estas gasolinas pueden ser vendidas directamente a un precio de $24,83 por barril o bien mezcladas para obtener gasolinas de aviacin (avgas A y avgas B). La calidad de estas dos ltimas junto con sus precios de venta son:
NP
RV
26,45 25,91
El NP y RVP de cada mezcla es un promedio de los respectivos NP y RVP de las gasolinas empleadas. Se desea obtener un plan de venta de las distintas gasolinas que maximice los retornos.
41
MAESTRIA PYGE
xjA: cant. de gas j usado en avgas A. xjB: cantidad de gas j usado en avgas B.
Funcin objetivo: Max 24,83 (x1 + x2 + x3 + x4) + 26,45xA + 25,91xB Restricciones: x1 + x1A + x1B = 3814
42
MAESTRIA PYGE
NP, avgas B:
RVP, avgas A:
5x1A 8x 2 A 4x 3 A 21x 4 A 7 xA
RVP, avgas B:
5 x 1B 8 x 2 B 4 x 3 B 21x 4 B 6 xB
43
MAESTRIA PYGE
MAESTRIA PYGE
yt : cantidad de MW expandidos en planta a gas al inicio del ao t, con t = 1, 2, ..., T. zt : cantidad total de MW disponible en plantas nuevas a petrleo al inicio del ao t. wt : cantidad total de MW disponible en plantas nuevas a gas al inicio del ao t.
45
p t x t gt yt
T t 1
MAESTRIA PYGE
c t z t w t dt z t xk
k 1 t t
w t yk
k 1 t
t 15 t 15 t 1...T
t 20 t 20
wt
k t 14
zt
k t 19
wt 0,30 ct zt w t x t , yt , z t , w t 0
46
MAESTRIA PYGE
Las reglas sindicales sealan que cada empleado debe trabajar por 5 das consecutivos y despus descansar 2 das. La oficina de correos quiere cumplir con sus requerimientos diarios y utilizar solamente empleados de tiempo completo. Formule un PL para minimizar los empleados de tiempo completo
47
MAESTRIA PYGE
Funcin Objetivo:
Ing. Rodrigo Semprtegui lvarez
Min z= x1 + x2 + x3 + x4 + x5 + x6 + x7 x1 x1 + x2 x1 + x2 + x3 x1 + x2 + x3 + x4 x1 + x2 + x3 + x4 + x5 x2 + x3 + x4 + x5 + x6 x3 + x4 + x5 + x6 + x7 + x4 + x5 + x6 + x7 + x5 + x6 + x7 + x6 + x7 + x7 17 13 15 19 14 16 11
Restricciones:
xi 0 (i=1,2,3,7)
48
MAESTRIA PYGE
49
MAESTRIA PYGE
50
Consideremos el siguiente problema a resolver grficamente: Max sa: z = 3x1 + 5x2 x1 2x2 3x1 + 2x2 x1 , x2 0 4 12 18
51
MAESTRIA PYGE
x2 9
6 4
x1
52
MAESTRIA PYGE
x2 9
6 4
x1
53
MAESTRIA PYGE
x2 9 x* 6 4
6
54
MAESTRIA PYGE
x2 9 x* 6 4
6
55
MAESTRIA PYGE
x2 9 x* 6 4
6
56
MAESTRIA PYGE
x2 9 x* 6 4 x*
6
57
MAESTRIA PYGE
En primer lugar, se debe obtener la regin de puntos factibles en el plano, obtenida por medio de la interseccin de todos los semi-espacios que determinan cada una de las inecuaciones presentes en las restricciones del problema. A continuacin, con el desplazamiento de las curvas de nivel de la funcin objetivo en la direccin de crecimiento de la funcin (que corresponde a la direccin del vector gradiente de la funcin, z(x1,x2) = (3,5)T, se obtiene la solucin ptima del problema en la interseccin de las rectas: 2x2 = 12 y 3x1+2x2 = 18 (restricciones activas). Esto es: x1 * = 2 x2 * = 6
z* = 3 x1* + 5 x2* = 36
58
MAESTRIA PYGE
Hemos tomado una situacin real y hemos construido sus equivalentes matemticos: MODELO MATEMTICO Durante la formulacin de los modelos matemticos, hemos considerado el mtodo cuantitativo que nos permitir resolver el modelo numricamente ALGORITMO El algoritmo es un conjunto de instrucciones que siguiendo de manera ordenada producen una solucin numrica
NUEVA DEFINICION
Ciencia para la representacin de problemas reales mediante modelos matemticos que junto con mtodos cuantitativos nos permiten obtener una solucin numrica a los mismos.
59
Notar que se pueden dar otras situaciones en la bsqueda de una solucin ptima para esta clase de problemas: 1) La solucin ptima exista pero haya ms de una. En el ejemplo, considerese la nueva funcin objetivo: z = 6x1+4x2. 2) El problema no tenga solucin, dada una regin de puntos factibles no - acotada. En el ejemplo, reemplace cada desigualdad por una . 3) El problema no tenga solucin, porque no existen puntos factibles. En el ejemplo, suponga que agregamos la restriccin: x1 5.
MAESTRIA PYGE
60
61
MAESTRIA PYGE
x 2y 6 x, y 0
6
62
63
MAESTRIA PYGE
x 2y 6 x, y 0
3 2
6
64
MAESTRIA PYGE
x 2y 6 x, y 0
6
65
MAESTRIA PYGE
x 2y 6 x, y 0
6
66
MAESTRIA PYGE
x 2y 6 x, y 0
6
67
MAESTRIA PYGE
x 2y 6 x, y 0
6
68
MAESTRIA PYGE
x 2y 6 x, y 0
6
69
MAESTRIA PYGE
x 2y 6 x, y 0
6
70
MAESTRIA PYGE
x 2y 6 x, y 0
6
71
Sea
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
z = c1x1+c2x2
c1 1 1 c2 2
72
MAESTRIA PYGE
Tambin podemos estudiar el intervalo de un slo coeficiente, dejando el resto de los parmetros fijos:
Ing. Rodrigo Semprtegui lvarez
c1 1 1 20 2
10 c1 20
15 1 1 c2 2
15 c 2 30
73
MAESTRIA PYGE
2) Cul es el intervalo de variacin de los coeficientes del lado derecho (trminos libres) de las restricciones, de modo que la actual solucin siga siendo la ptima? Estudiaremos por separado las variaciones de cada uno de los coeficientes del lado derecho de las restricciones, de modo preservar la geometra del problema, esto es, que se conserven las mismas restricciones activas de la solucin ptima inicial.
74
MAESTRIA PYGE
x 2y 6 x, y 0
8
75
MAESTRIA PYGE
x 2y 6 x, y 0
8
76
Primera restriccin.
MAESTRIA PYGE
La mayor variacin del coeficiente del lado derecho se alcanza en x=0 y y=4, de donde se obtiene: z(0,4) = 15 0 + 20 4 = 80 y b1* = 0 + 2 4 = 8 La menor variacin del coeficiente del lado derecho se alcanza en: x=4 ; y=0, de donde se obtiene: z(4,0) = 15 4 + 20 0 = 60 y b1 = 4 + 2 0 = 4 De aqu, se calcula el precio sombra P1, que indica la razn o tasa de cambio de la funcin objetivo con respecto al cambio en una unidad del lado derecho:
P1 z(0,4) z(4,0) 80 60 5 * b1 b1 84
77
MAESTRIA PYGE
x 2y 6 x, y 0
6
78
MAESTRIA PYGE
x 2y 6 x, y 0
6
79
MAESTRIA PYGE
z(0,4) = 15 x 6 + 20 x 0 = 90 y b1*= 2 x 6 + 2x0 = 12 La menor variacin del coeficiente del lado derecho se alcanza en: x=0 ; y= 3, de donde se obtiene: z(4,0) = 15 x 0 + 20 x 3 = 60 y b1= 2 x 0 + 2 x 3 = 6 De aqu, se calcula el precio sombra P2, que indica la razn o tasa de cambio de la funcin objetivo con respecto al cambio en una unidad del lado derecho: z(6,0) z(0,3) 90 60 P2 5 * b2 b2 12 6
80
MAESTRIA PYGE
81
MAESTRIA PYGE
82
MAESTRIA PYGE
83
x=20 y=60
84
85
MAESTRIA PYGE
86
Sea
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
z = c1x1+c2x2
87
MAESTRIA PYGE
Tambin podemos estudiar el intervalo de un slo coeficiente, dejando el resto de los parmetros fijos:
Ing. Rodrigo Semprtegui lvarez
Para C1:
c1 2 1 2
Para C2:
2 c1 4
3 2 1 c2
3 c2 3 2
88
MAESTRIA PYGE
2) Cul es el intervalo de variacin de los coeficientes del lado derecho (trminos libres) de las restricciones, de modo que la actual solucin siga siendo la ptima? Estudiaremos por separado las variaciones de cada uno de los coeficientes del lado derecho de las restricciones, de modo preservar la geometra del problema, esto es, que se conserven las mismas restricciones activas de la solucin ptima inicial.
89
s.a. 80 40 0
90
MAESTRIA PYGE
z(60,0) = 3 60 + 2 0 = 180 y b1* = 2 60 + 1 0 = 120 La menor variacin del coeficiente del lado derecho se alcanza en: x= 4 ; x2 = 0, de donde se obtiene:
z(40,0) = 3 40 + 2 0 = 120 y b1 = 2 40 + 1 0 = 80
Obsrvese que, aunque para 80b1120 la base actual es ptima los valores de las variables de decisin y de la funcin objetivo cambian. Por ejemplo si 80b1100 la solucin ptima cambiar del punto B a algn otro punto en el segmento AB. Similarmente si 100b1120 la solucin ptima cambiar del punto B a algn otro punto en el segmento AD.
91
MAESTRIA PYGE
Para ilustrar la idea, sea b1 el nmero de horas de acabado disponibles. Si cambiamos b1 a 100+ sabemos que la base ser ptima para -20 20 la solucin para el PL ser todava e punto en el que las restricciones de acabado y carpintera son obligatorias. Por lo tanto si cambiamos b1 = 100+ se puede encontrar los nuevos valores de las variables al resolver 2x + y = 100 + y x + y =80
Esto produce que x=20+, y y=60- lo que significa que si aumentamos el nmero de horas de acabado da como resultado un aumento de nmero de soldados producidos y una disminucin de trenes producidos.
92
MAESTRIA PYGE
93
MAESTRIA PYGE
94
MAESTRIA PYGE
z(0,100) = 30 + 2100 = 200 y b2*= 10 + 1100 = 100 La menor variacin del coeficiente del lado derecho se alcanza en: x=0 ; y=60, de donde se obtiene:
95
MAESTRIA PYGE
Ahora, sea b2 el nmero de horas de carpintera disponibles. Si cambiamos b2 a 80+ sabemos que la base ser ptima para -20 20 la solucin para el PL ser todava e punto en el que las restricciones de acabado y carpintera son obligatorias. Por lo tanto si cambiamos b2 = 80+ se puede encontrar los nuevos valores de las variables al resolver 2x + y = 100 y x + y =80 +
Esto produce que x=20-, y y=60+2 lo que significa que si aumentamos el nmero de horas de carpintera da como resultado una disminucin del nmero de soldados producidos y un disminucin de trenes producidos.
96
MAESTRIA PYGE
De aqu, se calcula el precio sombra P2, que indica la razn o tasa de cambio de la funcin objetivo con respecto al cambio en una unidad del lado derecho:
z(0,100 ) z(0,60 ) 200 120 P2 2 * b2 b2 100 60
97
Supongamos que cambiamos b3 = 40+ se puede encontrar los nuevos valores de las variables al resolver
MAESTRIA PYGE
2x + y = 100
x + y =80
Esto produce que la solucin original x=20, y y=60. Entonces se puede demostrar que la base actual es ptima para un -20, lo que significa que si se cambia el lado derecho de esta restriccin en el intervalo en la cual la base es ptima, la solucin del PL no cambia.
98
Un PL puede tener restricciones en forma de igualdad o de desigualdad. Tambin pueden tener variables que tienen que ser no negativas o no tener restriccin de signo. Para usar el Algoritmo Simplex hay que transformar el PL en un problema equivalente, en el cual:
Todas las restricciones son ecuaciones Todas las variables son no negativas
MAESTRIA PYGE
Una Empresa produce dos tipos de cinturones: El modelo de lujo y el modelo regular. Cada tipo requiere un metro de cuero. El cinturn regular requiere de una hora de trabajo especializado y el de lujo necesita de dos horas. Se dispone semanalmente de 60 horas de mano de obra especializada y 40 metros de cuero. Si cada cinturn regular y cada cinturn de lujo contribuyen a las ganancias con 3 y 4 dlares cada uno; cual es el plan de produccin para generar la mxima utilidad?
100
Si definimos:
MAESTRIA PYGE
X1 = Nmero de cinturones de lujo producidos X2 = Nmero de cinturones regulares producidos El modelo sera: Max Z= 4X1+3X2 s.a. X1+ X2 40 (Restriccin del cuero) 2X1+ X2 60 (Restriccin de mano de obra) X1, X2 0 (Restriccin de no negatividad)
101
MAESTRIA PYGE
102
II. Modelos de Programacin Matemtica Programacin Lineal Ejemplo 2 de transformacin en su forma Estandar MAESTRIA PYGE
Min Z= 50X1+20X2 +30X3 +80X4 s.a. 400X1+ 200X2 + 150X2 + 500X2 500 3X1 + 2X2 6 2X1 + 2X2 + 4X3 + 4X4 10 2X1 + 4X2 + X3 + 5X4 8 X1, X2 , X3 , X4 0 En la Forma Estadar: Min Z= 50X1+20X2 +30X3 +80X4 s.a. 400X1+ 200X2 + 150X2 + 500X2 E1 = 500 3X1 + 2X2 E2 =6 2X1 + 2X2 + 4X3 + 4X4 E3 = 10 2X1 + 4X2 + X3 + 5X4 E4 = 8 X1, X2 , X3 , X4 ,E1 ,E2 ,E3 ,E4 0
103
II. Modelos de Programacin Matemtica Programacin Lineal Ejemplo 3 de transformacin en su forma Estandar MAESTRIA PYGE
Max Z= 20X1+15X2 s.a. X1
Ing. Rodrigo Semprtegui lvarez
En la Forma Estadar: Max Z= 20X1+15X2 s.a. X1 +S1 = 100 X2 +S2 = 100 50X1 + 35X2 +S3 = 4000 20X1 + 15X2 E4 = 2000 X1, X2 ,S1 ,S2 ,S3 ,E4 0
104
MAESTRIA PYGE
105
MAESTRIA PYGE
En resumen, es posible reformular de manera equivalente el problema usando las siguientes justificaciones:
Ing. Rodrigo Semprtegui lvarez
1) Siempre es posible llevar un problema de maximizacin a uno de minimizacin. Si f(x) es la funcin objetivo a maximizar y x* es la solucin ptima: f(x*) f(x) , x factible - f(x*) - f(x) , x factible
106
MAESTRIA PYGE
3) De igual modo, cada restriccin del tipo puede ser llevada a una ecuacin de igualdad usando una variable de exceso no negativa. 4) Siempre es posible escribir una variable libre de signo como la diferencia de dos variables no negativas.
107
MAESTRIA PYGE
DEFINICION: Se obtiene una solucin bsica de Ax = b, haciendo n-m variables (variables no bsicas VNB) iguales a cero y resolviendo el sistema resultante de m variables que quedan (variables bsicas VB).
Naturalmente, las selecciones diferentes de variables no bsicas VNB llevaran a soluciones bsicas diferentes. La bsqueda de la solucin ptima se restringe a encontrar un vrtice ptimo y cada vrtice del conjunto de las restricciones del problema, llamado regin de puntos factibles, corresponde a una solucin bsica factible del sistema Ax = b. Esta solucin bsica factible, corresponde a su vez a aquellas soluciones que resultan de resolver el sistema para exactamente m variables, fijando las restantes n-m en cero, llamadas respectivamente variables bsicas y no-bsicas, que adems deben satisfacer condiciones de no-negatividad
108
Programacin Lineal
MAESTRIA PYGE
max s.a.
= 40 = 60
Solucin Bsica Factible s1=s2= 0 x2=s2= 0 x2=s1= 0 X1=s2=0 x1=s1= 0 x1=x2= 0 x1= 20 x1= 30 x1= 40 s1= -20 x2= 40 s1= 40
109
Programacin Lineal
MAESTRIA PYGE
110
Programacin Lineal
MAESTRIA PYGE
Punto
MAESTRIA PYGE
112
MAESTRIA PYGE
113
MAESTRIA PYGE
n
Ing. Rodrigo Semprtegui lvarez
B A=
D m
x1 xB x x 2 xn xD
nm
c B c c D
nm
n-m
cB :costos bsicos.
cD :costos no bsicos.
114
MAESTRIA PYGE
Max z= 60 x1 + 30 x2 +
Ing. Rodrigo Semprtegui lvarez
20 x3 + 0 s1 + 0 s2 + 0s3
s1 =48 + s2 = 20 + s3 = 20
s.a:
MAESTRIA PYGE
s.a. 1
8
4 2 0 0 0
S1
X3 X1 +
6
2 1.5
0
1 0 x2 s2 s3
0
0 1 0 0 0
x2
s2 s3
0 1.5 0 0.5 S1 X3 X1
116
MAESTRIA PYGE
Donde:
1 B-1= 0 0 2 2 -0.5 -8 -4 1.5
117
MAESTRIA PYGE
1 + 0
2 2
-8 -4
6 2
0 1
0 0
x2 s2 =
1 0
2 2
-8 -4
48 20
x3
x1
-0.5
1.5
1.5
s3
-1
1.5
O bien:
s1
x3 x1 +
-2
-2 1.25
2
2 -0.5
-8
-4 1.5
x2
s2 s3 =
24
8 2
118
MAESTRIA PYGE
CB xB + CB B-1DxD= CBB-1 b z- CB xB - CD xD =0 Sumando las dos ltimas relaciones: z+(CB B-1D- CD)xD= CBB-1 b
119
MAESTRIA PYGE
Criterio de Optimalidad:
Ing. Rodrigo Semprtegui lvarez
T T c T x cB x B cD xD T B
T c B b B D x B cD xD 1
T T T cB B b cD cD B D xB xD
120
MAESTRIA PYGE
La ecuacin que define cada uno de los costos reducidos es: T 1 rj c j cB B Aj Donde j es el ndice de variable no-bsica y Aj la respectiva columna en A de esa variable.
La actual solucin bsica factible es ptima ssi rj j, existe una variable no bsica xp con costo reducido negativo, que entra a la nueva base.
121
MAESTRIA PYGE
Para decidir quin deja la base, es necesario calcular el mayor valor que puede tomar la variable entrante que garantiza la factibilidad de la nueva solucin bsica, con:
y1 0 y1 p x y 2 0 B 1A 2 p B 1b j x y m0 mp
y se debe calcular:
MAESTRIA PYGE
123
MAESTRIA PYGE
x1
2 1 1 -40
x2
1 1 3 -60
s1
1 0 0 0
s2
0 1 0 0
s3
0 0 1 0 70 40 90 0
124
MAESTRIA PYGE
x1
2 1 1 -40
x2
1 1 3 -60
s1
1 0 0 0
s2
0 1 0 0
s3
0 0 1 0 70 40 90 0
Se calcula Min { 70/1, 40/1, 90/3 } = 30, por lo tanto sale s3.
125
MAESTRIA PYGE
Actualizando, queda la siguiente tabla (no ptima), donde la variable entrante a la base es x1 (pues r1<0).
x1
5/3 2/3 1/3 -20
x2
0 0 1 0
s1
1 0 0 0
s2
0 1 0 0
s3
-1/3 -1/3 1/3 20 40 10 30 1800
Se calcula Min { 40/(5/3), 10/(2/3), 30/(1/3) } = 15, por lo tanto s2 deja la base actual.
126
MAESTRIA PYGE
x1
0 1 0 0
x2
0 0 1 0
s1
1 0 0 0
s2
-5/2 -1/3 1/3 20
s3
- 10 15 15 25 2100
Como todos los costos reducidos son mayores o iguales que cero nos encontramos en la solucin ptima.
127
MAESTRIA PYGE
x 1 15 x B x 2 25 s1 15
z* = - 40 x 15 - 60 x 25 = - 2100
s2 0 xD s3 0
En la formulacin inicial, tenemos como solucin ptima x*=15, y *=25, con valor ptimo 2.100.
128
129
Mtodo 1
-Z
X1
X2 S1 S2 ld
Variable Bsica
Razn
MAESTRIA PYGE
1 0 0
2 1 1
-3 1 -1
0 1 0
0 0 1
-Z
1 0 0
X1
5 1 2
X2 S1 S2 ld
0 1 0 3 1 1 0 0 1
Variable
Bsica
Mtodo 2
X1
X2 S1 S2 ld
Variable Bsica
Razn
MAESTRIA PYGE
1 0 0
-2 1 1
3 1 -1
0 1 0
0 0 1
-Z
1 0 0
X1
-5 1 2
X2 S1 S2 ld
0 1 0 -3 1 1 0 0 1
Variable
Bsica
MAESTRIA PYGE
Paso 0: Escribir el problema de programacin lineal en su forma estndar. Paso 1: Escoger una solucin bsica factible inicial. Paso 2: Escoger una variable no - bsica con costo reducido negativo que determina la variable entrante y seguir al paso tres. Sin embargo, si todos los costos reducidos son mayores que cero , parar, ya que la actual solucin es la ptima. Paso 3: Calcular el criterio de factibilidad que determina que variable deja la base. Si todos los cuocientes son negativos: problema no - acotado, parar. Paso 4: Actualizar la tabla de modo de despejar el valor de las nuevas variables bsicas, los costos reducidos y el valor de la funcin objetivo. Volver al Paso 2.
132
II. Modelos de Programacin Matemtica Programacin Lineal II. Modelos de Programacin Matemtica Programacin Lineal
MAESTRIA PYGE
No siempre es fcil obtener una solucin bsica factible inicial, en las variables originales del modelo. Para conseguir esto existen varios procedimientos como son:
133
MAESTRIA PYGE
INTRODUCCIN Hasta el momento slo se han estudiado problemas en la forma estndar. Maximizar Z. Restricciones de la forma menor igual. Todas las variables no negativas
FORMA ESTANDAR Maximizar Z= Sujeto a: 3X1 X1 + 5X2 2X2 2X2
3X1 + X1, X2
4 12 18 0
134
MAESTRIA PYGE
135
MAESTRIA PYGE
136
MAESTRIA PYGE
137
MAESTRIA PYGE
Lo que se hace entonces es introducir variables artificiales Cambiemos la tercera restriccin de desiguladad de Wyndor Glas c.o. por una igualdad
3X1 X1
3X1 + X1, X2
4 12 18 0
138
MAESTRIA PYGE
139
MAESTRIA PYGE
140
MAESTRIA PYGE
141
MAESTRIA PYGE
142
MAESTRIA PYGE
143
MAESTRIA PYGE
144
MAESTRIA PYGE
PROBLEMA REAL
Ing. Rodrigo Semprtegui lvarez
PROBLEMA ARTIFICIAL
5 x2 Max z= s.a. 3 x1 + 5 x2 M x5
Max z= s.a.
x1
x1
x1
2
3 x1 + 2
x2
x2
12
18 3 x1 +
2
2
x2
x2
12
18
x1
x2
0
As: 3
x1
x1
,
+ 2
x2
x2 + x5
0
18
145
MAESTRIA PYGE
146
MAESTRIA PYGE
147
MAESTRIA PYGE
148
MAESTRIA PYGE
149
MAESTRIA PYGE
150
MAESTRIA PYGE
151
MAESTRIA PYGE
152
MAESTRIA PYGE
153
MAESTRIA PYGE
154
MAESTRIA PYGE
155
MAESTRIA PYGE
156
MAESTRIA PYGE
157
MAESTRIA PYGE
158
MAESTRIA PYGE
159
MAESTRIA PYGE
160
MAESTRIA PYGE
Se debe agregar una variable de holgura (x3) y una variable de exceso (x4), y llevarlo a su forma estndar. Min -2x1 - x2 sa: 10x1 + 10x2 +s3 =9 10x1 + 5x2 - e4 = 1 x1,x2, x3, x4 0
161
MAESTRIA PYGE
MAESTRIA PYGE
Fase 1:
Min sa:
x1
x2
s3
e4
a5
valor
MAESTRIA PYGE
xB=
xD=
x2
a5
e4
Luego se hace cero el costo reducido de la variable x5 de la tabla anterior, y queda la siguiente tabla inicial.
Variables
Bsicas w s3 a5 -10 10 10 -5 10 5 0 1 0 1 0 -1 0 0 1 -1 9 1
x1
x2
s3
e4
a5
valor
164
II. Modelos de Programacin Matemtica Programacin Lineal Mtodo Simplex de dos Fases. MAESTRIA PYGE
Variables Bsicas
x1
x2
s3
e4
a5
valor
w
s3 a5
-10
10 10
-5
10 5
0
1 0
1
0 -1
0
0 1
-1
9 1
MAESTRIA PYGE
Bsicas
Ing. Rodrigo Semprtegui lvarez
w s3 x1 s3
0 0 1
0 5 1/2 8
0 1 0
1 1 -1/10
0 -1 1/10
0 8 1/10 x2
xB=
xD=
e4
x1
1/10
a5
Donde, la anterior, corresponde a la solucin ptima del problema en la Fase 1, con valor ptimo 0. De aqu entonces tomamos x1 y x3 como variables bsicas.
166
MAESTRIA PYGE
z
Ing. Rodrigo Semprtegui lvarez
-2
-1
s3 x1
0 1 1/2
1 0 -1/10
1 1/10
Luego la variable entrante a la base es e4 (pues r4<0). Y calculando Min { 8/1, (-1/10)/(1/10) } = 8, se tiene que sale s3.
167
MAESTRIA PYGE
x1
x2
s3
e4
valor
x1
1/10
9/10
168
MAESTRIA PYGE
s.a.
Ing. Rodrigo Semprtegui lvarez
1/2 x1
1 x1
+
+
1/4 x2
3 x2
4
20
x1 x1
+ ,
x2 x2
10 0
2 x1 1/2 x1 1 x1 x1
+ + + +
3 x2 1/4 x2 3 x2 x2 + s1 e2 + a2 + a3 = = = 4 20 10
169
MAESTRIA PYGE
variable w s1 a2 a3
w 1
x1 0 1/2 1 1
x2 0 1/4 3 1
s1 0 1 0 0
e2 0 0 -1 0
a2 -1 0 1 0
a3 -1 0 0 1
valor 0 4 20 10
170
x2
4 0
s1
0 1
e2
-1 0
a2
0 0
a3
0 0
valor
30 4
prueba
MAESTRIA PYGE
w s1
16
a2
a3
0
0
1
1
3
1
0
0
-1
0
1
0
0
1
20
10
7
10
171
MAESTRIA PYGE
w
Ing. Rodrigo Semprtegui lvarez
s1 a2 a3
Realizando operaciones
VB w s1 x2 a3 w x1 1 0 0 0 2/3 5/12 1/3 2/3 x2 s1 e2 0 0 1 0 0 1 0 0 1/3 1/12 -1/3 1/3 a2 -4/3 -1/12 1/3 -1/3 a3 valor 0 0 0 1 10/3 7/3 20/3 10/3 28/5 20 1/5 prueba
172
VB
x2 0 0 1
s1 0 1 0
a3 0 0 0
MAESTRIA PYGE
w
Ing. Rodrigo Semprtegui lvarez
s1 x2
x1
1/2
-1/2
3/2
VB w
w x1 1 0
x2 0
s1 0
e2 0
a2 -1
a3 -1
valor 0
s1
x2 x1
0
0 0
0
0 1
0
1 0
1
0 0
-1/8
-1/2 1/2
1/8
1/2 -1/2
-5/8
-1/2 3/2
1/4
5 5
173
z x1
1 0 0 0 -2 0 0 1
x2
-3 0 1 0
s1
0 1 0 0
e2
0 -1/8 -1/2 1/2
valor
0 1/4 5 5
MAESTRIA PYGE
s1
Ing. Rodrigo Semprtegui lvarez
x2 x1
Como tanto x1 como x2 estn en la base ptima de la fase 1 debemos eliminarlas del rengln 0
VB z z x1 1 0 x2 0 s1 0 e2 -1/2 valor 25
s1
x2 x1
0
0 0
0
0 1
0
1 0
1
0 0
-1/8
-1/2 1/2
1/4
5 5
Debido a que no hay como mejorar la funcin objetivo original Esta es OPTIMA
174
MAESTRIA PYGE
s.a.
Ing. Rodrigo Semprtegui lvarez
1/2 x1
1 x1 x1
+
+ +
1/4 x2
3 x2 x2
4
36 10
x1
x2
= 36 = 10
175
+ a3
= = = 4 36 10
MAESTRIA PYGE
variable w s1 a2 a3
w x1 1 0 0 0 0 1/2 1 1
x2 0 1/4 3 1
s1 0 1 0 0
e2 0 0 -1 0
a2 -1 0 1 0
a3 -1 0 0 1
valor 0 4 36 10
176
MAESTRIA PYGE
VB
w x1
x2
s1
e2
a2
a3
valor
prueba
w
s1
1
0
2
1/2
4
1/4
0
1
-1
0
0
0
0
0
46
4 16
a2
a3
0
0
1
1
3
1
0
0
-1
0
1
0
0
1
36
10
12
10
177
MAESTRIA PYGE
w
Ing. Rodrigo Semprtegui lvarez
s1 a2 x2
Como ninguna variable en el rengln 0 tiene coeficiente positivo se trata de un cuadro ptimo de la fase 1. Debido a que el valor de w=6>0 el PL no tiene solucin factible
178
1) Problema Infactible. Esta situacin se detecta cuando el valor ptimo del problema de la Fase 1 da mayor que cero. 2) Mltiples soluciones ptimas. Esta situacin se detecta cuando existen costos reducidos iguales a cero en una o ms de las variables bsicas ptimas. 3) Problema no acotado. Esta situacin se detecta cuando al realizar el clculo de la variable que deja la base, todos los elementos ykj de la columna j en la tabla, son negativos para j el ndice de una variable no bsica con costo reducido negativo.
179
MAESTRIA PYGE
MAESTRIA PYGE
En lo que sigue, combinaremos las distintas restricciones del problema, ponderando por los valores 1, 2 y 3 cada una, respectivamente, de modo de obtener la mejor cota superior del valor ptimo del problema P). Vale decir: 1(2x1+2x2) + 2(x1+x2) + 3(x1+3x2) 70 1 + 40 2 + 90 3 Para garantizar que el lado derecho de esta ltima desigualdad sea una cota superior de la funcin objetivo se debe cumplir que : 2 1 + 2 + 3 40 21 + 2 + 3 3 60
181
MAESTRIA PYGE
70 1 + 40 2 + 90 3 2 1 + 2 + 3 40 21 + 2 + 3 3 60 1, 2, 3 0 Este problema se conoce como el problema Dual D) asociado al problema Primal P). D) Min sa: Tambin resulta que al formular el problema dual de D) se obtiene el problema primal (o uno equivalente). Cualquiera de los dos entrega la misma informacin y el valor ptimo alcanzado es el mismo.
182
Max
c jx j
j1 n j1
MAESTRIA PYGE
sa :
aijx j bi
xj 0
i 1,2,..., n
j 1,2,..., m
bi i
i 1 m
sa :
aiji c j
i 1
j 1,2,..., n
i 0
i 1,2,..., m
183
MAESTRIA PYGE
bT AT c 0 Si el problema primal corresponde a: P) Max -cTx sa: Ax b x0 Su dual resulta ser: D) Min -bT sa: AT c 0 Es decir, el dual del dual es el problema primal D) Min sa:
184
MAESTRIA PYGE
Teorema de dualidad dbil: Si x IRn, es una solucin factible del problema primal P) y IRm, una solucin factible del problema dual D), entonces:
c x c jx j bi i bT
T j 1 i 1
En particular, si ambas soluciones son los ptimos de sus respectivos problemas, sus valores ptimos cumplen que : v(P) v(D) Teorema de dualidad fuerte: Si x* = (x1*, x2*, ..., xn*)T, es una solucin ptima problema primal P), entonces el problema dual D) tiene solucin ptima * = (1*, 2*, ..., m*)T que satisface:
v(P) c x * c j x j * bi i * b T v(D)
T j1 i 1
185
MAESTRIA PYGE
Adems: i)Si P) es no-acotado entonces D) es infactible. ii)Si D) es no-acotado entonces P) es infactible. Ejemplo: P) Min sa: 3x1 + 4x2 + 5x3 x1+ 2x2 + 3x3 5 2x1 + 2x2 + x3 6 x1 , x2 , x3 0 5 1 + 6 2 1 + 22 3 21 + 22 4 31 + 2 5 1, 2 0
D)
Max sa:
186
1
MAESTRIA PYGE
1
Ing. Rodrigo Semprtegui lvarez
2
2 2 1 -6
3 4
1 0 0 0 0 1 0 0
5
0 0 1 0 3 4 5 0
2 3 -5
3 3 xB 4 4 5 5 1 0 xD 2 0
Luego la variable entrante a la base es 2 (pues r2<0). Y calculando Min { 3/2, 4/2, 5/1 } = 3/2, se tiene que sale 3
187
1 2 3 4 5
MAESTRIA PYGE
1 0 0
-1 -1/2
0 1 0
0 0 1
3/2 1 7/2
1 5/2
-2 0 3 0 0 9 Luego la variable entrante a la base es 1 (pues r2<0). Y calculando Min { (3/2)/(1/2), 1/1, (7/2)/(5/2)} = 1, se tiene que sale 4
2 3 / 2 xB 4 1 5 7 / 2 1 0 xD 3 0
188
MAESTRIA PYGE
1 2 3 4 5
0 1 0 0 1 0 0 0 1 -1 2 1 -1/2 1 -5/2 2 0 0 1 0 1 1 1 11
1 1 x B 2 1 5 1 3 0 xD 4 0
1* = 1; 2* = 1;
Sol. ptima de P): x1* = 1; x2* = 2; x3* = 0;
v(D) = 11
v(P) = 11
189
MAESTRIA PYGE
La idea de este mtodo consiste en resolver de alguna manera el problema dual asociado a P) en la tabla y variables del problema primal P), segn veremos en su aplicacin a un problema primal (ejercicio anterior). Min 3x1 + 4x2 + 5x3 sa: x1+ 2x2 + 3x3 5 2x1 + 2x2 + x3 6 x1 , x2 , x3 0
190
MAESTRIA PYGE
Min sa:
x(-1) x(-1)
x1
-1 -2 3
x2
-2 -2 4
x3
-3 -1 5
x4
1 0 0
x5
0 1 0 -5 -6 0
191
MAESTRIA PYGE
En la tabla anterior se toman dos variables de exceso x4 y x5 , y se multiplica por un nmero negativo con la finalidad de encontrar la matriz identidad IRn, adems es necesaria la condicin de que los costos reducidos de la tabla sean mayores que cero ( lo que en este caso se cumple). En la tabla anterior se escoge, usando el lado derecho, alguna variable con valor negativo.
Escogemos x5 , variable que dejar la base. Enseguida , se obtiene la variable entrante calculando: Min { (-3/-2) , (-4/-2),(-5/-1)} = 3/2. De donde resulta que x1 entra a la base.
192
MAESTRIA PYGE
x1
0 1 0
x2
-1 1 1
x3
-5/2
x4
1 0 0
x5
-1/2 -1/2 3/2 -2 3 -9
1/2 7/2
La tabla posee an un lado derecho negativo (costos reducidos negativos del problema dual), por lo cual no es factible en P).
193
MAESTRIA PYGE
x1
0 1 0
x2
1 0 0
x3 x4 x5
5/2 -1 -2 1 1 1 -1 1 2 1 -11
La tabla posee lados derechos no-negativos (costos reducidos positivos del problema dual) y tambin los costos reducidos de las variables no bsicas x3, x4 y x5 son no-negativos , por lo que tenemos una solucin factible en P) que es la solucin ptima del problema.
x 1 1 x x 2 2 x 3 0 v(P) 11
194
II. Modelos de Programacin Matemtica Programacin Lineal II.6. Anlisis de Sensibilidad o Post-Optimal MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
1) Qu ocurre con las actuales variables bsicas si se cambia algn coeficiente del lado derecho (b)? Si calculamos: xB B b y se cumple: xB 0 Las mismas variables bsicas lo son tambin de la nueva solucin ptima, calculada con el nuevo b . Si lo anterior no se cumple, se puede aplicar el Mtodo Simplex Dual. 2) Qu ocurre con la actual solucin ptima si se agrega una nueva variable al problema ?
1
Para decidir si la actual solucin bsica es ptima para el nuevo problema, calculamos el costo reducido de la nueva variable mediante la formula:
T 1 rk ck cB B Ak
195
MAESTRIA PYGE
P)
Min c x sa : Ax b x0
196
MAESTRIA PYGE
Siempre que los nuevos costos reducidos sean mayores o iguales a cero (notar que tambin cambia el valor de la funcin objetivo en la actual solucin ptima). Es decir se debe cumplir que:
T 1 rD cD cB B D0
o equivalentemente j
rj c j c B 1 Aj 0
T B
En caso contrario, se aplica el Simplex a partir de la tabla final de P) con los nuevos costos reducidos y nuevo valor de la actual solucin bsica.
197
MAESTRIA PYGE
Veamos los cambios que tienen lugar cuando slo vara un coeficiente del vector c de la funcin obj.
Ing. Rodrigo Semprtegui lvarez
a) Cambio de un coeficiente asociado a una variable nobsica xJ: Se conserva la misma solucin ptima del problema P) ssi. para esa variable xJ:
rj
T 1 c j cBB A j
198
c j c j j
Por lo tanto se conserva la misma solucin ssi:
j rj
c j c j rj
199
MAESTRIA PYGE
b) Cambio en un coeficiente de la funcin objetivo asociado a una variable bsica: En este caso para tener la misma solucin ptima, se debe cumplir que el costo reducido de todas las variables.a cero.
rj
T 1 c j cBB A j
ci ci i
0 cB cB i 1 cB ie i 0
200
MAESTRIA PYGE
MAESTRIA PYGE
Ejemplo:
Ing. Rodrigo Semprtegui lvarez
0,00 -0,03
MAESTRIA PYGE
Las xB del problema primal no cambian como base ptima, si los valores asociados a estas variables.
xB B1b y se cumple xB 0
Para calcular estos intervalos de recursos, se necesita la matriz inversa asociada a las variables bsicas del tabla final.
203
MAESTRIA PYGE
4 10 B 1 40
1/ 15 4 / 15 B 1 / 150 2 / 75
1
Intervalo recurso 1:
1 / 15 6000 b1 4 / 15 1 / 150 2 / 75 x 4000 0
20000 4 b1 0 15 15 10000 b1 0 150 150
b1 5000
b1 10000
204
MAESTRIA PYGE
Intervalo recurso 2:
1 / 15 6000 4 / 15 1 / 150 2 / 75 x 4000 b 0 2
205
MAESTRIA PYGE
Variable x1:
Ing. Rodrigo Semprtegui lvarez
0 D1 2
Variable x4:
10 C1* 12
MAESTRIA PYGE
Variable x2:
Ing. Rodrigo Semprtegui lvarez
C3* = C3 + 3 3 - r3
1. Linear Programming and Network Flow, M.Bazaraa, J.Jarvis and H.Sherali. John Wiley & Sons, Inc., New York, Second Edition 1990. 2. Introduction to Linear Optimization, D.Bertsekas and J.Tsitsiklis. Athena Scientific USA, 1997. 3. Linear Programming, V.Chvtal. W.H. Freeman and Company, New York, 1983. 4. Linear Programming and Extensions, G. Dantzig. Princeton University Press, New Jersey, tenth printing, 1993. 5. Introduccin a la Programacin Lineal y No Lineal, D.Luenberger. Adisson Wesley Iberoamericana, 1989. 6. Linear and Combinatorial Programming, K. Murty. John Wiley & Sons, Inc., New York, Second Edition 1976. 7. Model Building in Mathematical Programming, H.P. Williams. John Wiley & Sons, Inc., New York, 4rd Edition 1999.
208
MAESTRIA PYGE
MAESTRIA PYGE
209
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
III.1. Introduccin y ejemplos de modelamiento. III.2. Resolucin de problemas de P. E. III.3. Mtodo de Branch and Bound.
211
MAESTRIA PYGE
Seat considera la fabricacin de 3 tipos de automviles: compacto, mediano y largo. En la tabla se presentan los recursos requeridos y las ganancias proporcionadas por cada tipo de automvil. En la actualidad se cuenta con 6000 ton de acero y 60000 horas de trabajo. Para que la produccin de un tipo de automvil sea econmicamente factible hay que fabricar al menos 1000 unidades de este tipo. Formular un modelo matemtico para optimizar las ganancias.
LARGO 5 40
ganancia (dolares)
2000
3000
4000
212
MAESTRIA PYGE
213
Se dan dos restricciones de la forma: f(x1, x2, xn) 0 g(x1, x2, xn) 0 Se desea asegurar que se cumpla por lo menos una de las dos restricciones. Deber agregarse las siguientes dos restricciones: f(x1, x2, xn) My g(x1, x2, xn)M(1-y) Donde: y es una variable 0-1 y M es un nmero muy grande que asegure que se satisfagan: f(x1, x2, xn)M g(x1, x2, xn)M Para todos los valores de x1, x2, xn que satisfagan las otras restricciones.
214
MAESTRIA PYGE
x1 M1 y1 1000- x1 M1(1-y1) R2 x2 M2 y2 1000- x2 M2(1-y2) R3 x3 M3 y3 1000- x3 M3(1-y3) R4 1.5 x1 + 3 x2 + 5 x3 6000 R4 30 x1 +25 x2 + 40 x3 60000 x1, x2, x3 0 y1, y2, y3 = 0 1 R1
215
MAESTRIA PYGE
Se desea asegurar que se cumpla g(x1,x2,xn)0 si se satisface f(x1,x2,xn)>0, mientras que si no se satisface f(x1,x2,xn)>0 entonces g(x1,x2,xn)0 no debe satisfacerse. Deber agregarse las siguientes dos restricciones: -g(x1, x2, xn)My f(x1, x2, xn) M(1-y) Donde: y es una variable 0-1 y M es un nmero muy grande que asegure que se satisfagan: f(x1, x2, xn)M g(x1, x2, xn)M Para todos los valores de x1, x2, xn que satisfagan las otras restricciones.
217
II. Modelos de Programacin Matemtica Programacin Entera Funciones lineales por partes
Supongamos que una funcin lineal por partes tiene los puntos de ruptura b1,b2,bn. Para algn k (k=1,2,3n-1), bk x bk+1. Entonces, para algn nmero zk (0 zk 1) se puede escribir x como: x= zk bk+(1-zk)bk+1 Ya que f(x) es lineal para bk x bk+1 , podemos escribir como: f(x)=zk f(bk)+(1-zk) f(bk+1 ) Paso 1: donde se presenta f(x) remplazar por z1f(b)+z2f(b2)+z3f(b3)..znf(bn) Paso 2: aadir las siguientes restricciones: z1y1, z2y1+y2 , z3y1+y2. zn-1 yn-2+yn-1, zn 1 y1+y2+y3.+yn-1=1 z1+z2+z3.+zn-1=1 X=z1b1+z2b2+z3b3+.znbn
MAESTRIA PYGE
218
MAESTRIA PYGE
Supngase que se produce gasolina a partir de petoleo. Al comprar el petrleo de nuestro proveedor se recibe un descuento por cantidad. Los 500 primeros galones de petrleo cuestan 0.25 dlares el galn; los siguientes 500 galones a 0.20 dlares; y los siguientes 500 galones a 0.15 dlares Se pueden comprar como mximo 1500 galones
219
MAESTRIA PYGE
220
MAESTRIA PYGE
Disp. Recursos
30 15 12
Interesa determinar en cules proyectos invertir de modo de conseguir el mayor V.P.N. de la inversin.
221
MAESTRIA PYGE
con i 1,2,3,4
222
MAESTRIA PYGE
= 30 + s2 = 15 + s1 12 + s2
2) Sin invertir el dinero no utilizado en un perodo, pero utilizando el retorno de los proyectos concludos: Ao1: 10x1 + 8x2 + 6x3 + 12x4 30 Ao2: 8x1 + 15x2 + 4x3 Ao3: 18x1 xi {0,1} + 16x3 i = 1,2,3,4 15 + 16x4 12 + 18x2
223
MAESTRIA PYGE
Ao1: 10x1+ 8x2+ 6x3+ 12x4+ s1 Ao2: 8x1+ 15x2+ 4x3 Ao3: 18x1 + 16x3
= 30 12 + s2 + 18x2
+ s2 = 15 + s1 + 16x4
xi {0,1}
i = 1,2,3,4
224
MAESTRIA PYGE
225
MAESTRIA PYGE
x1 + x2 + x3 1
- El proyecto 2 no puede ser tomado a menos que el proyecto 3 si sea tomado: x2 x3 - Se puede tomar el proyecto 3 o 4 pero no ambos: x3 + x4 1 - No se puede invertir en ms de dos proyectos: x1 + x2 + x3 + x4 2
226
MAESTRIA PYGE
12x1 + 24x2 + 18x3 2400 15x1 + 32x2 + 12x3 1800 20x1 + 15x2 + 20x3 2000
Supongamos adems, que nos basta con obtener alguna solucion ptima que verifique el cumplimiento de al menos 2 de las 3 restricciones anteriores.
Variables de decisin:
227
MAESTRIA PYGE
15x1 + 32x2 + 12x3 1800 + M2 (1- y2) 20x1 + 15x2 + 20x3 2000 + M3 (1- y3) Adems, debemos agregar la restriccin que permita que a lo ms una de las restricciones no se cumpla: y1 + y2 + y3 2 Mi = constante lo suf. grande
228
MAESTRIA PYGE
Supongamos que se desea tener lotes de compra de un producto dado, para satisfacer demandas que fluctan en el tiempo sobre un horizonte de planificacin dividido en T perodos. Asumimos conocidos: una estimacin de la demanda dt, con t = 1, 2, ..., T, los costos fijos asociados a la compra de una unidad pt, los costos asociados al mantenimiento de una unidad en inventario de cada perodo ht y los costos fijos asociados a la gestin de compra en el perodo t, st. Observacin: no se permite unidades de faltante.
229
MAESTRIA PYGE
x t: It :
nmero de unidades compradas en t. nivel de inventario al final del perodo t. con t: 1, 2, ..., T
230
MAESTRIA PYGE
Min
Restricciones xt + It-1 - It = dt xt Mt yt
st
t 1
y t p t x t h t It
231
MAESTRIA PYGE
Se conoce la informacin relativa a que zonas pueden ser atendidas por cada uno de los n potenciales servidores, es decir, se conoce la matriz de incidencia A = (aij) donde :
1, si la zona i puede ser atendida por el servidor j aij 0, sin o con i 1,2,..., m y j 1,2,..., n
232
MAESTRIA PYGE
Variables de decisin:
Ing. Rodrigo Semprtegui lvarez
Min
cj xj
j 1
j1
aij x j 1
Se agrega la siguiente restriccin, si adicionalmente, hay algn lmite en el nmero de servidores que se pueden instalar (digamos k) : m
xj k
j 1
233
MAESTRIA PYGE
Si se tiene un conjunto de m clientes que demandan di unidades de un producto determinado. Una compaa desea satisfacer esas demandas desde un cierto conjunto de plantas elegidas de n potenciales lugares donde se instalarn. Sean cj los costos asociados a la instalacin de la planta j , vj el costo unitario de produccin de la planta j y tij el costo de transporte de una unidad desde la planta j al cliente i . Se desea decidir cules plantas abrir y el tamao de cada una de modo de satisfacer las demandas estimadas.
234
Variables de decisin:
MAESTRIA PYGE
235
MAESTRIA PYGE
Funcin objetivo:
m Min c j y j v j xij i 1 j 1 j 1 n n
tijxij
j 1 i 1
Costo de Instalacin
Costo de Produccin
Costo de Transporte
236
1) Demanda cliente i:
Ing. Rodrigo Semprtegui lvarez
x ij di
i 1
2) Relacionar variables de produccin con las asociadas a la apertura de plantas (variables binarias):
xij Mj y j
j 1
donde Mj es una constante grande (por ejemplo, capacidad mxima de produccin de la planta j), con xij 0 e yj {0,1}.
237
Supongamos que tenemos el siguiente problema de programacin lineal: PL) Max s.a. cTx Ax=b x0 Pero todas o una parte de las variables deben restringir su valor a nmeros enteros, dando origen a un problema de Programacin Entera (puro) o de Programacin Entera- Mixta, respectivamente.
238
Por ejemplo:
MAESTRIA PYGE
PLE)
Ing. Rodrigo Semprtegui lvarez
Max s.a.
El problema PL) corresponde a la relajacin continua del problema PLE), que resulta de eliminar las condiciones de integralidad de las variables de decisin en PLE). El valor ptimo de PL) provee slo una cota superior del valor ptimo de PLE). Notar sin embargo, que si la solucin ptima de PL) cumple con la integralidad de los valores requiridos, entonces esta solucin es tambin solucin ptima de PLE).
239
MAESTRIA PYGE
240
MAESTRIA PYGE
241
MAESTRIA PYGE
242
enteros
243
MAESTRIA PYGE
2x1 + x2 7
7
x2
- 2x1 + 2x2 1
. . . .
. . . . .
3.5
x1
244
MAESTRIA PYGE
x1 * = 1
x1 * = 2
x2 * = 1
x2 * = 1
Esta alternativa de enumeracin queda naturalmente restringida a problemas muy pequeos. Alternativamente, podemos resolver la relajacin continua asociada al problema PLE). Si la solucin ptima de la relajacin continua da una solucin entera, esa es la solucin ptima no solo del problema lineal sino que tambin lo es del problema lineal entero. En el ejemplo, la solucin de la relajacin continua es: x1 = 3/2 x2 = 2
245
MAESTRIA PYGE
MAESTRIA PYGE
x2 = 9/10
Ing. Rodrigo Semprtegui lvarez
x1 = 1
x2 = 1 x1 = 0;
infactible
x1 = 1
x2 = 0
f(1,0)=1
247
MAESTRIA PYGE
PLE)
Max s.a.
Consideremos inicialmente la resolucin de la relajacin continua de PLE), que consiste en eliminar las condiciones de integralidad.
248
x2 = 3
3/2 2 x2 = 2
21x1+11x2=39
1 x2 = 1 13/7 sol. relajada
x1 = 1
x1
x1 = 2
21x1+11x2
7x1+4x2=13
249
MAESTRIA PYGE
MAESTRIA PYGE
Solucin ptima de PLE), la solucin entera asociada a la actual cota inferior de v(PLE), si existe (si no existe entonces PLE) es infactible) Si el problema no est agotado pasar al paso 2. Paso 2 Seleccionar una variable xj= j, cuyo valor en la solucin ptima de Pi) no de entero. Eliminar la regin correspondiente a j < j < j + 1 Crear dos nuevos problemas de programacin lineal que incorporen a Pi) dos restricciones mutuamente excluyentes: xj j, xj j +1 una en cada problema y volver al paso 1.
251
MAESTRIA PYGE
P0
x1 = 1 x2 = 3/2 z = 37.5
x21 x1 = 1 x2 = 1 z = 32 x11
x1 = 13/7 x2 = 0 z = 39 x12
P1
x22
P2
infactible x1 = 5/7 x2 = 2 z = 37 x11
P11 P121
P12
P122
x24 infactible
P1212
infactible
P0) Relajacin continua -< z 39 P1) Max 21x1 + 11x2 s.a. 7x1 + 4x2 13 x1 1 x1 0 x2 0 P2) Max 21 x1 + 11x2 s.a. 7x1 + 4x2 13 x1 1 x2 1 x1 0 x2 0 De donde 32 z 39 Solucin ptima x1* = 0; x2* = 3; z = 33 252
MAESTRIA PYGE
1) Integer Programming, L.A.Wolsey. John Wiley & Sons, Inc., New York, 1998. 2) Combinatorial Optimization C.H.Papadimitriou and K.Steiglitz. Prentice Hall Inc., USA, 1982. 3) Linear and Combinatorial Programming, K. Murty. John Wiley & Sons, Inc., New York, Second Edition 1976. 4) Integer and Combinatorial Optimization, George L. Nemhauser and Laurence A. Wolsey. John Wiley & Sons, Inc., New York, 1999. 5) Model Building in Mathematical Programming, H.P. Williams. John Wiley & Sons, Inc., New York, 4rd Edition 1999.
253
MAESTRIA PYGE
254
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
PROGRAMACION SEPARABLE
max o min Z f j ( x j )
i 1
s.a.
g
j 1
ij
( x j ) bi
Ejemplo
max Z 30 x1 2 x 35 x 2 2 x
2 1
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
2 2
s.a. x 2 x 250
2 1 2 2
x1 x 2 20 x1 , x1 0
max Z f 1 ( x1 ) f 2 ( x 2 ) donde
MAESTRIA PYGE
f1 ( x1 ) 30 x1 2 x12
Ing. Rodrigo Semprtegui lvarez
f 2 ( x 2 ) 35x 2 2 x
2 2
g 22 ( x 2 ) x 2
i xi i
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
0 x1 16 0 x2 12
sup ongamos pi , j xi pi , j 1
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
generali z a ndo
j1 j 2 j 3 j 4 ...... jn 1
MAESTRIA PYGE
x j j1 p j1 j 2 p j 2 j 3 p j 3 j 4 p j 4 ...... jn p jn f j ( x j ) j1 f j ( p j1 ) j 2 f j ( p j 2 ) j 3 f j ( p j 3 ) j 4 f j ( p j 4 ) ...... jn f j ( p jn ) gij ( x j ) j1 gij ( p j1 ) j 2 gij ( p j 2 ) j 3 gij ( p j 3 ) j 4 gij ( p j 4 ) ...... jn gij ( p jn ) donde exiten solamente 2 con sec utivos que son diferentes de cero
para el ejemplo
MAESTRIA PYGE
f1 ( x1 ) 11 f1 ( p11 ) 12 f1 ( p12 ) 13 f1 ( p13 ) 14 f1 ( p14 ) 15 f1 ( p15 ) f1 ( x1 ) 11 0 12 88 13 112 14 72 15 ( 32) f 2 ( x2 ) 21 f 2 ( p21 ) 22 f 2 ( p22 ) 23 f 2 ( p23 ) 24 f 2 ( p24 ) 25 f 2 ( p25 ) f 2 ( x2 ) 21 0 22 87 23 138 24 153) 25132 g11 ( x1 ) 11 g11 (0) 12 g11 ( 4) 13 g11 (8) 14 g11 (12) 15 g11 (16) g11 ( x1 ) 11 0 1216 13 64 14144 15 256 g12 ( x2 ) 12 g12 (0) 22 g12 (3) 23 g12 (6) 24 g12 (9) 25 g12 (12) g12 ( x2 ) 12 0 2218 23 72 24162 25 288 g 21 ( x1 ) 12 g12 (0) 22 g12 ( 4) 23 g12 (8) 24 g12 (12) 25 g12 (16) g 21 ( x1 ) 12 0 22 4 23 8 2412 2516 g 22 ( x1 ) 21 g 22 ( p21 ) 22 g 22 ( p22 ) 23 g 22 ( p23 ) 24 g 22 ( p24 ) 25 g 22 ( p25 ) g 22 ( x1 ) 21 0 22 6 23 6 24 9 2512
MAESTRIA PYGE
11 y1 12 y1 y2 13 y2 y3 14 y3 y4 15 y4 21 z1 22 z1 z2 23 z2 z3 24 z3 z4 25 z4
y1 y2 y3 y4 y5 1 z1 z2 z3 z4 z5 1 0 ij 1 y {0,1} z {0,1}
263
MAESTRIA PYGE
PROGRAMACIN NO
MAESTRIA PYGE
265
MAESTRIA PYGE
Contenidos
Ing. Rodrigo Semprtegui lvarez
III. Programacin
Dinmica
266
como no lineales. La programacin dinmica es til para resolver un problema donde se deben tomar una serie de decisiones interrelacionadas. A diferencia de la P.L., la programacin dinmica no tiene formulacin matemtica estndar. Se trata de un enfoque de tipo general para la solucin de problemas, y las ecuaciones se derivan de las condiciones individuales de los mismos.
267
El problema de la diligencia
Un cazafortunas desea ir de Missouri a California en una diligencia, y quiere viajar de la forma ms segura posible. Tiene los puntos de salida y destino conocidos, pero tiene mltiples opciones para viajar a travs del territorio. Se entera de la posibilidad de adquirir seguro de vida como pasajero de la diligencia. El costo de la pliza estndar (Cij) se muestra en la tabla de la siguiente pgina. Cul es la ruta que minimiza el costo de la pliza de seguro?
MAESTRIA PYGE
268
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
269
2.
3.
Enumeracin exhaustiva: enumerar todas las rutas posibles, calcular su costo y elegir la de menor valor. En total son 18. Elegir la ruta ms barata en cada etapa. Esta solucin no necesariamente conduce al ptimo global. Un pequeo sacrificio en una etapa puede permitir mayores ahorros ms adelante Programacin Dinmica:
Estrategia de solucin: un problema complejo es desagregado en problemas ms simples que se resuelven etapa por etapa En el caso de la diligencia un problema simple sera pensar qu pasara si al viajero slo le faltara una etapa del viaje
270
MAESTRIA PYGE
271
MAESTRIA PYGE
fn(S, Xn)
= =
272
MAESTRIA PYGE
f5*(J) = 0 El objetivo es hallar f1*(A) y su ruta correspondiente. Cuando el cazafortunas tiene slo una etapa por recorrer (n=4) , su ruta de ah en adelante, estar determinada por el estado actual (H o I) y su destino final X4=J La ruta ser: S J donde S= H o I Luego f4*(S) = cS,J + f5*(J) = cS,J f4(H) = cH,J = 3 f4(I) = cI,J = 4
273
MAESTRIA PYGE
274
Etapa n=2
En la segunda etapa, el cazafortunas tiene 3 jornadas por recorrer (n=2). Suponga que sale de C.
MAESTRIA PYGE
275
Etapa n=1
En la primera etapa, el cazafortunas tiene todas las jornadas por recorrer (n=1). Necesriamente debe salir de A.
MAESTRIA PYGE
276
277
MAESTRIA PYGE
278
MAESTRIA PYGE
6. El procedimiento de la solucin termina cuando se obtiene la poltica ptima de la ltima etapa (por lo general la solucin en esta etapa es trivial) 7. Siempre se dispone de una relacin recursiva (esto es lo que permite trabajar las decisiones interrelacionadas). La relacin recursiva ser: fn* (Sn)= Max Xn { fn (Sn, Xn) } fn* (Sn)= Min Xn { fn (Sn, Xn) }
Donde: N: Nmero de etapas n: etiqueta para la etapa actual (n=1,2,3,N) Sn: Estado actual para la etapa n Xn: Variable de decisin para la etapa n
279
MAESTRIA PYGE
8.Cuando se tiene una relacin recursiva como la de la funcin, el procedimiento de solucin hacia atrs inicia en la ltima etapa y se mueve hacia la primera, etapa por etapa
Ing. Rodrigo Semprtegui lvarez
280
MAESTRIA PYGE
ALGORITMO DE P.D. HACIA ATRS: Para cada valor probable de la variable de estado al inicio de la etapa, determinar el mejor estado final
ALGORITMO DE P.D. HACIA ADELANTE: Para cada valor probable de la variable de estado al final de la etapa, determinar el mejor estado inicial
281
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
INTRODUCCIN
Esta es una nueva herramienta que nos ayuda en la
beneficios.
Estando dados de la siguiente forma:
INTRODUCCIN
Un algoritmo eficiente que se puede utilizar es el simplex.
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
Para el caso de tener muchas variables este ya no resulta eficiente, por el tiempo que requiere para realizar los clculos.
INTRODUCCIN
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
Existen los llamados Problemas de Optimizacin Combinatorios que contienen variables de decisin enteras y el espacio de soluciones est formado por ordenaciones o subconjuntos de nmeros naturales.
Los ms conocidos son:
Problema de la mochila Problema del viajante
Consiste en: Seleccionar de entre un conjunto de n productos, cada uno con un valor ci y un volumen vi. Determinando aquellos que quepan en un recipiente con volumen V y que tengan el mayor valor posible. As se determina un subconjunto para el cual:
|* 1..n
c
i
max
i 1..n
c
i
con la restricci n
v
i
xi toma el valor de 1 cuando el item se introduce en la mochila xi toma el valor de 0 en caso contrario. As:
max ci xi
i 1
v x
i 1 i
V i 1..n
xi 0,1
ciudades de forma que se recorra el menor nmero de kilmetros y solo se visite una sola vez cada ciudad. Solucin Optima es la ruta
<1, 6, 3, 4, 2, 5>
6
15 20
1
12 16
16
3
20 7 8 21 10
5
9
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
Complejidad Computacional
Problema de la mochila
Para el problema de la mochila el nmero de
Problemas P
Aquellos para los cuales se conocen algoritmos que
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
Problemas NP
Aquellos para los cuales no se conoce un algoritmo
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
P NP
Problemas NP-completos
La mayora de problemas de inters empresarial son
MAESTRIA PYGE
NP-completos.
Ing. Rodrigo Semprtegui lvarez
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
Heursticas
HEURSTICAS
Del griego heuriskein que significa encontrar.
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
Son algoritmos que ofrecen soluciones factibles, que aunque no optimicen la funcin objetivo, se supone que al menos se acercan al valor ptimo.
Las heursticas son utilizadas como una herramienta til que da soluciones a problemas reales.
HEURSTICAS
Definicin:
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
Procedimientos simples, a menudo basados en el sentido comn, que se supone ofrecern una buena solucin (aunque no necesariamente la ptima) a problemas difciles, de un modo fcil y rpido Zanakis, Evans 1981
HEURSTICAS
Factores para la utilizacin de mtodos heursticos:
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
Cuando no existe un mtodo exacto de resolucin o ste requiere mucho tiempo de clculo o memoria. b. Cuando no se necesita la solucin ptima. c. Cuando los datos son poco fiables. d. Cuando las hay limitaciones de tiempo, espacio, etc. e. Como paso intermedio en la aplicacin de otro algoritmo.
a.
HEURSTICAS
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
Ventajas: Permiten mayor flexibilidad para el manejo de las caractersticas del problema.
Ofrecen ms de una solucin.
HEURSTICAS
Inconvenientes:
No es posible conocer la calidad de la solucin.
MAESTRIA PYGE
procedimiento es el de relajar el problema, de manera que sea ms fcil de resolver. utilizar mtodos que indican que la heurstica no es buena. heurstica.
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
TIPOS DE HEURSTICAS
TIPOS DE HEURSTICAS
MTODOS CONSTRUCTIVOS
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
MTODOS DE DESCOMPOSICIN MTODOS DE REDUCCIN MANIPULACIN DEL MODELO MTODOS DE BSQUEDA POR ENTORNOS
MTODOS CONSTRUCTIVOS
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
Aaden paulatinamente componentes individuales a la solucin, hasta que se obtiene una solucin factible.
ALGORITMOS
GOLOSOS O DEVORA-DORES
(GREEDY)
MTODOS DE DESCOMPOSICIN
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
MTODOS DE REDUCCIN
MAESTRIA PYGE
Identifican
Ing. Rodrigo Semprtegui lvarez
Modifican la estructura del modelo, hacindolo ms sencillo de resolver, y deduciendo de su solucin, la solucin del problema original.
REDUCIR EL ESPACIO DE SOLUCIONES AUMENTAR EL ESPACIO DE SOLUCIONES
Partiendo de una solucin factible inicial, y mediante alteraciones de esa solucin, pasan de forma iterativa a otras soluciones factibles de su entorno.
CRITERIO DE PARADA
(NEIGHBORHOOD)
En cada iteracin se pasa de una solucin actual a una de su entorno que sea mejor que ella,
Solucin inicial
ptimo local
Xi-2 Xi-1
Xi
Xi+1
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
Metaheursticas
Recocido Simulado
MAESTRIA PYGE
Recocer:
Ing. Rodrigo Semprtegui lvarez
Para cada temperatura durante el proceso de recocido, el slido puede alcanzar el equilibrio trmico solo si el enfriamiento se produce lentamente.
Recocido Simulado
MAESTRIA PYGE
La evolucin de un slido en el bao trmico puede
ser simulado mediante el algoritmo de la Metrpolis, basado en las tcnicas Monte Carlo. Realiza en paso de un estado a otro segn: Si el estado generado posee una energa menor que el estado que actualmente se tiene, se acepta el estado generado como actual. Caso contrario el estado se aceptar con una determinada probabilidad, la cual est en funcin de:
Algoritmos genticos
MAESTRIA PYGE
en un amplio nmero de
Algoritmos genticos
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
En los algoritmos biolgicos, la informacin hereditaria es pasada a travs de los cromosomas o genes, los cuales a su vez estn formados de un determinado nmero de valores (alelos).
Los organismos pueden agruparse formando poblaciones, los que mejor se adaptan tienen mayor probabilidad de sobrevivir y reproducirse.
Algoritmos genticos
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
Los alelos pueden representar valores de las variables de decisin, que se correspondern con los genes.
Los cromosomas representan las soluciones. Los algoritmos genticos, trabajan sobre una poblacin de soluciones generando una nueva en cada iteracin.
Bsqueda Tab
Los orgenes de la bsqueda tab se ubican a
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
heurstica para tratar de resolver problemas de cubierta no lineal, aunque varios de sus principios tambin fueron delineados independientemente por P. Hansen .
Bsqueda Tab
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
diseadas para permitir una mejor explotacin de los criterios de evaluacin y la informacin histrica de la bsqueda que se conseguira con estructuras rgidas de memoria o con sistemas carentes de memoria condiciones que limitan y hacen ms flexible el proceso de bsqueda. Este mecanismo se encuentra inmerso en la tcnica en la forma de restricciones y criterios de aspiracin plazo), para implementar estrategias que intensifiquen y diversifiquen la bsqueda. Las estrategias de intensificacin refuerzan las propiedades de las combinaciones de movimientos que han demostrado ser buenas, mientras que las estrategias de diversificacin dirigen la bsqueda hacia nuevas regiones del espacio de soluciones factibles.
GRASP
Una entre las metaheursticas ms exitosas que
MAESTRIA PYGE
problemas difciles en optimizacin combinatoria. En su versin bsica cada iteracin consiste en dos fases: una fase constructiva cuyo producto es una solucin factible buena, aunque no necesariamente un ptimo local, y una bsqueda local, durante la cual se examinan vecindades de la solucin, al llegar a un ptimo local la iteracin termina. encontrada en cada una de ellas, hasta que se alcanza un criterio de terminacin.
Redes Neuronales
Una de las misiones en una red neuronal consiste en
MAESTRIA PYGE
simular las propiedades observadas en los sistemas neuronales biolgicos a travs de modelos matemticos recreados mediante mecanismos artificiales
Una red neuronal se compone de unidades llamadas neuronas. Cada neurona recibe una serie de entradas a travs de interconexiones y emite una salida. Esta salida viene dadas por 3 enlaces:
Enlace Sinptico Enlace de Activacin Enlace de Transferencia
Redes Neuronales
Enlace Sinptico: que por lo general consiste en la
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
sumatoria de cada entrada multiplicada por el peso de su interconexin (valor neto). Si el peso es positivo, la conexin se denomina excitatoria; si es negativo, se denomina inhibitoria. Puede no existir, siendo en este caso la salida la misma funcin de propagacin. devuelto por la funcin de activacin. Se utiliza para acotar la salida de la neurona y generalmente viene dada por la interpretacin que queramos darle a dichas salidas
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
costos en el control de inventarios. FINANZAS: Costos de oportunidad o intereses por concepto de la inversin en inventarios. SISTEMAS DE INFORMACIN: Desarrollo y mantenimiento de sistemas para la administracin de inventarios. MARKETING Y VENTAS: Se depende de los inventarios disponibles para atender a los clientes. OPERACIONES: Nos permite ejercer un control total del inventario en la empresa.
MAESTRIA PYGE
CONCEPTOS DE INVENTARIO
MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
se tiene debido a que las partes o bienes terminados que se recibe, es mayor que el volumen de los mismos que se distribuye.
la entrega y mejora la puntualidad en el reparto de mercanca; reduce la posibilidad de que se presenten faltantes o/y ordenes atrasadas. COSTO DE HACER PEDIDOS: Generalmente el costo de hacer pedidos es el mismo, independientemente del tamao del pedido, influyen factores como: 1.Tiempo para seleccionar un proveedor y negociar las condiciones de la operacin. 2.Tiempo para preparar la documentacin. 3.Tiempo para realizar el seguimiento y recibir la mercanca solicitada.
MAESTRIA PYGE
reajustar una mquina para que fabrique un componente o artculo diferente del que estuvo fabricando anteriormente. Esto incluye la mano de obra y el tiempo requerido para efectuar las modificaciones, la limpieza y la instalacin de nuevas herramientas o aparatos. Adems se debe considerar los costos de material desperdiciado. El costo de preparacin es tambin independiente del tamao del pedido.
MAESTRIA PYGE
de salida de la planta puede reducir aumentando los niveles de inventario. Tener inventario disponible permite realizar ms embarques con cargas completas y minimizar la necesidad de acelerar los embarques utilizando otras modalidades de transporte ms costosas. PAGOS A PROVEEDORES: Se puede reducir el total de los pagos a los proveedores obteniendo descuentos unitarios por volmenes o cantidades grandes de pedidos.
MAESTRIA PYGE
MAESTRIA PYGE
MAESTRIA PYGE
inventario, las empresas tienen que conseguir un prstamo o perder la oportunidad de hacer un a inversin que nos puede devolver un rdito atractivo.
COSTOS DE ALMACENAMIENTO Y MANEJO:
El inventario requiere espacio y tiene que ser transportado para entrar o salir del almacn. Ese espacio tiene un costo, al igual que el personal y equipo (ejemplo montacargas), utilizados para manejar el inventario.
IMPUESTOS, SEGUROS Y MERMAS: Al final del ao, se pagan
ms impuestos cuando los inventarios son ms altos, y el seguro sobre los activos es ms caro cuando los elementos por asegurar son ms numerosos.
por:
1.
2. 3.
Robo o sustraccin de elementos del inventario por clientes o empleados Obsolescencia o caducidad del producto Deterioro por desperdicios o por daos fsicos (golpes accidentales).
TIPOS DE INVENTARIO
Los tipos de inventario no pueden identificarse por sus rasgos fsicos, sin embargo en trminos conceptuales cada uno tiene una gestacin diferente. Los tipos son: Inventario de Ciclo Inventario de Seguridad Inventario de Previsin Inventario en Trnsito
MAESTRIA PYGE
TIPOS DE INVENTARIO
INVENTARIO DE CICLO: La porcin del inventario
total que vara en forma directamente proporcional al tamao del lote se conoce como Inventario de Ciclo o Inventario del Ciclo. El tiempo transcurrido entre los pedidos se lo conoce como CICLO. La cantidad de inventario pedido o comprado, se lo conoce como TAMAO DE LOTE (Q). El tamao del lote, Q, vara en forma directamente proporcional al tiempo transcurrido (o ciclo) entre los pedidos. Cuanto ms tiempo transcurra entre dos pedidos sucesivos de un artculo determinado, tanto mayor tendr que ser el inventario de ciclo.
MAESTRIA PYGE
TIPOS DE INVENTARIO
Al inicio del ciclo el inventario es el mximo es decir Q Al final se encuentra en su mnimo valor, es decir 0 (cero) Por lo tanto, el inventario promedio de ciclo ser: Inventario Promedio de Ciclo = (Q+0)/2 = Q/2 Esta frmula es exacta, cuando la tasa de demanda es constante y uniforme. Sin embargo cuando la demanda por ciclo no es contante, proporciona una estimacin razonable y satisfactoria.
MAESTRIA PYGE
TIPOS DE INVENTARIO
MAESTRIA PYGE
INVENTARIO DE SEGURIDAD: Para evitar problemas en el servicio al cliente, se suele mantener un inventario extra, llamado inventario de seguridad, que es una proteccin contra la incertidumbre de la demanda, las cantidades suministradas y los tiempos de entrega. El inventario de seguridad garantiza que las operaciones no se interrumpan cuando esos problemas se presenten.
TIPOS DE INVENTARIO
INVENTARIO DE PREVISIN: Inventario que se
utiliza en las empresas para absorber las irregularidades que se presentan a menudo en la tasa de demanda o en el suministro. ejemplo. Los fabricantes de aparatos de aire acondicionado, suelen recibir hasta el 90% de su demanda anual durante solo tres meses al ao. Esa irregularidad en la demanda provoca que un fabricante acumule un inventario de previsin en los perodos de baja demanda, a fin de no tener que incrementar demasiado sus niveles de produccin cuando la demanda alcance sus puntos mximos.
MAESTRIA PYGE
TIPOS DE INVENTARIO
INVENTARIO EN TRNSITO: El inventario que se
mueve de un punto a otro recibe el nombre de inventario en trnsito. El inventario en trnsito esta constituido por los pedidos que los clientes han hecho, pero que todava no han sido entregado. El inventario en trnsito entre dos puntos puede medirse como la demanda promedio durante el tiempo de entrega, que es la demanda promedio del artculo por perodo (d) multiplicada por el nmero de perodos comprendidos dentro del tiempo de entrega del artculo (L), para trasladarse entre los dos puntos:
MAESTRIA PYGE
Inv. en Trnsito = dL
TIPOS DE INVENTARIO
EJEMPLO INVENTARIO DE CICLO E INVENTARIO
EN TRNSITO: Una planta enva mensualmente llantas a un mayorista, en partidas cuya lote promedio es de 2800 llantas. La demanda promedio del mayorista es de 700 llantas por semana y el tiempo de entrega desde la planta es de 3 semanas. Cuanto inventario del ciclo e inventario en trnsito maneja el mayorista? Inv. Del Ciclo = Q/2=2800/2= 1400 llantas Inv. En trnsito= DL= dL=(700 llantas/sem.)*3 sem. = 2100 llantas
MAESTRIA PYGE
336
MAESTRIA PYGE
1.
Ing. Rodrigo Semprtegui lvarez
Reducir el tamao del lote Q; sin embargo el hecho de efectuar tales reducciones en Q, sin realizar ningn otro cambio, nos puede complicar la situacin, por lo que es necesario tambin considerar adicionalmente aspectos como: Perfeccionar los mtodos para hacer pedidos y ajustes iniciales, lo cual abaratar los costos de hacer pedidos y de preparacin, y as permitir que el lote Q, pueda efectivamente ser menor. El incremento de la repetibilidad para suprimir o eliminar la necesidad de realizar cambios o alteraciones importantes.
MAESTRIA PYGE
MAESTRIA PYGE
100
Porcentaje del valor monetario
CLASE C CLASE B
90 80 70 60 50 40 30 20 10 10 20 30 40 50 60 70 80 90 100
Porcentaje de los artculos CLASE A
MAESTRIA PYGE
MAESTRIA PYGE
El objetivo del anlisis ABC, es identificar los niveles de inventario de los artculos clase A y permitir que la gerencia los controle cuidadosamente, utilizando adecuadamente los criterios analizados anteriormente
La cantidad econmica de pedido EOQ (del ingls economic order quantity) permite determinar el mejor ciclo del nivel de inventario para un artculo, y por tanto, minimiza el total de los costos anuales de hacer pedidos y de manejo de inventarios.
MAESTRIA PYGE
CLCULO DE LA EOQ
Cuando las suposiciones de la EOQ han sido satisfechas, el inventario del ciclo se comporta como se muestra en la figura: MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
Q/2
1 ciclo
CLCULO DE LA EOQ
Un ciclo comienza con Q unidades en inventario (pedido), durante el ciclo es utilizado a una tasa constante. Se pide un nuevo lote calculando que el nuevo pedido se reciba precisamente cuando el inventario llegue a 0 (porque la demanda es constante). Ya que el inventario vara entre Q y 0, el inventario del ciclo promedio ser igual a la mitad del tamao del lote.
MAESTRIA PYGE
CLCULO DE LA EOQ
Costo total = CM+CP
Costo Anual (dolares) Costo Anual (dolares)
Costo de pedidos
CM
MAESTRIA PYGE
Costo de manejo
Tamao del lote (Q) (a) Costo anual de manejo de inventario
CP
Costo anual manejo de inventarios= (inventario del ciclo promedio Q/2)(costo de manejo unitario H) Costo anual de hacer pedidos= (nm. pedidos/ao: D/Q)(costo de hacer pedidos o de preparacin S) Costo total (C) = costo anual de manejo + costo anual de hacer pedidos
CLCULO DE LA EOQ
En el siguiente ejemplo grfico observamos que el mejor tamao de lote, o EOQ, est en el punto ms bajo de la curva total (entre 50 y 100 unidades) MAESTRIA PYGE
Ing. Rodrigo Semprtegui lvarez
Costo actual
3000
50
100
150
200
250
300
350
400
El mejor Q (EOQ)
Actual Q
CLCULO DE LA EOQ
Esta frmula la obtenemos de las siguientes formas: - Con la derivada de la funcin del costo total con respecto a Q. - Igualando las frmulas correspondientes al costo anual de hacer pedidos y el costo anual de manejo de inventario y resolviendo para Q.
MAESTRIA PYGE
El tiempo entre pedidos TBO (time between orders) para un tamao de lote es el tiempo promedio que transcurre entre la recepcin o solicitud de dos pedidos de reabastecimiento constituidos por Q unidades.
- Un anlisis de sensibilidad de la EOQ, permite modificar sistemticamente los parmetros de importancia crucial a fin de determinar los efectos del cambio: - Un cambio en la tasa de demanda. A medida que aumenta la demanda, el tamao de lote tambin debe aumentar, pero ms lentamente que la demanda real.
S, aumenta la EOQ, y en consecuencia, aumenta el inventario del ciclo promedio. - Un cambio en los costos de manejo de inventario. La EOQ disminuye a medida que H aumenta y viceversa. - Errores en la estimacin de D, H y S. El costo total es muy poco sensible a los errores, aun en el caso de que las estimaciones sean errneas por un amplio margen.
MAESTRIA PYGE
Los sistemas de control de inventarios especficamente son: - Sistema de revisin continua (Q) - Sistema de revisin peridica (P) - Sistemas hbridos:
-
El punto de reorden es igual a la demanda durante el tiempo de entrega. R En la figura anterior la O con pendiente Tiempo lnea descendente es el inventario disponible. Cuando aquel llega a R, se presenta un nuevo pedido. El inventario sigue descendiendo en el perodo de entrega L. El nuevo pedido llega justo cuando el inventario desciende a 0.
IP IP IP
Pedido recibido Pedido recibido Pedido recibido Pedido recibido
MAESTRIA PYGE
Inventario disponible
OH
OH
OH
Pedido prestado
Pedido prestado
Pedido prestado
TBO
TBO
TBO
MAESTRIA PYGE
IP
Pedido recibido Pedido recibido
IP
Pedido recibido
Pedido recibido
R
Pedido prestado Pedido prestado Pedido prestado
O
L1 TBO 1 TBO 2 L2 TBO 3 L
3
Tiempo
Cuando la demanda es incierta: Punto de reorden = demanda promedio durante el perodo de entrega + inventario de seguridad.
Inventario disponible
La lnea ondulada con pendiente descendente indica que la demanda vara de un da a otro. La tasa de demanda cambiante denota que el tiempo entre perodos es variable, de modo que TBO1 TBO2 TBO3
Debido a la incertidumbre de la demanda, las ventas en el perodo de entrega son imprevisibles por lo que se aade un inventario de seguridad como medida de proteccin contra posibles perdidas de ventas. Cuanto ms grande sea el inventario de seguridad, y por ende ms alto el punto de reorden, tanto menos probable ser que se presenten faltantes.
zs
Fig. Cmo encontrar el inventario de seguridad con una distribucin de probabilidad normal para un ciclo del nivel de servicio del 85%
Conocido a veces como sistema de reorden a intervalos fijos o sistema de reorden peridico, e el cual la posicin de inventario de un artculo se revisa peridicamente y no en forma continua. Los nuevos pedidos se colocan siempre al final de cada revisin y el TBO tiene un valor fijo de P. La demanda total entre revisiones es variable.
- Que no existan restricciones en el tamao del lote. - Que los costos pertinentes sean los de manejo de
inventario y pedidos. - Que las decisiones de un artculo sean independientes de las decisiones del resto. - Que no exista incertidumbre en los tiempos de entrega ni en el suministro.
MAESTRIA PYGE
Inventario de seguridad Se calcula de forma similar al sistema Q. Se multiplica las desviaciones estndar deseadas, z, por la desviacin estndar de la demanda en el curso de intervalo de proteccin, P+L Inventario de seguridad = z P+L