• Categories
• Programación lineal
• Optimización matemática
• Matemática
• Ciencia
• Informática y tecnología de la información

Citation preview

UNIVERSIDAD TECNOLÓGICA DEL PERÚ Vicerrectorado de Investigación

OPTIMIZACIÓN DE SISTEMAS I TINS Básicos INGENIERÍA INDUSTRIAL, INGENIERÍA DE SISTEMAS

TEXTOS DE INSTRUCCIÓN BÁSICOS (TINS) / UTP

Lima - Perú

OPTIMIZACIÓN DE SISTEMAS I

© OPTIMIZACIÓN DE SISTEMAS I Desarrollo y Edición:

Vicerrectorado de Investigación

Elaboración del TINS:

• Ing. César Canelo Sotelo • Ing. Luis Medina Aquino

Diseño y Diagramación:

Julia Saldaña Balandra

Soporte académico:

Instituto de Investigación

Producción:

Imprenta Grupo IDAT

Queda prohibida cualquier forma de reproducción, venta, comunicación pública y transformación de esta obra.

2

OPTIMIZACIÓN DE SISTEMAS I

“El presente material contiene una compilación de obras de Optimización de Sistemas publicadas lícitamente, resúmenes de los temas a cargo del profesor; constituye un material auxiliar de enseñanza para ser empleado en el desarrollo de las clases en nuestra institución. Éste material es de uso exclusivo de los alumnos y docentes de la Universidad Tecnológica del Perú, preparado para fines didácticos en aplicación del Artículo 41 inc. C y el Art. 43 inc. A., del Decreto Legislativo 822, Ley sobre Derechos de Autor”.

3

OPTIMIZACIÓN DE SISTEMAS I

4

OPTIMIZACIÓN DE SISTEMAS I

PRESENTACIÓN

El asunto del presente texto está engastado en el espacio ilimitado de la matemática aplicada, surgida en la segunda mitad del siglo XX, como producto de la ciencia y la tecnología desarrollada en el fragor de la segunda guerra mundial. Producto que tomó el nombre de Investigación de Operaciones, compuesto por los temas de Programación Lineal, Programación Dinámica, Programación Geométrica, entre otros. Técnicas que han venido aplicándose intensamente como soporte de la maximización de indicadores concurrentes a la excelencia de calidad con altos niveles de productividad y bajos costos de producción. En este contexto, la educación superior universitaria ha integrado en los currícula de diferentes Carreras de Ingeniería algunas Asignaturas de Investigación Operativa, con enfoque metodológico generado en el espacio de la Teoría General de Sistemas. Con tal método, mediante un acucioso trabajo académico de recopilación y selección, ha venido en prepararse para estudios de la Carrera de Ingeniería de Sistemas el Curso de Optimización de Sistemas I, aplicado en el VI ciclo de estudios, con temas de Programación Lineal. El acopio y selección de materias pertinentes a la Asignatura, arriba mencionada, ha sido realizada por el profesor Ing. Luis Medina Aquino , en unión con el profesor Ing. César Canelo Sotelo, en congruencia al sillabus correspondiente, en el nivel de calidad académica requerido; comprende los siguientes temas: El capítulo 1 trata sobre la descripción del modelo de programación lineal, que consta de variables de decisión, función objetivo, restricciones y condición de no negatividad. El objetivo del modelo es maximizar (utilidades) o minimizar (costos) de una función lineal Z. A través de un problema que tiene dos variables de decisión se explica la formulación del programa lineal y cómo hallar una solución gráfica. El capítulo 2 complementa el capítulo uno con la construcción de modelos de programación lineal formulados en base a diferentes problemas que se presentan en la industria. Existen diferentes problemas tipo con su respectiva formulación, con más de dos variables.

5

OPTIMIZACIÓN DE SISTEMAS I

El capítulo 3 explica el método general para hallar una o varias soluciones de un modelo de programación lineal. Para ello es necesario colocar las restricciones del modo estándar, convirtiendo las desigualdades en igualdades, añadiendo variables de holgura, de exceso y/o artificiales. Luego se coloca todos los valores en una tabla y se empieza a realizar iteraciones hasta hallar el valor óptimo de la función objetivo. En el capítulo 4 se explica la dualidad de un programa lineal. Si la función objetivo original (primal) es de maximización entonces la función objetivo del programa dual será de minimización y viceversa. La cantidad de variables de decisión en el programa primal será la misma cantidad de restricciones del dual. El número de restricciones del primal será la misma del número total de variables duales. Los valores del lado derecho del programa primal serán los coeficientes de la función objetivo del dual. Y los coeficientes de las variables primales serán los mismos valores del lado derecho de las restricciones del programa dual. El capítulo 5 trata acerca del análisis de sensibilidad que se puede hacer al modelo de programación lineal. El análisis de sensibilidad responde a la pregunta: ¿Qué pasa con el valor de la función objetivo si existen cambios en los coeficientes de la función objetivo o el lado derecho de las restricciones? ¿Se mantiene o no la solución óptima si existen cambios en los coeficientes de las variables de la función objetivo? ¿Si se agrega una nueva variable a la función objetivo y a las restricciones, formará parte de la nueva solución? ¿Si se agrega una nueva restricción, se modificará la solución óptima? En el capítulo 6 se aborda el modelo de programación lineal entera. La programación lineal resuelta por el método simplex nos da un resultado con valores continuos, pero existen problemas de programación lineal que exigen una solución con valores enteros. En realidad, para este tipo de problemas se añaden restricciones de variable entera a través de un algoritmo de ramificación que va restringiendo las variables a soluciones enteras. En este campo se incluye problemas de programación lineal cuyas variables de decisión son valores binarios {0, 1}. El contenido del texto representa el esfuerzo académico de los profesores Ing. Luis Medina Aquino e Ing. César Canelo Sotelo, a quienes la Institución agradece de manera especial en el camino de contribución, tendiente a la mejora continua de la calidad académica. LUCIO HERACLIO HUAMÁN URETA VICERRECTOR DE INVESTIGACIÓN

6

OPTIMIZACIÓN DE SISTEMAS I

ÍNDICE GENERAL

CAPÍTULO 1 PROGRAMACION LINEAL: EL METODO GRAFICO......................

11

CAPÍTULO 2 PROGRAMACION LINEAL: FORMULACIÓN DE PROBLEMAS.....

43

CAPÍTULO 3 EL MÉTODO SIMPLEX....................................................................

67

CAPÍTULO 4 EL PROBLEMA DUAL .....................................................................

83

CAPÍTULO 5 ANALISIS DE SENSIBILIDAD .........................................................

97

CAPÍTULO 6 PROGRAMACION LINEAL ENTERA...............................................

109

BIBLIOGRAFÍA ...............................................................................

135

7

OPTIMIZACIÓN DE SISTEMAS I

8

OPTIMIZACIÓN DE SISTEMAS I

DISTRIBUCIÓN TEMÁTICA Clase N°

1

2

Tema

Introducción a la Investigación de Operaciones: Origen, Definición, Modelo, Tipos de modelo. Metodología de la Investigación de Operaciones. Programación Lineal: Definición, Presentación del modelo de P.L., Suposiciones del modelo de P.L., Interpretación económica del modelo de PL., Propiedades del modelo de PL., Formas de mostrar el modelo de PL., Variable de holgura, Transformaciones en el modelo de PL.

Semana

Horas

1

04

2

04

3

Formulación de problemas con P.L.

3

04

4

Solución del modelo P.L.: Región Factible y Soluciones, Conjuntos convexos. Métodos de solución de modelos lineales: Gráfico y Simplex.

4

04

5

Casos especiales de soluciones de PL usando el método gráfico.

5

04

6

Método simplex: Teoremas. Algoritmo simplex primal.

6

04

7

Casos especiales usando el algoritmo simplex: Solución única, Solución con región factible no acotada, Problema no factible, Soluciones Múltiples. Solución degenerada. Prevención del ciclado

7

04

8

Técnica de las variables artificiales. Método de penalización. Método de las dos fases.

8

04

9

Revisión – Nivelación

9

04

10

EXAMEN

10

02

PARCIAL

9

OPTIMIZACIÓN DE SISTEMAS I

Clase N°

Tema

Semana

Horas

11

Simplex revisado: Inversión explícita de la base, Método de las dos fases.

11

04

12

Problema dual: El programa Dual, Relaciones Primal – Dual. Teoremas. Propiedades Primaldual.

12

04

13

Algoritmo Simplex Dual. Interpretación Económica del problema Dual.

13

04

14

Análisis de sensibilidad. Rangos de sensibilidad: Cambios en el vector de costos y recursos. Cambios en la matriz de coeficientes tecnológicos.

14

04

15

Análisis de sensibilidad: Adición de una nueva variable. Adición de una nueva restricción.

15

04

16

Uso de software computacional para PL; Linear Interactive Discrete Optimizer (LINDO).

16

04

17

Programación entera: Formulación.

17

04

18

Métodos de solución de problemas de Programación Lineal Entera: Método gráfico, Método de ramificar y acotar.

18

04

19

02

19

10

EXAMEN FINAL

OPTIMIZACIÓN DE SISTEMAS I

CAPÍTULO 1 PROGRAMACION LINEAL: EL METODO GRAFICO

1.1. INTRODUCCIÓN Existen problemas de decisión administrativos que pueden ser resueltos a través de un modelo matemático llamado programación lineal. Por ejemplo el fabricante desea elaborar un programa de producción de costo mínimo; exigido por la demanda a atender y limitado por su capacidad de producción. Un modelo de programación lineal busca el objetivo de maximizar (minimizar) una función lineal, sujeta a un conjunto de restricciones lineales.1 Un modelo de programación lineal esta compuesto de lo siguiente: 1. Un conjunto de variables de decisión 2. Una función objetivo 3. Un conjunto de restricciones Para formular un modelo de programación lineal primero se debe entender el problema y responder a las siguientes preguntas: ¿Cuál es nuestro objetivo económico? ¿Maximizar utilidades o minimizar costos? ¿Qué limitaciones de recursos existen? ¿Qué requerimientos mínimos se necesitan? Con esto podemos: (i) Identificar las variables de decisión del modelo, a las que llamaremos X1, X2, X3,...., Xn. (ii) Expresar la función objetivo como: Maximizar (o Minimizar) Z = C1 X1 + C2 X2 + C3 X3 +....+ Cn Xn (iii) Determinar las restricciones del modelo que son funciones lineales de las variables de decisión. Estas restricciones pueden ser igualdades (=) o desigualdades de la forma (>, bi ò ai1 X1 + ai2 X2 + ai3 X3 +....+ ain Xn = bi

11

OPTIMIZACIÓN DE SISTEMAS I

También incluye restricciones de no negatividad de las variables: X1, X2, X3,...., Xn > 0. 1.2.

FORMULACIÓN Y SOLUCIÓN GRÁFICA DE UN EJERCICIO DE PROGRAMACION LINEAL La solución gráfica se emplea para resolver modelos de programación lineal con dos variables, ya que resulta bastante difícil dibujar planos de tres variables, e imposible hacerlo para cuatro o más variables. Ejercicio Juan es un prospero negociante que se dedica a la compra y venta de naranja y papaya. Él tiene su cartera de clientes que son aquellos comerciantes que tienen su puesto de frutas en los diferentes mercados del distrito de Jesús María. Todos los días temprano en la mañana visita a su proveedor de frutas en el mercado mayorista y hace las compras del día. El día anterior recibe los pedidos de sus clientes y esta suma 600 kilos de papaya y 1200 kilos de naranja. Juan lleva su camión para el transporte cuya capacidad de carga es de 1600 kilos. Entonces ¿Cuántos kilos de cada fruta debe comprar Juan para maximizar los beneficios? Para resolver esta pregunta se tienen los siguientes precios y costos por kilo de fruta:

Papaya Naranja

Precio de compra Precio de venta al por mayor al minorista S/. 1.30 S/. 1.60 S/. 1.00 S/. 1.20

Procedimiento de Solución (Método Gráfico) 1) Establecer la formulación del problema 2) Graficar en el plano cartesiano las restricciones del tipo >, < ó =, como si fueran rectas. 3) Ubicar el espacio de la solución factible (región factible), el cual está dado por el área común a todas las restricciones. 4) Obtener la solución óptima.

12

OPTIMIZACIÓN DE SISTEMAS I

Formulación del Ejercicio Definición de las Variables de Decisión X1 = Cantidad, en kilos, de papaya que se debe comprar. X2 = Cantidad, en kilos, de naranja que se debe comprar. Función Objetivo Maximizar la utilidad total de los dos productos Restricciones Cantidad máxima de Papaya < 600 kilos. Cantidad máxima de Naranja < 1200 kilos. Carga máxima del camión < 1600 kilos. Condición de No Negatividad X1, X2 > 0 El Modelo Maximizar Z = 0.30 X1 + 0.20 X2 (Beneficio Total) Sujeto a: R1 X1 < 600 (Cantidad máxima de Papaya) R2 X2 < 1200 (Cantidad máxima de Naranja) R3 X1 + X2 < 1600 (Carga máxima del camión) X1, X2 > 0 (Condición de no negatividad) Gráfica en el Plano Cartesiano Primero graficar la igualdad de la restricción, luego escoger un punto de ensayo (por ejemplo el punto 0,0) y se sustituye este punto en la desigualdad para comprobar si cumple esta restricción. Si lo cumple entonces sombrear el área que cubre este punto de ensayo y si no lo cumple sombrear el área que no lo cubre. X2

X2

(0,1200)

R2

X2 < 1200

R1 X1 < 600

(0,0)

(600,0)

X1

(0,0)

(600,0)

X1

Gráfica 1

13

OPTIMIZACIÓN DE SISTEMAS I

En el plano cartesiano de la izquierda de la gráfica 1, primero se dibuja la recta X1 = 600 y luego se escoge un punto de ensayo, para nuestro caso (0,0), y se verifica que cumpla con la desigualdad, por tanto se sombrea los puntos que cumplen con todos los puntos de X1 que sean menores o iguales a 600 Kg. (restricción R1). Con el mismo procedimiento se dibuja la recta X2=1200, que se muestra en el plano cartesiano de la derecha de la gráfica 1, y se acota más el área de los puntos factibles con los puntos de X2 que sean menores o iguales a 1200 Kg. (restricción R2). Por último graficamos la tercera restricción el cual restringe aún más el área de puntos factibles, como se muestra la gráfica 2. Región factible es el conjunto de puntos que satisface todas las restricciones

simultáneamente.

Existen

infinitos

puntos

factibles

(soluciones). Se llaman puntos extremos a los vértices de la región de factibilidad.

14

OPTIMIZACIÓN DE SISTEMAS I

X2

Los valores que optimizan la función objetivo siempre se encuentran en uno de los

R1

(0,1600)

puntos extremos, en este caso A, B, C, D ó E.

(0,1200)

(400,1200) B

R2

C D

(600,1000)

R3

A (0,0)

(1600,0)

E (600,0)

X1

Gráfica 2

Maximizar Z = 0.30 X1 + 0.20 X2 en la región factible

Se debe dibujar el contorno de la función objetivo (línea iso-beneficio) mediante rectas paralelas, en cada vértice, según la relación: X2 = – 1.5 X1 + K, ver gráfica 3.

15

OPTIMIZACIÓN DE SISTEMAS I

Pendiente de la función objetivo

0.30

Gráfica 3

X2

0.20 Z1 = 0.30 Z3

Z2 = 0.30 (600) + 0.20

(0,1200) Z2

(0) + 0.20

(400,1200) B C D

Z3 = 0.30

(0) = 0 (0) = 180

(0) + 0.20 (1200) = 240

Z4 = 0.30 (400) + 0.20 (1200) = 360 (600,1000) Z5

Z5 = 0.30 (600) + 0.20 (1000) = 380 (óptimo) Se debe adquirir 600 Kg. de papaya

Z1

Z4

y 1000 Kg. de naranja para obtener S/. 380 de utilidad

A (0,0)

1.3.

E (600,0)

X1

PROBLEMAS RESUELTOS: FORMULADOS Y CON SOLUCION GRAFICA DE PROGRAMACION LINEAL, UTILIZANDO WINQSB2.

1.3.1 CREDIFONDO, una empresa que administra fondos mutuos, tiene $50,000 de un fondo de pensiones, y desea invertir en bonos tipo A y bonos tipo B que producen una rentabilidad de 6% y 10% anual respectivamente. Por motivos de liquidez no puede invertir más del 25% en bonos tipo A, y lo mínimo a depositar en bonos tipo B es $10,000. Determinar un plan óptimo de inversiones VARIABLES DE DECISIÓN: X1 = Cantidad, en dólares, que se debe invertir en bonos tipo A. X2 = Cantidad, en dólares, que se debe invertir en bonos tipo B. 2 WINQSB es un software creado por el Dr. Yih-Long Chang que se puede bajar de internet. Buscando en google hallamos varias direcciones, una de ellas es http://www.investigacionoperaciones.com/Metodos_computacionales.htm

16

OPTIMIZACIÓN DE SISTEMAS I

FUNCION OBJETIVO: Se debe maximizar la rentabilidad total de la inversión en los dos tipos de bonos. Maximizar Z = 0.06 X1 + 0.10 X2 RESTRICCIONES: R1 = Fondo máximo a depositar: X1 + X2 ≤ 50,000 R2 = Máximo a invertir en bonos tipo A: X1 ≤ 0.25 (X1 + X2) 0.75 X1 – 0.25 X2 ≤ 0 R3 = Mínimo a invertir en bonos tipo B: X2 ≥ 10,000 CONDICION DE NO NEGATIVIDAD: X1 ≥ 0, X2 ≥ 0 SOLUCION CON WINQSB:

1.3.2. INTERBANK tiene un total de $20 millones asignados a préstamos para adquisición de casas y automóviles. En promedio, los préstamos hipotecarios tienen una tasa anual de recuperación del 10%, y los préstamos para autos una tasa anual de recuperación del 12%. La gerencia ha estipulado que la cantidad total de préstamos hipotecarios debe ser mayor o igual cuatro veces la cantidad total de préstamos para autos. Determine la 17

OPTIMIZACIÓN DE SISTEMAS I

cantidad total de los préstamos de cada tipo que debe realizar INTERBANK para maximizar el monto de recuperación. VARIABLES DE DECISIÓN: X1= Cantidad, en dólares, que se debe asignar para créditos hipotecarios X2 = Cantidad, en dólares, que se debe asignar para créditos de autos. FUNCION OBJETIVO: Se debe maximizar la recuperación total de los préstamos Maximizar Z = 0.10 X1 + 0.12 X2 RESTRICCIONES: R1 = Fondo máximo para asignar créditos: X1 + X2 ≤ 20,000,000 R2 = Relación de préstamos X1 ≥ 4 X2 X1 – 4 X2 ≥ 0 CONDICION DE NO NEGATIVIDAD: X1 ≥ 0, X2 ≥ 0 SOLUCION CON WINQSB:

18

OPTIMIZACIÓN DE SISTEMAS I

1.3.3. MAQUINASA es una pequeña fábrica situada en los alrededores de una gran ciudad. Su producción se limita a dos productos industriales: Alfa y Beta. El departamento de contabilidad de la empresa ha calculado las contribuciones de cada producto: 10 dólares para el producto Alfa y 12 dólares para el Beta. Cada producto pasa por tres departamentos de la fábrica. Los requerimientos de tiempo para cada producto y el total del tiempo disponible en cada departamento son los siguientes: Horas Requeridas Horas Producto Producto Disponibles Departamento ALFA BETA este mes Determine la cantidad 1 2.0 3.0 1,500 de productos Alfa y Beta 2 3.0 2.0 1,500 de tal forma que maximice 3 1.0 1.0 600 la contribución total. VARIABLES DE DECISIÓN: X1 = Cantidad, en unidades, del producto Alfa que se debe producir por mes. X2 = Cantidad, en unidades, del producto Beta que se debe producir por mes. FUNCION OBJETIVO: Se debe maximizar la utilidad total de los dos productos Maximizar Z = 10 X1 + 12 X2 RESTRICCIONES: R1 = Horas disponibles del Departamento 1: 2X1 + 3X2 ≤ 1500 R2 = Horas disponibles del Departamento 2: 3X1 + 2X2 ≤ 1500 R3 = Horas disponibles del Departamento 3: 1X1 + 1X2 ≤ 600 CONDICION DE NO NEGATIVIDAD: X1 ≥ 0, X2 ≥ 0

19

OPTIMIZACIÓN DE SISTEMAS I

SOLUCION CON WINQSB:

1.3.4. Una dietista del hospital Rebagliati es responsable de la planeación y administración de los requerimientos alimenticios de los pacientes. La especialista examina en estos momentos el caso de un paciente que se le ha restringido a una dieta especial que consta de dos fuentes alimenticias. Al paciente no se le ha restringido la cantidad de los dos elementos que se puede consumir; sin embargo, se deben satisfacer los siguientes requerimientos nutritivos mínimos por día: - 1000 unidades del nutriente A. - 2000 unidades del nutriente B; y - 1500 unidades del nutriente C. Cada onza de la fuente alimenticia #1, contiene 100 unidades del nutriente A, 400 unidades del nutriente B y 200 unidades del nutriente C.

20

OPTIMIZACIÓN DE SISTEMAS I

Cada onza de la fuente alimenticia #2, contiene 200 unidades del nutriente A, 250 unidades del nutriente B y 200 unidades del nutriente C. Ambas fuentes alimenticias son algo costosas: La fuente alimenticia #1 cuesta $6 por libra y la fuente #2 $8 por libra. La dietista desea determinar la combinación de fuentes alimenticias que arroje el menor costo y que satisfaga todos los requerimientos nutritivos. Nota: 1 libra = 16 onzas VARIABLES DE DECISIÓN: X1 = Cantidad, en onzas, de la fuente alimenticia #1 que se debe asignar a la dieta por día. X2 = Cantidad, en onzas, de la fuente alimenticia #2 que se debe asignar a la dieta por día. FUNCION OBJETIVO: Se debe minimizar el costo total de la dieta. Minimizar Z = 6/16 X1 + 8/16 X2 = 0.375 X1 + 0.5 X2 RESTRICCIONES: R1 = Cantidad mínima de nutriente A: 100X1 + 200X2 ≥ 1000 R2 = Cantidad mínima de nutriente B: 400X1 + 250X2 ≥ 2000 R3 = Cantidad mínima de nutriente C: 200X1 + 200X2 ≥ 1500

CONDICION DE NO NEGATIVIDAD: X1 ≥ 0, X2 ≥ 0

21

OPTIMIZACIÓN DE SISTEMAS I

SOLUCION CON WINQSB:

1.3.5. La fábrica ABC vende dos tipos de bombas hidráulicas: (1) normal y (2) extra grande. El proceso de manufactura asociado con la fabricación de las bombas implica tres procesos: ensamblado, pintura y pruebas de control de calidad. Los requerimientos de recursos para ensamble, pintura y prueba de las bombas se muestran en la siguiente tabla: Tabla de Requerimientos de Manufactura Tiempo de Tiempo de Tiempo de Tipo Ensamble Pintado Prueba Normal 3.6 1.6 0.6 Extra Grande 4.8 1.8 0.6 La contribución a las utilidades por la venta de una bomba normal es $50, en tanto que la utilidad por una bomba extra grande es $75. Existen disponibles por semana 4,800 horas en tiempo de ensamble, 1,980 horas en tiempo de pintura y 900 horas en tiempo de prueba. Las experiencias anteriores de renta señalan que la compañía puede esperar vender cuando menos 300 bombas normales y 180 de los extra grandes por semana. A la fábrica ABC le gustaría determinar la cantidad de cada tipo de bomba que debe fabricar semanalmente con el objeto de maximizar sus utilidades. 22

OPTIMIZACIÓN DE SISTEMAS I

VARIABLES DE DECISIÓN: X1 = Cantidad, en unidades, de bombas hidráulicas normales que se debe producir por semana X2 = Cantidad, en unidades, de bombas hidráulicas extragrandes que se debe producir por semana. FUNCION OBJETIVO: Se debe maximizar la utilidad total de los dos productos Maximizar Z = 50 X1 + 75 X2 RESTRICCIONES: R1 = Horas disponibles de ensamble: 3.6 X1 + 4.8 X2 ≤ 4800 R2 = Horas disponibles de pintado: 1.6 X1 + 1.8 X2 ≤ 1980 R3 = Horas disponibles de prueba: 0.6 X1 + 0.6 X2 ≤ 900 R4 = Demanda mínima de X1: X1 ≥ 300 R5 = Demanda mínima de X2: X2 ≥ 180 CONDICION DE NO NEGATIVIDAD: X1 ≥ 0, X2 ≥ 0 SOLUCION CON WINQSB:

23

OPTIMIZACIÓN DE SISTEMAS I

1.3.6. PARAMONGA tiene dos tipos de papel, para libros y para revistas. Cada tonelada de papel para libros requiere 2 toneladas de abeto y 3 ton. de pino. Cada tonelada de papel para revistas requiere 2 toneladas de abeto y 2 toneladas de pino. La empresa debe proveer al menos 25000 tons de papel para libros y 10000 tons de papel para revistas por año. La disponibilidad anual de materiales es de 300000 tons de abeto y 450000 de pino. Por razón de mercado la cantidad de papel fabricado para revistas debe ser al menos 1.5 veces a la cantidad de papel fabricado para libros. Cada tonelada de papel para libros da una utilidad de $215 y de revistas de $270. Determine un plan óptimo de producción VARIABLES DE DECISIÓN: X1 = Cantidad, en toneladas, de papel para libros que se debe producir por año. X2 = Cantidad, en toneladas, de papel para revistas que se debe producir por año. FUNCION OBJETIVO: Se debe maximizar la utilidad total de los dos productos Maximizar Z = 215 X1 + 270 X2 RESTRICCIONES: R1 = Disponibilidad de abeto: R2 = Disponibilidad de pino: R3 = Razón de mercado: R4 = Demanda mínima de X1: R5 = Demanda mínima de X2:

2 X1 + 2 X2 ≤ 300000 3 X1 + 2 X2 ≤ 450000 X2 ≥ 1.5 X1 1.5 X1 – X2 ≤ 0 X1 ≥ 25000 X2 ≥ 10000

CONDICION DE NO NEGATIVIDAD: X1 ≥ 0, X2 ≥ 0

24

OPTIMIZACIÓN DE SISTEMAS I

SOLUCION CON WINQSB:

1.3.7. CATS es un nuevo producto alimenticio para mascotas. Cada lata de 16 onzas de Cats es una mezcla, o combinación, de dos ingredientes alimenticios para mascotas. Sean X1 = número de onzas del ingrediente A en lata de 16 onzas. X2 = número de onzas del ingrediente B en lata de 16 onzas. Cada onza del ingrediente A contiene 1/2 onzas de proteínas y 1/8 de onza de grasas. Cada onza del ingrediente B contiene 1/10 de onza de proteínas y 1/3 de onza de grasas. Las restricciones implican que una lata de 16 onzas de Cats debe contener cuando menos 4 onzas de proteínas y no más de 2.5 onzas de grasas. Si el ingrediente A cuesta $0.04 por onza y el ingrediente B cuesta $0.03 la onza. a) b) c)

Formule el problema de programación lineal. ¿cuál es la mezcla de costo mínimo de los ingredientes A y B para cada lata de 16 onzas? Identifique e interprete los valores de las variables de excedente para este problema.

25

OPTIMIZACIÓN DE SISTEMAS I

VARIABLES DE DECISIÓN: X1 = Cantidad, en onzas, del ingrediente A en la lata de 16 onzas. X2 = Cantidad, en onzas, del ingrediente B en la lata de 16 onzas. FUNCION OBJETIVO: Se debe minimizar el costo total de los ingredientes en la lata de 16 onzas. Minimizar Z = 0.04 X1 + 0.03 X2 RESTRICCIONES: R1 = Cantidad de los ingredientes A y B en la lata de 16 oz.: X1+X2=16 R2 = Cantidad mínima de proteínas: 0.5 X1 + 0.10 X2 ≥ 4 R3 = Cantidad máxima de grasas C: 0.125 X1 + 0.333 X2 ≤ 2.5 0.375 X1 + X2 ≤ 7.5 CONDICION DE NO NEGATIVIDAD: X1 ≥ 0, X2 ≥ 0 SOLUCION CON WINQSB:

26

OPTIMIZACIÓN DE SISTEMAS I

1.3.8. Un expendio de carnes de la ciudad acostumbra preparar la carne para albondigón con una combinación de carne molida de res y carne molida de cerdo. La carne de res contiene 80% de carne y 20% de grasa, y le cuesta a la tienda 80 ctvs por libra; la carne de cerdo contiene 68% de carne y 32% de grasa, y cuesta 60 ctvs por libra. ¿Qué cantidad de cada tipo de carne debe emplear la tienda en cada libra de albondigón, si se desea minimizar el costo y mantener el contenido de grasa no mayor de 25%? VARIABLES DE DECISIÓN: X1 = Cantidad, en libras, de carne molida de res contenida en una libra de albondigón. X2 = Cantidad, en libras, de carne molida de cerdo contenida en una libra de albondigón. FUNCION OBJETIVO: Se debe minimizar el costo total de los ingredientes en una libra de albodigón. Minimizar Z = 0.80 X1 + 0.60 X2 RESTRICCIONES: R1 = Cantidad de ingredientes en una libra de albondigón: X1 + X2 = 1 R2 = Cantidad máxima de grasa 0.25 libras: 0.20 X1 + 0.32 X2 ≤ 0.25 27

OPTIMIZACIÓN DE SISTEMAS I

CONDICION DE NO NEGATIVIDAD: X1 ≥ 0, X2 ≥ 0 SOLUCION CON WINQSB:

1.3.9. Una fábrica de automóviles y camiones consta departamentos que a continuación se enumeran: 1. 2. 3. 4.

de

los

Estampado de planchas metálicas Armado de motores Montaje de automóviles Montaje de camiones

El Departamento 1 puede estampar, por mes, las planchas necesarias para 25,000 automóviles 35,000 camiones, o las correspondientes combinaciones de automóviles y camiones. El Departamento 2 puede armar, por mes, 33,333 motores de automóviles o 16,667 motores de camión, o las correspondientes combinaciones de motores de automóvil y camión. El Departamento 3 puede montar y terminar 22,500 automóviles, y el Departamento 4 puede montar 15,000 camiones. Si cada automóvil deja una utilidad de 300 dólares y cada camión de 250, ¿qué cantidades de automóviles y camiones deben producirse, de 28

OPTIMIZACIÓN DE SISTEMAS I

manera que las utilidades que se obtengan sean las máximas posibles? VARIABLES DE DECISIÓN: X1 = Cantidad, en unidades, de automóviles que se debe producir por mes. X2 = Cantidad, en unidades, de camiones que se debe producir por mes. FUNCION OBJETIVO: Se debe maximizar la utilidad total de los dos productos Maximizar Z = 300 X1 + 250 X2 RESTRICCIONES: R1 = El departamento 1 puede estampar, por mes, planchas metálicas para 25000 automóviles o 35000 camiones. Supongamos que los primeros 15 días (1/2 mes) se producen 12500 automóviles, entonces los últimos 15 días se deben producir 17500 camiones. El análisis parte del tiempo de producción para cada producto, cuyo límite máximo es un mes. De esta forma el tiempo para producir automóviles es X1/25000 (fracción de mes) y el tiempo para producir camiones es X2/35000, de tal forma que la suma de tiempos sea menor o igual a un mes: X1/25000 + X2/35000 ≤ 1 7 X1 + 5 X2 ≤ 175000 R2 = Similar análisis para el Departamento 2: X1 + 2X2 ≤ 33333 R3 = Cantidad máxima en el departamento de montaje para automóviles: X1 ≤ 22500 R4 = Cantidad máxima en el departamento de montaje para camiones: X2 ≤ 15000 CONDICION DE NO NEGATIVIDAD: X1 ≥ 0, X2 ≥ 0

29

OPTIMIZACIÓN DE SISTEMAS I

SOLUCION CON WINQSB:

1.3.10. Un fabricante de gasolina para aviación vende dos clases de combustible, A y B. El combustible de clase A tiene 25 % de gasolina grado 1, 25% de gasolina grado 2, y 50 % de gasolina grado 3. El combustible de clase B tiene 50% de gasolina grado 2 y 50% de gasolina grado 3. Disponibles para producción hay 75 galones/hora de grado 1, 150 galones/hora de grado 2, y 200 galones/hora de grado 3. Los costos son 30 centavos por galón de grado 1, 60 centavos por galón de grado 2, y 50 centavos por galón de grado 3. Las clases A y B, pueden venderse a 75 y 90 centavos por galón, respectivamente. ¿Qué cantidad por hora debe fabricarse de cada combustible? VARIABLES DE DECISIÓN: X1 = Cantidad, en galones, del combustible A que se debe producir por hora. X2 = Cantidad, en galones, del combustible B que se debe producir por hora. FUNCION OBJETIVO: Se debe maximizar la utilidad total (ingreso menos costo) de los dos productos 30

OPTIMIZACIÓN DE SISTEMAS I

Ingreso en centavos = 75 X1 + 90 X2 Costo en centavos = 30 (0.25 X1) + 60 (0.25 X1 + 0.50 X2) + 50 (0.50 X1 + 0.50 X2) Maximizar Z = Ingreso – Costo = 27.5 X1 + 35 X2 RESTRICCIONES: R1 = Gasolina grado 1 disponible: 0.25 X1 ≤ 75 R2 = Gasolina grado 2 disponible: 0.25 X1 + 0.50 X2 ≤ 150 R3 = Gasolina grado 3 disponible: 0.50 X1 + 0.50 X2 ≤ 200 CONDICION DE NO NEGATIVIDAD: X1 ≥ 0, X2 ≥ 0 SOLUCION CON WINQSB:

1.4.

PROBLEMAS FORMULADOS PERO SIN SOLUCION GRAFICA

1.4.1. Una compañía petrolera que tiene dos refinerías, necesita al menos 800, 1400 y 500 barriles de petróleo de grados bajo, medio y alto, respectivamente. Cada día, la refinería I produce 200 barriles de grado bajo, 300 de medio y 100 de alto grado, mientras que la refinería II produce 100 barriles de grado alto, 100 de bajo y 31

OPTIMIZACIÓN DE SISTEMAS I

200 de grado medio. Si los costos diarios son de $2,500 para operar la refinería I y de $2,000 para la refinería II, ¿cuántos días debe ser operada cada refinería para satisfacer los requerimientos de producción a un costo mínimo?¿cuál es el costo mínimo? VARIABLES DE DECISIÓN: X1 = Cantidad de días que debe ser operada la refinería I para cumplir con los requerimientos de producción. X2 = Cantidad de días que debe ser operada la refinería II para cumplir con los requerimientos de producción. FUNCION OBJETIVO: Se debe minimizar el costo total de operación de las dos refinerías. Minimizar Z = 2500 X1 + 2000 X2 RESTRICCIONES: R1 = Cantidad mínima de barriles de petróleo de grado bajo requerido: 200X1 + 100X2 ≥ 800 R2 = Cantidad mínima de barriles de petróleo de grado medio requerido: 300X1 +200X2 ≥ 1400 R3 = Cantidad mínima de barriles de petróleo de grado alto requerido: 100X1 + 100X2 ≥ 500 CONDICION DE NO NEGATIVIDAD: X1 ≥ 0, X2 ≥ 0 1.4.2. A causa de reglamentaciones gubernamentales nuevas sobre la contaminación, una compañía química ha introducido en sus plantas un nuevo y más caro proceso para complementar o reemplazar un proceso anterior en la producción de un químico en particular. El proceso anterior descarga 15 gramos de dióxido de azufre y 40 gramos de partículas a la atmósfera por cada litro de químico producido. El nuevo proceso descarga 5 gramos de dióxido de azufre y 40 gramos de partículas a la atmósfera por cada litro de químico producido. La compañía obtiene una utilidad de 30 y 20 centavos por litro en los procesos anterior y nuevo, respectivamente. Si el gobierno permite a la planta descargar no más de 10,500 gramos de dióxido de azufre y no más de 30,000 gramos de partículas a la atmósfera cada día, ¿cuántos litros de químico deben ser producidos diariamente, por cada uno de los procesos, para maximizar la utilidad diaria? ¿Cuál es la utilidad diaria?

32

OPTIMIZACIÓN DE SISTEMAS I

VARIABLES DE DECISIÓN: X1 = Cantidad, en litros, del producto químico que se debe producir con el proceso anterior. X2 = Cantidad, en litros, del producto químico que se debe producir con el proceso nuevo. FUNCION OBJETIVO: Se debe maximizar la utilidad total del producto químico en los dos procesos Maximizar Z = 30 X1 + 20 X2 RESTRICCIONES: R1 = Descarga máxima de dióxido de azufre: 15 X1 + 5 X2 ≤ 10500 R1 = Descarga máxima de partículas: 40 X1 + 40 X2 ≤ 30000 CONDICION DE NO NEGATIVIDAD: X1 ≥ 0, X2 ≥ 0 1.4.3. INDUMIL, un fabricante progresista de mecanismos civiles y militares, fabrica actualmente una línea de armas para civiles, con una producción actual diaria de 30 unidades del modelo Z-1200 y de 120 unidades del modelo Z-1500. El gerente de manufactura quiere saber si podrían aumentarse las ganancias cambiando la mezcla de productos entre los dos modelos. Se compiló la siguiente información sobre las horas requeridas para la fabricación de cada modelo y las capacidades de los departamentos de la fábrica.

Departamento 1 2 3 4 Contribución por unidad a) b)

Horas-Hombre Requeridas Modelo Modelo Z-1200 Z-1500 2.0 0.0 0.0 3.0 2.0 2.0 1.2 1.5 $50

Horas Disponibles por día 300 540 440 300

$40

Determínese la mezcla óptima de productos suponiendo que pueden venderse las cantidades. Use el método gráfico. ¿Cuánto aumentaría la mezcla óptima la contribución a los costos fijos y a las ganancias? 33

OPTIMIZACIÓN DE SISTEMAS I

c)

Suponga que el precio del modelo Z-1200 se reduzca a $10 ¿Cuál será la mezcla óptima de productos? Use el método gráfico

VARIABLES DE DECISIÓN: X1 = Cantidad, en unidades, del producto Z-1200 que se debe producir por día. X2 = Cantidad, en unidades, del producto Z-1500 que se debe producir por día. FUNCION OBJETIVO: Se debe maximizar la utilidad total de los dos productos Maximizar Z = 50 X1 + 40 X2 RESTRICCIONES: R1 = Horas disponibles del Departamento 1: 2 X1 ≤ 300 R2 = Horas disponibles del Departamento 2: 3 X2 ≤ 540 R3 = Horas disponibles del Departamento 3: 2 X1 + 2 X2 ≤ 440 R4 = Horas disponibles del Departamento 4: 1.2 X1 + 1.5 X2 ≤ 300 CONDICION DE NO NEGATIVIDAD: X1 ≥ 0, X2 ≥ 0 1.4.4. JUGUETES SAC fabrica dos tipos de juguetes de madera: soldados y trenes. Se vende un soldado a $27 y se usan $10 de materia prima. Cada soldado que se produce aumenta los costos variables de mano de obra y los costos generales en $14. Se vende un tren a $21 y se usan $9 de materia prima. Cada tren producido aumenta los costos variables de mano de obra y los costos generales en $10. La producción de soldados y trenes de madera necesita dos tipos de trabajo especializado: carpintería y acabado. Un soldado requiere 2 horas de acabado y una hora de carpintería. Un tren requiere 1 hora de acabado y 1 hora de carpintería. Cada semana, la empresa puede conseguir toda la materia prima que necesita, pero solamente dispone de 100 horas de acabado y 80 horas de carpintería. La demanda de los trenes no tiene límite, pero se venden a lo más 40 soldados semanalmente. La firma quiere maximizar su ganancia semanal (ingresos - costos). a) Formule un modelo matemático para la situación de JUGUETES SAC que se pueda utilizar para maximizar su ganancia semanal. b) Determine gráficamente la región factible del problema y el punto donde se hace máxima la ganancia de la compañía. 34

OPTIMIZACIÓN DE SISTEMAS I

VARIABLES DE DECISIÓN: X1 = Cantidad de soldados de madera que se debe producir por semana. X2 = Cantidad de trenes de madera que se debe producir por semana. FUNCION OBJETIVO: Se debe maximizar la utilidad total (ingreso menos costo) de los dos productos: Ingreso por ventas = 27 X1 + 21 X2 Costo de materia prima = 10 X1 + 9 X2 Costo de mano de obra y costos generales = 14 X1 + 10 X2 Maximizar Z = Ingresos – Costos = 3 X1 + 2 X2 RESTRICCIONES: R1 = Horas disponibles de carpintería: 1 X1 + 1 X2 ≤ 80 R2 = Horas disponibles de acabado: 2 X1 + 1 X2 ≤ 100 R3 = Demanda máxima de soldados: X1 ≤ 40 CONDICION DE NO NEGATIVIDAD: X1 ≥ 0, X2 ≥ 0 1.4.5. Financiera Solución administra fondos de empresas y clientes pudientes. La estrategia de inversión se adecua a las necesidades de cada cliente. Para un cliente nuevo, a Financiera Solución (FS) se le ha autorizado invertir 1.2 millones de dólares en dos fondos de inversión: un fondo de acciones y un fondo de bonos. Cada unidad del fondo de acciones cuesta 50 dólares, con una tasa de rendimiento anual de 10%; cada unidad del fondo de bonos cuesta 100 dólares, con una tasa de rendimiento anual de 4%. El cliente desea minimizar el riesgo, pero quiere tener un ingreso anual sobre la inversión de por lo menos 60,000 dólares. De acuerdo con el sistema de medición de riesgo de FS, cada unidad adquirida en el fondo de acciones tiene un índice de riesgo de 8, y cada unidad adquirida en el fondo de bonos tiene un índice de riesgo de 3. El índice de riesgo más elevado asociado con el fondo de acciones indica, simplemente, que se trata de la inversión más riesgosa. El cliente de la financiera también ha especificado que se inviertan por lo menos 300,000 dólares en el fondo de bonos.

35

OPTIMIZACIÓN DE SISTEMAS I

¿Cuántas unidades de cada uno de los fondos deberá adquirir la financiera para el cliente, si el objetivo es minimizar el índice de riesgo total para esta cartera? ¿Cuál es el valor del riesgo total? VARIABLES DE DECISIÓN: X1 = Cantidad, en unidades, que se debe adquirir del fondo de acciones X2 = Cantidad, en unidades, que se debe adquirir del fondo de bonos FUNCION OBJETIVO: Se debe minimizar el riesgo total de la inversión. Minimizar Z = 8 X1 + 3 X2 RESTRICCIONES: R1 = Cantidad máxima a invertir de dólares: X1 + X2 ≤ 1200000 R2 = Rendimiento anual mínimo requerido: 5 X1 + 4 X2 ≥ 60000 R3 = Inversión mínima en bonos: 100 X2 ≥ 300000 CONDICION DE NO NEGATIVIDAD: X1 ≥ 0, X2 ≥ 0

1.5.

PROBLEMAS PROPUESTOS

1.5.1. Un frutero necesita 16 cajas de naranjas, 5 de plátanos y 20 de manzanas. Dos mayoristas pueden suministrarle para satisfacer sus necesidades, pero sólo venden la fruta en contenedores completos. El mayorista A envía en cada contenedor 8 cajas de naranjas, 1 de plátanos y 2 de manzanas. El mayorista B envía en cada contenedor 2 cajas de naranjas, una de plátanos y 7 de manzanas. Sabiendo que el mayorista A se encuentra a 150 km de distancia y el mayorista B a 300 km, calcular cuántos contenedores habrá de comprar a cada mayorista, con objeto de ahorrar tiempo y dinero, reduciendo al mínimo la distancia de lo solicitado. 1.5.2. Una compañía tiene dos minas: la mina A produce diariamente 1 tonelada de carbón de antracita de alta calidad, 2 toneladas de carbón de calidad media y 4 toneladas de carbón de baja calidad; la mina B produce 2 toneladas de cada una de las tres clases. La compañía necesita 70 toneladas de carbón de alta calidad, 130 de calidad media y 150 de baja calidad. Los gastos diarios de la mina A ascienden a 150 dólares y los de la mina B a 200 dólares.

36

OPTIMIZACIÓN DE SISTEMAS I

¿Cuántos días deberán trabajar en cada mina para que la función de coste sea mínima? 1.5.3. Imaginemos que las necesidades semanales mínimas de una persona en proteínas, hidratos de carbono y grasas son, respectivamente, 8, 12 y 9 unidades. Supongamos que debemos obtener un preparado con esa composición mínima mezclando dos productos A y B, cuyos contenidos por Kg son los que se indican en la siguiente tabla:

a) b)

Proteínas

Hidratos

Grasas

Costo (pts/kg)

A

2

6

1

600

B

1

1

3

400

¿Cuántos Kg de cada producto deberán comprarse semanalmente para que el costo de preparar la dieta sea mínimo? ¿Cuántos Kg de cada producto deberíamos comprar si el precio de A subiera a 1000 pts/Kg ?

1.5.4. En la elaboración de un producto A se necesita una sustancia B. La cantidad de A obtenida es menor o igual que el doble de B utilizada, y la diferencia entre las cantidades del producto B y A no supera los 2g mientras que la suma no debe sobrepasar los 5g. Además se utiliza por lo menos 1g de B y se requiere 1 g de A. La sustancia A se vende a 5 millones y la B cuesta 4 millones el gramo. Calcular la cantidad de sustancia B necesaria para que el beneficio sea máximo. 1.5.5. En una encuesta realizada por una televisión local se ha detectado que un programa con 20 minutos de variedades y un minuto de publicidad capta 30.000 espectadores, mientras que otro programa con 10 minutos de variedades y 1 minuto de publicidad capta 10.000 espectadores. Para un determinado período, la dirección de la red decide dedicar 80 minutos de variedades y los anunciantes 6 minutos de publicidad. ¿Cuántas veces deberá aparecer cada programa con objeto de captar el máximo número de espectadores?

37

OPTIMIZACIÓN DE SISTEMAS I

1.5.6. Una empresa tiene dos factorías A y B. En ellas fabrica un determinado producto, a razón de 500 y 400 unidades por día respectivamente. El producto ha de ser distribuido posteriormente a tres centros I, II y III, que requieren, respectivamente, 200, 300 y 400 unidades. Los costos de transportar cada unidad del producto desde cada factoría a cada centro distribuidor son los indicados en la tabla siguiente: I II III FABRICACIÓN A 50 60 10 500 u. B 25 40 20 400 u. DEMANDA 200 300 400 ¿De qué manera deben organizar el transporte a fin de que los gastos sean mínimos? 1.5.7. Una empresa fabrica dos tipos de tarjetas gráficas, de 16Mb y 32Mb de memoria, respectivamente. Se utilizan dos máquinas que emplean 2 min. en fabricar las de 16Mb y 3 min. en fabricar las de 32Mb. La cadena de montaje sólo puede funcionar, como máximo, 300 minutos diarios. Además cada máquina tiene una capacidad máxima de fabricación diaria de 125 unidades, entre las cuales no puede haber más de 90 tarjetas de 16Mb ni más de 80 tarjetas de 32Mb, siendo el beneficio neto de las primeras de 45$ y el de las segundas de 60$. ¿Cuántas tarjetas de 16Mb y 32Mb debe fabricar diariamente cada máquina para que el beneficio sea máximo?. 1.5.8. Una multinacional farmacéutica desea fabricar un compuesto nutritivo a base de dos productos A y B. El producto A contiene 30% de proteínas, un 1% de grasas y un 10% de azúcares. El producto B contiene un 5% de proteínas, un 7% de grasas y un 10% de azúcares. El compuesto tiene que tener, al menos, 25g. de proteínas, 6g. de grasas y 30g. de azúcares. El coste del producto A es de 0.6 pts/g. y el de B es de 0.2 pts/g. ¿Cuántos gramos de cada producto debe tener el compuesto para que el coste total sea mínimo?

38

OPTIMIZACIÓN DE SISTEMAS I

1.5.9. Una compañía minera tiene abiertas dos minas M1 y M2, desde las cuales transporta carbón a dos grupos G1 y G2 de una central térmica. De la mina M1 salen diariamente para la central 800T de antracita y de la mina M2 300T. De las 1100T, 500 tienen que ir hasta el grupo G1 y 600T hasta el grupo G2. El coste de cada tonelada transportada de M1 a G1 es de 60$, el de A1 a G2 de 80$, el de M2 a G1 de 40$ y el de M2 a G2 de 50$. ¿Cuántas toneladas hay que transportar desde cada mina hasta cada grupo para que el coste total sea mínimo?. 1.5.10. Una asociación agrícola tiene de dos parcelas: la parcela P1 tiene 400Ha de tierra utilizable y dispone de 500m3 de agua, mientras la parcela P2 tiene 900Ha de tierra utilizable y dispone de 1200m3 de agua. Los cultivos aconsejados son: remolacha y algodón. La remolacha consume 3m3 de agua por Ha y tiene un beneficio de 700$ por Ha y el algodón consume 2m3 de agua por Ha y tiene un beneficio de 500$ por Ha. Se ha establecido una cuota máxima por Ha para cada cultivo: 800 para la remolacha y 600 para el algodón, siendo el porcentaje total de terreno cultivado el mismo en cada parcela. Plantear el problema de programación lineal. 1.5.11. Una empresa constructora dispone de dos tipos de camiones C1 y C2 y quiere transportar 100T de arena a una obra. Sabiendo que dispone de 6 camiones tipo C1 con capacidad para 15T y con un coste de 4000pts por viaje y de 10 camiones tipo C2 con una capacidad de 5T y con un coste de 3000pts por viaje. a) b) c)

¿Cuál es el número posible de camiones que puede usar (gráficamente)? ¿Cuál es el número posible de camiones que debe usar para que el coste sea mínimo? ¿Cuál es el valor de dicho coste?.

1.5.12. Un quiosco de prensa vende bolígrafos a 20pts y cuadernos a 30pts. Llevamos 240pts y pretendemos comprar los mismos cuadernos que bolígrafos por lo menos. ¿Cuál será el número máximo de piezas que podemos comprar?

39

OPTIMIZACIÓN DE SISTEMAS I

1.5.13. Una compañía aérea dispone de dos tipos de aviones A1 y A2 para cubrir un determinado trayecto. El avión A1 debe hacer más veces el trayecto que el avión A2 pero no puede sobrepasar 120 viajes. Entre los dos aviones deben hacer más de 60 vuelos, pero menos de 200. En cada vuelo, A1 consume 900 litros de combustible y A2 700 litros. En cada viaje del avión A1 la empresa gana 30.000$ y 20.000$ por cada viaje del avión A2. a) b)

¿Cuántos viajes debe hacer cada avión para obtener el máximo de ganancias?. ¿Cuántos vuelos debe hacer cada avión para que el consumo de combustible sea mínimo?

1.5.14. Un joyero fabrica dos tipos de anillos: los anillos A1 precisan 1g. de oro y 5g. de plata vendiéndolos a $40 cada uno. Para los anillos tipo A2 emplea 1,5g. de oro y 1g. de plata y los vende a $50. El joyero dispone en su taller de 750g. de cada metal. ¿Calcular cuántos anillos debe fabricar de cada clase para obtener el máximo beneficio? 1.5.15. Electrón S.A. produce 2 tipos de monitores para PC; de 17” y 15” (conocidos como M17 y M15). Los pronósticos de mercado indican que será posible vender todos los monitores que se puedan producir para el siguiente mes. Cada monitor pasa por un proceso en el departamento electrónico (DE) y otro en el departamento mecánico (DM) y además es sometido a verificación de calidad en el dpto. CC. En el DE se disponen de 150 hrs. de operación, en el DM de 160 hrs. Por acuerdo con los trabajadores deben utilizarse al menos el 90% del total de una meta de 150 hrs. en el departamento de verificación de calidad CC. El monitor M17 tiene un costo de producción de $ 1,200 y se venderá a $ 1,700 y requiere de 10 hrs. de operación en el DE, 20 hrs en el DM y 30 hrs de control de calidad.

40

OPTIMIZACIÓN DE SISTEMAS I

El monitor M15 tiene un costo de producción de $ 1,000 y se venderá a $ 1,400 y requiere de 15 hrs. de operación en el DE, 10 hrs en el DM y 10 hrs de control de calidad. La gerencia de ventas por estrategia de ventas exige que se produzca al menos un monitor M15 por cada 2 monitores M17. E informa que se debe de cumplir con el pedido ya recibido de un cliente de por lo menos 4 monitores (en cualquier combinación M17 y M15). 1.5.16. La Smith Motors, Inc., vende automóviles normales y vagonetas. La compañía obtiene $300 de utilidad sobre cada automóvil que vende y $400 por cada vagoneta. El fabricante no puede proveer más de 300 automóviles ni más de 200 vagonetas por mes. El tiempo de preparación para los distribuidores es de 2 horas para cada automóvil y 3 horas para cada vagoneta. La compañía cuenta con 900 horas de tiempo de taller disponible cada mes para la preparación de automóviles nuevos. Plantee un problema de PL para determinar cuántos automóviles y cuántas vagonetas deben ordenarse para maximizar las utilidades. 1.5.17. La empresa La Preferida, del Valle de Ica, cultiva brócoli y coliflor en 500 acres de terreno en el valle. Un acre de brócoli produce $500 de contribución a las utilidades y la contribución de un acre de coliflor es de $1000. Debido a reglamentos gubernamentales, no pueden cultivarse más de 200 acres de brócoli. Durante la temporada de plantación, habrá disponibles 1200 horas-hombre de tiempo de plantadores. Cada acre de brócoli requiere 2.5 horas-hombre y cada acre de coliflor requiere 5.5 horas-hombre. Plantee un problema de PL para determinar cuántos acres de brócoli y cuántos de coliflor deben plantarse para maximizar la contribución a las utilidades. 1.5.18. Los supervisores de la producción de una refinería deben programar dos procesos de mezclado. Cuando se realiza el proceso 1 durante una hora se consumen 100 barriles de petróleo nacional y 300 barriles de petróleo importado. De manera similar, cuando se efectúa el proceso 2 durante una hora, se consumen 100 barriles de petróleo nacional y 200 barriles de petróleo importado, Con respecto a la producción, el proceso 1 genera 4,000 galones de gasolina y 1,750 galones de petróleo para uso doméstico por hora de operación. El proceso 2 genera 3,500 galones de gasolina y 2,250 galones de petróleo para uso 41

OPTIMIZACIÓN DE SISTEMAS I

doméstico, por hora. Para la siguiente corrida de producción, existen disponibles 1,200 barriles de petróleo nacional y 1 ,800 barriles de petróleo importado. Los contratos de ventas exigen que se fabriquen 28,000 galones de gasolina y 12,000 galones de petróleo para consumo doméstico. Las contribuciones a las utilidades por hora de operación son $1,000 y $1,100 para los procesos 1 y 2, respectivamente. a.

b.

42

Plantee un modelo de programación lineal para determinar el programa de producción que maximice la contribución total. Asegúrese de indicar las unidades de medición para sus variables de decisión y las unidades en las que se mide cada restricción. El Ministerio de Energía y Minas puede emitir un dictamen que limite la producción total de gasolina a no más de la mitad del petróleo que se fabrique para uso doméstico. ¿Qué restricción debe añadirse al modelo para plantear esta condición?

OPTIMIZACIÓN DE SISTEMAS I

CAPÍTULO 2 PROGRAMACION LINEAL: FORMULACIÓN DE PROBLEMAS

2.1.

EJERCICIOS FORMULADOS DE PROGRAMACION LINEAL

Ejercicio 2.1.1: Problema de transporte Un fabricante de jabón y detergentes tiene tres plantas, localizadas en Cincinnati, Denver y Atlanta. Los almacenes principales se encuentran en Nueva York, Boston, Chicago, Los Angeles y Dallas. En las tablas siguientes se proporcionan los requerimientos de ventas del próximo año para cada almacén y los costos de envío desde cada planta a cada almacén. Requerimientos de los almacenes Ubicación del almacén Nueva York Boston Chicago Los Angeles Dallas Total

Ventas anuales (miles de cajas) 50 10 60 30 20 ---170

Costos de envío de 1000 cajas de jabón De / A Cincinnati Denver Atlanta

Nueva York $ 120 210 150

Boston

Chicago

$ 150 220 170

$ 80 150 150

Los Angeles $ 250 100 240

Dallas $ 180 110 200

La capacidad de producción para las plantas de Cincinnatti, Denver y Atlanta son 100 mil, 60 mil y 50 mil cajas, respectivamente. La compañía quiere determinar un programa de entregas que minimice los costos totales de transporte de la compañía.

43

OPTIMIZACIÓN DE SISTEMAS I

FORMULACIÓN: Variables de Decisión: Sea: X11 = Número de cajas enviadas de la primera fábrica (Cincinnati) al primer almacén (Nueva York), en miles de cajas. Análogamente: X12, X13, X14, X15

=

Número de cajas enviadas de la primera fábrica (Cincinnati) al segundo, tercer, etc., almacén (Boston, Chicago, etcétera)

X21, X22, X23, X24, X25 =

Número de cajas enviadas de la segunda fábrica (Denver) al primer, segundo, etc., almacén.

X31, X32, X33, X34, X35 =

Número de cajas enviadas de la tercera fábrica (Atlanta) al primer, segundo, etc., almacén.

Función Objetivo: El objetivo es minimizar costos de transporte. Minimizar: Z =

120X11 + 150X12 + 80X13 + 250X14 + 180X15 + 210X21 + 220X22 + 150X23 + 100X24 + 110X25 + 150X31 + 170X32 + 150X33 + 240X34 + 200X35

El costo total es la suma de los productos, de cada fábrica a cada almacén, del costo de envío de la tabla multiplicado por el número de millares de cajas que se envían. Restricciones: Hay dos conjuntos de restricciones para este problema. El primero garantiza que se cumplirán las necesidades del almacén, entonces, para Nueva York: X11 + X21 + X31 = 50 Lo anterior estipula que la suma de las cajas que se envían a Nueva York de la primera fábrica (Cincinnati), la segunda (Denver) y la tercera (Atlanta) debe ser 50 mil cajas, el requerimiento de ventas de Nueva York. Para los otros almacenes se tiene: Boston: Chicago: 44

X12 + X22 + X32 = 10 X13 + X23 + X33 = 60

OPTIMIZACIÓN DE SISTEMAS I

Los Angeles: Dallas:

X14 + X24 + X34 = 30 X15 + X25 + X35 = 20

El segundo conjunto de restricciones garantiza que las fábricas no excedan sus capacidades de producción. De esta manera, para la fábrica de Cincinnati: X11 + X12 + X13 + X14 + X15 < 100 Esta expresión indica que la cantidad que se envía de la primera fábrica al primer almacén, al segundo, al tercero, etc., no debe exceder la capacidad de 100,000 cajas de la fábrica. En forma similar: Denver: Atlanta:

X21+ X22 + X23 + X24 + X25 < 60 X31 + X32+ X33 + X34 + X35 < 50

Por último, todas las X deben ser mayores o iguales que cero. La solución de este problema de programación lineal ofrecerá el programa óptimo de envíos (es decir, el de menor costo) para la compañía. Es un ejemplo de un tipo especial de problema, conocido, de manera bastante natural, como el problema de transporte. En resumen, la formulación completa de este problema es: Minimizar Z = 120X11 + 150X12 + 80X13 + 250X14 + 180X15 + 210X21 + 220X22 + 150X23 + 100X24 + 110X25 + 150X31 + 170X32 + 150X33 + 240X34 + 200X35 Sujeto a:

X11 + X21 + X31 = 50 X12 + X22 + X32 = 10 Restricciones de X13 + X23 + X33 = 60 requerimientos de X14 + X24 + X34 = 30 almacenes X15 + X25 + X35 = 20 X11 + X12 + X13 + X14 + X15 < 100 Capacidad X21 + X22 + X23 + X24 + X25 < 60 disponible de X31 + X32 + X33 + X34 + X35 < 50 fábricas

Condición de no negatividad X11, X12, …., ……, X35 > 0

45

OPTIMIZACIÓN DE SISTEMAS I

Ejercicio 2.1.2: Problema de mezcla Se obtienen distintos tipos de gasolina mezclando ciertas gasolinas que se obtienen directamente de las operaciones de refinería. En un proceso de refinamiento real hay varias gasolinas para mezcla, varias gasolinas que son productos finales (por ejemplo, distintos grados de gasolina para aviación y para motores) y varias características de importancia para la composición química de los diversos grados de gasolina (por ejemplo, octanaje, presión de vapor, contenido de azufre, contenido de goma). En este ejemplo simplificando se supondrá que la refinería sólo tiene dos tipos de gasolina para mezcla, con las características que se presentan en la siguiente tabla: Tabla 2.1. Características de las gasolinas para mezcla Mezclas disponibles Gasolina para mezcla, tipo 1 Gasolina para mezcla, tipo 2

Octanaje

Presión de vapor

Cantidad disponible

104

5

30 000 barriles

94

9

70 000 barriles

Estas gasolinas para mezcla pueden combinarse para obtener dos productos finales: gasolina para aviación y gasolina para motores. En la siguiente tabla se presentan las características que requieren estos productos finales. Tabla 2.2. Características de las gasolinas finales Productos finales Gasolina para aviación Gasolina para motores

Octanaje mínimo

Presión de vapor máxima

102

6

96

8

Ventas máximas 20 000 barriles Cualquier cantidad

Precio de venta (por barril) $45.10 $32.40

FORMULACIÓN: Variables de Decisión: Sea: X1 = Número de barriles de gasolina para la mezcla 1 utilizados en gasolina para aviación X2 = Número de barriles de gasolina para la mezcla 2 utilizados en gasolina para aviación 46

OPTIMIZACIÓN DE SISTEMAS I

X3 = X4 =

Número de barriles de gasolina para la mezcla 1 utilizados en gasolina para motores Número de barriles de gasolina para la mezcla 2 utilizados en gasolina para motores

Función Objetivo: La función objetivo es maximizar Z = Ingresos totales: Maximizar: Z = 45.10 (X1 + X2) + 32.40 (X3 + X4) = 45.10 X1 + 45.10 X2 + 32.40 X3 + 32.40 X4 Observe que X1 + X2 es la cantidad total de gasolina para aviación mezclada (en barriles); como se vende a 45.10 dólares por barril, los ingresos por este producto son 45.10 (X1 + X2). De manera análoga, los ingresos por la gasolina para motor son de 32.40 (X3 + X4) y la suma de estos términos representa los ingresos totales. Restricciones: Hay varios tipos de restricciones que afectan la forma en que la refinería mezclará su gasolina. La primera es el nivel de ventas o tamaño de la demanda, el hecho de que no pueden venderse más de 20 000 barriles de gasolina para aviación. Lo anterior puede expresarse así: X1 + X2 < 20 000 Otro conjunto de restricciones se refiere a las cantidades disponibles de las gasolinas para mezcla. Entonces: X1 + X3 < 30 000 Observe que X1 + X3 representa la cantidad total de gasolina para mezcla 1 (la suma de la cantidad utilizada en gasolina para aviación, X1, y la cantidad usada en gasolina para motores, X3). La ecuación anterior establece que la cantidad de gasolina para mezcla 1 no debe exceder la cantidad disponible, 30 000 barriles. Hay una restricción similar para la gasolina para la mezcla 2: X2 + X4 < 70 000 Otro conjunto de restricciones tiene que ver con el octanaje de las gasolinas finales. Recuerde que la cantidad total de la gasolina para aviación es X1 + X2 y su octanaje estará definido por las cantidades relativas de X1 + X2, de acuerdo con la siguiente fórmula:

47

OPTIMIZACIÓN DE SISTEMAS I

Octanaje de la gasolina para aviación

104 * X1 + 94 * X2 -------------------------X1 + X2

Las cifras 104 y 94 provienen de la tabla 2.1 y son los octanajes de las gasolinas para mezcla 1 y 2, respectivamente. En la tabla 2.2 se observa que el octanaje de la gasolina para aviación debe ser por lo menos 102, por lo cual se tiene la siguiente restricción: 104 X1 + 94 X2 ------------------- > 102 X1 + X2 Al acomodar la expresión para convertirla en restricción lineal, se tiene: 104 X1 + 94 X2 > 102 X1 + 102 X2 o (104 X1 – 102 X1) + (94 X2 – 102 X2) > 0, ó 2 X1 – 8 X2 > 0 Análogamente, para la gasolina para motores, se tiene: 104 X3 + 94 X4 > 96(X3 + X4) ó 8 X3 – 2 X4 > 0 El último conjunto de restricciones tiene que ver con los requisitos de presión de vapor de las gasolinas finales. En el caso de la gasolina para aviación, la restricción es: 5 X1 + 9 X2 < 6(X1 + X2) ó – X1 + 3X2 < 0 y el requisito de presión de vapor de la gasolina para motores es: 5 X3 + 9 X4 < 8(X3 + X4) ó – 3X3 + X4 < 0 En resumen, la formulación total del modelo de programación lineal es: Maximizar Z = 45.10X1 + 45.10X2 + 32.40X3 + 32.40X4 Sujeto a:

48

X1 + X2 < 20 000

Restricción de la demanda

X1 + X3 < 30 000 X2 + X4 < 70 000

Restricciones de disponibilidad de gasolina para mezclas

2 X1 – 8 X2 > 0 8 X3 – 2 X4 > 0

Restricciones de octanaje

OPTIMIZACIÓN DE SISTEMAS I

– X1 + 3X2 < 0 – 3X3 + X4 < 0

Restricciones de presión de vapor

Condición de no negatividad X1, X2, X3, X4 > 0 Ejercicio 2.1.3: Problema de programación de actividades Una empresa tiene un programa estricto de compromisos de entrega de un producto para los próximos seis meses. El costo de producción varía por mes, por los cambios anticipados en costos de materiales. La capacidad de producción de la compañía es de 100 unidades por mes con tiempo normal y hasta 15 unidades adicionales por mes con tiempo extra. La siguiente tabla muestra los requerimientos de entrega y los costos de producción por mes. Requerimientos y costos Mes Compromisos de entrega (unidades) Costo por unidad en tiempo normal Costo por unidad en tiempo extra

1

2

3

4

5

6

95

85

110

115

90

105

$30

30

32

32

31

32

$35

35

37

37

36

37

El costo de almacenar en inventario una unidad que no se vende es de dos dólares por mes. El problema de la compañía es determinar el número de unidades que debe producir cada mes en tiempo normal y tiempo extra para cubrir los requerimientos con el menor costo. La empresa no tiene unidades disponibles al iniciar el mes 1 y no quiere que sobren unidades al terminar el mes 6. FORMULACIÓN: Variables de Decisión: Sea: X1, X2, X3, X4, X5, X6 = Y1, Y2, Y3, Y4 , Y5 , Y6 =

Número de unidades producidas en tiempo normal cada mes Número de unidades producidas en tiempo extra cada mes

49

OPTIMIZACIÓN DE SISTEMAS I

I1, I2, I3, I4, I5 , I6

=

Número de unidades en almacén vendidas) al final de cada mes

(no

Función Objetivo: El objetivo es minimizar costos: Minimizar: Z = 30X1 + 30X2 + 32X3 + 32X4 + 31X5 + 32X6 + 35Y1 + 35Y2 + 37Y3 + 36Y4 + 37Y5 + 37Y6 + 2I1 + 2I2 + 2I3 + 2I4 + 2I5 + 2I6 La primera parte de esta expresión es el costo de la producción en tiempo normal multiplicado por las cantidades producidas en tiempo normal cada mes. La segunda parte representa el costo de producción en tiempo extra multiplicado por las cantidades que se producen en tiempo extra cada mes. La tercera parte es el costo de almacenamiento de las unidades que no se vendan, multiplicados por el número de unidades no vendidas cada mes. Restricciones: Las restricciones de la producción en tiempo normal son: Por último, se requiere un grupo de restricciones de enlace o de equilibrio para unir los períodos y asegurar que se cumplan los compromisos de entrega. Estas restricciones tienen la siguiente forma: (Fuentes de las unidades) = (Usos de las unidades) ó ⎡ Inventario ⎤ ⎡ Pr oducción en ⎤ ⎡ Pr oducción en ⎤ ⎡Compromisos ⎤ ⎡ Inventario ⎤ ⎢ Inicial ⎥ + ⎢tiempo normal ⎥ + ⎢ tiempo extra ⎥ = ⎢ de entrega ⎥ + ⎢ final ⎥ ⎣ ⎦ ⎣ ⎦ ⎣ ⎦ ⎣ ⎦ ⎣ ⎦ Para el mes 1, esto es: 0 + X1 + Y1 = 95 + I1 puesto que no hay inventario inicial. Al acomodar la ecuación: X1 + Y1 – I1 = 95 Para el mes 2: I1 + X2 + Y2 = 85 + I2, ó I1 + X2 + Y2 - I2 = 85 Para los otros meses: Mes 3: I2 + X3 + Y3 - I3 = 110 Mes 4: I3 + X4 + Y4 - I4 = 115 50

OPTIMIZACIÓN DE SISTEMAS I

Mes 5: Mes 6:

I4 + X5 + Y5 - I5 = 90 I5 + X6 + Y6 - I6 = 105

Puesto que el inventario final debe ser cero, la última restricción es: I6 = 0 En resumen, la formulación es: Minimizar: Z = 30 X1 + 30 X2 + 32 X3 + 32 X4 + 31 X5 + 32 X6 + 35 Y1 + 35 Y2 + 37 Y3 + 36 Y4 + 37 Y5 + 37 Y6 + 2 I 1 + 2 I 2 + 2 I 3 + 2 I 4 + 2 I 5 + 2 I6 Sujeto a:

X1 + Y1 I1 + X2 + Y2 I2 + X3 + Y3 I3 + X4 + Y4 I4 + X5 + Y5 I5 + X6 + Y6

- I1 = 95 - I2 = 85 - I3 = 110 - I4 = 115 - I5 = 90 - I6 = 105 X1 < 100 X2 < 100 X3 < 100 X4 < 100 X5 < 100 X6 < 100 Y1 < 15 Y2 < 15 Y3 < 15 Y4 < 15 Y5 < 15 Y6 < 15 I6 = 0

Restricciones de balance de inventario

Restricciones de producción en tiempo normal

Restricciones de producción en tiempo extra Restricción de inventario final

Condición de no negatividad X1, X2, X3, X4, X5, X6 > 0 Y1, Y2, Y3, Y4, Y5, Y6 > 0 I1, I2, I3, I4, I5, I6 > 0 Ejercicio 2.1.4: Problema de producción Una empresa fabricante de tostadoras eléctricas, debe tomar una decisión sobre la producción de un nuevo modelo. La empresa tiene la posibilidad de emplear 3 técnicas alternativas de producción: manual, semi-automática y mediante el empleo de robots.

51

OPTIMIZACIÓN DE SISTEMAS I

Los requerimientos de cada técnica se resumen en el siguiente cuadro:

Mano de Obra Especializada Mano de Obra no-especializada Tiempo de Taller de Ensamblado

TECNICA DE ENSAMBLADO Manual Semi-automática Robotizada 1 min 4 min 8 min 40 min 30 min 20 min 3 min 2 min 4 min

La disponibilidad de recursos para este producto son los siguientes: 4500 minutos de mano de obra especializada, 36000 minutos de mano de obra no-especializada y 2700 minutos de tiempo disponible de taller de ensamblado. El costo total de producción manual es de $7 por tostadora, de $8 por tostadora para la producción semi-automática, y de $8.5 por tostadora para la producción robotizada. Se necesita producir 1000 tostadoras. Calcular la cantidad de tostadoras que se debe producir con cada técnica de ensamblado. FORMULACIÓN: Variables de Decisión: Xi: Cantidad de tostadoras ensambladas con la técnica i: (= M, S, R) Función Objetivo: Se debe minimizar Costos Minimizar Z = 7 XM + 8XS + 8.5XR Restricciones: Disponibilidad de Recursos: Ensamblado Manual: 1 XM + 4 XS + 8 XR < 4,500 Ensamblado Semiautomático: 40 XM + 30 XS + 20 XR < 36,000 Ensamblado Robotizado 3XM + 2XS + 4XR < 2,700 Condición de Producción: XM + XS + XR = 1000 Condición de No Negatividad: XM > 0 , XS > 0 , XR > 0 52

OPTIMIZACIÓN DE SISTEMAS I

Ejercicio 2.1.5: Problema de asignación de personal El gerente de personal de “La Tortuga Veloz, S.A.”, está analizando la necesidad de mano de obra semi calificada durante los próximos seis meses. Se lleva 1 mes adiestrar a una persona nueva. Durante este período de entrenamiento un trabajador regular, junto con uno en adiestramiento (aprendiz), producen el equivalente a lo que producen 1.2 trabajadores regulares. Se paga $500 mensuales a quien está en entrenamiento, mientras que los trabajadores regulares ganan $800 mensuales. La rotación de personal entre los trabajadores regulares es bastante alta, del 10% mensual. El gerente de personal debe decidir cuántas personas necesita contratar cada mes para adiestramiento. En seguida se da el número de meses-hombre necesarios. También se desea tener una fuerza de trabajo regular de 110 al principio de julio. En cuanto al 1º de enero, hay 58 empleados regulares. Mes Enero Febrero Marzo

Meses-hombre requeridos 60 50 60

Mes Abril Mayo Junio

Meses-hombre requeridos 80 70 100

FORMULACIÓN: Este problema tiene un aspecto dinámico, ya que la fuerza de trabajo en cualquier mes depende de la fuerza de trabajo regular y en adiestramiento del mes anterior. Para cualquier mes, el número total de meses-hombre disponibles se puede expresar como sigue: Variables de Decisión Meses-hombre disponibles: Ri + 0.2Ai en donde: Ri = número de trabajadores regulares al principio del mes Ai = número de aprendices contratados en el mes. Función Objetivo El objetivo global del gerente de personal es minimizar el costo de personal Minimizar Z = 800 (R1 + R2 + R3 + R4 + R5 + R6) + 500 (A1 + A2 + A3 + A4 + A5 + A6)

53

OPTIMIZACIÓN DE SISTEMAS I

Restricciones Entonces los requerimientos de cada mes pueden expresarse por las restricciones: enero

R1 + 0.2A1

≥

60

febrero

R2 + 0.2A2

≥

50

marzo

R3 + 0.2A3

≥

60

abril

R4 + 0.2A4

≥

80

mayo

R5 + 0.2A5

≥

70

junio

R6 + 0.2A6

≥

100

julio (principio)

R7

≥

110

Debido a la rotación, el 10% de los trabajadores regulares se van cada mes. Así, el número de trabajadores regulares disponibles, por ejemplo, al principio de febrero sería: R2 = 0.9R1 + A1 En la misma forma, pueden escribirse las ecuaciones para el número de trabajadores disponibles al principio de cada mes: Enero febrero marzo abril mayo junio julio

R1 R2 R3 R4 R5 R6 R7

= = = = = = =

58 (dado) 0.9 R1 + A1 0.9 R2 + A2 0.9 R3 + A3 0.9 R4 + A4 0.9 R5 + A5 0.9 R6 + A6

Condición de No Negatividad: R1, R2, R3, R4, R5, R6, R7, A1, A2, A3, A4, A5, A6 > 0 Ejercicio 2.1.6: Problema de presupuesto militar Un gobierno ha dispuesto 7500 millones de dólares de su presupuesto general para fines militares. El 60% del presupuesto militar se usara para comprar tanques, aviones y misiles. Éstos pueden adquirirse a un costo por unidad de 3, 10 y 4 millones de dólares, respectivamente. Se ha 54

OPTIMIZACIÓN DE SISTEMAS I

decidido que se deben adquirir al menos 200 tanques y 200 aviones. Debido a la escasez de pilotos experimentados, también se ha decidido no comprar más de 300 aviones. Por razones estratégicas, la proporción de misiles respecto a aviones debe estar en el rango de ¼ a ½ .El objetivo es maximizar la utilidad total de esta armas, siendo las utilidades individuales 1, 3 y 2. Encuentre la solución óptima (ignore el hecho de que la respuesta deba darse en unidades enteras). FORMULACIÓN: Variables de Decisión: X1 = Cantidad de unidades de adquisición de tanques. X2 = Cantidad de unidades de adquisición de aviones. X3 = Cantidad de unidades de adquisición de misiles. Función Objetivo: Se debe maximizar la utilidad total de las armas: Maximizar: Z = 1 X1 + 3 X2 + 2 X3 Restricciones: Esta primera restricción nos indica la cantidad mínima que se deberá comprar de tanques es 200: X1 > 200 En esta segunda restricción nos indica la cantidad mínima que se deberá comprar de aviones es 200: X2 > 200 Según los datos mencionados no se pueden comprar más de 300 aviones debido a que no hay pilotos experimentados por tanto la siguiente restricción será: X2 < 300 También nos indica que por medidas estratégicas la proporción de misiles con respecto a aviones esta en el rango de ¼ a ½. 1 < 4

X3 < 1 X2 2

X3 > 1 ……...4X3 – X2 > 0 X2 4

X3 < 1 …………..2X3 – X2 < 0 X2 2

55

OPTIMIZACIÓN DE SISTEMAS I

La ultima restricción será:

3 X1 + 10X2 + 4X3 < 4500

Ya que el 60% del presupuesto general (7500) es destinado para fines militares y esta a su vez será igual a la suma de la cantidad de armas multiplicada por sus respectivos costos. En resumen la formulación completa de este problema es: Maximizar:

Z = 1X1 + 3X2 + 2X3

Sujeto a: Cantidad mínima de tanques: Cantidad mínima de aviones: Cantidad máxima de aviones: Relación de misiles con aviones: Presupuesto militar:

X1 > 200 X2 > 200 X2 < 300 4X3 – X2 > 0 2X3 – X2 < 0 3 X1 + 10X2 + 4X3 = 4500

Condición de No Negatividad: X1, X2, X3 > 0. Ejercicio 2.1.7: Problema de Mezclas. Un comerciante compra azúcar a granel y vende al detalle. Para venderla tiene dos alternativas: envases de 1 kg y envases de 5 kg. El precio de venta es $300 y $250 por kg respectivamente, y en el mercado del azúcar al detalle se pueden vender 20.000 kg en envases de 1 kg y 17.000 en envases de 5 kg. Debido a un contrato anterior se deben entregar 5.000 kg en envases de 5 kg a un determinado cliente. El comerciante se puede abastecer de azúcar desde dos proveedores. El primero le puede vender hasta 15.000 kg a un precio de $90 por kg, y el segundo le ofrece la cantidad de azúcar que el comerciante desee, pero a un precio de $110 por Kg. y debido a requerimientos de sus distribuidores el comerciante debe vender menos del tercio del azúcar en envases de 1 kg. Además, suponga que el precio de los envases y el proceso de envasado son nulos, y que el comerciante no tiene azúcar almacenada y vende todo el azúcar que compra.

56

OPTIMIZACIÓN DE SISTEMAS I

Formule un problema de programación lineal que permita al comerciante decidir cual es el plan de abastecimiento y ventas de modo de obtener el mayor beneficio en su negocio. FORMULACIÓN: Variables de Decisión X i j : Cantidad , en kilogramos, de azúcar que van en envases de i (=A, B) suministrada por los proveedores j(= 1,2) Donde A: envases de 1 Kg. y B: envases de 5 Kg. 1: Proveedor uno y 2: Proveedor 2 Función Objetivo Maximizar Z = $300(XA1 + XA2) + $250(XB1 + XB2) - $90(XA1 + XB1) - $110(XA2 + XB2) Maximizar Z = 210XA1 + 190 XA2 + 160 XB1 + 140XB2 Esta ecuación establece que el mayor beneficio en el negocio del comerciante se forma con el ingreso obtenido de las ventas del azúcar en envases de 1Kg y envases de 5Kg (Trescientos dólares multiplicados por la sumatoria de los kilogramos del azúcar de 1 kg. del proveedor 1 2) menos el costo por la adquisición del azúcar del proveedor 1 (Noventa dólares multiplicados por la sumatoria del azúcar del proveedor 1 en envases de 1 Kg. Y 5 Kg.) menos el costo por la adquisición del azúcar del proveedor 2 (Ciento diez dólares multiplicados por la sumatoria del azúcar del proveedor 2 en envases de 1 Kg. Y 5 Kg.). Restricciones Ventas en el mercado de azúcar Contrato anterior Abastecimiento del azúcar Requerimiento de distribuidores

R1

(XA1 + XA2) ≤ 20,000 Kg.

R2 R3

(XB1 + XB2) ≤ 17,000 Kg. (XB1 + XB2) ≥ 5,000 Kg

R4

(XA1 + XB1) ≤ 15,000 Kg.

R5

(XA2 + XB2) ≥ 0.

R6

(XA1 + XA2) ≤ ⅓(XA1 + XA2 + XB1 + XB2) ⅔XA1 + ⅔XA2 - ⅓XB1 + ⅓XB2 ≤ 0

57

OPTIMIZACIÓN DE SISTEMAS I

Condición de No Negatividad Xij ≥ 0 ∀i = A, B ∀ j = 1. 2 Ejercicio 2.1.8: Planificación de Personal Las enfermeras de un hospital llegan cada 4 horas y trabajan en turnos de 8 horas continuas. La administración ha decidido la idea de definir 6 cambios de turno al día para minimizar las distracciones y los problemas de comunicación que ocurren en los cambios de turno. El hospital ha realizado un análisis del trabajo requerido durante cada uno de los seis períodos del día. Las características de cada período son las siguientes: Número Período Mínimo Enfermeras 2 AM - 6 AM 1 25 6 AM - 10 AM 2 60 10 AM - 2 PM 3 50 2 PM - 6 PM 4 35 6 PM - 10 PM 5 55 10 PM - 2 AM 6 40 HORA DEL DIA

Las enfermeras que empiezan a trabajar en los períodos 2, 3 y 4 ganan US$40 al día, y aquellas que comienzan en los períodos 1, 5 y 6 ganan US$50 al día. ¿Cuántas enfermeras deben empezar a trabajar en cada turno para minimizar los costos por salarios? FORMULACIÓN: Variables de Decisión En este caso podemos identificar como variable de decisión el número de enfermeras Ni que comienza a trabajar en el turno "i" (i = 1 : : : 6). Función Objetivo El objetivo será disminuir costos Minimizar Z = 50 N1 + 40 N2 + 40 N3 + 40 N4 + 50 N5 + 50 N6

58

OPTIMIZACIÓN DE SISTEMAS I

Restricciones Para construir las restricciones es conveniente recurrir a una representación gráfica de los turnos:

De la gráfica anterior se observa que en cada turno trabajan las enfermeras que comenzaron en dicho turno, pero también las que empezaron en el turno anterior. Por lo tanto, las restricciones de personal mínimo por turno quedan: N1 + N2 ≥ N2 + N3 ≥ N3 + N4 ≥ N4 + N5 ≥ N5 + N6 ≥ N6 + N1 ≥

60 50 35 55 40 25

Condición de No Negatividad Ni ≥ 0 ∀i = 1, 2, …., 6 Ejercicio 2.1.9: Programación de la producción Aceros Arequipa produce un acero especial usado en la industria de la construcción. El departamento de ventas de Aceros Arequipa ha recibido pedidos de 2400, 2200, 2700 y 2500 toneladas de acero para cada uno de los siguientes cuatro meses. La empresa puede satisfacer estas demandas produciendo el acero, extrayéndolo de su inventario, o usando cualquier combinación de las dos alternativas. Se proyecta que los costos de producción por tonelada de acero durante cada uno de los siguientes cuatro meses sean de $7400, $7500, $7600 y $7650. Como los costos suben cada mes, debido a las presiones inflacionarias, tal vez sea mejor que Aceros Arequipa produzca más acero del que necesita en un mes determinado y que almacene el exceso. La capacidad de 59

OPTIMIZACIÓN DE SISTEMAS I

producción, sin embargo, no puede exceder las 4000 toneladas en ningún mes. La producción mensual se termina al final del mes, cuando la demanda se satisface. Cualquier acero remanente se almacena en inventario a un costo de $120 por tonelada por cada mes que permanece allí. Estos datos se resumen en la siguiente tabla:

Demanda (tons) Costo de producción ($/ton) Costo de inventario ($/ton/mes)

Mes 1 2400 7400 120

Mes 2 2200 7500 120

Mes 3 2700 7600 120

Mes 4 2500 7650 120

Si el nivel de producción se incrementa de un mes al siguiente, entonces la compañía incurre en un costo de $50 por tonelada de producción incrementada para cubrir la mano de obra adicional y/o el tiempo extra. Cada tonelada de producción disminuida incurre en un costo de $30 para cubrir los beneficios de empleados no utilizados. El nivel de producción durante el mes anterior fue de 1800 toneladas, y el inventario que comienza es de 1000 toneladas. El inventario al final del cuarto mes debe ser de al menos 1500 toneladas para cubrir la demanda anticipada. Formule un plan de producción para Aceros Arequipa que minimice los costos totales en los siguientes cuatro meses. FORMULACIÓN: Variables de Decisión Pi = Cantidad, en toneladas, de acero que se produce en el mes i(=1,2,3,4). Ii = Cantidad, en toneladas, de acero que se quedan en el inventario al final del mes i(=1,2,3,4). Ai = Aumento, en toneladas, de la producción de acero en el mes i(=1,2,3,4) Di = Disminución, en toneladas, de la producción de acero en el mes i(=1,2,3,4) Función Objetivo Minimizar Z: costo de Producción + costo de Inventario + costo por Aumento de producción. + costo por Disminución de producción 7400*P1 + 7500*P2 + 7600*P3 + 7650* P4 +120*( I1 + I2 + I3 + I4) + 50* (A1 + A2 + A3 + A4) + 30*(D1 + D2 + D3 + D4)

60

OPTIMIZACIÓN DE SISTEMAS I

Restricciones: P0 = 1800 Producción del mes anterior al mes 1. P1 < 4000 Cantidad máxima de producción en el mes 1 P2 < 4000 Cantidad máxima de producción en el mes 2 P3 < 4000 Cantidad máxima de producción en el mes 3 P4 < 4000 Cantidad máxima de producción en el mes 4 I0 = 1000 Inventario inicial del mes 1 es 1000 toneladas I4>= 1500 Inventario para el mes 4 es al menos 1500 toneladas Relación Contable Mes Inventario Inicial + Producción – Inventario Final = Demanda 1 I0 + P1 I1 = 2400 2 I1 + P2 I2 = 2200 + P3 I3 = 2700 3 I2 4 I3 + P4 I4 = 2500 Relación de Aumento y Disminución de Producción Mes 1 P1 – P0 – A1 + D1 = 0 Mes 2 P2 – P1 – A2 + D2 = 0 Mes 3 P3 – P2 – A3 + D3 = 0 Mes 4 P4 – P3 – A4 + D4 = 0 El artificio para esta restricción es darle un valor positivo a Ai si la diferencia de producción de un mes a otro es positivo o darle un valor positivo a Di si la diferencia es negativa. Los dos valores (Ai y Di) no pueden ser positivas, uno de ellos es positivo y el otro cero, dado que la función objetivo determinará la más conveniente para minimizarlo. Condición de No Negatividad Ii > 0; para todo i(=1,2,3,4) Pi > 0; para todo i(=1,2,3,4) Ai > 0; para todo i(=1,2,3,4) Di > 0; para todo i(=1,2,3,4)

61

OPTIMIZACIÓN DE SISTEMAS I

2.2.

PROBLEMAS PROPUESTOS PARA FORMULAR DE PROGRAMACION LINEAL

2.2.1. PescaPerú S.A. desea determinar la mejor forma de cargar las bodegas de su bolichera “Sol” que le proporcione la mejor utilidad posible. “Sol” tiene una capacidad de carga total de 300 tn, divididas en 3 bodegas. La bodega-1 está designada para carga de camarones, la bodega-2 para mariscos y la bodega-3 para pescado. Por disposiciones de mercado a lo más deben de cargarse 130 tn de pescado y la carga de mariscos mas la de camarones no debe superar las 170 tn. Por medidas de seguridad, la carga en la bodega-2, no debe exceder a la tercera parte de la carga realizada en la bodega-3. La carga en la bodega-1, a lo más debe ser las dos quintas partes de lo cargado en la bodega-3. Las utilidades por ton. de camarones es de $1000, por ton de mariscos $700 y por pescado $ 600. 2.2.2. Bata para su campaña escolar, va a fabricar 3 modelos de calzado E1, E2 y E3. El proceso de manufactura del calzado hace que se requieran 2 operaciones de producción. Cada par de calzado de cualquier modelo requiere una hora de tiempo en la operación-1. En la operación-2 el modelo E1 requiere de 1 hora, el modelo E2 de 2 horas y el E3 de 2.5 horas. Bata cuenta con mano de obra suficiente para operar hasta 300 hrs. de tiempo semanal en la operación-1 y hasta 360 hrs. para la operación-2. Se ha calculado una utilidad de S/.14, S/16 y S/17 por cada par de calzado E1, E2 y E3 respectivamente. Se tiene un pronóstico de que la demanda del modelo E1 no será de mas de 500 pares por semana, y que l demanda combinada de los modelos E2 y E3 será como mínimo 300 pero no mas de 600. ¿Cuántos pares de calzado de cada modelo deben fabricarse para obtener la mayor utilidad?. 2.2.3. Dos aleaciones A y B se hacen de cuatro metales diferentes I, II, III, IV de acuerdo con las especificaciones siguientes.

62

OPTIMIZACIÓN DE SISTEMAS I

Aleación A B

Especificaciones A lo más 80% de I A lo más 30% de II Al menos 50% de IV Entre 40% y 60% de II Al menos 30% de III A lo más 70% IV

Los cuatro metales se extraen de diferentes minerales cuyos constituyentes en porcentaje de estos metales, cantidad máxima disponible y costo por tonelada se tabulan como sigue.

Constituyentes (%)

Mineral

Cantidad max. (tons)

I

II

III

IV

OTROS

Precio ($ / ton)

1

1000

20

10

30

30

10

30

2

2000

10

20

30

30

10

40

3

3000

5

5

70

20

0

50

Los precios de venta de las aleaciones A y B son 200 y 300 $/tonelada. Formule el problema como un modelo de programación lineal eligiendo la función objetivo apropiada que hará el mejor uso de la información dada [Sugerencia: sea xijk la cantidad (en toneladas) del metal i(i = I, II, III, IV) obtenida del mineral j(j = 1,2,3) y asignada en la aleación k-ésima (k = A, B).] 2.2.4. Una fábrica de zapatos ha pronosticado la demanda para los siguientes seis meses: 5.000, 6.000, 5.000, 9.000, 6.000 y 5.000 pares de zapatos respectivamente. A principios del mes 1, la empresa tiene 13 empleados. Cada empleado de la fábrica utiliza 15 minutos para producir un par de zapatos y trabaja al mes 150 horas y hasta 40 horas de tiempo extra. Asimismo, cada empleado recibe un salario de 2.000 dólares al mes más 50 dólares por cada hora de tiempo extra que trabaje. A principios de cada mes la empresa puede contratar o despedir empleados. A la empresa le cuesta 1.500 dólares contratar un trabajador y 1.900 dólares despedir un empleado. El costo mensual de almacenamiento de cada par de zapatos es igual al 3% del costo de producir un par de zapatos en tiempo normal. La materia prima necesaria para producir un par de zapatos cuesta 10 dólares.

63

OPTIMIZACIÓN DE SISTEMAS I

Formule un modelo de programación lineal, y resuélvalo utilizando el software SOLVER, de modo tal que minimice el costo de satisfacer la demanda de los siguientes seis meses. 2.2.5. La Corporación Financiera de Desarrollo (COFIDE) está estudiando un plan inversiones para los próximos dos años. Actualmente, COFIDE tiene tres millones de dólares para invertir. COFIDE espera recibir en 6, 12 y 18 meses un flujo de ingresos de las inversiones previas. En la siguiente tabla se presentan los datos: Ingresos de inversiones previas 6 MESES 12 MESES 18 MESES Ingreso $300,000 $400,000 $180,000 Hay tres proyectos de desarrollo (a los que llamaremos proyecto 1, proyecto 2 y proyecto 3 en los que COFIDE está planeando participar. En la siguiente tabla se muestra el flujo de caja que se tendría si COFIDE participara a un nivel del 100% en el proyecto 1 (los números negativos son inversiones y los positivos son ingresos). Así, para participar en el proyecto 1 a un nivel de 100% COFIDE tendría que desembolsar de inmediato $1,000,000. A los seis meses erogaría otros $700,000, etc.

Ingreso

TABLA 1: Flujo de Caja del Proyecto 1 INICIAL 6 MESES 12 MESES 18 MESES 24 MESES -$1,000,000 -$700,000 $1,800,000 $400,000 $600,000

En las siguientes tablas se muestran los flujos de caja de los proyectos 2 y 3: TABLA 2: Flujo de Caja del Proyecto 2 INICIAL 6 MESES 12 MESES 18 MESES 24 MESES Ingreso -$800,000 $500,000 -$200,000 -$700,000 $2,000,000 TABLA 3: Flujo de Caja del Proyecto 3 Ingreso

64

INICIAL 6 MESES 12 MESES 18 MESES 24 MESES -$500,000 $100,000 -$200,000 $300,000 $1,000,000

OPTIMIZACIÓN DE SISTEMAS I

Debido a la política de COFIDE la inversión mínima en cada uno de los tres proyectos tiene que ser de $100,000. A COFIDE no se le permite pedir prestado dinero. Sin embargo, al comienzo de cada período de 6 meses todos los fondos excedentes (esto es, los que no sean colocados en los proyectos 1, 2 y 3) se invierten con un interés del 7% para este período de 6 meses. COFIDE puede participar en cualquiera de los proyectos a un nivel menor que el 100%, en cuyo caso todos los flujos de efectivo de ese proyecto se reducirán en forma proporcional. Por ejemplo, si COFIDE opta por participar en el proyecto 1, a un nivel de 30%, el flujo de caja asociado con esta decisión sería 0.3 veces los datos de la Tabla 1. El problema que actualmente encara COFIDE es decir qué parte de los tres millones en efectivo debe invertirse en cada proyecto y cuánto debe colocarse simplemente por la renta del 7% semestral. La meta del administrador consiste en maximizar el efectivo que habrá al final de los 24 meses. Formule este problema como modelo de programación lineal. 2.2.6. Gracias a una adecuada estrategia de marketing y a la calidad del producto, cierta pequeña fábrica de canastos de mimbre ha recibido pedidos que superan su actual capacidad de producción. Durante las próximas cuatro semanas debe entregar 52, 65, 70 y 85 canastos, respectivamente. Actualmente cuenta con seis artesanos. La gerencia general de la fábrica ha decidido contratar personal nuevo para poder cumplir sus compromisos comerciales. Dada la escasez de artesanos, se deberá contratar personal sin experiencia. Un novato puede ser entrenado para llegar a ser aprendiz durante una semana. La segunda semana trabaja como aprendiz para ganar experiencia. Comenzando la tercera semana (después de dos semanas de trabajo) se transforma en artesano. La producción estimada y sueldos de los empleados es la siguiente:

65

OPTIMIZACIÓN DE SISTEMAS I

PRODUCCIÓN Canastos/sem ana Artesano dedicado sólo a la producción 10 Artesano dedicado a prod. y entrenamiento 5 Aprendiz 5 Novato 1

SALARIOS $/semana 30.000 40.000 15.000 5.000

Cada artesano puede entrenar hasta dos novatos por semana (el entrenamiento de un novato sólo dura una semana). Todo excedente de producción semanal puede ser guardado para cumplir los siguientes compromisos comerciales. Los analistas de la empresa estiman que la demanda semanal de canastos difícilmente superará los noventa canastos, por lo que han decidido terminar el período sin novatos y aprendices, pero con al menos nueve artesanos. Los reglamentos sindicales de la empresa prohíben los despidos por reducción de personal. Formule un modelo de programación lineal que permita definir las contrataciones a realizar, de modo de cumplir los compromisos comerciales a costo mínimo.

66

OPTIMIZACIÓN DE SISTEMAS I

CAPÍTULO 3 EL MÉTODO SIMPLEX

El algoritmo del simplex provee una metodología rápida y efectiva para resolver problemas de programación lineal. El método simplex es un método que llega a la solución óptima a través de iteraciones sucesivas. Este método utiliza los conceptos básicos de álgebra matricial para determinar la solución de un problema lineal. Comienza una solución básica factible inicial y sucesivamente obtiene soluciones en las intersecciones que ofrecen mejores valores para la función objetivo. En cada iteración el método proporciona un indicador que evalúa la optimalidad de la solución encontrada, este indicador es el que nos permite identificar la solución óptima. DEFINICIONES Todo modelo de programación lineal, luego de habérsele agregado las variables de holgura y/o exceso, se convierte en un sistema de ecuaciones con n variables y m ecuaciones, siendo n>m, en donde las m restricciones del modelo dan origen a las m ecuaciones del sistema. Una solución de tal sistema es un vector n-dimensional que satisface la relación Ax = b. Solución básica.- Es aquella solución en la que (n – m) variables se han igualado a cero y los valores de las m variables restantes se han determinado resolviendo las m ecuaciones con m variables. Una solución de este tipo puede tener como máximo m componentes no nulos. Solución básica no degenerada.- Es una solución básica que tiene exactamente m componentes no nulos. Solución básica degenerada.- Es una solución básica que tiene menos de m componentes distintos de cero. PROPIEDADES DE LAS SOLUCIONES El problema de programación lineal consiste en hallar el vector columna X que es solución de: a11x1 + a12x2 + . . . + a1nxn = b1 (1) a21x1 + a22x2 + . . . + a2nxn = b2 ... ... ... ... .. am1x1 + am2x2 + . . . + amnxn = bm 67

OPTIMIZACIÓN DE SISTEMAS I

que haga máximo a: z = c1x1 + c2x2 + . . . + cnxn

(2)

y de tal manera que las variables Xj estén sujetas a las condiciones: xj >= 0 ; j= 1, 2, ..., n (3) En las ecuaciones (1), se supone que: 1. Algunas de las restricciones pueden haber sido desigualdades antes de que les fueran sumadas o restadas nuevas variables para convertirlas en ecuaciones. 2. Todas las bi >= 0, lo cual puede requerir que algunas de las ecuaciones deba multiplicarse por – 1. 3. n > m, y A es de orden m x n. La función (2) es la función objetivo, se observará que se ha definido el problema general de programación lineal como un problema de maximización, que es el problema que se presenta más frecuentemente en la práctica. TEOREMA 1. - Dado un problema de programación lineal, en el cual no puedan existir soluciones básicas factibles degeneradas y en el que se ha formado una solución básica factible en función de las m primeras variables, puede formarse una nueva solución básica factible introduciendo la variable xk, tal que k > m, si al menos un elemento de la k-ésima columna de la matriz reducida es positivo. LA SOLUCIÓN ÓPTIMA Cada solución tiene un valor y la función objetivo controla cuál de las muchas soluciones es la óptima. Si aplicamos el teorema 1 podemos, efectivamente, encontrar dicha solución óptima. TEOREMA 2.- Dado un problema de programación lineal en el cual son imposibles las soluciones básicas factibles degeneradas, en el que la solución óptima es única y donde pueden formarse soluciones básicas factibles adicionales, la solución máxima debe ser una solución factible. ELECCIÓN DE LA VARIABLE DE ENTRADA Si se sabe de qué tipo es la solución óptima y cómo generar una solución de tal índole a partir de otra, se necesita una regla que permita decidir qué variable se debe introducir en la nueva solución.

68

OPTIMIZACIÓN DE SISTEMAS I

TEOREMA 3.- Supongamos que tenemos un problema de programación lineal y una solución básica factible del mismo. Si existe una variable Xk para la cual se puede realizar el cálculo de θ mediante su propia regla, se puede generar otra solución que aumente el valor de la función objetivo en un valor determinado de antemano. Por el teorema 1, tenemos que: Θ = min bi / aik , para todo aik >0 Y por el teorema 2 tenemos que: z’ = z + Θ(ck – zk) Puesto que estamos tratando de encontrar la solución máxima, es lógico exigir que z’ > z. Ya que z’ depende de la variable xk que elegimos, es lógico también requerir que k se elija de forma que el cambio del valor de la solución, z’ – z, resulte tan grande como sea posible, como quiera que esto no sucede normalmente, es necesario el cálculo de todas las θ posibles. En vez de hacer esto, el método usual consiste en determinar como variable que formará parte de la nueva solución, la que satisface ck – zk = max ( cj –zj) Esta es la regla llamada del “ascenso más rápido”. EL ALGORITMO DEL SIMPLEX El algoritmo del Simplex presupone que ya se ha hallado una solución básica factible en el momento de iniciar las iteraciones. El procedimiento detallado para hallar la solución máxima del problema general de programación lineal es el siguiente: 1. a) b) c) 2.

Calculamos cj –zj para cada variable que no está en la presente solución. Si para al menos un j, cj –zj es positivo y si al menos un aij para este j es positivo, entonces existe una mejor solución factible. Si para un j, cj –zj es positivo, pero los aij para este j son no positivos, entonces la función objetivo no está acotada. Si cj –zj es no positivo para todo j, entonces la solución óptima se ha encontrado. Si estamos en el caso 1(a), identificamos la variable que da el mayor cj –zj como xk ( es la columna pivote). Llamamos xr a la variable que se reduciría a cero al aplicar la regla del θ (es la fila pivote). El elemento ark se llama elemento pivote.

69

OPTIMIZACIÓN DE SISTEMAS I

3.

4.

Dividimos la r-ésima fila por ark (elemento pivote), para reducir a 1 el correspondiente elemento de ark en la tabla siguiente. Efectuamos luego las operaciones de fila que reducirán a cero todos los otros aik . Repetimos los pasos 1, 2 y 3 hasta que en alguna tabla se cumpla la condición 1(c). Entonces se ha obtenido la solución óptima.

La teoría del método simplex asegura que la solución hallada en cada paso tiene un valor mayor, o al menos igual, que el de la solución anterior. Puesto que el número de soluciones básicas del sistema de ecuaciones es finito, el algoritmo debe converger hacia la solución óptima en un número finito de iteraciones. La experiencia ha demostrado que el número de iteraciones que deben efectuarse en la mayoría de los problemas que se encuentran en la práctica, oscila entre m y 2m . UN PROBLEMA DE MAXIMIZACIÓN El tipo más sencillo de problema con el que se trabaja en programación lineal es el de maximización, en el que todas las restricciones son de tipo “menor o igual que”. Problema.- Una empresa cuenta con 1000 toneladas del mineral b1, 2000 toneladas del mineral b2 y 500 toneladas del b3. A partir de dichos minerales pueden extraerse y fundirse los productos metálicos 1, 2 y 3. Los requerimientos de fabricación sobre los productos son los siguientes: Una tonelada del producto 1 requiere 5 toneladas de b1, 10 de b2 y 10 de b3. Una tonelada del producto 2 requiere 5 toneladas de b1, 8 de b2 y 5 de b3. Una tonelada del producto 3 requiere 10 toneladas de b1, 5 de b2 y nada de b3. El fabricante obtendrá $100 de utilidades por tonelada del producto 1, $200 por tonelada del producto 2 y $50 por tonelada del producto 3. La empresa desea determinar cuántas toneladas de cada producto que debe fabricar, a partir de los minerales aprovechables, para obtener las máximas utilidades de la operación.

70

OPTIMIZACIÓN DE SISTEMAS I

En este problema de programación lineal, las variables funcionales o de decisión son X1, X2 y X3, que representan las unidades a fabricar de los productos o niveles de actividad. El problema es encontrar el vector X que hace máxima a la función objetivo: Z = 100x1 + 200x2 + 50x3 y que satisface las restricciones: 5x1 + 5x2 + 10x3 = 0

Puesto que cada una de las restricciones es del tipo “menor o igual”, debemos sumar nuevas variables no negativa, variables de holgura, para obtener: 5x1 + 5x2 + 10x3 + s1 10x1 + 8x2 + 5x3 10x1 + 5x2

+ s2 + s3

= 1000 = 2000 = 500

Vemos inmediatamente que este problema tiene una solución básica inicial: s1 = 1000; s2 = 2000; s3= 500 Que no sólo es una solución básica inicial sino que también es factible. El significado físico de las nuevas variables puede ilustrarse razonando sobre s1. Esta representa, tanto en la primera solución como en todas las posteriores, la cantidad del mineral b1 sobrante en el programa. Si s1 resulta igual a cero en la última iteración, será debido a que el programa emplea todo el material b1 para fabricar todos los productos del mismo. Puesto que x1=x2=x3=0 en la primera solución, s1 es igual a la cantidad total del mineral útil b1. Puesto que de no usar el mineral, no resulta de ningún provecho para la empresa, el coeficiente de s1 en la función objetivo es cero. Cuando se ha hallado una solución básica factible inicial, se está en condiciones de iniciar las iteraciones, para lo cual disponemos toda la información en la tabla simplex inicial:

71

OPTIMIZACIÓN DE SISTEMAS I

cj ck 0 0 0

xk s1 s2 s3

b 1000 2000 500 zj 0 cj – zj

100 x1 5 10 10 0 100

200 50 x3 x2 5 10 8 5 0 5 0 0 200 50 u x2: variable entrante

0 s1 1 0 0 0 0

0 s2 0 1 0 0 0

0 s3 0 0 1 0 0

θ 200 250 100t s3: variable saliente

En la tabla simplex, la primera fila cj, son los coeficientes de las variables en la función objetivo. La primera columna ck, es la de los coeficientes de las variables básicas en la primera solución, la columna xk es la “solución” y contiene a las variables básicas de la presente solución, la siguiente columna b contiene el vector solución, es decir los valores que toman las variables básicas en la solución actual. La sub-matriz bajo las columnas x1, x2 y x3 son los vectores estructurales, es decir los coeficientes de las variables en las restricciones. La sub-matriz bajo las columnas s1, s2 y s3, es la matriz de vectores unitarios. Nótese que cada variable básica está colocada en la fila en la que aparece el 1 de su vector unitario. En la última columna se encuentra el cociente θ. En la penúltima fila de la tabla se encuentra zj que se calcula a partir de la expresión: Zj = Σ ci aij Y la última fila de la tabla, cj – zj, se obtiene restando de cada elemento de la primera fila cj el correspondiente elemento zj. En la solución básica factible inicial, las variables básicas son: s1=1000; s2 = 2000; s3 = 500. Para hacer la prueba de optimalidad a esta solución, evaluamos la última fila de la tabla, cj–zj, vemos que existen varios valores positivos, por lo tanto la solución puede ser mejorada. Por la regla del ascenso más rápido, escogemos como variable entrante en la siguiente solución a la variable x2 por tener el cj–zj más positivo. Para determinar la variable saliente, calculamos el cociente θ = min bi / aik para todo los aik >0, la variable saliente 72

OPTIMIZACIÓN DE SISTEMAS I

resulta ser aquella variable que tiene el θ mínimo. La variable saliente es s3. En la siguiente tabla se ingresa la variable x2 en lugar de s3, el elemento pivote 5 se reduce a la unidad dividiendo toda la tercera fila entre 5, y se realizan operaciones elementales de fila para reducir a cero los elementos sobre y bajo el elemento pivote y conseguir formar el vector unitario bajo la columna correspondiente a la variable x2, variable que en esta nueva solución, es variable básica. Se continúan con el algoritmo del Simplex. Son necesarias dos iteraciones para encontrar la solución óptima. A continuación se presentan las dos tablas sucesivas.

ck 0 0 200

cj xk b 500 s1 s2 1200 x2 100 zj 20000 cj – zj

100 x1 -5 -6 2 400 -300

200 x2 0 0 1 200 0

50 x3 10 5 0 0 50 u x2: variable entrante

0 s1 1 0 0 0 0

0 s2 0 1 0 0 0

0 s3 -1 -1.6 0.2 40 -40

0 s1 0.1 -0.5 0 5 -5

0 s2 0 1 0 0 0

0 s3 -0.1 -1.1 0.2 35 -35

θ 50t 240 sal --

Tabla de 1ra Iteración

Tabla Óptima: cj ck xk b 50 50 x3 s2 950 0 x2 100 200 zj 20500 cj – zj

100 x1 -0.5 -3.5 2 375 -275

200 x2 0 0 1 200 0

50 x3 1 0 0 50 0

Tabla de la 2da Iterac.

La solución óptima encontrada en la segunda iteración, se interpreta como sigue: x1 = 0 x2 = 100 x3 = 50 s1 = 0

No se fabrica el producto 1. Fabricar 100 tons. del producto 2. Fabricar 50 tons. del producto 3. Se consume todo el mineral b1.

73

OPTIMIZACIÓN DE SISTEMAS I

s2 = 950 Hay un sobrante de 950 tons. del mineral b2. s3 = 0 Se consume todo el mineral b3. Z = $22,500 Máximo beneficio. VARIABLES ARTIFICIALES MAXIMIZACIÓN CON RESTRICCIONES DADAS POR “IGUALDADES” Consideremos el mismo ejemplo de maximización, una de cuyas restricciones sea una igualdad. Podemos cambiar una restricción del problema, exigiendo que todo el mineral b2 sea utilizado en la producción. Esto se hace incluyendo una igualdad en la segunda restricción. Considerando las variables de holgura en la primera y tercera restricción, tenemos: 5x1 + 5x2 + 10x3 + s1 10x1 + 8x2 + 5x3 10x1 + 5x2

+ s2

= 1000 = 2000 = 500

Este sistema no tiene una solución básica inicial, y por lo tanto no se cumplen las condiciones para empezar el algoritmo del simplex. Para corregir esta situación, se introduce una “variable artificial” en la segunda restricción, para obtener: 5x1 + 5x2 + 10x3 + s1 10x1 + 8x2 + 5x3 10x1 + 5x2

+ a1 + s2

= 1000 = 2000 = 500

Ahora se tiene la solución básica inicial s1=1000; a1=2000; s2=500. Entonces, una “variable artificial” en una restricción cumple el papel de variable básica en su restricción. A pesar de que se ha corregido el problema para que tenga la forma correcta para empezar el algoritmo, esta no es una solución factible. Las variables artificiales representan “productos imaginarios”, por lo que cualquier valor que tome a1 no satisface la segunda restricción. Dejamos a1 en el sistema de forma que se pueda empezar el algoritmo, pero intentaremos sacarla de éste asignándole un coeficiente en la función objetivo igual a –M. La función objetivo se convierte entonces en: Z = 100x1 + 200x2 + 50x3 – M a1 Si Z se debe maximizar, a1 debe dejar la solución y por lo tanto debe hacerse igual a cero, puesto que M se considera un número muy 74

OPTIMIZACIÓN DE SISTEMAS I

grande, tan grande que puede dominar a los demás números que aparecen en el problema. El hecho de penalizar a la variable artificial en la función objetivo con un coeficiente M, se conoce como el método de la penalización. Cuando se realiza el algoritmo simplex en un problema que contiene variables artificiales, existen tres casos posibles: 1)

2) 3)

Antes de obtener la tabla en la cual todos los cj – zj = 0 Sumando s1 y a1 respectivamente a la primera y segunda restricciones, restando s2 y sumando a2 a la tercera restricción, se tiene la siguiente solución por el método simplex: ck 0 M M

cj xk s1 a1 a2 zj cj-zj

0 M 1

s1 a1 x3 zj cj-zj

b 5 2 1 3M

11/2 1 ½ M+ ½

2 x1 1 2 -1 M 2-M

4 x2 2 -1 2 M 4-M

1 x3 -1 2 2 4M 1-4M

0 s1 1 0 0 0 0

0 s2 0 0 -1 -M M

M a1 0 1 0 M 0

M a2 0 0 1 M 0

½ 3 -½ - ½ +3M 5/2 -3M

3 -3 1 1-3M 3+3M

0 0 1 1 0

1 0 0 0 0

-½ 1 -½ -M ½ -M

0 1 0 M 0

½ -1 ½ M - ½ +2M

77

OPTIMIZACIÓN DE SISTEMAS I

0 2 1

0 0 1

s1 x1 x3

16/3 1/3 2/3

0 1 0

7/2 -1 ½

0 0 1

zj cj-zj

4/3

2 0

-3/2 11/2

1 0

6 1 1 1

2 3 1 1 1

3/2 -3 -1/2 -1/2 9/2

0 0 1 1 0

s1 s2 x3 zj cj-zj

1 0 0

-2/3 1/3 -1/3

-1/6 1/3 1/6

2/3 -1/3 1/3

0 0

1/3 -1/3

5/6 -5/6+M

-1/3 1/3+M

1 0 0 0 0

0 1 0 0 0

½ 1 ½ ½ -½+M

0 -1 0 0 M

Método Simplex: Problemas resueltos La empresa Industria Alfa importa componentes electrónicos para ensamblar dos modelos de computadoras personales. Uno de los modelos se denomina Deskpro Computer, y al otro, Portable Computer. Los administradores están interesados en elaborar el programa semanal de producción de ambos productos. El modelo Deskpro genera una contribución a las utilidades de $50 por unidad y el Portable genera una contribución de $40 por unidad. Se tendrán disponibles un máximo de 150 horas de ensamble para la producción de la semana siguiente. Cada unidad de Deskpro requiere de 3 horas de tiempo de ensamblado, y cada una de Portable requiere de 5 horas. Además, en estos momentos la empresa tiene en inventario sólo 20 monitores de los que se emplean en la Portable; por ello, no es posible ensamblar más de 20 unidades de este tipo. En el almacén sólo se dispone de 300 pies3 de espacio para el almacenamiento de estos nuevos productos. Cada unidad de Deskpro requiere de 8 pies3 de almacén y cada unidad de Portable requiere de 5 pies3. Determine el plan de producción que maximiza las utilidades. El modelo de programación lineal de este problema es el siguiente: x1 = Unidades del modelo Deskpro que se ensamblan. x2 = Unidades del modelo Portable que se ensamblan. Max Z = 50x1 + 40x2 s. a. : 3x1 + 5x2 1x2 8x + 5x2

=

125

+ a2 >=

350

= ct (5) y >= 0 (6) que se denomina el problema dual. En estos modelos se tiene que: c: Es un vector fila de n componentes. b : Es un vector columna de m componentes. A : Es una matriz de orden m x n. At es la transpuesta de A. x: Es un vector columna de n componentes, cuyos valores deben ser hallados para maximizar la función Z sujeta a las restricciones (2) y (3). y: Es un vector columna de m componentes, cuyos valores deben ser hallados para minimizar la función W sujeta a las restricciones (5) y (6). El problema dual tiene los mismos parámetros que el problema primal, pero en diferentes lugares. RELACIONES PRIMAL – DUAL: a) A todo programa lineal Primal, le corresponde un programa Dual. b) La función objetivo Z del programa Primal está dada por: Z = c x y las restricciones del Dual están dadas por: 83

OPTIMIZACIÓN DE SISTEMAS I

c)

d) e)

At y >= ct Como c es un vector fila de n componentes, entonces el Primal tiene n variables, mientras que el Dual tiene n restricciones. Por lo tanto cada variable del Primal corresponde a una restricción en el Dual. Como b es un vector columna de m componentes, entonces el primal tiene m restricciones ( A x = ct y sin restricción de signo (SRS).

El problema primal no siempre se presenta en su forma canónica, la siguiente tabla resume las reglas de transformación del problema primal a su correspondiente problema dual, sea cual sea la forma en que se presente el problema primal: Problema Primal Máx Mín Restricción1 i = bi Variable1 i >= 0 84

Problema Dual Mín Máx variable i >= 0 variable i SRS variable i = ci

OPTIMIZACIÓN DE SISTEMAS I

Variable1 i SRS Variable1 i = 0 b) Halle la solución óptima del primal o del Dual y determine qué máquina o máquinas trabajan a su máxima capacidad. Si se tratase de aumentar la producción ¿Que máquina o máquinas elegiría para contratar más horas? El valor de las variables duales en la solución óptima Dual también se puede determinar hallando la solución óptima del Primal. Se resolverá el Primal:

94

OPTIMIZACIÓN DE SISTEMAS I

Cj Ck

Xk

4

6

3

1

0

0

0

b

x1

x2

x3

x4

s1

s2

s3

0

s1

550

3/2

2

4

3

1

0

0

0

s2

700

4

1

2

1

0

1

0

0

s3

200

2

3

1

2

0

0

1

0

0

0

0

0

0

0

0

4

6

3

1

0

0

0

Zj Cj - Zj

0

s1 4162/3

1/6

0

10/3

5/3

1

0

-2/3

0

s2 6331/3

10/3

0

5/3

1/3

0

1

-1/3

6

x2

662/3

2/3

1

1/3

2/3

0

0

1/3

Zj

400

4

6

2

4

0

0

2

0

0

1

-3

0

0

-2

Cj - Zj

3

x3

125

3/60

0

1

15/30

3/10

0

-6/30

0

s2

425

195/60

0

0

-15/30 -5/10

1

0

6

x2

25

39/60

1

0

15/30 -1/10

0

12/30

1980

81/20

6

3

9/2

3/10 0

54/30

-1/20

0

0

-7/2

Zj

Cj - Zj

-3/10 0 -54/30

De la solución óptima del problema Primal, se tiene que los valores que toman las variables duales en la solución óptima son: y1 = 3/10 = $0.30 Valor de la hora adicional de la máquina A. y2 = 0 y3 = 54/30 = $1.80 Valor de la hora adicional de la máquina C. Las máquinas que trabajan a su máxima capacidad son las máquinas A y C, sus precios sombra son respectivamente $0.30 y $1.80. Si se tratase de elegir una de las dos máquinas para contratar horas adicionales y aumentar la producción, se elegiría la máquina C, por tener un mayor precio dual.

95

OPTIMIZACIÓN DE SISTEMAS I

96

OPTIMIZACIÓN DE SISTEMAS I

CAPÍTULO 5 ANALISIS DE SENSIBILIDAD Una vez que se ha resuelto un problema de programación lineal, puede darse el caso de que uno o varios parámetros de la formulación del problema original, tales como los precios unitarios o la disponibilidad de ciertos recursos cambien, dando origen a un nuevo problema. Si esto sucede, ¿Es necesario volver a resolver el problema desde el principio? La respuesta afortunadamente es no, porque existen métodos, llamados de análisis de sensibilidad, que permiten ahorrar muchas iteraciones, al resolver el nuevo problema partiendo de la solución óptima del problema original. El análisis de sensibilidad permite determinar el impacto del cambio en el modelo sin repetir el proceso de solución completo. Al estudiar y analizar el impacto que tiene cambiar el parámetro de un modelo, es deseable calcular los límites de un cambio, es decir, es deseable saber por ejemplo cuánto puede cambiarse un coeficiente específico de la función objetivo o un valor dado del segundo término de una restricción, sin cambiar la solución óptima que se tiene. Si el cambio propuesto de algún parámetro cae dentro de los límites permitidos de cambio, entonces no ocurre ningún cambio en la solución óptima y no es necesario resolver un problema nuevo. Si el cambio que se propone cae fuera de los límites, entonces la solución óptima existente ya no será óptima y deberá calcularse de nuevo la solución óptima del problema. Los procedimientos del análisis de sensibilidad permiten calcular estos límites del cambio. Se estudiaran varios casos de análisis de sensibilidad. Supongamos que tenemos el siguiente problema de programación lineal referente a la producción de tres tipos de fertilizantes: x1 = Toneladas de fertilizante 5-5-10 que se deben fabricar. x2 = Toneladas “ “ 5-10-5 “ “ “ “ x3 = Toneladas “ “ 5-5-5 “ “ “ “ Se trata de determinar qué combinación de producción de fertilizantes proporciona la máxima utilidad. Max Z = 18.5x1 + 20x2 + 14.5x3

97

OPTIMIZACIÓN DE SISTEMAS I

s.a.: 0. 05x1 + 0.05x2 + 0.05x3 = 4.0 se obtiene c3 >= 18.5

Cj ˆ Xk x1 x2 s3

Ci 18.5 20.0 0 Zj

Cj-Zj

b 8,000 14,000 500 428,000

18.5 x1 1 0 0 18.5 0

20 x2 0 1 0 20 0

14.5+Δ3 x3 1 0 -0.05 18.5 Δ3-4.0

0 s1 40 -20 -3 340 -340

0 s2 -20 20 1 30 -30

0 s3 0 0 1 0 0

Tabla 2. Tabla modificada para un cambio en c3.

Esto indica que si el precio de x3 se elevara un poco más de los $4.00, es decir, si su contribución a las utilidades fuera mayor que $18.50, entonces la producción de x3 se volvería más rentable que la mezcla actual de producción de 8,000 tons. de x1 y 14,000 tons. de x2. Si el precio se aumentara exactamente $4.00, se llegaría a un punto de decisión en el que podría fabricarse x3, pero no se obtendrían utilidades adicionales. Se obtendrían los mismos $428,000 de utilidades para esta solución óptima alternativa. 99

OPTIMIZACIÓN DE SISTEMAS I

En general, para cualquier variable no básica el cambio necesario para dicha variable se convierta en básica será: Δj >= | cj – zj | Si la utilidad de la variable básica disminuye, no hay cambio en la solución óptima; o si la contribución a las utilidades aumenta en un cantidad inferior a |cj – zj | para la variable, no habrá cambio en la solución óptima. Sólo si la contribución a las utilidades aumenta en una cantidad que sea mayor que el valor actual de |cj – zj | cambiará la solución óptima. Una variable no básica no se encuentra en la solución óptima porque las utilidades que se obtienen al fabricar ese producto son inferiores a lo que se perdería por hacerlo. CAMBIO EN EL COEFICIENTE DE LA FUNCIÓN OBJETIVO DE UNA VARIABLE BÁSICA (cj) Ahora se considerará el cambio en el coeficiente de utilidades de una variable básica, digamos x1, que representa el nivel de producción del fertilizante 5-5-10. Es deseable saber cuál es el valor máximo en el que puede cambiar este coeficiente básico de utilidades antes de cambiar las variables básicas restantes en la solución óptima. Si cambia la contribución a las utilidades de una variable básica, entonces puede producirse uno de dos resultados. Si el coeficiente de la variable básica disminuye, entonces es posible que la variable tuviera que dejar la base puesto que tal vez no fuera suficientemente rentable para seguir siendo básica. Por otro lado, si la contribución a las utilidades aumenta, podría obtenerse un mayor nivel de producción para la variable que se considera. A diferencia de los cambios en los coeficientes de contribución para las variables básicas, en el caso de las variables que sí son básicas deben considerarse tanto aumentos como disminuciones en los coeficientes cj. En este caso, los cambios en los coeficientes de contribución a las utilidades, tendrán de alguna manera algún impacto sobre la solución existente. Para analizar el efecto que tienen los cambios en la contribución a las utilidades para una variable básica, es posible añadir un coeficiente Δj al coeficiente cj que ya se tiene. De nuevo denotaremos la nueva contribución a las utilidades como cj = cj + Δj. En el caso de la variable no básica, la adición de Δj afectó sólo una columna de la tabla, sin embargo el caso de una variable básica puede resultar afectada más de una columna. Por ello, para determinar los límites de Δj, debemos examinar todos los valores (cj – zj) que se ven afectados por Δj. 100

OPTIMIZACIÓN DE SISTEMAS I

Para x1 la contribución a las utilidades es 18.5, se ha añadido un coeficiente Δ1 para utilizarlo en el análisis de cambios en las utilidades para x1. El resultado se muestra en la tabla 3. Cj ˆ Ck 18.5+Δ1 20.0 0 Zj

Xk x1 x2 s3

B 8,000 14,000 500 428,000+8,000 Δ1 Cj-Zj

18.5+ Δ1 x1 1 0 0 18.5+ Δ1 0

20 x2 0 1 0 20 0

14.5 x3 1 0 -0.05 8.5+ Δ1 -4.0+ Δ1

0 s1 40 -20 -3 340-40 Δ1 -340-40 Δ1

0 s2 -20 20 1 30-20 Δ1 -30+20 Δ1

Tabla 3. Tabla modificada para un cambio en c1.

Para que la solución actual siga siendo óptima, debe asegurarse que ningún valor (cj – zj) de la tabla 3 se vuelva positivo. La pregunta es ¿Cuánto puede cambiar c1 en una dirección positiva o negativa de manera que mantengan las condiciones de optimalidad? Puede determinarse la magnitud de estos cambios, Δj, despejando una desigualdad para cada uno de los valores no básico, es decir: cj – zj 0, entonces 2000 Y1 > 0 puede ocurrir solamente si Y1 = 1. De manera similar se expresaría la dependencia entre las variables asociadas a los otros dos procesos. El modelo completo quedaría: PROBLEMA DE PROGRAMACION ENTERA Minimizar Z = 5X1 + 4X2 + 3X3 + 100Y1 + 200Y2 + 300Y3 Sujeto a: X1 + X2 + X3 = 3500 X1 ≤ 2000 Y1 X2 ≤ 3000 Y2 X3 ≤ 4000 Y3 Yi es binaria para i = 1, 2, 3 Xi ≥ 0 para i = 1, 2, 3 Ejercicio 6.3.6. Problema de presupuesto de capital Cierta compañía tiene la oportunidad de invertir el próximo de año en cinco proyectos diferentes, P1, P2 , P3 , P4 y P5 , cada uno con un beneficio neto estimado como se muestra en la tabla . Proyecto número

Beneficio neto esperado ( 000 s )

1

$ 100

2

80

3

70

4

60

5

90

127

OPTIMIZACIÓN DE SISTEMAS I

Ya que de los diferentes requerimientos de cada proyecto (mano de obra, equipo, etc.), los costos varían de proyecto a proyecto. Además las obligaciones de requerimientos de flujo de caja hacen que la Cia. no pueda invertir en todos los cinco proyectos. En la tabla se listan los costos totales o salidas de caja requeridos para invertir en cada proyecto. Proyecto número

Costo esperado ( miles)

1

$ 60

2

40

3

20

4

40

5

50

La Cia. estima que tendrá una disponibilidad de caja en la cantidad de $150.000 para el próximo año. a)

¿En cuales proyectos podrid invertir la Cia. el próximo año? Variables de decisión Yj = ⎡ 1 si el proyecto j es seleccionado j = 1, 2, 3, 4, 5. ⎣ 0 si el proyecto j no es seleccionado Y1, Y2, Y3 , Y4 y Y5 Son variables de decisión binarias. Función objetivo. La meta de la Cia. es seleccionar los proyectos que maximicen la utilidad total esperada. Sea Z = utilidad total esperada, entonces la función objetivo es Maximizar Z = 100 Y1 + 80Y2 + 70Y3 + 60Y4 + 90Y5. Restricciones sobre las variables de decisión. Cantidad disponible para inversión: 60Y1 + 40Y2 + 20Y3 + 40Y4 + 50Y5 ≤ 150

128

OPTIMIZACIÓN DE SISTEMAS I

RESUMEN DEL MODELO DE PROGRAMACION ENTERA Maximizar Z = 100Y1 + 80Y2 + 70Y3 + 60 Y4+ 90Y5 Sujeto a 60Y1 + 40Y2 + 20Y3 + 40Y4 + 50Y5 ≤ 150 Yi es binaria para i = 1, 2, 3, 4, 5. b)

Suponga que la gerencia ha decidido que exactamente un proyecto puede ser seleccionado del conjunto de proyectos P1, P3 y P5. Pero, los proyectos P2 y P4 pueden ser seleccionados sujetos a la restricción de presupuesto. ¿Cuál restricción (o restricciones) necesita ser agregada al modelo original? Solución.- Ya que uno de P1, P3, ó P5 y solo uno, puede ser seleccionado, exactamente una de las tres variables Y1, Y3 y Y5 puede tomar el valor de uno. Y1 + Y2 + Y3 = 1, restricción a ser agregada

c)

Suponga que la Cia. ha decidido que no más de uno de los dos proyectos, P2 y P4, puede ser seleccionado. ¿Cuál restricción adicional necesita ser agregada al modelo original? Y2 + Y4 ≤ 1. Restricción a ser agregada

d)

Suponga que la Cia. ha decidido que si P3 es seleccionado, entonces P4 tiene que ser seleccionada. ¿Cuál restricción adicional necesita ser agregada al modelo original? Solución.Y3 ≤ Y4 restricción a ser agregada RESUMEN DEL MODELO DE PROGRAMACION ENTERA Maximizar Z = 100Y1 + 80Y2 + 70Y3 + 60 Y4+ 90Y5 Sujeto a 60Y1 + 40Y2 + 20Y3 + 40Y4 + 50Y5 ≤ 150 Y1 + Y2 + Y3 = 1 Y2 + Y4 ≤ 1 Y3 - Y4 ≤ 0 Yi es binaria para i = 1, 2, 3, 4, 5.

129

OPTIMIZACIÓN DE SISTEMAS I

EJERCICIOS PROPUESTOS 1.

Presupuesto de capital Se esta evaluando el capital de cinco proyectos a lo largo de un horizonte de planificación de tres años. La siguiente tabla proporciona las utilidades para cada proyecto, y los egresos anuales asociados Egresos (millones US$)/ anuales Proyecto

1

2

3

1 2 3 4 5 Fondos disponibles Mill.US.$

5 4 3 7 8

1 7 9 4 6

8 10 2 1 10

25

25

25

Utilidades Mill .US.$ 20 40 20 15 30

Determinar los proyectos que se van a ejecutar durante un horizonte de tres años 2.

Sea p1, p2, y p3 variables cuyos valores son 1, si se va a abrir una planta en particular y 0 en cualquier otro caso. Escriba una restricción lineal separada para cada caso abcdef-

3.

130

si se abre la planta 1, entonces la planta 2 no debería abrirse si se abre la planta 1, entonces la planta 2 debería abrirse al menos una de las tres plantas debería abrirse no más de dos de las tres plantas deberían abrirse si ni la planta 2 ni la planta 3 se abre, entonces la planta 1, no debería abrirse si se abre la planta 1 o la planta 3 no se abre, entonces la planta 2 debe abrirse

Inversiones La junta de directores estudia el conjunto de inversiones, donde Ri y Ci representan el rendimiento total y el costo de la inversión i. Se quiere maximizar no más de M dólares en total. Determinar un plan óptimo de inversión.

OPTIMIZACIÓN DE SISTEMAS I

inversión 1 2 3 4 5 6 7 4.

Condiciones Ninguna Sólo si 1 Sólo si 2 Se hará si 1 y 2 No si 1 o 2 No si 2 y 3 Solo si 2 y no 3

La ciudad de Lima piensa reubicar subestaciones de policía para obtener una mejor vigilancia en áreas de alta criminalidad. Las ubicaciones bajo consideración junto con las áreas que pueden ser cubiertas a partir de dichas ubicaciones son: Ubicación potencia de las subestaciones A B C D E F G

Áreas cubiertas 1.5.7 1.2.5.7 1.3.5 2.4.5 3.4.6 4.5.6 1.5.6.7

Formular un modelo de programación de enteros que se pudiera utilizar para encontrar el número mínimo de localizaciones necesarias a fin de proporcionar cobertura para todas las áreas. 5.

El Congreso está en proceso de formar un comité para que maneje las quejas de los estudiantes. Las instrucciones recibidas de la administración son incluir por lo menos una mujer, un hombre, un estudiante, un administrador y un miembro del profesorado. Diez personas ha sido nominadas identificadas para simplificar, por las letras A a la J. La mezcla de estas personas en las diferentes categorías se da así: Categoría Mujeres Hombres Estudiantes Administradores Profesorado

Personas A,b,c.d,e F,g,h,i,j A,b,c,j E,f D,g,hi

El Congreso desea formar el comité mas reducido posible y, al mismo tiempo, garantizar la representación de cada una de las cinco categorías. Formular el modelo. 131

OPTIMIZACIÓN DE SISTEMAS I

5.

MAYORISTA CORP, fabrica 2 productos A y B. Cada unidad de producto A requiere 1 hora de servicio de ingeniería, 8 horas de mano de obra directa y 4 libras de material. Una unidad de B necesita 3 horas de ingeniería, 3 horas de mano de obra directa y 2 libras de material. Se dispone de 800 horas de ingeniería, 880 horas de manos de obra directa y 3000 libras de material cada mes. Utilidades

Producto A Venta de unidades 0-40 40-100 100-150 150-más

Utilidad unitaria $ 10 9 8 6

Producto A Venta de unidades

Utilidad unitaria $

0-50 50-100 100-más

6 4 3

Determinar un plan óptimo para maximizar la utilidad mensual 6.

Programación en una aerolínea Lan Perú. Desea programar no más de un vuelo desde Lima hasta cada una de las siguientes ciudades: Arequipa, Cusco, Trujillo y Iquitos, Los horarios de salida disponibles son 8, 10, 12 de la mañana. Lan Perú arrienda los aviones al costo de $5000 hasta la 10am, y de $3000 después de las 10 y está en posibilidad de arrendar cuando mucho 2 por horario de salida. La tabla muestra las utilidades esperadas por vuelo antes de los costos de arrendamientos. Determinar un plan óptimo de vuelos Tabla de utilidades en millares de dólar Arequipa Cusco Trujillo Iquitos

7.

132

8 am 10 9 14 18

10 am 6 10 11 15

12 am 6 9 10 10

Instalación Cada día un electricista debe decidir que generadores conectar. Tiene 3 generadores. Hay dos períodos en el día. En el primer período se necesitan 2900 MEGAWATTS. En el segundo, 3900 MW. Un generador que se conecte para el primer período puede ser usado en el segundo sin causar un nuevo gasto de conexión.

OPTIMIZACIÓN DE SISTEMAS I

Todos los generadores principales como lo son A, B y C son apagados al término del día. Formular el modelo Costo fijo de conexión $3.000 $2.000 $1.000

Generador A B C

8.

Costo por período por MW usado $5 $4 $7

Capacidad máxima por período MW 2.100 1.800 3.000

Natación. El entrenador de natación debe asignar competidores para la prueba de 200 metros combinados por equipos para mandarlos a las olimpiadas juveniles. En la siguiente tabla se muestra los cinco mejores nadadores y sus mejores tiempos en segundos

Tipo nado Dorso Pecho Mariposa Libre

Carl 37.7 43.4 33.3 29.2

Cristina 32.9 33.1 28.5 26.4

David 33.8 42.2 38.9 29.6

Antonio 37.0 34.7 30.4 28.5

José 35.4 41.8 33.6 31.1

Determinar la mejor asignación para minimizar la suma de los mejores tiempos

133

OPTIMIZACIÓN DE SISTEMAS I

134

OPTIMIZACIÓN DE SISTEMAS I

BIBLIOGRAFÍA

G.D.EPEN, F.JGOULD, C.P SCHMIDT, JEFFREY H MOORE, LARRY R. WEATHERFORD. (2000). Investigación de Operaciones. En la ciencia administrativa. Quinta edición. Prentice Hall, México. HILLIER FREDERIC, LIEBERMAN GERALD. (2004). Investigación de Operaciones. Mc Graw-Hill Interamericana. Editores, S.A. de CV.

RIOS INSUA SIXTO, RIOS INSUA DAVID, CABALLERO ALFONSO, JIMENEZ JACINTO. (1997). Programación Lineal y Aplicaciones. Ejercicios resueltos. Alfaomega Grupo Editor, España. SWEENEY, ANDERSON, WIILIAMS (1998). Métodos Cuantitativos para los negocios. Séptima Edición. Internacional Thomson Editores, México TAHA, HAMDY A. (1998). Investigaciones de Operaciones. Una introducción. Sexta Edición. Prentice Hall. S.A. México

WINSTON WAYNE. (1994). Investigación de Operaciones. Aplicaciones y Algoritmos. Grupo Editorial Iberoamericana. México.

135