ORT Simulacion Practica PDF
ORT Simulacion Practica PDF
ORT Simulacion Practica PDF
de Computación
1er.Cuatrimestre
2006
MATERIA:
Simulación
V1.0 - 11 -
PLANIFICACION de la MATERIA
1) OBJETIVOS GENERALES:
Se tratará que el alumno al finalizar de cursar la materia logre:
• Ensamblar los elementos aprendidos en las materias de matemáticas y los
conocimientos de programación a efectos de una mejor comprensión de los temas desarrollados
con especial énfasis en la técnica de simulación.
3) PRESUPUESTO de TIEMPO:
Horas
planificadas
Unidad
1 10
2 10
3 6
TOTAL 26
UNIDAD 1 –SIMULACIÓN.
Concepto. Simulación de juegos de azar. Simulación de juegos con variable
uniforme (random). Generación de números al azar. Generación de valores
de variable discreta. Generación de valores de variable continua. Casos
particulares: Uniforme – Normal – Exponencial – Binomial - Poisson. Uso de
planilla de cálculo (MsExcel TM) para simulación. Modelos sencillos de
simulación. Problemas de aplicación.
UNIDAD 2 –COLAS.
Una cola con un canal. Una cola con dos canales. Dos colas con dos
canales. Colas con prioridad. Colas limitadas. Colas sucesivas.
Trabajo de aplicación en MsExcel TM presentado por los alumnos.
UNIDAD 3 - STOCK
Stock: Modelo simple, con plazo de reposición, con costo
diferenciado, con limitaciones. Stock con perecederos.Problemas de
aplicación.
Trabajo de aplicación en MsExcel TM presentado por los alumnos.
5) RECURSOS DIDACTICOS:
Se alentará la modalidad de trabajo en equipo. Se fomentará el espíritu crítico sobre las
soluciones propuestas y los hábitos de orden, claridad y optimización en la resolución
de los problemas propuestos.
7) BIBLIOGRAFIA:
a) De lectura obligatoria.
• Apuntes y prácticos dictados por el curso.
• Guía de sitios en Internet.
• Administración de Producción y operaciones
• Chase; Aquino; Sacors Mc Graw Hill
• Administración de Operaciones. Estrategia y análisis
• Krajewski-Ritzman Prentice Hall
b) De consulta facultativa.
• Métodos de simulación.
•
Definiciones:
La simulación es una técnica que permite generar valores aleatorios imitando una situación
real en términos matemáticos.
Simulación: es una técnica numérica que a través de relaciones matemáticas y lógicas,
permiten describir el comportamiento sistemas complejos del mundo real
Simulación: es el proceso de diseñar y desarrollar un modelo de un sistema o proceso y
conducir experimentos con este modelo con el propósito de entender el comportamiento del
sistema o evaluar varias estrategias con las cuales se puede operar el sistema.
Conceptos necesarios
Variable aleatoria discreta y continua funciones de probabilidad y densidad. media variancia
intervalo de confianza. distribuciones particulares : Binomial. Geométrica Pascal, Binomial
Negativa. Poisson Exponencial Gamma. Normal . Uniforme. Función distribución
F(x)=P{X<=x}; siendo F(x) la función de distribución acumulada de la v.a.x.
X= 0 1 2 3
P(x)= 0,1 0,2 0,5 0,2
F(x)= 0,1 0,3 0,8 1
R x
0,2034 1
0,8849 3
0,6989 2
0,9808 3
0,0644 0
0,5093 2
Esta es la forma en que se puede simular una v.a.discreta por el método de la transformada
inversa.
Se desean simular una muestra de 30 varones y mujeres, solamente se debe copiar esa
fórmula en 30 casillas.
Lo interesante del comando ALEATORIO() es que cada vez que se modifica algún dato
de alguna celda en la planilla, se vuelve a generar otra vez toda la simulación, cambiando
los valores de la matriz.
Glosario :Los modelos de colas o espera en fila son aquellos donde se producen llegadas
( personas autos, objetos etc. ) de manera aleatoria a un lugar de atención o ventanilla donde
deben ser atendidos también con un tiempo aleatorio . El estudio del comportamiento general
del sistema se denomina modelo de colas.
La ventaja de aplicar la teoría de simulación para los fenómenos de espera o colas reside en
que se pueden abordar problemas muy sofisticados donde los modelos clásicos matemáticos
no pueden resolver o lo hacen de manera muy complicada.
En cualquier cola sencilla, ciertos entes (clientes) llegan a punto de servicio, se ponen en
una fila, son atendidos en cierto orden cuando el servicio está disponible, y después salen del
sistema. Los fenómenos básicos que se necesitan para diseñar un modelo del fenómeno de
espera, son:
a)Forma como los clientes llegan al punto de servicio.
b)Forma en como se realiza el servicio.
c)Modo de elegir los clientes de la fila de espera para el servicio.
¿Qué características serán objeto de especial interés cuando se simule un fenómeno de
espera o de colas?
• Longitud de la cola en los diversos tiempos
• El tiempo que el cliente se pasa esperando en el sistema. Es decir, el tiempo que pasa
haciendo la cola más el tiempo que está recibiendo el servicio.
• El tiempo que el sistema está inactivo.
Glosario:
En general, se admite que existe un problema de inventario cuando es necesario tener
almacenado un conjunto de recursos útiles, en algún período de tiempo, con el propósito de
satisfacer unas determinadas necesidades de consumo aleatorio.
Los factores a controlar en un modelo de INVENTARIO O STOCKS son:
• Momento en que se decide realizar un nuevo pedido (punto de pedido).
• Cantidad a pedir en cada reabastecimiento.
El modelo de simulación de un modelo de inventario consistirá en, fijado un " punto de pedido"
y una "cantidad pedida en cada reabastecimiento"; simular el modelo para un tiempo
determinado de tiempo y calcular el costo total.
Puede ser que se requiera la optimización de la operatoria que se logra haciendo variar el
punto de pedido y la cantidad pedida hasta minimizar costos. Mediante la simulación se
obtendrá información acerca del punto de pedido y de la cantidad pedida que hacen mínimo el
costo total.
PRÁCTICA 1
13.Una empresa vende un producto a un precio que es variable según una N(10,1,2)
mientras que las cantidades vendidas son Gamma(k=1; Lambda=1/100). Encontrar el
valor de venta. Y analizar los resultados.
PROBLEMA RESUELTO
Aplicación de generación de variables aleatorias en un problema:
El responsable de un negocio de electrodomésticos está cuestionando el procedimiento que sigue
actualmente a la hora de hacer pedidos a fábrica de una determinada marca de televisores. Este
procedimiento consiste en lo siguiente: se hace un pedido siempre que en la tienda haya 4 o menos
televisores y el número de televisores de los que se piden cada vez es de 10 – y, donde 10 es el
número de televisores con que empieza el estudio e y el número de televisores que hay en la tienda en
el momento de hacer el pedido. Es decir, cada vez que se hace un pedido se solicitan tantos
televisores como sean necesarios hasta completar la cifra 10 . Se supone que si acude un cliente a
comprar un televisor y en ese momento no hay stock en la tienda, se ha perdido la venta. Además se
ha estimado que la demanda diaria de televisores es una v.a. binomial de parámetros n=4, p= 0.45, el
tiempo que transcurre, en días, entre la fecha en que se hace el pedido y la fecha en que la fábrica lo
recibe es también una v.a. binomial de parámetros n=4 y p=0.45. .
4
5 0.1235 1 1 3
6 0.8889 3 3 0
7 0.987 4 0 4 0
8 0.768 3 0 3 0
9 6 6(10-4)
9 0.2456 1 1 5
10 0.876 3 3 2
10 Pedido 3
y ≤ 4 =>
3
11 0.865 3 2 1 0
12 0.654 2 0 2 0
13 0.875 3 0 3 0
14 8 8 (10-2)
14 0.2875 1 1 7
15 0.089 0 0 7
16 0.876 3 3 4
16 0
17 10 6(10-4)
Obsérvese que la simulación de tres semanas da como resultado que se han dejado de vender 13
televisores.
Naturalmente, el dueño del establecimiento puede preguntarse si este valor se mantendrá de esta
forma estable, por este motivo, es lógico que se repita varias veces este procedimiento.
PRÁCTICA 2
MODELOS DE COLAS
1. Supóngase que se quiere simular un periodo de 15 unidades de tiempo y que tanto las
llegadas como los servicios son constantes e iguales a 2,5 y 3 unidades de tiempo
respectivamente. Se supone que el sistema está inicialmente vacío y que la primera
llegada se produce en el instante t=0.Realizar las simulaciones necesarias hasta las 15
unidades. Calcular el tiempo medio de espera en el sistema.
2. Con los datos del ejercicio 1- realizar la misma simulación, suponiendo que en lugar de
haber un punto de servicio hay dos puntos de servicio (por ej, en un supermercado en
lugar de haber una cajera hay dos). Calcular el tiempo que permanece ocioso cada
servicio. Realice los análisis y conclusiones pertinentes.
4. Considere un fenómeno de espera en el que los tiempos entre llegadas de los clientes
al sistema, están distribuidos normalmente con media=1,5 minutos y var.=1/9 min 2 y
que el; tiempo de servicio es constante e igual a 1,5 minutos por cliente. Realizar 10
simulaciones y calcular: a) el tiempo medio de espera de los clientes en el sistema.
b)el número medio de clientes que esperan en el sistema. c) el % de tiempo de ocio.
6. Realice este ejercicio simulando con la computadora en Excel . Intente primero con
100 ejecuciones y luego con 500. Conteste las mismas preguntas del ej3.
El servicio es N(2;05). El 60% de los clientes compra y paga en una caja en medio minuto,
en caso contrario demora el tiempo de servicio. Simular para seis llegadas, en caso
contrario . Hallar: a) la longitud de la cola b) el tiempo medio de espera en el sistema c)
hallar el costo total. Costo del empleado 3$/hora. Costo cliente perdido 5$.
8.En una oficina la secretaria atiende los e-mail que le llegan a su computadora con una
distribución uniforme (2,4). La respuesta a los correos sigue una distribución exponencial
con media 5 ( es decir α=1/5). Calcule el tiempo de espera de cada recibidor en caso de
no tener inconveniente en las comunicaciones.
9.Un empleado atiende los clientes que llegan a una estación de servicio. El tiempo de
servicio está distribuido exponencialmente con una media de 6 minutos. Cuando hay más
de un automóvil en espera de servicio, otro mecánico llega a ayudar, siendo su tiempo de
servicio uniforme (2,4) minutos.
La tasa de llegada de los clientes al sistema es una v.a. Poisson con parámetro λ=1 .
Realice 6 simulaciones y calcule :a)La longitud promedio de la cola. b)El costo del
empleado adicional si éste cobra $2 por hora. c)El % de tiempo de ocio del sistema. d)El
tiempo promedio de espera en el sistema.
10.(Cola simple) Los vehículos llegan a un peaje de la ruta a razón de 3,5 por minuto. Si el
tiempo de pago es una variable exponencial con lambda=6, describir el proceso ( tiempo
de espera, longitud de la cola, tiempo inactivo de ventanilla etc.) y sacar el máximo de
conclusiones. Analizar la situación para valores cambiantes de la tasa de llegadas ( hasta
6 por minuto.
11.(Cola múltiple) Comparar el funcionamiento de una oficina de atención al publico con 2
o con 3 ventanillas si la llegada de personas se produce según una exponencial de media
4 minutos y el tiempo de tramite es N( 11; 3 ) minutos. Realice un análisis con las
conclusiones correspondientes. Realizar las conclusiones pertinentes.
12.(Cola secuencial) La revisión de camiones de transporte debe pasar 1° por la
verificación de carrocería 2° por la de motor y transmisión. Los tiempos respectivos de
atención son U(40,80)min. y N(65,15)min.. Si los vehículos tienen horario asignado de 1
cada hora. Analizar el proceso.
13-(Cola con impaciencia) Los tiempos entre llegadas de personas a una cola de colectivo
es una variable exponencial de media 30 segundos si la cola es de 4 personas la P de que
la quinta se quede es del 50%. En cambio si hay 5 la P de que la sexta se queda es del
PRÁCTICA 3
MODELOS DE INVENTARIOS
Los factores estacionales para cada uno de los meses del año son como se muestra a
continuación:
Demanda (D) 1 2 3 4 5 6 7 8 9
Prob(D=i) 0,135 0,271 0,271 0,18 0,09 0,036 0,012 0,004 0,001
4-La demanda diaria de un cierto artículo está regida por una distibución binomial con
parámetros n=6 y p=1/2. El tiempo de entrega en días es una v.a. Poisson con λ=3. El
costo de mantener una unidad en inventario es de 1$ por día, el costo faltante es de 10$
por unidad, y el costo de ordenar es de 50$ por orden. Se ordena c/8 días hasta 30
artículos. El inventario inicial es 10. Cuál es el costo total.
5-La demanda diaria de un cierto artículo está regida por una distribución binomial con
párametros n=6 ;p=1/2. El tiempo de distrib. en días, es decir el tiempo en que tarda en
llegar el pedido, es una v.a.uniforme (2,4),tiempo en que tarde en llegar el pedido. El stock
inicial es de 10 unidades y el p. de pedido se realiza cuando el st<=4 y en este caso, el
stock solicitado debe completar las 10 unidades. El costo de mantenimiento es de 100$
por día, el costo faltante es de 1000$ por unidad y el costo de ordenar es de 500$ por
orden. Calcular el costo de operación durante una semana.
6-En un comercio el costo de compra es de 200$ por compra, el de almacenamiento es de
0,50 $ / unidad,dia. Si la venta diaria es U(200,300). La reposición es inmediata.
Determinar el tamaño optimo de la partida a comprar y la frecuencia media de compra.
6- (stock con estacionalidad) Si la venta responde a N (1000*sent+500; 100) donde t es el
numero de mes. (sent, es la función:CONTAR.SI).Determinar cual debe ser la producción
mensual constante para cumplir con la entrega. Cual debe ser el stock mínimo para que el
riesgo de no poder entregar sea menor al 1%.
7-(stock perecedero) Una pescadería debe diseñar el tamaño de sus heladeras. La
operatoria consiste en comprar cuando llega el camión y llenar las heladeras. La venta
diaria es N(200;30) . Los días entre llegadas de camión proveedor se produce según una
variable geométrica de media 2 . Los pescados duran como máximo 3 días ( al cuarto
deben tirarse). El costo de almacenamiento es de 0,3$/día Kg . Se compra el pescado a 2
$ y se vende a 3. Analizar el proceso una vez definido el diseño.
Otros modelos
1) Una playa de estacionamiento tiene 30 autos de capacidad ( abre a las 8 y cierra a las 20
hs.), Y llegan autos a razón de 8 por hora, El tiempo de permanencia es v.a. Exponencial
a)
1-En días buenos los clientes llegan a una cola con un servidor según una distribución exponencial con
promedio 4 clientes por minuto, mientras que en otros días llegan según una distribución exponencial
con promedio 12 por minuto. Los tiempos de servicio, de todos los días, se distribuye uniformemente
en el intervalo (0,3) minutos. Suponga que de manera independiente, un día bueno tiene igual
probabilidad que un día no bueno.
Realice 5 simulaciones corrida para cada tipo de día y calcule y analice para cada caso:
a)La longitud de la cola.
b)El tiempo promedio de espera hasta ser atendido.
c)El % de ocio del servidor, si es que existe.
d)Cuál es el tiempo real de servicio en cada caso, es decir debe descontar el tiempo de ocio.
e) Indique qué otras conclusiones podría obtener con este modelo.
f)Compare los resultados obtenidos y realice las interpretaciones correspondientes.
2-La demanda diaria de cierto artículo sigue una distribución uniforme (20,40) El stock inicial es de 100
unidades. Se realiza un pedido cuando el stock final es inferior a 30. Se piden 60 unidades. Los
tiempos de entrega siguen la siguiente distribución de probabilidad:
Tiempo de entrega 1 2 3 4
(en días)
Probabilidad 0.2 0.3 0.25 0.25
Nota: En cada ejercicio indique en la forma más precisa posible como se ha realizado el cálculo
y de dónde obtiene los valores.
b)
1-Dos clubes A y B, con cuatro corredores cada uno, compiten una carrera con postas: Del club A, los
tres primeros corredores tienen una velocidad que es una v.a. N(6;1)m/seg ; en cambio el último
corredor tiene una velocidad que es una v.a. Uniforme con valor mínimo 6 y máximo 7, m/seg.
El club B tiene los cuatro corredores parejos cuya velocidad es una v.a. N(6,3; 0,6) m/seg.
La distancia de cada posta es de 100m.
Simule cinco carreras. ¿Cuál es la probabilidad de que gane el equipo A?
2-Los colectivos pasan exactamente cada cinco minutos. Suponiendo que la cola comienza vacía, las
personas llegan a la parada con intervalos de tiempo exponencial de media un minuto. La cantidad de
asientos vacíos responde a la siguiente función de probabilidad:
Ejercicios de auto-evaluación.
1-Una fábrica de tomate en lata coloca 5 tomate en cada una. Cada tomate pesa según una distribución N(50,5) y
completa el contenido con salsa hasta llegar a los 300 gr.
Complete la tabla. Indique claramente como debe ser completada y llene 5 filas para calcular:
a) El peso neto en tomate (sin salsa) b) El peso neto (con salsa) c) La cantidad de salsa
d) La proporción de latas con menos de 240 gr en tomate sólo.
e) Calcule la ganancia media por envase si se vende a 1,2 $, el tomate se paga 0,70 $/ Kg la salsa cuesta 1,5$/Kg y el
envase 0,20 $ el costo de la mano de obra es de 0,30 $ por lata.
2-Una empresa comercial recibe diariamente 300 docenas de facturas La demanda diaria una distribución uniforme
(250-310) docenas en días normales y Poisson de media 100 en días de lluvia. Si se considera que hay un 30 % de
lluvia. La mercadería no vendida se entrega a un productor agropecuario a razón de 0,2$/doc. La información con
respecto a los costos relevantes es la siguiente:
Precio de venta 2,30 $/docena. Costo 1,4 $/unidad. Gastos fijos diarios 150 $. Costo por demanda insatisfecha 1,5
$/docena.
Complete la tabla. Indique claramente como debe ser completada y realice 6 simulaciones y calcule el costo total de la
operación.
Nota: En cada ejercicio indique en la forma más precisa posible como se ha realizado el
cálculo y de dónde obtiene los valores. En cada ejercicio simule un valor de cada variable
indicada en el mismo. Analice los resultados obtenidos en cada ejercicio.
Ejercicios de autoevaluación:
1-Supongamos que el número de horas que le lleva a una persona aprender a operar cierta máquina es una
variable aleatoria uniforme en el intervalo (6,10) horas. Supongamos que para esta operación se necesitan
dos personas. Simule el tiempo de aprender a operar la máquina; realizando 8 simulaciones. El costo para
aprender una persona es de $25 la hora.
a)Calcule el costo total de las dos personas para la empresa.
b)¿Cuál es la probabilidad de que la primera le gane a la segunda y en qué tiempo?
2.-Una peluquería cuenta con dos empleadas, una para el lavado y otra para el corte y peinado. Los clientes
arriban los días sábado, de acuerdo a un proceso exponencial con promedio 15 por minuto. El tiempo que
dura la el lavado sigue una distribución uniforme entre 0 y 10 minutos. El corte y peinado sigue una
Nota: En cada ejercicio indique en la forma más precisa posible como se ha realizado el cálculo y de
dónde obtiene los valores. En cada ejercicio simule un valor de cada variable indicada en el mismo.
Analice los resultados obtenidos en cada ejercicio.
Herr---An.Datos---Gen
BINOMIAL cant. de éxitos r ( n ) pr (1-p)n-r n p r = 0,1,2,.......,n Σ (Si (Aleatorio()<p,1,0) Nos Aleat.---Binom. n.p n.p (1-p)
r
r =Σ(Si (Aleatorio()<p,1,0)
PASCAL
cant. de ensayos hasta q' la suma = r
hasta r-esimo n ( n-1 )pr(1-p)n-r r p n = r,r+1,......., no esta r/p r (1 - p)
r-1 2
éxito p
GEOMETRICA cant.de ensayos n (1-p)n-1p p ---- n = 1,2,3,......., igual a la anterior hasta que no esta 1/p 1-p
hasta 1erexito la suma sea igual a 1 p2
BINOMIAL
NEGATIVA cant.de fracasos
hasta r-esimo s (s+r-1) pr(1-p)s r p s=n-r igual a Pascal - r NEGBINOMDIST (r / p) - r r (1 - p)
r-1
éxito p2
DISTR.BETA.INV.
BETA (*) exitos p (n+1)! pr(1-p)n-r r n 0 <=p<=1 ___________________
α = r+1 ; β = n-r-1
r+1
r!(n-r)! Herr---Anal.Datos--- n+2
POISSON cant. de exitos k (λ t)ke- λ t λ t k=0,1,2,........, __________________ Gen Aleat.-----Poisson λt λt
k! DISTR.GAMMA INV.
k k-1 - t
GAMMA (*) lapso hasta t λ t eλ k λ t >=0 Σ -1 / λ ln(Aleatorio()) α = k ; β= 1 / λ k/λ k / λ2
k-esimo éxito (k-1)! DISTR.GAMMA INV.
-
EXPONENCIAL tiempo entre t λe λt λ --- t >=0 t = - µ ln (Aleatorio()) α=1 ; β = 1/ λ 1/λ 1 / λ2
éxito y éxito F(t) = 1-e-λt
NORMAL (*) X 1 e-½ ( (x- m)/s)^2 µ σ <= x <= {[ Σ (Aleatorio())] - 6} x σx + µx DISTR.NORMAL.INV. µ σ
2
πσ µ,σ
UNIFORME u F(u)= u - a min max a <= u <= b Aleatorio() x (b-a) + a Aleatorio() x (b-a) + a a+b (b - a)2
b-a a b 2 12