Practicas Simulacion Msse
Practicas Simulacion Msse
Practicas Simulacion Msse
Modelado y Simulacin de Sistemas de Eventos Discretos Xabier Basogain / Miguel ngel Olabe
Prcticas de LABORATORIO
PROGRAMA PRCTICO
SESIN 0.- CONCEPTOS BSICOS DE ESTADSTICA
1.- once!tos y "erminolog#a
$.- %!licaciones &r'cticas de Modelos Estad#sticos
(.- Distribuciones Discretas
).- Distribuciones ontinuas
SESIN .- INTROD!CCIN AL S" DE SIM!LACIN ARENA
1.- Descri!cin %rena
$.- &rimer E*em!lo
(.- Estudio de Mdulos
).- E+E,--O - &roceso de Solicitud de una .i!oteca
SESIN #.- CONSTR!CCIN DE MODELOS CON ARENA STANDARD
1.- Me*oras en la /isuali0acin de la Simulacin del Modelo
$.- Modi1icaciones en el &roceso de Solicitudes de una .i!oteca
(.- Modelos +er'r2uicos3 submodelos
).- E+E,--O - &roceso de Solicitud de &r4stamo !ara %utomviles
SESIN $.- MODELOS ESTADSTICOS EN SIM!LACIN
1.- .erramienta -n!ut %naly0er
$.- Datos a %nali0ar
(.- Datos y /entanas
).- 5eneracin de Datos
6.- %*uste de los Datos a una Distribucin
7.- Modi1icacin de &ar'metros
8.- E*ercicios
SESIONES % & '.- SIM!LACIN de SISTEMAS de COLAS
1.- -ntroduccin
$.- "eor#a de olas y %rena
(.- E*ercicios
SESIN (.- ANLISIS DE DATOS DE ENTRADA
1.- -ntroduccin
$.- E*ercicios
(.- %!4ndice3 Distribuciones
3
Xabier Basogain / Miguel ngel Olabe Modelado y Simulacin de Sistemas de Eventos Discretos
Sesin 9.- once!tos B'sicos de Estad#stica
1.$.- /ariable %leatoria ontinua.-
1
Modelado y Simulacin de Sistemas de Eventos Discretos Xabier Basogain / Miguel ngel Olabe
SESIN 0.- CONCEPTOS BSICOS DE ESTADSTICA
Ob*etivo3 -ntroducir y re!asar los conce!tos b'sicos de estad#stica utili0ados en las sesiones
!r'cticas de la asignatura Modelado y Simulacin --. :os modelos !robabil#sticos utili0ados
en el modelado y simulacin de sistemas de eventos discretos re2uieren el conocimiento de
los t4rminos y conce!tos elementales de la estad#stica b'sica.
;ndice3
1.- once!tos y "erminolog#a
/ariable %leatoria Discreta< /ariable %leatoria ontinua< =uncin de
Distribucin %cumulativa< /alor Es!erado< Moda.
$.- %!licaciones &r'cticas de Modelos Estad#sticos
Sistemas de olas< Sistemas -nventario< Mantenimiento y =iabilidad< Datos
:imitados o -ncom!letos< Otras Distribuciones.
(.- Distribuciones Discretas
Bernoulli< Binomial< 5eom4trica y &oisson.
).- Distribuciones ontinuas
>ni1orme< E?!onencial< 5amma< Erlang< @ormal< Aeibull y "riangular
.- CONCEPTOS ) TERMINOLOGA
1.1.- /ariable %leatoria Discreta.-
El nBmero de !osibles valores de la variable es 1inito o in1inito !ero contable.
&ara cada !osible valor ?i de la variable X se tiene 2ue !C?iD E !CX E ?iD es la
!robabilidad de 2ue la variable X tome el valor ?i.
Se cum!len las siguientes condiciones3
aD !C?iD 9 &ara todo ?i
bD
pC xiD = 1
i =1
Distribucin de &robabilidad o =uncin Masa de &robabilidad C!m1D de X es el
con*unto de !ares C?i<!C?iDD con iE1< $ ...
E*e m !lo :an0amiento del dado trucado
$
Modelado y Simulacin de Sistemas de Eventos Discretos Xabier Basogain / Miguel ngel Olabe
El es!acio de valores de la variable X C,?D es un intervalo o un con*unto
de intervalos.
:a !robabilidad de 2ue el valor de X se encuentre en un intervalo Fa<bG viene
dada !or la e?!resin3
b
PCa x bD =
f C xDdx
a
1C?D se llama 1uncin densidad de !robabilidad C!d1D de la variable X
!d1 satis1ace las siguientes condiciones3
aD 1C?D 9 &ara todo ?
bD
f C xDdx = 1
Rx
cD 1C?D E 9 Si x Rx
E*e m !lo =uncionamiento de una bombilla
1.(.- =uncin de Distribucin %cumulativa.-
:a 1uncin de distribucin acumulativa Ccd1D< denotada !or =C?D< mide la
!robabilidad de 2ue la variable X tenga un valor menor o igual 2ue ?H es decir =C?D
E & CX ?D.
Si X es discreta
Si X es continua
&ro!iedades de
=C?D3
F C xD =
pC
xiD
xi x
x
F C xD =
f Ct D
dt
x f C xD dx
Si X es discreta
(
Modelado y Simulacin de Sistemas de Eventos Discretos Xabier Basogain / Miguel ngel Olabe
E
C
X
D
=
x
i
p
C
x
i
D
todoi
:a media es una medida de la tendencia central de la variable aleatoria
1.6.- /arian0a.-
/CXD
$
se denomina varian0a y se de1ine como3
/CXD E E FCX - ECXDD
$
G e2uivalentemente /CXD E ECX
$
D - FECXDG
$
:a varian0a de X mide la variacin de los valores de ? res!ecto de la media
:a desviacin est'ndar se de1ine como la ra#0 cuadrada de la varian0a de X
1.7.- Moda.-
Se de1ine de la siguiente 1orma3
/ariable Discreta3 :a moda es el valor de la variable 2ue a!arece m's 1recuentemente.
/ariable ontinua3 :a moda es el valor m'?imo de !d1
:a Moda !uede no ser Bnica.
Si el valor de la moda ocurre en dos valores la distribucin es bimodal
#.- APLICACIONES PRCTICAS DE MODELOS ESTADSTICOS
$.1.- Sistemas de olas.-
:a distribucin del tiem!o entre llegadas y la distribucin del nBmero de llegadas !or
!eriodo de tiem!o son im!ortantes a la Kora de simular los sistemas de colas.
El tiem!o de servicio !uede ser constante o !robabil#stico.
Distribuciones utili0adas3
E?!onencia l. Si los tiem!os de servicio son com!letamente aleatorios.
@o r m al . Si los tiem!os de servicio son constantes !ero e?iste una variabilidad 2ue
!roduce 1luctuaciones !ositivas y negativas.
@o r m al "runcada . Si e?isten valores de la variable 2ue deben ser mayores o
menores 2ue un cierto valor y el resto siguen una distribucin normal.
5 a m m a y Aeibul l . Se utili0an !ara modelar "iem!os de Servicio.
$.$.- Sistemas -nventario.-
E?isten tres variables aleatorias3
aD @Bmero de items solicitados en cada !edido o en cada !eriodo de
tiem!o
bD "iem!o entre !edidos
cD "iem!o entre reali0ar un !edido y recibir dicKo !edido
:a variable nBmero de items o tamaLo del !edido suele ser re!resentada !or
las siguientes distribuciones3
5eo m 4trica . Se caracteri0a !or tener la moda centrada en la unidad
Bino m i al @egativa. Se caracteri0a !or tener una cola larga
&oisson. Est' tabulada y se conoce con !ro1undidad. "iene una cola
m's corta 2ue la binomial negativa
$.(.- Mantenimientos y =iabilidad.-
El tiem!o de 1allo !uede ser modelado !or varias 1unciones3
Distribucin E?!onencial. Si solamente ocurren 1allos aleatorios
Distribucin 5 a m m a . Surge del modelado !or redundancia.
Distribucin A eibul l. uando Kay mucKos com!onentes en un sistema
y el 1allo se debe al de1ecto m's serio del con*unto de de1ectos
Distribucin @o r m al . uando los 1allos se deben al desgaste normal
Distribucin :ogonor m a l. Se utili0a !ara describir el tiem!o de 1allo de
algBn ti!o de com!onentes
$.).- Datos :imitados o -ncom!letos.-
@o se dis!one del su1iciente nBmero de datos
:as 1unciones 2ue se utili0an son3
Distribucin uni1or m e . Se utili0a cuando se sabe 2ue el tiem!o entre llagadas o
de servicio es aleatorio !ero no se dis!one de m's in1ormacin
Distribucin "riangular. Se !uede utili0ar cuando se Kacen su!osiciones sobre
el m'?imo< el m#nimo y la moda
Distribucin Beta. &ro!orciona una gran variedad de 1ormas en su distribucin
M
C
od
D
elado
p
$
y Simulacin de Sistemas de Eventos
Discretos
Xabier Basogain / Miguel ngel
Ola
p
be
$.- DISTRIB!CIONES DISCRETAS
(.1.- Distribucin de Bernoulli
Es una distribucin de !robabilidad con dos !untos de !robabilidad discreta
de1inida como3
!C9D E 2
!C1D E ! siendo !M2 E 1 con !<2 N9
Sea un e?!erimento consistente en n ensayos y cada uno de ellos !uede tener
4?ito o 1racaso C1 9D. Sea !or e*em!lo ?*E1 4?ito y ?*E9 1racaso.
Se tiene3
!C?1< ?$ ... ?nD E !1C?1D. !$C?$D ... !nC?nD
p
?* E 1 * E 1< $ ... n
pjC xjD = pC xjD =
1 p =
q
9
?* E 9 * E 1< $ ... n
otros
(.$.- Distribucin Binomial
:a variable X indica el nBmero de 4?itos en n ensayos de Bernoulli y tiene una
distribucin binomial dada !or3
p
x
q
n x
pC xD =
9
? E 9< 1< $ ... n
resto
:a media y la varian0a son3
E C?D E n! /C?D E n!2
(.(.- Distribucin 5eom4trica C,elacionada con la secuencia de ensayos de BernoulliD
:a variable X indica el nBmero de ensayos !ara obtener el !rimer 4?ito. :a
distribucin de esta variable es3
q
x 1
p
pC xD =
9
? E 1< $ ...
otros
El evento OX E ?P ocurre cuando Kay ?-1 1allos seguidos de un 4?ito
ada uno de los 1allos tiene asignada una !robabilidad de 2E1-! y cada uno de
los 4?itos tiene !robabilidad !. %s#3
!C===....=SD E 2
?-1
.!
:a Media y la /arian0a vienen dados !or3
EC xD =
1
M
C
od
D
elado
p
$
y Simulacin de Sistemas de Eventos
Discretos
Xabier Basogain / Miguel ngel
Ola
p
be
V x =
q
5
7
Modelado y Simulacin de Sistemas de Eventos Discretos Xabier Basogain / Miguel ngel Olabe
(.).- Distribucin de &oisson
Se utili0a !ara modelar tiem!os entre eventos aleatorios ocurridos en un
intervalo de tiem!o 1i*o.
:a 1uncin masa de !robabilidad C!m1D est' dada !or3
e
? E 9< 1 ...
otros con N 9
pC xD =
xQ
9
>na !ro!iedad im!ortante de la distribucin de &oisson es EC?D E /C?D E
:a 1uncin de distribucin acumulativa es3
x
F C xD =
i
i
Esta 1uncin est' tabulada
E*e m !lo3 Servicio de =ontaner#a
%.- DISTRIB!CIONES CONTIN!AS
).1.- Distribucin >ni1orme
i = 9
Q
Se utili0a cuando todos los valores en un rango 1inito se !ueden considerar
iguales
:a variable aleatoria X est' uni1ormemente distribuida en el intervalo Ca<bD si
!d1 est' dada !or3
cd1 est' dada !or3
1
f C xD =
b
a
9
a ? Ib
otros
F C x D =
9
x
a
b
a
1
? I a
a ? I b
? b
Se tiene 2ue !C?1 I ? I ?$D E =C?$D-=C?1
x
D
$
E
x1
b a
es !ro!orcional a la
longitud del intervalo !ara todo ?1 y ?$ 2ue satis1aga a ?1 I ?$ b
:a media y la /arian0a son3
EC X D =
a + b
$
Cb aD
$
V C X D =
1$
9
? 9
resto
1
e
x
F C xD =
9
? 9
resto
:a media y la /arian0a son3
E C xD =
1
V C xD =
1
$
).(.- Distribucin 5amma
Se utili0a !ara re!resentar el tiem!o re2uerido !ara 1inali0ar una tarea
>na variable aleatoria X tiene una distribucin 5amma con !ar'metros y
si !d1 es3
f C xD =
C D
C xD
1
e
x
9
? N 9
otros
donde es un !ar'metro de 1orma y es un !ar'metro de escala
y adem's3
C D =
x
R S1
e
S ?
dx
9
!ara entero C D = C 1DQ
:a Media y la /arian0a son3
E C xD =
1
V C xD =
1
$
uando es entero< la distribucin 5amma est' relacionada con la
E?!onencial
&ara E 1 se obtiene una distribucin E?!onencial
).).- Distribucin Erlang
:a e?!resin !d1 de 5amma< !ara E T< con T entero se denomina
distribucin de Erlang de orden T
:a Media y la /arian0a son3
E C x D =
1
V C xD =
1
K
$
Se veri1ica 2ue3
F C xD =
k 1
e
k x
CK
xD
i
i = 9
? N 9
? 9
).6.- Distribucin @ormal
>na variable aleatoria X con media C- I I D y varian0a
$
tiene una distribucin normal si !d1 es3
1 C ?D E
1
e?!
S
1 x
!
$ " #
$
#
$
%
%
&
- I ? I
Se utili0a la notacin @C< D
:a cd1 de la Distribucin @ormal es3
x
F C xD = pC X xD =
e?!
#
1 t
$
!
%
dt
$ " $
1
e
$ "
t
$
$
dt
).7.- Distribucin Aeibull
Se utili0a en modelos de 1iabilidad !ara re!resentar tiem!os de vida de
dis!ositivos. >n sistema 1ormado !or mucKas !artes inde!endientes y el
sistema 1alla cuando Kay un 1allo.
>na variable aleatoria tiene una distribucin Aeibull si !d1 es
x (
1
x (
!
f C xD =
9
e?!
#
%
? (
otros
:os tres !ar'metros de una distribucin Aeibull son ( C- I ( I D 2ue
es el !ar'metro de locali0acinH C N 9D 2ue es el !ar'metro de escala y
CN9D 2ue es el !ar'metro de 1orma. &ara ( E 9 se tiene 2ue la !d1 es3
x
1
x
!
f C x D =
9
e?!
#
%
? 9
otros
9
iQ
$
$)c x *
f C xD =
)c b*)c a *
a ?
b b I ?
c otros
:a cd1 de una distribucin "riangular es3
9
) x a *
$
)b a *)c a
*
F C x D =
1
)c x *
? a
a I ? b
1
)c b *)c a
*
b I ? c
? N c
"%B:% - ,ES>ME@
Distribuci
n
1C?D !C?D =C?D EC?D /C?D
F C xD =
pC x
i
D
x
i
x
x
F C xD =
f
EC X D =
x
i
pC x
i
D
+i
E C X D =
x f C xD
/CXD E
ECX
$
D -
FECXDG
$
P*iss*+
e
x
pCxD =
? =
9<1<$
?Q
9 resto
EC?D E /C?D E
$
f CxD =
$
E,-*+e+c
ia.
e
x
? 9
x
F C xD =
1 e
9
E C xD =
1
1
V C xD =
$
Er.a+/-0
1
x
C xD
e f C xD =
C D
9
k 1
e
k x
CK
xD
i
= C ?D E
iQ
1
i = 9
9
E C x D =
1
V C xD =
1
K
$
"ei12..
1
!
x
e?! #
x
%
f C x D =
#
9
9
F C xD =
!
e?!
x (
1
#
%
#
$
x
=
? = 9< 1< ..
D
E 9
pC xD
Q?
9
F
resto
EC xD =
V C xD =
Distribucin >ni1orme.-
f C xD
=
b -
a
D
a G ? G b
F
9 resto
E C xD =
a + b
$
V C xD =
Cb a D
1$
Distribucin E?!onencial.-
e
-?
f C xD =
? 9
D
F
9
E C xD =
1
resto
H
V C xD =
1
$
61
Modelado y Simulacin de Sistemas de Eventos Discretos Xabier Basogain / Miguel ngel Olabe
Distribucin 5amma.-
f C xD
=
C
xD
CxD
1
e
x
D
x E 9
F
resto
V C xD =
1
$
Distribucin Aeibull.-
x /
1
C
x /
D
D
x v
f C xD =
9
F
resto
!ar'metro de =ormaC E
9D
!ar'metro de Escala C E
9D
C-n!ut %naly0er %rena EscalaD
C-n!ut %naly0er %rena =ormaD
v !ar'metro de locali0acin C - G v G D
E C xD =
V C xD =
Distribucin "riangular.-
$ C x a D
Cb aDCc
aD
a x b
D
f C xD =
$ C c
x D
Cc bDCc
aD
b x c
9 resto
H
E C xD =
a + b + c
(
V C xD =
a + b + c
(
9
H
C D
e
F
ES"-M%DO,ES de M'?i m a /erosi m ilitud
Distribucin 5amma 3 &ar'metros R y j
0 = ln X
n
ln X
i
1
/e Tab%as 1 / 0
= 1 / X
Distribucin Aeibull 3 &ar'metros R y %l1aCkD de M'?ima /erosimilitud
n
n
X
i
nX
l
= C
1
X
l
D
1/
l
-
$
=
i =1
n
i =1
n
n
n
n
X
i
ln X
i
n 1
i n
f C D =
+
i
=1
ln X
i =1
i
i =1
n n
f VC D =
n X
Cln X D
$
i =1
nC
+
i
=1
X
ln X D
$
i
i =1
l l
f C
l
D
C
i =1
X
D
$
j
=
j 1
j 1
f VC
l
D
1
n
$ $
i
i
X
j
1