Source: http://docplayer.es/954559-1-1-52-1-8675-0-035-6-5-6.html
Timestamp: 2018-04-23 06:14:42+00:00

Document:
,1*(1,(52,1'8675,$/(0%$ (035(6$5,$/(6 - PDF
,1*(1,(52,1'8675,$/(0%$ (035(6$5,$/(6
Download ",1*(1,(52,1'8675,$/(0%$ (035(6$5,$/(6"
Patricia Contreras Gallego
1 7(6,6'2&725$/ 3/$1,),&$&,Ï1<*(67,Ï1'(,19(17$5,26 (181(172512'()$%5,&$&,Ï1&21 '(6(&+2<5(352&(62(672&È67,&2-26e0$5Ë$%(5*(6$6Ë1,1*(1,(52,1'8675,$/(0%$ '(3$57$0(172'(25*$1,=$&,Ï1'( (035(6$6 )$&8/7$''(&,(1&,$6(&21Ï0,&$6< (035(6$5,$/(6
3 25*$1,=$&,Ï1'((035(6$6 )$&8/7$''(&,(1&,$6 (&21Ï0,&$6<(035(6$5,$/(6 3/$1,),&$&,Ï1<*(67,Ï1'(,19(17$5,26 (181(172512'()$%5,&$&,Ï1&21 '(6(&+2<5(352&(62(672&È67,&2 Director: Dr. Javier Conde Collado Autor: D. José María Berges Asín Ingeniero Industrial / E.M.B.A.
5 'HGLFDRULD: A mi amigo Antonio, que me animó hasta el fin de sus días.
7 $*5$'(&,0,(1726: A todos aquellos profesores y personal adjunto de las enseñanzas nocturnas y a distancia, que permiten que miles de personas podamos seguir creciendo.
9 Ë1',&( ÍNDICE DE CONTENIDOS ANEXOS SÍMBOLOS, ABREVIATURAS Y SIGLAS ÍNDICE DE FIGURAS ÍNDICE DE GRÁFICOS Ë1',&('(&217(1,'26 &$3,78/2,1752'8&&,Ï INTRODUCCIÓN DIFERENCIA ENTRE CERTIDUMBRE, INCERTIDUMBRE Y RIESGO OBJETIVO DECISIONES BAJO INCERTIDUMBRE INTRODUCCIÓN OPTIMIZACIÓN DETERMINISTA OPTIMIZACIÓN ESTOCÁSTICA REPRESENTACIÓN EN ÁRBOL Y MATRICIAL MÉTODO DE MONTE CARLO BREVE RESEÑA HISTÓRICA ANTECEDENTES PLANIFICACIÓN Y GESTIÓN DE INVENTARIOS PLANEACIÓN AGREGADA Y PROGRAMACIÓN CONCLUSIONES METODOLOGÍA DE TRABAJO MODELOS 32
10 CASOS INDICADORES 42 &$3Ë78/2'(),1,&,21(635(9,$ INCERTIDUMBRE EN PRODUCCIÓN INVENTARIOS EL PLANIFICADOR POSICIÓN DE INVENTARIO INVENTARIO DE PIEZAS DE REPROCESO 50 &$3,78/202'(/26(5,(3$5$/(/ MODELO SERIE PARALELO PROCESOS SERIE CASO PROCESO SERIE PROCESO PARALELO CASO PROCESO PARALELO ACONDICIONAMIENTO FLUJO TEMPORAL REPROCESO EN PARALELO FLUJO DE SALIDA RETARDO TRANSITORIO DE SALIDA CASO FLUJO TEMPORAL PARALELO 64 &$3,78/2'(),1,&,21'(275$69$5,$%/( DEFINICIÓN DE OTRAS VARIABLES TIEMPO TOTAL DE PROCESO O LEAD TIME 70
11 4.3.- DEFINICIÓN DE INVENTARIO AJUSTADO ROTACIÓN SISTEMAS CON DESECHO Y REPROCESO TIEMPO TOTAL DE PROCESO CONSTANCIA INVENTARIO AJUSTADO POR ROTACIÓN CASO DESARROLLO APLICACIÓN INVENTARIO AJUSTADO 73 &$3Ë78/2$1È/,6,63/$1,),&$&,Ï1'(/$352'8&&,Ï INTRODUCCIÓN CASO PLANIFICACIÓN ANÁLISIS CAUSAL CORTO, MEDIO Y LARGO PLAZO MÉTRICAS PARA EL PLANIFICADOR BALANCE DE INVENTARIOS RATIOS CUENTA DE INVENTARIOS CAPACIDAD FUTURA Y AMORTIZACIÓN CUENTA DE RESULTADOS BALANCE INVENTARIO DE MANIOBRA ÁNALISIS DE SENSIBILIDAD VARIACIÓN DEMANDA DESARROLLO CASO ANÁLISIS DE SENSIBILIDAD EVALUACIÓN DE LA INCERTIDUMBRE TRADICIONAL INTRODUCCIÓN ESTIMACIÓN POR INTERVALOS DE CONFIANZA CASO 95
12 LÍNEA DE TENDENCIA 96 &$3Ë78/23/$1,),&$&,Ï1%$-2,1&(57,'80%5( DEFINICIÓN CAPACIDAD FUTURA DISPONIBLE GRÁFICA DE RENDIMIENTOS GRÁFICA DE RENDIMIENTOS MONÓTONOS DEFINICIÓN DE VARIABLES MODELO REPARTO CAPACIDAD EN EL PROCESO INICIAL FUNCIÓN OBJETIVO RESTRICCIONES SOFTARE CASO OPTIMIZACIÓN PROGRAMACIÓN LINEAL ESCENARIO ESCENARIO ESCENARIO OPTIMIZACIÓN ESTOCÁSTICA VALOR ESPERADO DE LA INFORMACIÓN PERFECTA RENDIMIENTOS MEDIOS VALOR DE LA SOLUCIÓN ESTOCÁSTICA CASO DOS PROGRAMACIÓN EN SOLVER PROGRAMACIÓN EN GAMS 135
13 ÁRBOL MATRIZ MODIFICACIÓN DE RENDIMIENTO PLANIFICACIÓN SIN SERIES HISTÓRICAS, MÉTODO DE MONTE CARLO DESARROLLO MEDIANTE RISK SOLVER PLATFORM 148 &$3Ë78/2(5525'((67,0$&,Ï INTRODUCCIÓN MINIMIZACIÓN DEL ERROR DE ESTIMACIÓN FUNCIÓN OBJETIVO Y RESTRICCIONES CASO: DESARROLLO MEDIANTE RISK SOLVER PLATFORM 157 &$3Ë78/25(680(1<&21&/86,21(6 &$3Ë78/2)8785$6/Ë1($6'(,19(67,*$&,Ï1 %,%/,2*5$)Ë$ 5()(5(1&,$6(1,17(51(7 š
14 ANEXO 1 $1(;26 INDICADORES LOGÍSTICOS MÁS UTILIZADOS 195 ANEXO 2 PROGRAMAS SOLVER 199 PROGRAMAS GAMS 270 š
15 /,67$'(6Ë0%2/26$%5(9,$785$6<6,*/$6 &DSLXOR, : posición de inventario / : tiempo de proceso / : tiempo de proceso línea principal / : tiempo de proceso línea reproceso 6 : stock de productos comercializables al inicio del periodo t. 6 : el inventario al final del periodo anterior ) : demanda insatisfecha en el periodo t ) : demanda en el periodo t-1 4 : flujo de productos del proceso del periodo t 4 : flujo de productos de la línea principal del periodo 5 : flujo de la línea de reproceso del periodo 6 : inventario de reproceso en t 6 : inventario de reproceso en t 1 $5 : reproceso generado por la línea principal en el periodo t G : desecho de la línea de reproceso en el periodo t &DSLXOR 0L: i=1,2, M : talleres GL : desecho del taller i 3 : producción del taller i $5 : piezas enviadas a reproceso generado por el taller 2 $5 : piezas enviadas a reproceso generado por el taller 2 š
16 DU : porcentaje de reproceso del taller 2 DU : porcentaje de reproceso del taller 2 4: Piezas que entran al bucle de reproceso desde la línea principal $5 : Tanto por uno de reproceso de la línea paralela G : Tanto por uno de desecho de la línea paralela M: número de ciclos ;: tanto por uno de piezas que salen del sistema en el ciclo j 3 : número de piezas buenas que salen de la línea de reproceso en el ciclo j Q: número de periodos elementales., 2 : inventario proceso paralelo 7: número de periodos elementales en procesar el lote Q &DSLXOR % ', 3 piezas netas que el taller i puede producir en el periodo n. : número de órdenes recibidas al taller i en el tiempo n. : inventario a la salida de cada taller i al principio del tiempo n. :producción del taller i en el periodo n. /7 : Tiempo total de proceso de la planta /7 :Tiempo de proceso del taller i Q: periodo considerado 5 : rotación total del proceso o factoría 5 : rotación del inventario del taller i :Rendimiento del taller i 3 : Producción neta que sale del proceso en el periodo n,d : inventario ajustado en la sección i š
17 # : Inventario total ajustado &DSLXOR <: variación porcentual demanda respecto producción periodo anterior. Q: tamaño de la muestra o periodos considerados. *P: gap o desviación media de la muestra 2L: ordenes fabricación solicitadas al taller P1, i = 1, 2,3... n 3L: material real fabricado en cada periodo i, i = 1, 2,3... n. : desviación típica muestra distribución normal 6: cuasidesviación típica muestral, Student &DSLXOR 35 "! : Probabilidad del escenario sc L : Tipo de molde que identifica el modelo N: Fila o filas de moldes o máquinas que son susceptibles de montar moldes de tipo i Q: número máximo de moldes susceptibles de montarse del tipo k en nuestras instalaciones. : Periodo considerado : Rendimiento de la pieza # & : Proporción de moldes del producto i en la máquina k en el periodo t. : máxima rotación o número de piezas suministradas del molde i que puede ser producido en la máquina k en el periodo t. / :máxima capacidad de la máquina k en el periodo t. # D :capacidad necesaria por unidad del molde i de la máquina k en el periodo t. Es el tanto por uno de utilización del molde en función de las necesidades. š
18 E # : es la pérdida de capacidad de la máquina k, si es ajustada para colocar el molde i en el periodo t. # \ : variable binaria que toma el valor uno si la máquina k está disponible para producir la pieza i en el periodo t, igual a cero en caso contrario. # [ :es la producción del molde i en la máquina k en el periodo t. V : es el inventario en el almacén del producto i, al final del periodo t. L :es el inventario en curso ajustado del producto i, al final del periodo t. G : variable que representa retrasos, pedidos a corto y medio plazo del producto i en el momento t. S variable que representa pedidos largo plazo y previsiones de ventas en el periodo. t. F : costes de fabricación del modelo i # [ : producción del molde i en la máquina k en el periodo t. 35 : Probabilidad del escenario s, de la pieza i 6 : número máximo de escenarios S : penalización no servicio pieza i [$ : piezas no servidas del modelo i en el escenario s [%& : piezas almacenadas del modelo i en el escenario s D : coste de almacenamiento del modelo i S'( : precio de venta del modelo i ['( : piezas vendidas del modelo i en el escenario s : Rendimiento de la pieza i en el escenario s [$ : piezas no servidas del modelo i, en el escenario s S$ : piezas necesarias del modelo i. ['( : piezas vendidas del modelo i en el escenario s [%& : piezas almacenadas del modelo i en el escenario s ) : rendimientos medios š
19 (93,: Valor esperado de la información perfecta o arrepentimiento esperado, (([SHFHGYDOXHRISHUIHFLQIRUPDLRQ :6: valor esperado de la solución óptima, (ZDLDQGVHH) espera y observa. (9:3, : valor esperado con información perfecta (([SHFHGYDOXHZLKSHUIHF LQIRUPDLRQ). Otra denominación de S 53: problema estocástico con recurso (Recourse Problem) que tiene una solución óptima. 966 Valor de la solución estocástica (9DOXHRIKHVRFKDVLFVROXLRQ (9 solución del problema determinista del valor medio (H[SHFHGYDOXH(9). (99 Valor esperado de la solución del valor medio de los parámetros estocásticos. 966 Valor de la solución estocástica (YDOXHRIVRFKDVLFVROXLRQ966) 0 : punto de decisión &DSLXOR(UURUGHHVLPDFLyQ V denominación de escenario. V número total de escenarios. 35 : probabilidad de ocurrencia del escenario s. Q: número de elementos de la serie o tabla Q " : número de elementos del escenario s Q : número de orden en la serie de la pieza i Rendimiento de la pieza i en un periodo elemental t Rendimiento medio ) : Rendimiento medio de la pieza i en el escenario s HUURU error de estimación š
20 HUURUV error del escenario s (6&(1$5,26L variable binaria que vale 1 si el valor con número de orden L, pertenece al escenario 6. En caso contrario toma el valor 0. 6HOHFRUVL variable auxiliar que tomará el valor uno si el número de orden L de la tabla ordenada es menor o igual que el número de elementos del escenario V más el número de elementos de todos los escenarios precedentes. Si no es así tomara el valor cero. : rendimiento medio de un sub periodo. š
21 Ë1',&('(),*85$6 Figura 1.1 Árbol. Figura 1.2 Orden económica de pedido Figura 1.3 Probabilidad en la gestión de inventarios Figura 1.4 Talleres serie. Figura 1.5 Tamaño del lote Figura 1.6 Células de inspección Figura 1.7 Lote con reproceso Figura 1.8 Lote con reproceso, coste de arranque y de inventarios Figura 1.9 Planificación de la producción, clásico Figura 3.1 Diagrama proceso serie paralelo con desecho y reproceso. Figura 3.2 Diagrama proceso serie Figura 3.3 Flujo serie Figura 3.4 Diagrama proceso paralelo Figura 3.5 Flujo proceso paralelo Figura 5.1 Análisis causal Figura 5.2 Balance financiero. Figura 5.3 Balance inventarios pedidos Figura 5.4 Ejemplo Balance inventario, capacidad Pedidos Figura 5.5 Ejemplo Balance en t+1 inventario, capacidad Pedidos Figura 5.6 Fondo de maniobra financiero Figura 5.6 Inventario de maniobra Figura 6.1 Capacidad futura disponible š
22 Figura 6.2 Parámetros Solver Figura 6.3 Opciones de Solver Figura 6.4 Árbol de escenarios Figura 6.5 Matriz de escenarios Figura 6.6 Matriz de probabilidades Figura 6.9 Ejemplo distribución de probabilidades uniforme 1 Figura 6.10 Ejemplo distribución de probabilidades uniforme 2 Figura 6.11 Generación de escenarios de Monte Carlo 1 Figura 6.12 Generación de escenarios de Monte Carlo 2 Figura 7.1 Programación minimización error, variable binaria selector Figura 7.2 Minimización error, variable binaria escenario y rendimiento Figura 7.3 Programación minimización error, formulación error por escenario Figura 7.4 Caso programación minimización error, función objetivo y restricciones Figura 7.5 Caso programación minimización error, función objetivo y restricciones antes de procesar Figura 7.6 Caso programación minimización error: resultados Figura 7.7 Programación minimización error en Risk Solver cuatro escenarios, resultados š
23 Ë1',&('(*5È),&26 Gráfico 3.1 Flujo proceso serie Gráfico 3.2 Flujo proceso paralelo Gráfico 3.3 Reproceso paralelo primer lote Grafico 3.4 Gráfico flujo según rotaciones desde ciclo paralelo Gráfico 3.5 Gráfico flujo según ciclos elementales ciclo paralelo Gráfico 3.6 Ejemplo flujo de entregas de la línea serie paralelo por ciclos Gráfico 5.1 Balance pedidos-inventarios, órdenes Gráfico 5.2 Relación demanda órdenes producción-rotación Gráfico 5.3 Ejemplo desviaciones Gráfico 5.4 Caso línea de tendencia Gráfico 5.5 Ejemplo línea de desviaciones y línea de tendencia Gráfico 6.1 Caso temporal de rendimientos, fuente propia Gráfica 6.2 Caso de rendimientos monótonos Gráfico 6.3 Resultados escenarios programación lineal Gráfico 6.4 Resultados escenarios programación lineal estocástica Gráfico 6.5 Rendimientos medios š š
25 &$3,78/2,1752'8&&,Ï1,1752'8&&,Ï1 La planificación de la producción en procesos estocásticos es un problema difícil de resolver, cuyas desviaciones afectan a todas las áreas de la empresa: fallos en las entregas, desviaciones financieras, operaciones de marketing truncadas, ventas frustradas, afectando directamente a la imagen de la empresa... En definitiva, es clave para el éxito o supervivencia del negocio. La mayoría de las plantas de fabricación con procesos estocásticos se diseñan con modelos deterministas, dentro de un concepto estático, influyendo en las decisiones estratégicas, tácticas y operativas. Aunque facilita el diseño, la realidad posterior le aparta de las condiciones de diseño creando cuellos de botella o sobrecapacidades con altos costes financieros y dificultando las operaciones. No sólo el diseño de la planta sigue este patrón, la gestión de órdenes de fabricación, el cálculo de stocks o el cálculo de los plazos de entrega se realizan generalmente utilizando modelos deterministas. Sin embargo, muchas industrias, tienen procesos o parámetros no controlables y difícilmente cuantificables. El diseñador de plantas duda en el momento de prever capacidad de instalaciones o espacios requeridos, el planificador tiene enormes dificultades a la hora de gestionar los recursos para alcanzar o mantener los inventarios requeridos y al industrializador le es imposible determinar cuántos recursos tiene que emplear o cuándo va a alcanzar el nivel de inventario solicitado para el lanzamiento de un nuevo producto. Por ejemplo en la industria cerámica, Domenico Fortuna (2000), se pueden tener rechazos de hasta el 50% en piezas de complejidad elevada y niveles de reproceso que pueden alcanzar e incluso superar el 80%. Además, para gestionar, necesariamente tenemos que pronosticar y planificar, Jay Heizer, Barry Render (2003) define previsión como: es el arte y la ciencia de predecir acontecimientos futuros. Puede suponer la toma de datos históricos y su proyección hacia el futuro con algún tipo de modelo matemático, puede ser una predicción
26 subjetiva o intuitiva del futuro, o puede englobar una combinación de estás; es decir, un modelo matemático ajustado al buen juicio del directivo. La naturaleza de los procesos es dinámica y cualquier decisión se extiende en el tiempo condicionando la toma de decisiones posteriores. Esta naturaleza dinámica de los acontecimientos se extiende incluso en el quehacer más cercano, Steven Nahmias anota (2007): En la práctica, la programación de trabajos es un problema dinámico y define el término dinámico: implica que los trabajos llegan aleatoriamente con el tiempo y deben tomarse decisiones sobre la marcha en cuanto a cómo programar estos trabajos. Muchos autores han realizado estudios sobre planificación estocástica, aunque ninguno sobre el tema referido en esta tesis, sin embargo como se verá un poco más adelante autores de reconocido prestigio reconocen que los modelos actuales de planificación bajo incertidumbre son inapropiados o muy limitados con horizontes de planificación de varios periodos, por lo que se van a desarrollar planteamientos que sortean estas dificultades minimizando el problema de la falta de capacidad de computación disminuyendo el número de periodos y buscando indicadores óptimos que representen los datos seleccionados. La falta de precisión en los resultados mediante optimización estocástica a largo plazo, según pruebas de campo, se corrige mediante la formulación de indicadores causales como herramientas en soluciones heurísticas. Puesto que estamos en procesos con inventarios variables, como consecuencia de desechos y reproceso, a lo largo de la cadena de valor se definirán nuevas definiciones de inventario o tiempo de proceso que nos permita aplicar los conceptos utilizados habitualmente. ',)(5(1&,$(175(&(57,'80%5(,1&(57,'80%5(<5,(6*2 Certidumbre sería el caso en que cada acción da un resultado conocido e invariable. La incertidumbre es una característica intrínseca de una parte de la naturaleza, siendo común a todos los actores, el riesgo desde un punto de vista subjetivo es específico para cada actor. Así, como ejemplo, el rendimiento de una cosecha local en un país
27 subdesarrollado tiene un gran grado de incertidumbre observado por cualquier actor, pues depende de la climatología, plagas etc., pero tiene un gran riesgo para el pequeño agricultor autónomo de ese país subdesarrollado y ningún riesgo para un funcionario de un país desarrollado lejano. Desde un punto de vista matemático se definiría como riesgo en el caso de que cada acción tenga asociado un posible resultado y este tiene una cierta probabilidad de ocurrencia. La certidumbre sería una situación de riesgo con probabilidades cero y uno. Podríamos definir incertidumbre como la expresión del grado de desconocimiento de una condición futura es decir se desconocen las probabilidades de los resultados. La gran mayoría de las veces, las decisiones implican riesgos. El agricultor del ejemplo mostrado previamente empleará su conocimiento, es decir, tendrá sus propias reglas heurísticas, valorando y calculando el riesgo antes de tomar una decisión. Sin embargo, estas reglas, que en el ejemplo del agricultor estarán basadas en la experiencia de varias generaciones, no existen o son difícilmente aplicables en muchos contextos empresariales, por su complejidad y por los sesgos cognitivos que influyen en la toma de decisiones. En este documento se eliminan estos sesgos simplificando a la vez el modelo para que con los medios actuales, las decisiones sean óptimas. El objetivo perseguido es optimizar el beneficio, pero el modelo tiene datos de entrada que tienen un cierto grado de incertidumbre. En este caso en concreto, la variable con incertidumbre es el rendimiento, de tal forma que dependiendo del dato de rendimiento que se tome, afectará sobre el beneficio. Uno de los primeros pasos que deberán plantearse es cuantificar la incertidumbre. Para ello, la serie histórica será ordenada en una serie monótona. Hacer notar que podría tomarse la media de un conjunto de valores de toda la muestra como único parámetro de rendimiento, tomando solo una parte de la población y efectuando un análisis de riesgos, el cual ayudaría en la decisión final a tomar. El estudio va a dar un paso más cuantificando la incertidumbre, tomando la totalidad de la población mediante parámetros que le representen, minimizando el error de esta
28 elección y realizando un proceso de matemático de optimización que busque la mejor solución posible teniendo en cuenta unos costes de computación aceptables. 2%-(7,92 Conscientes de que siempre se necesita el buen saber del gestor, nuestro objetivo es desarrollar las herramientas necesarias que ayuden al planificador a cumplir con su misión, en una planta con varios talleres serie - paralelo con desecho y reproceso, teniendo en cuenta la estocasticidad del proceso. '(&,6,21(6%$-2,1&(57,'80%5(,1752'8&&,Ï1 En cualquier actividad humana se da un proceso continuo de toma de decisiones, algunas de pequeña importancia y otras muchas, de gran transcendencia. El hombre siempre ha utilizado las herramientas que ha tenido a su alcance para intentar asegurar el buen juicio de sus elecciones mediante todo tipo de herramientas, formación específica, decisiones colegiadas, consejo de expertos, crítica constructiva, incluso muchas tradiciones, algunas convertidas en leyes, intentaban trasladar a generaciones futuras el buen saber hacer. Actualmente y bajo ciertas circunstancias somos capaces de trasladar algunos tipos de cuestiones al lenguaje matemático, que junto con el incremento de la capacidad de cálculo gracias al desarrollo de los computadores, facilitan la toma de decisiones en algunos casos bajo el paraguas de la optimización. Sin embargo, la realidad en muchas ocasiones es demasiado compleja y cambiante por lo que al intentar reflejarla en modelos matemáticos se convierte en problemas de grandes dimensiones teniendo parámetros desconocidos y estocásticos, además todo ello puedo ir acompañado de incertidumbre debido a la poca fiabilidad de algunos datos, enfatizado por la falta de datos referente a acontecimientos futuros.
29 Existe un excelente trabajo que recopila y sintetiza una gran cantidad de información sobre el tema de Optimización bajo incertidumbre, Ramos, A., Alonso-Ayuso, A., Perez, G. y otros (2008) en el libro del mismo título. En el prólogo dice actualmente no existen aplicaciones estándar o comerciales, potentes y fiables para resolver problemas estocásticos y en el capítulo uno cita: muchos de estos modelos presentan las siguientes características: Son de dimensiones enormes, son estocásticos, son enteros y son no lineales. Cada una de estas tres últimas características dificultan enormemente la resolución del problema, considerándose problemas de complejidad computacional elevada. Sin embargo, desde mediados del siglo XX (Dantzig 1, 1955, 1960, 1963, 1990, 1991, 1997) se reconoce el hecho de que la optimización determinista tradicional no siempre es la más adecuada para capturar el verdadero comportamiento aleatorio de la mayoría de las aplicaciones de la vida real. La principal razón estriba en que las aplicaciones están involucrados datos inciertos; que aparecen porque la información necesaria en etapas siguientes aún no está disponible en el momento en que la decisión ha de ser tomada. Para atacar este problema de optimización bajo parámetros no deterministas surge el desarrollo de la Optimización bajo incertidumbre o Programación Estocástica. Prékopa, (1995) define la programación estocástica como la ciencia que proporciona formulaciones y soluciones de optimización de los sistemas estocásticos, utilizando la técnica de programación matemática. La optimización trata de minimizar o maximizar una función que denominamos objetivo, la cual está compuesta por una o varias variables. Estas variables generalmente están restringidas por ecuaciones o desigualdades. 1 George Bernard Dantzig ( ) se le considera el padre de la programación lineal, durante su carrera académica mostró cómo explotar la generalidad del modelo de programación lineal: la formulación y solución de problemas de gestión en la mayoría de las principales industrias, la resolución de problemas estratégicos y tácticos en la defensa; la evaluación de planes y soluciones operativas de energía y otros problemas de recursos limitados, los avances teóricos en matemáticas, estadística, economía, la teoría de juegos, la informática y la adaptación del modelo de programación lineal a un conjunto más general de los problemas matemáticos de programación. También trabajo como consultor para una serie de grandes empresas siendo muy prolífico en artículos de investigación y libros de texto.
30 Los matemáticos dividen este tipo de problemas en dos ramas, los que denominan modelos deterministas y los modelos estocásticos. En el primer caso los parámetros asociados al modelo son conocidos con certeza absoluta, en los modelos estocásticos la totalidad o una parte de los parámetros son una distribución de probabilidad. 237,0,=$&,Ï1'(7(50,1,67$ Newton, Bernouilli, y Lagrange, ocupados en el desarrollo del cálculo infinitesimal, empezaron a estudiar hace más de 300 años la forma de obtener máximos y mínimos de determinadas funciones. Sin embargo, no fue hasta el siglo XIX cuando el matemático francés Fourier empezó a indagar en los métodos que ahora llamamos programación lineal. En el año 1939 el matemático ruso Kantarovith define por primera vez de una forma precisa lo que hoy llamamos programación lineal en su obra Métodos matemáticos de organización y planificación de la producción. En los años posteriores a la Segunda Guerra Mundial, EEUU llegó a la conclusión de que la coordinación de todos los recursos era un problema de tal complejidad que solo podría resolverse por métodos de optimización cuyo modelo resuelve la programación lineal. Este concepto junto con el desarrollo de los ordenadores hizo que un grupo de investigadores, Richard. Cottle, B. Curtis Eaves and Michael A. Saunders (2006), liderados por Dantzing en 1947 formaran un grupo, Scientific Computation of Optimum Programs, cuya primera aplicación fue el puente aéreo de Berlín tras su aislamiento por parte del bloque soviético. Después de este grupo de trabajo, se crean otros donde se desarrollan distintas aplicaciones y en 1951 el grupo dirigido nuevamente por Dantzig ayudado por modelos de ordenador se desarrolla el método simplex, siendo este el método más utilizado en la resolución de ecuaciones lineales. Podría definirse como programación lineal la maximización o minimización de algunas funciones, las cuales se encuentran sujetas a ciertas restricciones. Se caracteriza por funciones lineales de las incógnitas, el objetivo es lineal en las incógnitas y las
31 restricciones son igualdades o desigualdades lineales en las incógnitas. David E. Luenberger (1989). Se puede expresar como: A. Alonso-Ayudo y otros (2008): Z = minimizar c 1 x 1 + c 2 x c n x n Sujeto a a 11 x 1 + a 12 x a 1n x n = b 1 a 21 x 1 + a 22 x a 2n x n = b A m1 x 1 + a m2 x a mn x n = b m x 1, x 2,..., x n Y en forma matricial: Z = minimizar cx Sujeto a Ax = b x Donde x es un vector columna de decisiones de dimensión n de términos no negativos, c que representa los costos es un vector fila de dimensión n, A es la matriz de restricciones de dimensión m x n y b es el vector columna m dimensional de términos independientes, siendo c,a, b datos reales conocidos. Un valor optimo x * es una solución factible tal que cx F[ * para cualquier otra solución factible de x. Los programas lineales obtienen soluciones óptimas de mínimo coste bajo restricciones lineales de demanda o máximo beneficio bajo una situación de recursos limitados. Las restricciones definen los límites del sistema que se intenta optimizar. Para definir la restricción, primero deberá calcularse un valor de las variables de decisión, y después será definido un límite mediante los operadores <=, =, >=, sobre este valor calculado. Restricciones generales: Por ejemplo, en un reparto en porcentajes de un conjunto total la suma de estos porcentajes tiene que ser cien. Límites en las variables: las variables de decisión están limitadas a un rango de valores.
32 Política de restricciones: definidas por las políticas que se establezcan en la definición del problema. Por ejemplo, podría establecerse un límite mínimo al porcentaje de capacidad de producción reservado para un determinado producto. Limitaciones físicas: están determinadas por la naturaleza del problema. Por ejemplo, el número de productos a producir tiene que ser positivo (en este caso números negativos no tienen sentido). Limitaciones de entero: en algunos casos la variable de decisión tiene que ser un entero, por ejemplo en un problema cuya variable de salida sea el número de personas a ser operadas o volar en un determinado avión. Restricción variables binarias de decisión: serán utilizadas como variables de decisión, por ejemplo si ponemos una máquina en marcha o no en función de variables como los costes fijos y variables. Diversos estudios señalan que si un país utilizase métodos de programación lineal, su producto interior bruto aumentaría entre un 10 y un 15%. 237,0,=$&,Ï1(672&È67,&$ Siguiendo con la notación empleada en el apartado anterior se definen los problemas lineales estocásticos como problemas de optimización lineal en los que algunos de los parámetros c,a,b se consideran inciertos. Z = mín c x s.a. A x = b x GRQGH HVFDGDXQRGHORVHVFHQDULRV\ Hl conjunto de escenarios. Generalmente la toma de decisiones es función del tiempo, por lo que se define etapa A. Alonso-Ayudo y otros (2008) como un conjunto de periodos de tiempo en los que tiene lugar la realización de parámetros inciertos.
33 Sin embargo, la resolución del sistema de ecuaciones puede dar lugar a confusión pues puede haber soluciones en un escenario y no en el otro. Haciendo un breve resumen de una parte del trabajo de Andrés Ramos y Santiago Cerisola (2009), podría explicarse cómo: El método que recoge adecuadamente la incertidumbre, es el denominado análisis de escenarios El primer proceso puede resumirse como la elección de un conjunto finito de valores de los parámetros estocásticos que representen el conjunto de todos los valores posibles de los mismos. Los parámetros estadísticos son una consecuencia inevitable del propósito esencial de la estadística: modelar la realidad. Tradicionalmente se ha atacado el problema de la incertidumbre bajo varios planteamientos que pueden resumirse como sigue: a) Análisis de escenarios: 1.- Un problema es resuelto para cada escenario posible, es decir, el problema se resuelve independientemente para cada escenario, por lo que a priori tendremos una decisión para cada escenario. En la bibliografía a este tipo de resolución se la denomina espera y observa, análisis de escenarios o qué pasaría si (ZDLDQGVHHVFHQDULR DQDO\VLV R ZKD LI DQDO\VLV). Un caso particular es la resolución del escenario medio (H[SHFHGYDOXH(9). Con este tipo de resoluciones la solución de un escenario puede ser irrealizable para otro. 2.- Con las soluciones obtenidas de forma determinista y con criterios heurísticos, se busca la mejor solución, que difícilmente será la solución óptima para todos e incluso alguno de los escenarios. Generalmente se elegirá la solución que aparezca en más escenarios deterministas buscando la flexibilidad en el cambio de las decisiones a lo largo del tiempo mientras se resuelve la incertidumbre. b) Optimización estocástica: Las decisiones dependen únicamente de la información disponible en ese momento, teniendo en cuenta la distribución estocástica. Las decisiones de la primera etapa, en
Sistemas de distribución: avances en la gestión de inventarios
Curso 2004/05 CIENCIAS Y TECNOLOGÍAS/24 I.S.B.N.: 84-7756-663-1 BEATRIZ ABDUL-JALBAR BETANCOR Sistemas de distribución: avances en la gestión de inventarios Director JOAQUÍN SICILIA RODRÍGUEZ SOPORTES

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