Source: https://ru.scribd.com/document/118588275/Modelos-Matematicos-de-Optimizacion
Timestamp: 2020-05-27 23:41:39+00:00

Document:
Modelos Matematicos de Optimizacion | La investigación de operaciones | Optimización Matemática
Modelos Matematicos de Optimizacion
сохранитьСохранить «Modelos Matematicos de Optimizacion» для последующего чтения
Álvaro Baíllo Pedro Linares Andrés Ramos Pedro Sánchez Ángel Sarabia Begoña Vitoriano
Alberto Aguilera 23 – E 28015 Madrid – Tel: 34 91 542 2800 – Fax: 34 91 541 1132 – www.doi.icai.upco.es
I.1. OPTIMIZACIÓN
I.1.1. Investigación operativa y optimización I.1.2. Referencias 9
I.2. MODELOS DE OPTIMIZACIÓN
I.2.1. Modelo y modelado
I.2.2. Etapas en el desarrollo de un modelo
I.2.3. Referencias
I.3. FORMULACIÓN DE PROBLEMAS DE OPTIMIZACIÓN
I.3.1. Modelos característicos de programación lineal
I.3.2. Modelos característicos de programación entera
I.3.3. Modelado de restricciones con variables binarias
I.3.3.1. Modelado de algunas restricciones especiales
I.3.3.2. Modelado de implicaciones lógicas
I.3.3.3. Modelado de proposiciones condicionales y/o compuestas
I.3.3.4. Modelado de productos con variables binarias
I.3.4. Modelos característicos de programación no lineal
I.3.5. Referencias
I.3.6. Biblioteca de problemas
I.3.7. Resultados de la biblioteca de problemas
I.4. CODIFICACIÓN DE PROBLEMAS DE OPTIMIZACIÓN
I.4.1. Lenguajes de modelado (OAE)
I.4.1.1. Lenguajes de modelado
I.4.1.2. Lenguajes algebraicos de modelado
I.4.1.3. Referencias
I.4.2. Modelado en GAMS (OAE)
I.4.2.1. Ejemplo de transporte
I.4.2.2. Ejemplo de planificación de la producción
I.4.2.3. Ejemplo de secuenciación de órdenes de trabajo
I.4.2.4. Ejemplo del viajante de comercio
I.4.2.5. Ejemplo de asignación de grupos térmicos
I.4.2.6. Ejemplo de flujo de cargas óptimo
I.4.3. Elementos de estilo de programación (DOCT)
I.4.3.1. Generales
I.4.3.2. Específicos de GAMS
I.4.3.3. Referencias
I.5. OPTIMIZACIÓN LINEAL
I.5.2. Solución gráfica
I.5.3. Geometría de la programación lineal
I.5.4.1. Problema de maximización
I.5.4.2. Múltiples óptimos
I.5.4.3. Convergencia del algoritmo
I.5.4.4. Variables acotadas superiormente
I.5.4.5. Forma tabular
I.5.4.6. Solución básica factible inicial
I.5.4.7. Método simplex revisado (DOCT)
I.5.4.8. Forma producto de la inversa (DOCT)
I.5.4.9. Factorización de la matriz base (DOCT)
I.5.4.10. Estrategias de cálculo de costes reducidos (DOCT)
I.5.5. Dualidad
I.5.6. Análisis de sensibilidad
I.5.6.1. Cambios en cotas de restricciones
I.5.6.2. Cambio en un coeficiente de una variable no básica
I.5.6.3. Introducción de una nueva variable
I.5.6.4. Cambio en un coeficiente de una variable básica
I.5.6.5. Introducción de nueva restricción
I.5.7. Método simplex dual
I.5.8. Programación lineal paramétrica
I.5.8.1. Cambios simultáneos en coeficientes de la función objetivo
I.5.8.2. Cambios simultáneos en cotas de las restricciones
I.5.9. Método de punto interior primal-dual (DOCT)
I.5.10. Referencias
I.5.11. Biblioteca de problemas
I.5.12. Resultados de la biblioteca de problemas
I OPTIMIZACIÓN
I.1.1. Investigación operativa y optimización
“In the last decade, new advances in algorithms have been as important as the impressive advances in computer technology” George L. Nemhauser
“The technology improvements in algorithms, modeling languages, software, and hardware have made the methodology accessible, easy to use, and fast. So the Age of Optimization has arrived” George L. Nemhauser (1994).
Definir el término investigación operativa no es una tarea fácil ya que su evolución permanente hace que sea difícil dar con precisión una definición. La investigación operativa se puede definir como la aplicación de métodos científicos en la mejora de la efectividad en las operaciones, decisiones y gestión, ver [Robinson, 1999]. Otra definición más extensa es la siguiente: la investigación operativa es la aplicación, por grupos interdisciplinarios, del método científico a los problemas complejos producidos en la dirección y gestión de grandes sistemas de hombres, máquinas, etc. La principal característica consiste en construir un modelo científico del sistema del cual se pueden predecir y comparar los resultados de diversas estrategias, decisiones, incorporando medidas del azar y del riesgo. El objetivo es ayudar a los responsables a determinar su política y actuaciones en forma científica. Los profesionales de la investigación operativa colaboran con los decisores en el diseño y mejora de las operaciones y decisiones, resuelven problemas y ayudan en las funciones de gestión, planificación o predicción, aportan conocimiento y ayuda en la toma de decisiones. Aplican las técnicas científicas más adecuadas seleccionadas de la matemática, ingeniería o cualquier ciencia social o de administración de empresas. Su trabajo normalmente consiste en recoger y analizar datos, desarrollar y probar modelos matemáticos, proponer soluciones o recomendaciones, interpretar la información y, en definitiva, ayudar a implantar acciones de mejora. Como resultado desarrollan e implantan aplicaciones informáticas, sistemas, servicios técnicos o productos. La investigación operativa tiene sus orígenes en la Segunda Guerra Mundial, debido a la necesidad urgente de asignación de recursos escasos en las operaciones militares, en problemas tácticos y estratégicos. Estas mismas técnicas se han extendido con posterioridad a las empresas.
I.1 OPTIMIZACIÓN
Disciplinas típicas de la investigación operativa son la optimización con sus múltiples sabores (lineal, no lineal, entera, estocástica, multiobjetivo), teoría de la decisión y de juegos, teoría de colas y simulación, teoría de grafos o flujos de redes. Otras disciplinas como algoritmos metaheurísticos y lógica borrosa, redes neuronales artificiales, reconocimiento de patrones y otras técnicas de inteligencia artificial, aunque conceptualmente se encuadran dentro de la investigación operativa, habitualmente se estudian dentro de otras disciplinas ligadas a la ingeniería informática como la inteligencia artificial. Los contenidos de algunas de estas últimas disciplinas también están muy ligados a la estadística. La optimización es una parte relevante dentro de la investigación operativa. Tuvo un progreso algorítmico inicial muy rápido. Muchas técnicas – programación lineal (linear programming) LP, programación dinámica (dynamic programming) DP– son anteriores a 1960. Por ejemplo, el método Simplex 1 de programación lineal debido a Dantzig 2 es de 1947, el principio de optimalidad de Bellman base de la programación dinámica se formuló en 1957. En la última década se han producido avances significativos generados por el desarrollo en 1984 por parte de Karmarkar de un método de punto interior para programación lineal. Por ejemplo, en una nota técnica de ILOG se presenta que desde su optimizador CPLEX 3.0 en 1994 a CPLEX 7.0 en 2000 la reducción de tiempo de resolución ha sido de 28 veces en el método simplex dual para un problema lineal concreto. Para otro caso se observa una mejora global, de software y algorítmica, de 10000 veces entre la versión de CPLEX 1.0 de 1988 y la 7.0 del 2000. Como referencia, se estima que la mejora en el rendimiento del hardware ha sido del mismo orden de magnitud. Si tomamos conjuntamente ambas mejoras hoy se pueden resolver problemas en segundos que habrían tardado años en ser resueltos hace una docena de años. Estos avances han sido tan importantes como los realizados en el campo de la informática, según la opinión de George L. Nemhauser uno de los expertos actuales en programación entera, y se han producido acompasadamente con ellos. Hoy es posible resolver un problema LP de 200000 ecuaciones con 200000 variables y 1000000 de elementos no nulos en la matriz de restricciones en un PC con suficiente memoria principal. Aproximadamente, para un problema LP se puede decir que se requiere 1 MB de memoria principal por cada 1000 ecuaciones.
1 En castellano la traducción de esta palabra es símplice pero no es habitual su uso para denominar este método de optimización lineal.
2 En http://www.e-optimization.com/directory/trailblazers/dantzig/ se puede encontrar un resumen de sus logros así como una entrevista sobre diversos temas, incluyendo imágenes en vídeo.
El estilo de este documento es eminentemente aplicado, práctico, ingenieril, a caballo entre una visión matemática de los problemas y de los algoritmos y la visión económica o de gestión empresarial de algunas de sus aplicaciones. Este documento trata de explicar suficientemente los fundamentos matemáticos como para permitir desarrollar aplicaciones de optimización de manera rigurosa y precisa. Al mismo tiempo, se presentan algunas aplicaciones a problemas concretos de ingeniería. Al final del capítulo se citan algunos libros generales o de referencia de investigación operativa que pueden servir de consulta o como texto para un nivel de pregrado y postgrado. Luego, en cada capítulo se indican además referencias específicas de los diferentes temas. Dentro de los libros generales, [Hillier y Lieberman, 2002] es un libro clásico de investigación operativa muy ampliamente utilizado que compendia numerosos temas y tiene una orientación ingenieril. [Taha, 1998] presenta los temas con una orientación más matemática mientras que [Winston, 1994] los presenta con una perspectiva más de administración de empresas. [Sarabia, 1996] da una base teórica suficiente para poder resolver una colección de problemas relacionados con el temario de investigación operativa. Entre las revistas principales que tratan sobre optimización se pueden incluir: Interfaces, Operations Research, Management Science, European Journal of Operational Research, Mathematics of Operations Research, OR/MS Today, Mathematical Programming, INFORMS Journal on Computing, Journal of the Operational Research Society, Omega, Journal of Optimization Theory and Applications, Transportation Science, Transportation Research. Existe una enciclopedia de investigación operativa que puede servir como consulta inicial y referencia de un tema específico, ver [Gass, 2001]. Además se puede encontrar información sobre los temas de investigación operativa en las direcciones de la Sociedad Española de Estadística e Investigación Operativa (SEIO) www.seio.es, de la Association of European Operational Research Societies (EURO) http://www.euro-online.org/, de la International Federation of Operational Research Societies (IFORS) www.ifors.org y del Institute for Operations Research and the Management Sciences (INFORMS) www.informs.org. La optimización consiste en la selección de una alternativa mejor, en algún sentido, que las demás alternativas posibles. Es un concepto inherente a toda la investigación operativa. Sin embargo, determinadas técnicas propias de la investigación operativa se recogen bajo el nombre de optimización o programación matemática. Los problemas de optimización se componen generalmente de estos tres ingredientes:
• función objetivo Es la medida cuantitativa del funcionamiento del sistema que se desea optimizar (maximizar o minimizar). Como ejemplo de funciones objetivo se pueden mencionar: la minimización de los costes variables de operación de un sistema eléctrico, la maximización de los beneficios netos de venta de ciertos productos, la minimización del cuadrado de las desviaciones con respecto a unos valores observados, la minimización del material utilizado en la fabricación de un producto, etc.
• variables Representan las decisiones que se pueden tomar para afectar el valor de la función objetivo. Desde un punto de vista funcional se pueden clasificar en variables independientes o principales o de control y variables dependientes o auxiliares o de estado, aunque matemáticamente todas son iguales. En el caso de un sistema eléctrico serán los valores de producción de los grupos de generación o los flujos por las líneas. En el caso de la venta, la cantidad de cada producto fabricado y vendido. En el caso de la fabricación de un producto, sus dimensiones físicas.
• restricciones Representan el conjunto de relaciones (expresadas mediante ecuaciones e inecuaciones) que ciertas variables están obligadas a satisfacer. Por ejemplo, las potencias máxima y mínima de operación de un grupo de generación, la capacidad de producción de la fábrica para los diferentes productos, las dimensiones del material bruto del producto, etc.
Resolver un problema de optimización consiste en encontrar el valor que deben tomar las variables para hacer óptima la función objetivo satisfaciendo el conjunto de restricciones. Los métodos de optimización los podemos clasificar en: métodos clásicos (que son los algoritmos que habitualmente se explican en los libros de optimización) y métodos metaheurísticos (que aparecieron ligados a lo que se denominó inteligencia artificial e imitan fenómenos sencillos observados en la naturaleza). Dentro de los primeros se encuentra la optimización lineal, lineal entera mixta, no lineal, estocástica, dinámica, etc. que se explican en el documento. En el segundo grupo se incluyen los algoritmos evolutivos (genéticos entre otros), el método del recocido simulado (simulated annealing), las búsquedas heurísticas (método tabú, búsqueda aleatoria, avariciosa, etc.) o los sistemas multiagente. De forma muy general y aproximada se puede decir que los métodos clásicos buscan y garantizan un óptimo local mientras que los métodos metaheurísticos
tienen mecanismos específicos para alcanzar un óptimo global aunque no garantizan su alcance. En la siguiente tabla se muestran las expresiones matemáticas generales de algunos tipos de problemas de optimización dentro de los métodos clásicos. Los problemas se distinguen por el carácter de las funciones que intervienen (lineales o no lineales) y de las variables (reales/continuas o enteras/discretas).
(mixed integer programming)
1 x Qx
(non linear programming)
Existen decisiones que no pueden ser representadas de forma adecuada mediante variables continuas. Por ejemplo, las decisiones de inversión son variables discretas (por ejemplo, planificación de la expansión de la generación o de la red, adquisición de equipos singulares, contratación de personas) o binarias (como localización de plantas o almacenes). Los problemas lineales con variables enteras se pueden clasificar en: programación entera pura PIP (pure integer
programming) si todas las variables son enteras, programación entera binaria BIP (binary integer programming) si todas son binarias o programación lineal entera mixta MIP (mixed integer programming) si algunas son enteras o binarias y el resto continuas. Un caso particular, pero muy frecuente, de variables enteras son las variables binarias (0/1), ya que permiten modelar condiciones de asignación o condiciones lógicas. Por otra parte, toda variable entera x se puede expresar como suma de
i siendo u una cota superior de x ,
. Existen algunos tipos de problemas de optimización que alteran ligeramente este esquema:
0 ≤ x ≤ u , y estando u comprendida en el intervalo
• sistemas de ecuaciones lineales – no lineales No existe una función objetivo como tal. Únicamente interesa encontrar una solución factible a un problema con un conjunto de restricciones.
• optimización sin restricciones Se trata de encontrar el conjunto de valores de las variables que determinan el mínimo/máximo de una función. Algunas de las técnicas que se verán en programación no lineal son para optimización sin restricciones.
• optimización multiobjetivo Existe más de una función objetivo. El problema que se plantea es cómo tratar varias funciones objetivo a la vez, teniendo en cuenta que el óptimo para un objetivo no lo es para otro, son objetivos en conflicto entre sí. Ésta se enmarca dentro de lo que se conoce de forma más general como decisión multicriterio (multicriteria decision making MCDM).
La formulación matemática de algunos problemas de optimización especiales por no incluir alguno de los componentes se presenta en la siguiente tabla.
Problema mixto complementario
(mixed complementarity problem)
Ajuste no lineal mínimo cuadrático
(multiobjective programming)
I.1.2. Referencias
Gass, S.L. and Harris, C.M. (eds.) (2001) Encyclopedia of Operations Research and Management Science. Centennial Edition. Kluwer Academic Publishers.
Hillier, F.S., Lieberman, G.J. (2002) Investigación de Operaciones. 7ª edición. McGraw Hill.
Robinson R. (1999) “Welcome to OR Territory” OR/MS Today pp. 40-43 August.
Sarabia, A. (1996) La Investigación Operativa. UPCO.
Taha, H.A. (1998) Investigación de operaciones. Una introducción. Prentice Hall.
Winston, W.L. (1994) Investigación de Operaciones. Aplicaciones y Algoritmos. Grupo Editorial Iberoamericana.
Modelo. Esquema teórico, generalmente en forma matemática, de un sistema o de una realidad compleja (por ejemplo, la evolución económica de un país), que se elabora para facilitar su comprensión y el estudio de su comportamiento. DICCIONARIO DE LA LENGUA ESPAÑOLA. REAL ACADEMIA ESPAÑOLA.
Un modelo es una representación matemática simplificada de una realidad compleja. Modelar es la acción de construir un modelo, de encorsetar la realidad. Implica la relación entre dos figuras (no necesariamente encarnadas por personas únicas sino por equipos): el modelador (encargado de la especificación y desarrollo del modelo) y el experto sobre la realidad (conocedor del problema real). La mayoría de las veces, el desarrollo de un modelo puede involucrar a un equipo multidisciplinar compuesto por matemáticos, estadísticos, ingenieros, economistas, psicólogos, etc. que aportan diferentes perspectivas y conocimiento en la representación de la realidad. Un modelo debe equilibrar la necesidad de contemplar todos los detalles con la factibilidad de encontrar técnicas de solución adecuadas. Un modelo es, en definitiva, una herramienta de ayuda a la toma de decisiones. Por esta razón, sus resultados deben ser inteligibles y útiles. Modelar se puede entender simultáneamente como ciencia y como arte. Es una ciencia pues se basa en un conjunto de procesos estructurados: análisis y detección de las relaciones entre los datos, establecimiento de suposiciones y aproximaciones en la representación de los problemas, desarrollo o uso de algoritmos específicos de solución. Es un arte porque materializa una visión o interpretación de la realidad no siempre de manera unívoca. Cada persona imprime su estilo en el modelo mismo y en la especificación, en el desarrollo y en la documentación. Características tales como elegancia o simplicidad pueden atribuirse a un modelo. El desarrollo de un modelo es una creación hecha con ayuda de ciencias básicas o herramientas de apoyo. Entre los beneficios explícitos o implícitos, tanto para el modelador como para el experto, derivados del proceso de modelado además del modelo en sí mismo, se pueden mencionar:
Ayuda a establecer un diálogo con intercambio de información entre el modelador y el experto
I.2 MODELOS DE OPTIMIZACIÓN
• Organiza los datos, la información disponible sobre el sistema
• Organiza, estructura y mejora la comprensión del sistema
• Internaliza la estructura organizativa de la empresa
• Permite compartir supuestos y resultados entre el modelador y el experto
• Proporciona un entorno ágil para el análisis y la sensibilidad
• Indica la dirección de mejora en las decisiones
En este capítulo se tratará exclusivamente de modelos de optimización, es decir, aquellos donde existe un conjunto de variables de decisión que deben maximizar/minimizar una función objetivo sometidas a un conjunto de restricciones. Los modelos de programación lineal son más utilizados que todos los otros tipos de optimización juntos y abarcan cualquier tipo de actividad humana como micro y macroeconomía, finanzas, marketing, economía de la energía, organización de la producción, planificación de la operación, selección de procesos, asignación de tareas, ingeniería química, forestal, agrónoma, comercio internacional, desarrollo económico, etc. Como referencias generales de modelado de problemas de optimización que se pueden utilizar en la enseñanza de pregrado o postgrado cabe citar a [Schrage, 1997] y [Williams, 1999].
Las etapas que componen el ciclo de vida de un modelo son las siguientes:
Consiste en la recolección y análisis de la información relevante para el problema, en el intercambio de información entre el modelador y el experto, en establecer una relación simbiótica y una estrecha coordinación entre ambos. Los problemas reales suelen estar definidos en términos vagos e imprecisos. Se debe hacer la tarea de traducción o interpretación en frases precisas, convertibles en ecuaciones matemáticas. En esta etapa se establecen y documentan los supuestos realizados que en etapas posteriores deberán ser validados. Esta etapa es fundamental para que las soluciones proporcionadas, las conclusiones obtenidas sean útiles, las decisiones adoptadas sean correctas. Los datos suelen ser vitales para conseguir un realismo o aplicabilidad en las soluciones. A menudo representan el cuello de botella del proceso de modelado.
Especificación matemática y formulación
Escritura matemática del problema de optimización, definiendo sus variables, sus ecuaciones, su función objetivo, sus parámetros. En esta etapa se analiza el tamaño del problema, la estructura de la matriz de restricciones, su tipo (LP,
MIP, NLP). Es una etapa de creación donde se debe prestar especial atención a
la precisión en la formulación y a la escritura de las ecuaciones que describen el
problema. En LP la elección de una formulación de un problema, aunque importante,
afecta de manera significativa la resolución del mismo. Sin embargo, en NLP
MIP la elección de la formulación es crucial. Pueden existir diversas
alternativas de modelado que afectan de manera fundamental en la resolución del mismo, existiendo un desarrollo cada vez mayor en la reformulación de problemas. En problemas MIP la calidad de una formulación se mide por la
cercanía entre la envoltura convexa del poliedro de soluciones enteras factibles y
la del poliedro del problema MIP relajado linealmente. En el apartado I.6.5.2 se
explica en más detalle algunas técnicas de reformulación de problemas MIP. La caracterización de un problema LP según su tamaño resulta difícil y ha sufrido un gran cambio desde los recientes desarrollos de algoritmos simplex
mejorados y, sobre todo, desde la aparición de los métodos de punto interior. En
la tabla 1.1 se propone una clasificación de tipos de problemas LP según su
tamaño. Esta clasificación debe ser tomada como guía o referencia relativa actual pero téngase en cuenta que los tamaños relativos de los problemas cambiarán conforme evolucionen los códigos de optimización. Actualmente se puede afirmar que los códigos de optimización lineal implantan algoritmos muy eficientes, son fiables y numéricamente robustos y están ampliamente disponibles.
Muy gran tamaño
Tabla 1.1 Tipos de problemas LP según su tamaño.
En lo referente a MIP o NLP ni siquiera se pueden dar criterios generales de tamaño ya que la dificultad de resolución no tiene por qué estar ligada al tamaño del problema, puede ser incluso preferible reformular un problema aunque aumenten las dimensiones, para lograr una resolución más eficiente.
Se trata de implantar un algoritmo de obtención de la solución numérica (muy próxima a la matemática) óptima o cuasióptima. El algoritmo puede ser de propósito general (método simplex) o específico. Puede haber diferentes
métodos de solución de un problema o diferentes implantaciones de un mismo método. El tiempo de resolución de un problema también puede depender drásticamente de cómo esté formulado. La solución óptima debe ser suficientemente satisfactoria, debe ser una guía de actuación para el experto.
Verificación, validación y refinamiento
Esta etapa conlleva la eliminación de los errores en la codificación, es decir, conseguir que el modelo haga lo que se ha especificado matemáticamente en la etapa anterior mediante su escritura en un lenguaje informático (depurar y verificar). Es necesario comprobar la validez de las simplificaciones realizadas a través de los resultados obtenidos, incluso contrastando éstos con situaciones reales ya transcurridas (validar) o comprobando que los resultados son coherentes con respecto a lo que sucedería en la realidad. Esta etapa de verificación, validación y comprobación da lugar a nuevas necesidades de refinamiento en el modelado para mejorar la capacidad de representación del sistema. Por ejemplo, eliminar la linealidad y hacer el modelo no lineal o hacer el modelo estocástico si la realidad lo fuera. Además, también se puede abordar el refinamiento matemático en la formulación del problema para hacerla más eficaz.
Esta etapa consiste en proponer soluciones. Permite conocer en detalle el comportamiento del modelo al hacer un análisis de sensibilidad en los parámetros de entrada, estudiar diferentes escenarios plausibles de los parámetros, detectar soluciones alternativas cuasióptimas pero suficientemente atractivas, comprobar la robustez de la solución óptima.
Implantación, documentación y mantenimiento
Ésta es una etapa fundamental del desarrollo de un modelo para garantizar su amplia difusión. La documentación ha de ser clara, precisa y completa. El manual de usuario debe incluir la especificación técnica funcional, matemática e informática. El propio código debe incluir una buena documentación para facilitar la tarea del mantenimiento. Piénsese que la mayor parte del ciclo de vida de un modelo no está en el desarrollo sino en la fase de uso y mantenimiento. En esta etapa se incluye también la tarea de formación para los usuarios del modelo.
Schrage, L. (1997) Optimization Modeling with LINDO. Duxbury Press.
Williams, H.P. (1999) Model Building in Mathematical Programming. 4th Edition. John Wiley and Sons.
A continuación se presentan algunos problemas característicos de
programación lineal y entera. Éstos se utilizan como referencia y clasificación para otros problemas. En particular, para los problemas enteros existen numerosas referencias de investigación dedicadas a la solución de los mismos. A pesar de la enorme atención que se ha dedicado a su solución su importancia práctica es limitada.
El problema por excelencia de programación lineal es el de asignación óptima de recursos. Un caso particular de éste es el denominado problema de la dieta. Consiste en determinar la composición de la dieta de mínimo coste que satisface las necesidades específicas de nutrientes. Pongamos un caso particular muy sencillo de alimentación de ganado bovino. Aprovechamos este ejemplo para seguir paso a paso las etapas en el desarrollo de un modelo.
Se ha determinado que las necesidades mínimas diarias en la alimentación de
ternera son de 700 g de proteínas, 28 g de calcio y 150 mg de vitaminas.
alimentos disponibles son pienso y forraje con un coste unitario de 0.30 y
0.35 €/kg respectivamente. La composición nutritiva por kg de alimento se muestra en la siguiente tabla.
I.3 FORMULACIÓN DE PROBLEMAS DE OPTIMIZACIÓN
Para ello analizamos y organizamos los datos del problema. Sean i los alimentos disponibles (pienso y forraje) y sean j los nutrientes (proteínas,
b la cantidad mínima diaria requerida de cada
calcio y vitaminas). Sea
nutriente. Sea
c el coste unitario de
x la cantidad
cada alimento. A continuación definimos las variables. Sea
diaria en kg de cada alimento i . Además indicamos la función objetivo y las restricciones del problema. La función objetivo es la minimización del coste diario de la dieta
correspondiente a los valores de la tabla dada. Sea
a la cantidad de nutriente por kg de alimento
• Después viene la resolución. Vamos a resolver gráficamente el problema. Para ello se dibujan las ecuaciones en forma de igualdad en el espacio de las variables y se indica la región factible del problema. Es decir, el conjunto de puntos que cumple todas las restricciones. Se traza la recta de la función objetivo para un valor cualquiera y se desplaza paralela a sí misma en el sentido de minimizar dicho valor hasta el último punto de la región factible. Dicho punto será el óptimo del problema.
• Seguidamente se realiza la interpretación y análisis de los resultados. Los resultados indican que la decisión óptima es comprar 18.83 kg de pienso y 8.33 kg de forraje cada día. Con estas decisiones el coste diario de los alimentos es de 6.1667 €. Al ganadero le ha llegado una oferta de otro fabricante de piensos a un precio de 0.25 €/kg pero con menor contenido en calcio, 1.5 g de calcio por kg de pienso, y tiene interés en analizar si le interesa comprar o no a dicho fabricante. Para ello planteamos este nuevo problema de optimización
x 1 ≥
La solución óptima para este nuevo problema es comprar 14.93 kg de pienso y 5.6 kg de forraje diariamente con un coste de 5.6933 €. Luego, esta oferta es atractiva económicamente.
Se trata de minimizar el coste total de transporte de un cierto producto desde los diferentes orígenes a los destinos, satisfaciendo la demanda de cada destino sin superar la oferta disponible en cada origen (ver ejemplo en I.4.2.1 Ejemplo de transporte). Se supone que todos los m orígenes están conectados
con los todos los n destinos. Sea
a i la oferta de producto en el origen i ,
demanda de producto en el destino j y c el coste unitario de transporte desde el origen i al destino j .
El problema de optimización consiste en determinar las unidades de
transportadas desde i hasta j , ∀i, j , que minimizan los costes
de transporte sujeto a las restricciones de oferta disponible en cada origen i (m restricciones de oferta) y demanda en cada destino j (n restricciones de demanda)
Implícitamente en esta formulación, se supone que la oferta del producto es
igual a la demanda del mismo
un sumidero universal con coste nulo. Si
universal conectada con todos los destinos con coste muy elevado. La estructura que presenta la matriz de restricciones del problema tiene el siguiente aspecto.
se añade una fuente
como las demandas de los productos
b son números
enteros, entonces el valor óptimo de
es entero por ser la matriz totalmente
unimodular 3 , por lo que no se necesita recurrir a métodos específicos de resolución de problemas de programación entera.
Consiste en determinar en una red con n nodos las cantidades óptimas para llevar unidades de un producto desde sus orígenes a sus destinos pasando por puntos de transbordo intermedios.
0 . El coste unitario de
cada transbordo ni genera ni consume unidades
transporte desde el origen i hasta el destino j en dicho sentido es c .
Hay que determinar las unidades de producto transportadas desde i a j ,
0 , ∀i, j , que minimizan los costes de transporte teniendo en cuenta la
restricción de balance o conservación del flujo en cada nudo i .
Implícitamente en esta formulación, se supone que la oferta es igual a la
demanda del producto, es decir,
Se trata de asignar la realización de n tareas a n personas (máquinas, etc.). Este problema es un caso particular del problema de transporte. Por consiguiente las variables toman valores enteros sin exigir esta condición en la formulación del problema. Consiste en minimizar el coste total de realizar las tareas sabiendo que cada tarea i debe ser hecha por una sola persona y cada persona j debe realizar una
única tarea, siendo
el coste de realizar la tarea i por la persona j . Las
variables del problema son x
1 si se asigna la tarea
3 Una matriz es totalmente unimodular si toda submatriz cuadrada tiene determinante 0, 1 ó –1. Si la matriz de un problema lineal es totalmente unimodular y las cotas de las restricciones son enteras, entonces todos los puntos extremos del poliedro tienen coordenadas enteras (se denomina politopo entero).
Los problemas de programación entera surgen en numerosos ámbitos de decisión.
• Por ejemplo, cuando se necesita representar número entero de productos o unidades enteras o indivisibles de recursos (aviones, personas, máquinas, etc.).
• También cuando se quieren imponer restricciones lógicas (si se fabrica el producto A también debe fabricarse el B y el C).
• Otro tercer ámbito son los denominados problemas combinatoriales, como por ejemplo los de secuenciación de tareas en máquinas, equilibrado de líneas de producción, asignación de tripulaciones, localización de almacenes o factorías, programación temporal de actividades.
• Así mismo surgen cuando hay funciones de naturaleza no lineal, como las de coste fijo o no convexas en general.
• Por último, también aparecen problemas enteros en teoría de grafos, por ejemplo el problema del tetracoloreado de un mapa.
En este apartado se presentan algunas técnicas de modelado que facilitan la formulación de problemas de optimización con variables enteras. No existe una manera sistemática de formular este tipo de problemas y plantear un “buen” modelado es, a menudo, un arte. Por esta razón, el método de aprendizaje se basa en la realización de ejemplos que muestren las posibles formulaciones de problemas característicos que aparecen con más frecuencia de lo que podría aparecer por su nombre. Algunos de los mencionados anteriormente se describen en este apartado y otro aparecen como problemas al final del capítulo. En los problemas que siguen, la envoltura convexa definida para el conjunto de soluciones ya no tiene soluciones enteras en los vértices como sucedía en los del apartado anterior, por lo que para su resolución hay que acudir a
procedimientos específicos de programación entera que se exponen en el apartado I.6.
Se trata de maximizar el valor total de la elección de un conjunto de n
coste de cada proyecto j respectivamente. El nombre procede de la decisión que toma un montañero que trata de maximizar el valor de lo que introduce en su mochila con una restricción de máximo peso admisible. Las variables del
proyectos sin sobrepasar el presupuesto b disponible, siendo
. Ésta es una utilización habitual
de las variables binarias como forma de seleccionar una alternativa, un proyecto
en este caso. La formulación del problema es la siguiente
Problema de recubrimiento (set covering)
características y n
combinaciones (subconjuntos) de dichas
características. La elección de una combinación implica realizar todas las
características de la misma. Se trata de
combinaciones elegidas de manera que se cubra o posea cada característica i al
menos una vez. Los datos son c el coste de elegir la combinación j y la matriz
minimizar el coste total de las
a cada combinación j ,
. El problema se formula de la siguiente
si se asigna la secuencia j
j = 1, … ,12
min 2xxxxxxxxxx+++++++++3 4 6 7 5 7 8 9 9 + 8x + 9x
+++ ≥
(SF-LA)
1 (SF-Denver)
(SF-Seatlle)
Las soluciones óptimas son xxx==
= 1 y el resto 0 ó xxx==
y el resto 0, ambas con coste 18 millones de €. Si no se permite que una tripulación fuera de servicio vuele en un avión las restricciones de cobertura de mayor o igual pasan a ser de igualdad. Luego, se trata de un problema de partición, cuya formulación se verá a continuación.
Problema de empaquetado (set packing)
Se tienen que realizar m proyectos divididos en n paquetes. La elección de un paquete implica realizar todos los proyectos del mismo. Se trata de maximizar el beneficio total de manera que cada proyecto i del conjunto de
todos los paquetes que lo incluyen no pueda ser elegido más de una vez.
beneficio de elegir el paquete j , la matriz de pertenencia de cada proyecto i a
cada paquete j
j 
. La formulación del problema es la siguiente
Problema de partición (set partitioning)
Problema del viajante de comercio (Traveling Salesman Problem TSP)
El problema consiste en hacer un recorrido que pase por n ciudades sin repetir ninguna y volviendo a la ciudad de partida de manera que la distancia (o tiempo o coste) total sea mínima. Es un problema de asignación pero con la condición de que la asignación sea un ciclo. Es uno de los problemas más importantes en la historia de la programación matemática por todas las investigaciones a las que ha dado lugar y por todas las aplicaciones que tiene, tanto directamente o apareciendo como subproblema dentro de otros más complejos. En una noticia de Optima (Mathematical Programming Society Newsletter) de junio de 1998, mencionaba que se había conseguido resolver un problema del viajante con 13509 ciudades. Los problemas de enrutamiento de vehículos (expedición o recogida de mercancías) pueden ser formulados de esta manera. Una de las características más interesantes de este problema es que existen muchas formulaciones conocidas para el mismo, ver [Williams, 1999] y [Nemhauser, 1999]. Dos de ellas se presentan a continuación, una tercera se
presenta en el apartado I.4.2.4 y se formulan en GAMS. Sea entre las ciudades i y j .
Formulación 1 (clásica): Se definen las variables
x ij ∈ {
Un ) ≤−
La primera restricción indica que a una ciudad j sólo se puede llegar una vez desde cualquier ciudad i . La segunda dice que desde una ciudad i sólo se puede salir una vez a cualquier otra ciudad j . Con esta formulación el problema se resuelve iterativamente. Primero se ignora el tercer tipo de restricciones. Se analiza la solución y se determina si contiene subciclos. En tal caso, se introducen restricciones para evitar subciclos de longitud U con la estructura del tercer bloque y se resuelve de nuevo el problema. Estos dos últimos pasos se repiten hasta encontrar la solución óptima. El número de restricciones del tercer tipo crece exponencialmente ya que el
número de subciclos de longitud U es
y siempre que aparece
un subciclo de longitud U aparece otro de longitud n −U . Por esta razón, su adición en la formulación se hace específicamente para los subciclos que van apareciendo y, usualmente, se necesita añadir solamente una pequeña fracción de las mismas.
Formulación 2: Se definen las variables
  1
El problema queda
ijk , ,
jrk + 1
∀ jk ,
La primera restricción indica que desde una ciudad i sólo se puede salir una vez a cualquier otra ciudad j en cualquier etapa k . La segunda dice que a una ciudad j sólo se puede llegar una vez desde cualquier ciudad i en cualquier etapa k . La tercera muestra que en la etapa k sólo se puede ir una vez de una ciudad cualquiera i a otra cualquiera j . La última indica que si en la etapa k estamos en la ciudad j habiendo venido de otra ciudad cualquiera i , en la siguiente etapa k + 1 saldremos de la ciudad j a otra ciudad cualquiera r . Esta última restricción es necesaria para establecer la continuidad espacial entre etapas sucesivas. La pregunta que cabe plantearse, es cuál de las dos formulaciones es mejor. En este caso, las formulaciones desde un aspecto teórico no son comparables,
aunque en la práctica resulta preferible la primera ya que tiene menos variables
(la primera tiene del orden de
n ), aunque tiene muchas más
restricciones. Las comparaciones teóricas de distintas formulaciones se presentan más adelante, dentro de lo que se denomina reformulación de problemas en programación lineal entera.
Problema de coste fijo
  
≤ My , j = 1, … , n , siendo M una
constante, cota superior de x , cuyo valor dependerá del problema, se distingue entre no realizar la actividad y realizarla al menos infinitesimalmente. El valor de la constante M debe ser el menor posible ya que esto es computacionalmente beneficioso. El problema lineal entero se formula como sigue
x . Introduciendo la condición

References: resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución