UNIDAD 2. Variables Aleatorias. Simulacion
UNIDAD 2. Variables Aleatorias. Simulacion
UNIDAD 2. Variables Aleatorias. Simulacion
MATERIA:
SIMULACIÓN
TITULAR:
ING. RAMON CORRES SHIBAYAMA
CARRERA:
ING. INDUSTRIAL
UNIDAD 2.
SIMULACIÓN DE VARIABLES ALEATORIAS.
SEMESTRE:
6°
GRUPO:
“Q”
ALUMNO:
DIANA MONSERRAT JIMENEZ LUNA
PERIODO ESCOLAR:
ENERO/JUNIO 2021
En simulación, para poder realizar las corridas de los sistemas o generar los
modelos se necesita de variables aleatorias, que a su vez requieren de
números pseudoaleatorios. Un número pseudoaleatorio es un número U (0,1)
producido por un algoritmo matemático. Es deseable que los números
pseudoaleatorios uniformes posean las siguientes características:
1. Uniformemente distribuidos.
2. Estadísticamente independientes.
3. Reproducibles.
4. Periodo largo.
5. Generados mediante un método rápido.
6. Generados mediante un método que no requiera mucha capacidad de
almacenamiento de la computadora.
Y otras características son que 𝜇 = 1/2, 𝛿 2 = 1/12 y debe existir 𝐼𝑛𝑑𝑒𝑝𝑒𝑛𝑑𝑒𝑛𝑐𝑖𝑎
𝑒𝑛𝑡𝑟𝑒 𝑣𝑎𝑟𝑖𝑎𝑏𝑙𝑒𝑠. A los conjuntos de números pseudoaleatorios también se les
conoce como ri y se pueden generar mediante distintos métodos, además
de que para ser representativos se les aplican distintas pruebas estadísticas
como:
Prueba de medias
Prueba de varianza
Pruebas de uniformidad
• Chi-cuadrada
• Kolmogorov-Smirnov
Pruebas de independencia
• Corridas arriba y abajo
• Corridas arriba y debajo de la media
• Prueba póker
• Prueba de series
• Prueba de huecos
Algoritmo lineal:
Este algoritmo genera una secuencia de números enteros por medio de la
ecuación 𝑋𝑖 + 1 = (𝑎𝑋𝑖 + 𝑐) 𝑚𝑜𝑑(𝑚) 𝑖 = 0,1,2,3, …, 𝑛
Donde:
• Xo>0, semilla
• a=1+4k, constante multiplicativa, k debe ser entero
• c primo a m, constante aditiva
• m=2^g es el módulo, donde g es un entero
A partir de los enteros generados por la ecuación anterior, para obtener los
números pseudoaleatorios (0,1) se usa la ecuación
𝑟𝑖 = 𝑋𝑖 𝑚 − 1 𝑖 = 0,1,2,3, …, n.
Algoritmo congruencial multiplicativo:
𝑋𝑖 + 1 = (𝑎𝑋𝑖)𝑚𝑜𝑑(𝑚) 𝑖 = 0,1,2,3, …, 𝑛
Donde:
• Xo debe ser un número impar
• a=3+8k, constante multiplicativa, k= 0, 1, 2, 3, …
• m=2^g es el módulo, donde g es un entero.
𝑟𝑖 = 𝑋𝑖 𝑚 − 1 𝑖 = 0,1,2,3, …, 𝑛
𝑥 = 𝐹 −1 (𝑅)
Método de convolución
En algunas distribuciones de probabilidad la variable a simular, Y, puede
generarse mediante la suma de otras variables aleatorias X de manera más
rápida que a través de otros métodos. Entonces el método de convolución se
puede expresar como:
𝑌 = 𝑋1 + 𝑋2 + ⋯ + 𝑋𝑘
Las variables aleatorias de Erlang, normal, binomial y Poisson pueden generarse
a través de este método.
Erlang
La variable aleatoria k. Eerlang con media 1/λ puede producirse a través de la
generación de k variables exponenciales con media 1/kλ:
𝑌 = 𝑋1 + 𝑋2 + ⋯ + 𝑋k
Distribución normal
La variable aleatoria normal con media µ y desviación estándar σ puede
generarse
mediante el teorema de límite central:
𝑌 = 𝑋1 + 𝑋2 + ⋯ + 𝑋𝑘
Al sustituir Xi por números pseudoaleatorios se obtiene:
Distribución binomial
Las variables pueden ser generadas a través de la suma de N variables
aleatorias con distribución de Bernoulli con parámetro p. 𝑌 = 𝐵𝑗 = 𝐵𝐸1 + 𝐵𝐸2 + ⋯
+ 𝐵𝐸N.
Método de composición
El método de composición permite generar variables aleatorias x cuando éstas
provienen de una función de densidad f(x) que puede expresarse como la
combinación convexa de m distribuciones de probabilidad fi(x). La combinación
convexa puede expresarse como:
𝑓(𝑥) = ∑𝑓𝑖(𝑥)𝐼𝑎(𝑥)
Donde
Ia (x)= 1 si x E A
Ia (x)= 0 si x no E A
Las distribuciones más usadas en este método son la triangular, la de Laplace y
la trapezoidal. El procedimiento general de la generación es:
1. Calcular la probabilidad de cada una de las funciones f(x).
2. Asegurarse que cada función f(x) se función de densidad.
3. Obtener mediante el método de la transformada inversa, expresiones para
generar variables aleatorias de cada una de las distribuciones.
4. Generar un número pseudoaleatorio ri.
5. Seleccionar la función generadora correspondiente a f(x).
6. Generar un segundo número pseudoaleatorio ri y sustituirlo en la función
generadora para obtener Y.
• Distribución triangular
A partir de la función triangular
• Sustituyendo