Simulacion de Sistemas
Simulacion de Sistemas
Simulacion de Sistemas
CEAIS HUARAZ
ESCUELA PROFESIONAL DE INGENIERA
INFORMTICA Y DE SISTEMAS
Sexta Edicin
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
CAPITULO I
MODELO DE SISTEMAS
1.1. MODELO
Un modelo es una representacin de un objeto, idea, o sistema en una forma
diferente a la entidad misma. En nuestro caso el modelo es un conjunto de
relaciones matemticas o lgicas derivadas de supuestos sobre el
comportamiento del sistema.
Es una Abstraccin de algn sistema real, que tiene la posibilidad de emplearse
para propsitos de prediccin y control.
Es una descripcin lgica de cmo un sistema, proceso o componente se
comporta. La simulacin incluye el diseo de un modelo de un sistema, proceso
o componente, y la realizacin de experimentos sobre el mismo. El propsito de
los experimentos es determinar cmo el sistema real se desempea, y
pronosticar el efecto de los cambios sobre el mismo en el tiempo.
Es una construccin intelectual y descriptiva de una entidad en la cual un
observador tiene inters.
Un modelo es un sistema desarrollado para entender la realidad y en supuestos
simples son usados para consecuencia para modificarla.
No es posible capturar modificar la realidad, en cierta direccin, si es que no
dispone de un modelo que la interprete.
Un modelo de un objeto puede ser una rplica exacta de ste (aunque en un
material diferente y a escala diferente), o puede ser una abstraccin de las
propiedades dominantes del objeto.
Pg. 2
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
Pg. 3
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
C1
C2
C3
Variables Exgenas:
TAi: Tiempo de Arribo de los clientes i.
TSij: Tiempo de Servicio de la caja j, para atender al cliente i.
Variables de Estado:
TEij: Tiempo de Espera del cliente i, para pasar a la caja j.
TOij: Tiempo de Ocio de la caja j, en espera del cliente i.
Variables Endgenas:
Ti: Tiempo Total que el isimo cliente est en el sistema.
Parmetros:
E(TAi): Media del tiempo de Arribo.
VAR(TAi): Varianza del tiempo de arribo.
E(TSj):Media del tiempo de Servicio.
VAR(TSj):Varianza del tiempo de servicio.
Caractersticas de Operacin:
f(TAi): Funcin de densidad de probabilidad para el tiempo de arribo.
F(TSij): Funcin de densidad de probabilidad para el tiempo de servicio.
Identidades:
Para los 3 primeros clientes:
TEij = 0
TO11 = TA1, TO22 = TA2, TO33 = TA3
T1 = TS11, T2 = TS22, T3 = TS33.
Si la diferencia D = TSij TAi, es:
: Existe espera en el sistema.
: Existe ocio en el sistema.
: No existe, ni espera ni ocio.
Pg. 4
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
Modelo de Tiempo Discreto: Cuando el estado del sistema est definido solo
para particulares instantes de tiempo. Una funcin, variable o sistema es discreto,
en contraposicin a continuo, si es divisible un nmero finito de veces. As, el
conjunto de los nmeros naturales es un conjunto discreto, as como la energa
de los estados cunticos.
Pg. 5
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
Modelo Dinmico
Pg. 6
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
Ejemplo
Modelos de sistema de computacin generalmente son de tiempo continuo, estado
discreto, probabilsticos, dinmicos y no-lineales. Algunos son abiertos y otros
cerrados. Tambin pueden ser estables o inestables.
1.5. CARACTERISTICAS DE UN SISTEMA
Estructura: Est dada por la secuencia de formas en el espacio. Ejemplo: auto,
hombre, pas, etc.
Funcionamiento: Est definido por las propiedades del sistema. Para el ejemplo
del auto, sera la velocidad, tipo de combustin, simple o doble traccin, etc.
Comportamiento: Consiste en analizar cul ser la reaccin del sistema
cuando recibe un estmulo externo ya sea de carcter natural o no natural.
Ejemplo: falta de gasolina, cambio de combustible, cambio de llantas, aceitar,
cambio de bujas, etc.
Estabilidad de sistemas: Habilidad para permanecer cercano a un punto de
equilibrio o para retornar a l despus de un disturbio.
Resistencia: Habilidad para resistir desplazamientos desde su estado inicial,
siguiendo un disturbio (cambio brusco de uno o ms variables)
Resiliencia: Habilidad para retornar al estado inicial despus de una disturbio
(recuperacin)
1.6. COMPONENTES DEL SISTEMA:
Entidad: Es un objeto de inters del Sistema: Unidades de: Entrada, Central de
Procesamiento y de Salida.
Atributo: Es una propiedad de la Entidad: Velocidad y almacenamiento.
Actividad: Es todo proceso que provoque cambios en el sistema: Clculos
Matemticos Lgicos y resultados.
Estado del Sistema: Es una descripcin de todas las Entidades, Atributos y
Actividades de acuerdo con su existencia en algn punto del tiempo.
1.7. MEDIO AMBIENTE:
Son los alrededores del Sistema que lo afectan. Existen en Simulacin dos
trminos relacionados con el medio ambiente que son de inters:
Endgeno: Describe las actividades que ocurren dentro del sistema.
Exgeno: Describe las actividades que ocurren en el medio ambiente que
afectan al sistema.
Pg. 7
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
Pg. 8
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
CAPITULO II
SIMULACION DE SISTEMAS
2.1. DEFINICIONES
Segn Robert Shannon, la simulacin es el disear y desarrollar un modelo
computarizado de un sistema o proceso y conducir experimentalmente con este
modelo con el propsito de entender el comportamiento del sistema del mundo real o
evaluar varias estrategias con los cuales puedan operar el sistema.
Para Thomas Taylor, Estos experimentos requieren de operaciones lgicas y
matemticas necesarias para descubrir el comportamiento y la estructura de sistemas
complejos del mundo real a travs de lago perodo de tiempo.
Para Shubik, la simulacin de un sistema o de un organismo es la operacin de un
modelo lo cual se va a llamar simulador el cual es una representacin del sistema.
Este modelo o simulador estar sujeto a diversas manipulaciones, las cuales seran
imposibles de realizar, demasiado costosas o imprcticas. La operacin de un modelo
puede estudiarse y con ello conocer las propiedades concernientes al comportamiento
del sistema o subsistema real - costoso.
Para el autor, la simulacin es la aplicacin de la ciencia y la tcnica en un modelo
abstrado en parte o en su totalidad del sistema real a simular con la finalidad de
analizar sus resultados y predecir cmo se comportar el sistema real simulado.
2.2. IMPORTANCIA DE LA SIMULACION
La simulacin es importante porque permite:
Evaluar cambios en modelos de un sistema existente, frecuentemente es el
mejor camino para reducir el riesgo de las principales decisiones.
Obtener un conocimiento preciso de la naturaleza del proceso.
Identificar problemas especficos o reas problemticas de un sistema.
Desarrollar planes o polticas especficas de un proceso.
Evaluar nuevos conceptos o sistemas antes de su implementacin.
Predecir el comportamiento de un nuevo sistema, sin necesidad de construirlo
fsicamente.
Predecir cambios de un sistema existente sin necesidad de afectar su
operacin.
Evaluar cualquier sistema ante un conjunto de condiciones experimentales
2.3. CUANDO SIMULAR
Sistema actual no existe o es dificultoso observarlo
El sistema actual es muy complejo para analizarlo.
El sistema actual no puede ser interrumpido.
Es costoso construir el sistema actual.
2.4. OBJETIVOS DE LA SIMULACION
Visualizacin: Ver lo que est pasando en el sistema
Clculos (Analizar/Optimizar) : Cuantificar lo que est pasando en el sistema
Pg. 9
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
MODEO DE
SIMULACION
Salidas (Respuestas)
EXPERIMENTACION
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
Pg. 11
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
Pg. 12
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
CAPITULO III
NUMEROS ALEATORIOS
3.1. DEFINICIN
Un nmero es aleatorio cuando no se conoce su generacin y tiene caracterstica
uniformes, independientes y corresponden a una cierta distribucin de probabilidad.
Si se desea desarrollar un Modelo eficiente de Simulacin, se deben de tomar en
cuenta 6 factores importantes
1. Uniformemente distribuidos.
2. Estadsticamente independientes.
3. Reproducibles, que se puedan copiar de un dispositivo a otro (CD, USB, etc.).
4. Con perodos amplios (que no se repitan frecuentemente).
5. Con grandes velocidades de generacin.
6. Con un mnimo de capacidad de almacenamiento.
3.2. PROPIEDADES DE NMEROS ALEATORIOS
Una secuencia de nmeros aleatorios R1, R2, ..., Rx debe tener dos importantes
propiedades estadsticas: uniformidad e independencia. Cada nmero aleatorio Ri es
una muestra independiente tomada de una distribucin continua uniforme entre cero y
uno. Esto es, la funcin de densidad de probabilidad es:
0 < X< 1
En otro caso
3.4. GENERACIN DE NMEROS PSEUDOS ALEATORIOS
La palabra pseudos se refiere a que los nmeros generados no son completamente
aleatorios puesto que se conoce el mtodo de generacin.
3.4.1. MTODO DE CUADRADOS CENTRALES
Fue uno de los primeros mtodos para generar nmeros aleatorios, creado por J.
Von Neumann en 1946. Consiste en el siguiente procedimiento:
1. Elegir un nmero X0 aleatorio de 2n dgitos, llamado semilla.
2. Elevarlo al cuadrado.
3. El nuevo nmero aleatorio se elige de la parte media (central), con 2n dgitos.
4. Repetir indefinidamente la operacin
Ejemplo:
X0 = 4122
x20 = 16 | 9908 | 84
X1 = 9908
x21 = 98 | 1684 | 64
X2 = 1684
x22 = 2 | 8358 | 56
X3 = 8358
x23 = 69 | 8561 | 64
Entre las desventajas de este mtodo est que tienen perodos breves (pocos
nmeros aleatorios), no satisfacen completamente las pruebas estadsticas de
aleatoriedad.
Pg. 13
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
BOTON EJECUTAR
Dim i, n As Integer
Dim x, k, w, y, z, v, t, nlimp As Double
Randomize()
k = Val(TextBox1.Text) REM nmero de simulaciones
x = Val(TextBox2.Text) REM primer nmero tomado al azar
de "N cifras
n = Val(TextBox3.Text) REM nmero de dgitos
nlimp = 0
i = 0
TextBox4.Text = TextBox4.Text & "X [" & i & "] = " & x &
Chr(13) & Chr(10)
For i = 1 To k
y = x * x
z = Int(y / 100)
v = Int(z / 10 ^ (2 * n))
t = v * (10 ^ (2 * n))
w = z - t
TextBox5.Text = TextBox5.Text & "Y [" & i & "] = " &
y & Chr(13) & Chr(10)
If i < nlimp * 10 + 10 Then
Pg. 14
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
Pg. 15
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
+1 = + ( )
Donde a, b y m son nmeros enteros convenientemente elegidos y x0 recibe el
nombre de semilla. Dada la semilla x0, la formula anterior permite construir una
sucesin de nmeros enteros, variando cada uno de ellos entre 0 y m-1,
mediante.
1 = 0 + ( )
2 = 1 + = 2 0 + ( + 1)( )
3 =
(3
(3 + 1)
( )
0 + + + 1) = 0 +
( 1)
2
Pg. 16
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
( + 1)
( ) ( 1)
= 0 +
( 1)
En donde los x son enteros entre 0 y m-1, y las constantes a y b son nonegativas. La seleccin de a, b, y m afectan el periodo y la autocorrelacin en la
secuencia. Entre los resultados de los estudios realizados con estos generadores
tenemos:
1. El modulo m debe ser grande. Dado que los x estn entre 0 y m-1, el periodo
nunca puede ser mayor que m.
2. Para que el computo de mod m sea eficiente, m debe ser una potencia de 2,
es decir, 2k. En este caso mod m puede ser obtenido truncando el resultado y
tomando en k bits a la derecha.
Pg. 17
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
EJEMPLO 02
Elaborar un programa en Visual Basic.Net o Visual Basic 6.0 para generar nmeros
pseudos aleatorios por el mtodo congruencial.
BOTON EJECUTAR
Private Sub Button1_Click(ByVal sender As System.Object, ByVal e As
System.EventArgs) Handles Button1.Click
Dim I, S, R, A, B, M As Integer
Dim X1 As Double
S = Val(TextBox1.Text) REM ingresa numero semilla
A = Val(TextBox2.Text) REM ingresa valor de a numero primo
B = Val(TextBox3.Text) REM ingresa valor de b numero primo <>b
M = Val(TextBox4.Text) REM igreso del modulo
R = Val(TextBox5.Text) REM ingrese veces a generar el n aleatorio
For I = 1 To R
X1 = (A ^ I * S + B * (A ^ I - 1) / (A - 1)) Mod M
TextBox6.Text = TextBox6.Text & "NA[" & I & "] = " & X1 & Chr(13)
& Chr(10)
Next
End Sub
BOTON NUEVO
Private Sub Button2_Click(ByVal sender As System.Object, ByVal e As
System.EventArgs) Handles Button2.Click
TextBox1.Text = ""
TextBox2.Text = ""
TextBox3.Text = ""
TextBox4.Text = ""
Pg. 18
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
TextBox5.Text = ""
TextBox6.Text = ""
TextBox1.Focus()
End Sub
BOTON SALIR
Private Sub Button3_Click(ByVal sender As System.Object, ByVal e As
System.EventArgs) Handles Button3.Click
End
End Sub
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
Hipotecario
Techo Propio
Vehculos
Libre disponibilidad
MONTO MINIMO
100,000.00
50,000.00
30,000.00
10,000.00
MONTO MAXIMO
10,000.00
10,000.00
5,000.00
5,000.00
De pide calcular las sumatorias por cada tipo de prstamo, promedio, varianza,
desviacin estndar, covarianza y hacer sus respectivas interpretaciones.
Pg. 20
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
BOTON SIMULAR
Dim np, i, j, m As Integer
Dim ph(0 To 100), ptp(0 To 100), pv(0 To 100), pld(0 To 100) As
Double
Dim sp(0 To 100) As Double
Dim t As Double
Randomize()
Dim sph, sptp, spv, splp, spld As Double
Randomize()
np = Val(TextBox1.Text)
t = Val(TextBox2.Text)
For i = 0 To np - 1
ph(i) = (Rnd() * (100000 - 10000) + 10000) * (1 + t / 100) ^
((Rnd() * (10 - 0.5) + 0.5))
sph = sph + ph(i)
ptp(i) = (Rnd() * (50000 - 10000) + 10000) * (1 + t / 100) ^
((Rnd() * (10 - 0.5) + 0.5))
sptp = sptp + ptp(i)
pv(i) = (Rnd() * (30000 - 5000) + 5000) * (1 + t / 100) ^ ((Rnd()
* (10 - 0.5) + 0.5))
spv = spv + pv(i)
pld(i) = (Rnd() * (10000 - 5000) + 5000) * (1 + t / 100) ^
((Rnd() * (10 - 0.5) + 0.5))
splp = spld + pld(i)
Next
REM pone sumas a un vector
sp(1) = sph : sp(2) = sptp : sp(3) = spv : sp(4) = splp
For i = 0 To np - 1
tabla.Rows.Add(np)
Next
For i = 0 To np - 1
tabla.Rows(i).Cells(0).Value = i + 1
tabla.Rows(i).Cells(1).Value = Format(ph(i), "#####.00")
tabla.Rows(i).Cells(2).Value = Format(ptp(i), "#####.00")
tabla.Rows(i).Cells(3).Value = Format(pv(i), "#####.00")
tabla.Rows(i).Cells(4).Value = Format(pld(i), "#####.00")
Next
BOTON NUEVO
Private Sub Command2_Click()
Dim k As Integer
k = tabla.Rows.Count
For i = 0 To k
tabla.Rows.Clear()
Next
TextBox1.Text = ""
TextBox2.Text = ""
TextBox1.Focus()
End Sub
BOTON SALIR
Private Sub Command3_Click()
Me.Close()
End Sub
PROGRAMAR TAMBIEN EVENTO FORM LOAD PARA CENTRAR EL FORMULARIO
Private Sub Form_Load()
Dim k As Integer
'Sita el formulario en el centro de la pantalla
Pg. 21
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
EJEMPLO 01
La demanda de un producto est en funcin del tiempo en aos, dicha demanda tiene
un rango de aleatoriedad en un lmite de Ls y Li y presenta una funcin Y=a+bX.
Elaborar un programa para determinar la ecuacin de ajuste y simular el pronstico de
la demanda para este producto para N aos.
Pg. 22
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
BOTON SIMULAR
Dim n, ns, ls, li, i As Integer
Dim y(0 To 100), x(0 To 100) As Double
Dim v, des, r As Double
Dim xy(0 To 100), x2(0 To 100), y2(0 To 100), sumas(0 To 100) As
Double
Randomize()
n = Val(Text1.Text)
ls = Val(Text2.Text)
li = Val(Text3.Text)
ns = Val(Text4.Text)
Dim sx, sy, sxy, sx2, sy2 As Double
sx = 0 : sy = 0 : sxy = 0 : sx2 = 0 : sy2 = 0
REM Dim x, y, x2, y2, xy As Double
For i = 0 To n - 1
x(i) = i + 1 'genera el nmero de aos
Pg. 23
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
BOTON NUEVO
Dim k As Integer
k = DGV1.Rows.Count
For i = 0 To k
Pg. 24
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
DGV1.Rows.Clear()
Next
Text1.Text = ""
Text2.Text = ""
Text3.Text = ""
Text4.Text = ""
Text5.Text = ""
Text6.Text = ""
Text7.Text = ""
Text8.Text = ""
Text9.Text = ""
Text10.Text = ""
Text11.Text = ""
Text12.Text = ""
Text13.Text = ""
Text14.Text = ""
Text15.Text = ""
Text1.Focus()
BOTON SALIR
Private Sub Command3_Click()
Me.Close()
End Sub
Pg. 25
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
EJEMPLO 4
Resolver el mismo problema del ejemplo 2 utilizando el modelo de regresin exponencial
BOTON PRONOSTICAR
Dim n, I, NS As Integer
Dim y(0 To 100), x(0 To 100) As Double
Dim v, des, r As Double
Dim logy(0 To 100), xlogy(0 To 100), x2(0 To 100), logy2(0 To 100),
sumas(0 To 100), prono(0 To 100) As Double
Randomize()
Dim m, sx, sy, sxy, sx2, sy2 As Double
n = Val(Text1.Text)
m = Val(Text2.Text)
sx = 0 : sy = 0 : sxy = 0 : sx2 = 0 : sy2 = 0
Dim Log(0 To 100) As Double
Dim slogy, sxlogy, slogy2 As Double
For I = 1 To n
x(I) = I 'ingrese nmero de datos
y(I) = Val(InputBox("INGRESE DATOS HISTORICOS"))
sy = sy + y(I)
logy(I) = Log(y(I))
xlogy(I) = x(I) * logy(I)
x2(I) = x(I) * x(I)
logy2(I) = logy(I) * logy(I)
sxlogy = sxlogy + xlogy(I)
sx = sx + x(I)
slogy = slogy + logy(I)
sx2 = sx2 + x2(I)
Pg. 26
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
BOTON NUEVO
Dim k As Integer
k = DGV1.Rows.Count
For i = 0 To k
DGV1.Rows.Clear()
Next
Dim N As Integer
k = DGV2.Rows.Count
For i = 0 To N
DGV2.Rows.Clear()
Next
Text1.Text = ""
Text2.Text = ""
Text3.Text = ""
Text4.Text = ""
Text5.Text = ""
Pg. 27
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
Text6.Text = ""
Text7.Text = ""
Text8.Text = ""
Text9.Text = ""
Text10.Text = ""
Text11.Text = ""
Text1.Focus()
BOTON SALIR
Me.Close()
EJEMPLO 5
Una empresa tiene una demanda aleatoria que oscila entre 200 y 100 unidades por da para
un determinado producto. El costo de hacer un pedido es de k soles, el costo de inventario es
de c soles y el costo de almacenamiento es de h soles por da. Tambin se sabe que el tiempo
de entrega para este producto es de L das. Si pide hacer un programa para que simule NS
simulaciones y de las seleccione el mejor dato, es decir el mejor registro simulado. El
programa debe reportar: Lote econmico, costo, ciclo, tiempo efectivo de entrega y punto de
reorden.
Pg. 28
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
BOTON SIMULAR
Dim i, ns, mvy, post As Integer
Dim k, h, l As Double
'dando valor a las variables de las cajas de texto
ns = Val(TextBox1.Text)
k = Val(TextBox2.Text)
h = Val(TextBox3.Text)
l = Val(TextBox4.Text)
Dim d(0 To 100), y(0 To 100), c(0 To 100), t(0 To 100), n(0 To 100),
le(0 To 100) As Double
For i = 0 To ns
d(i) = Rnd() * (200 - 100) + 100 'genera demanda de prodcuto
inventario
y(i) = (2 * k * d(i) / h) ^ 1 / 2 'genera lote economico
c(i) = k / (y(i) / d(i)) + h * (y(i) / 2) 'costo
t(i) = y(i) / d(i) 'coloca valores
n(i) = Int(l / t(i))
le(i) = l - n(i) * t(i) 'tiempo de entrega
Next
mvy = y(i)
For i = 0 To ns - 1
tabla.Rows.Add(ns)
Next
'utilizando un for para aser la simulasiones
For i = 0 To ns - 1
tabla.Rows(i).Cells(0).Value = i
tabla.Rows(i).Cells(1).Value = Format(d(i),
tabla.Rows(i).Cells(2).Value = Format(y(i),
tabla.Rows(i).Cells(3).Value = Format(c(i),
tabla.Rows(i).Cells(4).Value = Format(t(i),
Pg. 29
"###.00")
"###.00")
"###.00")
"###.00")
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
BOTON NUEVO
Dim k As Integer
k = tabla.Rows.Count
For i = 0 To k
tabla.Rows.Clear()
Next
TextBox1.Focus()
TextBox1.Text = ""
TextBox2.Text = ""
TextBox3.Text = ""
TextBox4.Text = ""
TextBox5.Text = ""
TextBox6.Text = ""
TextBox7.Text = ""
TextBox8.Text = ""
TextBox9.Text = ""
TextBox10.Text = ""
BOTON SALIR
Me.Close()
Pg. 30
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
Pg. 31
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
debe seleccionar solo una y no deben ser iguales en cada grupo). Recuerde que
debe modelar el sistema teniendo en cuenta una aproximacin real del sistema, es
decir tener en cuenta el mximo de variables que intervienen en dicho sistema.
11. Elaborar un programa para simular el comportamiento de las exportaciones de una
empresa de venta de flores.
12. Analice el cdigo fuente de un programa de simulacin en cualquier lenguaje de
programacin: de inventario, de produccin, de juegos, de colas, cualquier otro,
escrito en cualquier lenguaje, que no sea muy extenso y presntelo en su informe.
(solo presentar uno por grupo).
NOTA: La calidad de la presentacin, el uso de base de datos, grficos y animacin
permitir al alumno o al grupo obtener una mejor calificacin.
SUGERENCIA: Los programas pueden ser presentados en un solo archivo, para ello
debern usar men, en donde cada opcin de men pueda ejecutar a cada programa.
ADVERTENCIA: Si el alumno presenta cdigo fuente bloqueado para cualquier
ejercicio del laboratorio se considerar como que no ha presentado su laboratorio.
Pg. 32
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
TEORIA DE COLAS
3.6. MODELOS DE COLAS
3.6.1. DEFINICIN
Una Cola es una lnea de espera y la teora de colas es una coleccin de
modelos matemticos que describen sistemas de lneas de espera particulares o
de sistemas de colas. Los modelos sirven para encontrar un buen compromiso
entre costos del sistema y los tiempos promedio de la lnea de espera para un
sistema dado.
Problemas tpicos de Teora de Colas son:
Sistema
Llegadas
Cola
Aeropuerto
Dpto. de bomberos
Compaa telefnica
Panadera
Fbrica
Pasajeros
Alarmas de incendio
Nmeros marcados
Clientes
Piezas para ensamblar
Sala de espera
Incendios
Llamadas
Clientes con nmeros
Inventario en proceso
Mecanismo de
Servicio
Avin
Dpto. De Bomberos.
Conmutador
Vendedor
Estacin de trabajo
Pg. 33
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
Pg. 34
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
a. CARACTERISTICAS DE ARRIBO:
La fuente de ingreso que genera los arribos o clientes para el servicio tiene tres
caractersticas principales:
a. Tamao de la poblacin que arriba
b. Patrn de llegada a la cola
c. Comportamiento de las llegadas.
Tamao de la Poblacin:
El tamao de la poblacin puede ser: infinito (ilimitado) o limitado (finito). La
mayora de los modelos asume arribo infinito. Poblacin de arribo limitada o finita:
cuando se tienen muy pocos servidores y el servicio es restringido. Ejemplo: los
pacientes en un consultorio mdico.
Patrn de arribo al sistema:
Los clientes arriban a ser atendidos de una manera programada (un
paciente cada 15 minutos) o de una manera aleatoria.
Se consideran que los arribos son aleatorios cuando stos son
independientes de otros y su ocurrencia no puede ser precedida
exactamente.
Frecuentemente en problemas de colas, el nmero de arribos por unidad de
tiempo pueden ser estimados por medio de la Distribucin de Poisson que
es una distribucin discreta de probabilidad.
DISTRIBUCION DE POISSON:
x=
nmero de arribos por unidad de tiempo
l =
ratio promedio de arribo
e = 2.71828
Comportamiento de los arribos:
La mayora de los modelos de colas asume que los clientes son pacientes o sea
que esperan en la cola hasta ser servidos y no se pasan entre colas.
Desafortunadamente, la vida es complicada y la gente reniega. Aquellos que se
impacientan por la espera, se retiran de la cola sin completar su transaccin.
Pg. 35
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
Entrada
SERVIDOR
Salida
Servidor F2
Salida
Entrada
Servidor F1
CANAL1
CANAL2
Entradas
CANAL3
Pg. 36
Salidas
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
Fase 1
Canal 1
Fase 2
Canal 1
Salidas
Entradas
Fase 1
Canal 2
Fase 2
Canal 2
D = constante o determinstica
Ek = Erlang de orden k
P H = Tipo fase
H = Hiperexponencial
G = Arbitrario o general
GI = General independiente
K = Tamao de la poblacin.
Nota(1): A causa de las suposiciones de distribucin exponencial en los
procesos de arribo, estos modelos son llamados MARKOVIANOS
Pg. 37
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
Pg. 38
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
SIMBOLO/FORMULA
N
Ls = / ( )
=/
Ws = 1 / ( )
Lq = 2 / ( ( ))
Lq = Ls
Wq = / ( ( ))
Wq = Ws
Pn = (1 / )*( / )n
Pn = (1 ) * n
Po = 1 / = (1 )
Pn>k = ( / )k+1
SIMBOLO/FORMULA
Po
tasadepromedio
arribo
Probabilidad
que existan cerode
clientes
en el sistema
Po Probabilid
existan
CERO personas o unidades en el s
ad de que
n !
M ! 1 M
0
nP
para M
o
M
n M 1 1 n
1
M
M !
n 0 n ! de personas
Ls nmero promedio
unidades
en el sistema :
oM
M
Ls = nmero promedio
de personas o unidades en el sistema
LS
promedio
Ls nmero
Pode
personas o unidades en el sistema :
1.!.M M
LS
Po
2
M 1.!.M
Pg. 39
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
Ejemplo 01
Elaborar un programa que simule un sistema de colas M/M/1/1 dado un rango de
tiempo entre llegadas de los clientes y un rango de tiempo de servicio. Estos rangos
son aleatorios. Se pide calcular: promedio de clientes en el sistema, factor de
utilizacin del sistema, promedio de cliente en cola, servicio, tempo promedio en el
sistema, tiempo promedio en cola, tiempo en el sistema, tiempo en cola y
probabilidades de que haya N clientes en el sistema, ningn cliente en el sistema y
ms de K clientes en el sistema.
BOTON SIMULAR
Dim i, tsll, till, tsa, tia, ns, n, k, pos As Integer
Dim tll(0 To 50), ta(0 To 50), pcs(0 To 50), fu(0 To 50) As Double
Dim tpcs(0 To 50), npuc(0 To 50), tpcc(0 To 50), mvy, pn, po, pk As
Double
Randomize()
tsll = Val(TextBox1.Text)
till = Val(TextBox2.Text)
tsa = Val(TextBox3.Text)
tia = Val(TextBox4.Text)
ns = Val(TextBox5.Text)
n = Val(TextBox11.Text)
k = Val(TextBox12.Text)
For i = 0 To ns - 1
Pg. 40
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
Pg. 41
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
End Sub
BOTON NUEVO
Dim k As Integer
k = DGV.Rows.Count
For i = 0 To k
DGV.Rows.Clear()
Next
TextBox1.Text = ""
TextBox2.Text = ""
TextBox3.Text = ""
TextBox4.Text = ""
TextBox5.Text = ""
TextBox6.Text = ""
TextBox7.Text = ""
TextBox8.Text = ""
TextBox9.Text = ""
TextBox10.Text = ""
TextBox11.Text = ""
TextBox12.Text = ""
TextBox13.Text = ""
TextBox14.Text = ""
TextBox15.Text = ""
'TextBox1.SetFocus()
'CheckBox1.cheked = False
'CheckBox2.cheked = False
'CheckBox3.cheked = False
BOTON SALIR
End
Pg. 42
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
CAPITULO IV
Son aquellas medidas que nos ayudan a saber dnde estn los datos pero
sin indicar como se distribuyen.
Media aritmtica
La media aritmtica o simplemente media, que denotaremos por , es el
nmero obtenido al dividir la suma de todos los valores de la variable entre el
nmero total de observaciones, y se define por la siguiente expresin:
Pg. 43
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
=1
=1
=1
Media geomtrica
Sea una distribucin de frecuencias que se simboliza por G. se define como
la raz N-sima del producto de los N valores de la distribucin.
= 11 22 33 . .
Media armnica
Es el reciproco de la media aritmtica. No es aconsejable en distribuciones
de variables con valores pequeos. Se suele utilizar para promediar variables
tales como productividades, velocidades, tiempos, rendimientos, cambios,
etc.
=
1
=1
Varianza
Es una medida que cuantifica el grado de dispersin o de variacin de los
valores de una variable con respecto a su media aritmtica. Si los valores
tienden a concentrarse alrededor de su media, la varianza ser pequea, de
Pg. 44
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
Desviacin estndar
Es una medida de centralizacin o dispersin para variables de razn (ratio o
cociente) y de intervalo. Se define como la raz cuadrada de la varianza.
Junto con este valor, la desviacin tpica es una medida (cuadrtica) que
informa de la media de distancias que tienen los datos respecto de su media
aritmtica, expresada en las mismas unidades que la variable. Se simboliza
con s minscula.
Indica cunto tienden a alejarse los valores concretos del promedio en una
distribucin. De hecho, especficamente, la desviacin estndar es "el
promedio de la distancia de cada punto respecto del promedio".
= 2
4.3. CONCEPTOS BASICOS DE LA TEORA MUESTRAL
Pg. 45
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
Pg. 46
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
Pg. 47
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
CAPITULO V
Pg. 48
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
DERECHA
(50 , 100)
(100 , 50)
(0, 50)
(50, 0)
Pg. 49
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
CAPITULO VI
LA SIMULACIN EN LA PLANIFICACIN, PROGRAMACIN Y
CONTROL DE PROYECTOS
6.1. FUNDAMENTOS DE LA REPRESENTACIN GRFICA DE UN PROYECTO
6.1.1. PROYECTO
Es un conjunto de tareas u operaciones elementales bien diferenciables que se
ejecutan segn un orden determinado.
Los fundamentos de los sistemas PERT y CPM son las representaciones grficas del
proyecto mediante diagramas de flechas, o tambin lo podemos llamar red de flechas.
La red se crea segn el orden de realizacin de las tareas u operaciones, paso a paso,
hasta el final del proyecto. Originalmente estas tareas u operaciones se llaman
actividades. Un trabajo encargado a una persona responsable, bien lo realice
personalmente o bien lo hagan operarios a sus rdenes es lo que podemos definir
como actividades.
Una actividad puede comprender una sola tarea o bien una serie de ellas. Todo
depende de la designacin del responsable de los trabajos que se realizan bajo sus
rdenes segn la conveniencia de la realizacin del proyecto. Por tanto habr tantas
actividades como responsables.
Grficamente una actividad est compuesta de dos partes: la primera que es la
ejecucin del trabajo y est representada con una flecha con orientacin de izquierda
a derecha y la segunda se llama suceso que generalmente se dibuja con dos crculos
o dos rectngulos ponindolos en los dos extremos de la flecha:
i
El suceso que esta al final de la flecha se llama suceso inicial y el suceso que
conecta al comienzo de la flecha se le denomina suceso final. El suceso es un
instante de la actividad que sirve como el punto de control, describiendo el momento
de comienzo o terminacin de una actividad. La actividad es un smbolo de trabajo en
proceso. Por tanto, todas las actividades requieren tiempo y recursos.
La longitud de la flecha no representa la cantidad de tiempo como en los grficos de
GANTT. Por ejemplo, en la siguiente figura, la actividad A no es ms corta de duracin
que la B, aunque la longitud de las flechas lo sea:
B
A
1
Pg. 50
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
Tampoco es preciso que la flecha sea una lnea recta, sino que pueden dibujarse en
curva:
Esto depende de la facilidad que haya para representar las actividades en una red de
flechas que refleje el orden y secuencias de las relaciones del proyecto. Una actividad
debe estar terminada para que la subsiguiente pueda comenzar. Como todas las
actividades tienen sus sucesos iniciales y finales, el suceso final de la actividad
precedente es el mismo suceso inicial de la subsiguiente:
A
1
Sin embargo, hay una excepcin en los sucesos iniciales y finales. El primer suceso
inicial del proyecto no tiene una actividad que la preceda y el ltimo suceso final
tampoco tiene una actividad que la subsiga.
6.1.2. ENUMERACIN DE LOS SUCESOS
La enumeracin de los sucesos es otro es otro sistema para la identificacin de la
actividad. Pero para facilitar el clculo en el computador es conveniente asignar
nmeros naturales a los sucesos iniciales y finales. El siguiente grafico ser numerado
como sigue:
3
B
D
1
Pg. 51
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
Actividad B = (2,3).
Actividad C = (2,4).
Actividad D = (3,4).
Actividad E = (4,5).
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
Paso inicial. Determinar n = n para indicar que la primera y la ltima ocurrencia del
ultimo acontecimiento del proyecto son las mismas.
Paso general j. Dado que los nodos p,q,..., y v estn unidos directamente al nodo j
por las actividades de salida (j,p),(j,q),..., y ya se han calculado, entonces la ltima
fecha de ocurrencia del evento j se calcula como
1. i = i
2. j = j
3. j - i = j - i = Dij
Las tres condiciones plantean que la primera y la ltima fecha de ocurrencia de los
eventos finales i y j son iguales y que la duracin Dij ocupa exactamente su lapso. Por
tanto, una actividad que satisface las tres condiciones es no crtica. Las actividades
crticas de una red deben construir una ruta ininterrumpida que abarca toda la red de
principio a fin.
15
45
18
E - 12
B - 10
3
10
C - 10
K - 15
D - 15
A-6
0
F - 20
I - 30
30
G-7
30
10
J - 13
H - 26
21
47
10
36
60
60
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
2. Investigue sobre el mtodo PERT y elabore un programa para simular una red
de proyectos. Alcance toda la teora del dicho mtodo y desarrolle un ejemplo.
3. Analice y muestre el cdigo fuente de un software que simula el mtodo PERT
y mtodo CPM.
4. Nombre a cada actividad y determine la ruta crtica de la siguiente red de
proyecto, construya el diagrama de Gantt respectivo.
G-6
B-5
D-10
A-4
1
I-3
E-7
M-11
O-3
K-12
F-8
J-4
H-5
C-6
N-8
L-5
G - 15
B-8
E-4
C - 11
A-6
I-5
4
H-4
K-7
F-9
D - 12
L-3
J-8
Pg. 55
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
LIMPIAR
Y PINTAR
BASE
A
1
FICTICIA
C
2
SACAR
Y
DESARMAR
EL MOTOR
E
5
REBOBINAR LA
ARMADURA
FICTICIA
4
Reemplazar
los anillos
Pg. 56
6
ENSAMBLAR
E INSTALAR
EL MOTOR
EN LA BASE
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
[E]
[C]
ORD. Y REC.
2
[A]
DISEAR
PROD.
[ ]
FABR.
PROD.
[ ]
[G]
ENVASAR
PROD.
[I]
[ ]
8
[B]
[ ]
[H]
DISEAR
ENVASE
[D]
[ ]
[J]
ENTREGA
1
0
PRUEBA
PROD.
[F]
ORD. Y REC.
PRUEBA
MERCADO
FABR.
ENV.
[3]
[CODIGO]
DESCRIPCION
CLAVE:
DURACION
Dij
2.-
3.-
Pg. 57
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
to + 4tm + tp
te =
6
Veamos que ocurre con el tiempo con el caso de una empresa X en el cual se
proporcionan tres estimaciones de los tiempos que se requieren para terminar cada
una de las actividades del proyecto.
TABLA
Cdigo de
la actividad
A
B
C
D
E
F
G
H
I
J
Tiempo
optimista(to)
3.0
1.0
1.5
1.2
2.0
1.8
3.0
2.0
0.5
0.8
Tiempo mas
probable(tm)
5.5
1.5
3.0
3.2
3.5
2.8
6.5
4.2
0.8
2.1
Tiempo
pesimista(tp)
11.0
5.0
4.5
4.0
8.0
5.0
7.0
5.2
2.3
2.8
Pg. 58
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
Varianza de los
Tiempos de actividad
= t
(tp to) 2
36
Varianza
t2
1.78
0.44
0.56
0.22
1.00
0.28
0.44
0.28
0.09
0.11
Pg. 59
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
Z=
X-
Z=
26 22
=2
Pg. 60
SIMULACION DE SISTEMAS
UNIVERSIDAD SAN PEDRO
Pg. 61