Función de Densidad: Generar Con La Serie U (0,1) : 0,19 0,41 0,83 0,16 0,50
Función de Densidad: Generar Con La Serie U (0,1) : 0,19 0,41 0,83 0,16 0,50
Función de Densidad: Generar Con La Serie U (0,1) : 0,19 0,41 0,83 0,16 0,50
El método de la transformada inversa puede utilizarse para simular variables aleatorias continuas, lo cual se logra
mediante la obtención de la función acumulada F(X) y la generación de números pseudo aleatorios ri con distribución
uniforme entre 0 y 1.
Ejemplo:
función de densidad 1
1 𝑠𝑖 𝑎 ≤ 𝑥 ≤ 𝑏
𝑓(𝑥) = 𝑓(𝑥) = 𝑏 − 𝑎
𝑏−𝑎
0 𝑠𝑖 𝑎 > 𝑥 > 𝑏
1 1 1 𝑥 𝑥−𝑎
𝐹(𝑥) = 𝑑𝑥 = 𝑑𝑥 = [𝑥] 𝐹(𝑥) =
𝑏−𝑎 𝑏−𝑎 𝑏−𝑎 𝑎 𝑏−𝑎
𝑥−𝑎
𝐹(𝑥) = 𝑟 𝑟 = 𝑟 . (𝑏 − 𝑎) = 𝑥 − 𝑎 𝑥 = 𝑎 + 𝑟 . (𝑏 − 𝑎)
𝑏−𝑎
Ejercicios:
⎧ ; 𝑝𝑎𝑟𝑎 0 ≤ 𝑥 < 1
⎪
2. 𝑓(𝑥) = Generar con la serie U(0,1): 0,19 0,41 0,83 0,16 0,50
⎨
⎪ ; 𝑝𝑎𝑟𝑎 1 ≤ 𝑥 < 6
⎩
a) b)
Algunas soluciones:
. . 𝑥 . . 𝝀.𝒙
𝐹(𝑋) = 𝜆. ∫ 𝑒 𝑑𝑥 = −𝑒 = −𝑒 +𝑒 = −𝒆 +𝟏
0
Integral indefinida:
.
𝑑𝑢 𝜆 𝒙
𝜆. 𝑒 𝑑𝑥 = 𝜆. 𝑒 = 𝑒 𝑑𝑢 = −1. 𝑒 = −𝒆
−𝜆 −𝜆
𝑢 = −𝜆. 𝑥
𝑑𝑢 = −𝜆. 𝑑𝑥
𝑑𝑢
𝑑𝑥 =
−𝜆
𝑅𝑁𝐷 = .
−𝑒 +1
𝑅𝑁𝐷 − 1 = .
−𝑒
−𝑅𝑁𝐷 + 1 = .
𝑒
𝑙𝑛(−𝑅𝑁𝐷 + 1) = −𝜆. 𝑥
𝑥 = −1
. 𝑙𝑛(1 − 𝑅𝑁𝐷)
𝜆
Con esta última fórmula pueden obtenerse números pseudo-aleatorios que responden a una distribución
exponencial negativa.
⎧ ; 𝑝𝑎𝑟𝑎 0 ≤ 𝑥 < 1
⎪
2. 𝑓(𝑥) = Generar con la serie U(0,1): 0,19 0,41 0,83 0,16 0,50
⎨
⎪ ; 𝑝𝑎𝑟𝑎 1 ≤ 𝑥 < 6
⎩
𝑥 √𝑥 𝑥 𝑥 𝑥
𝐹(𝑋) = 𝑑𝑥 = 𝑑𝑥 = 𝑥 𝑑𝑥 = = . 𝑥 = 𝑥 =
16 4 0 0 6
𝑅𝑁𝐷 = 𝑥
6
𝑅𝑁𝐷. 6 = 𝑥
Para poder continuar se debe valuar la integral definida para el segmento en el que la F(X)1 es válida (entre cero y
uno).
1
𝐹(𝑋) = ∫ 𝑥 𝑑𝑥 = 𝑥 = −0 = Este valor de 1/6 pasará a ser el corte entre un generador y otro
0
Por lo tanto el área total de la primer fórmula F(X)1 es de 1/6. Este valor se debe agregar a la segunda función
acumulada F(X)2, de otro modo la misma no contemplaría el total de la superficie.
1 1 1 1 1 𝑥 1 1 1 1 1 1
𝐹(𝑋) = + 𝑑𝑥 = + 𝑑𝑥 = + [𝑥] = + (𝑥 − 1) = + . 𝑥 − = . 𝑥
6 6 6 6 6 1 6 6 6 6 6 6
𝑅𝑁𝐷 = 1
.𝑥
6
𝑅𝑁𝐷. 6 = 𝑥
Se procede a valuar la serie numérica 0,19 ; 0,41 ; 0,83 ; 0,16 ; 0,50, cada número con su correspondiente
generador. El primer generador es válido para valores de la serie hasta 1/6 (valor de corte).
Resultados:
( )
𝐹(𝑋) = − Valor de corte: 𝑋 = 18 − 32. (para 0 ≤ RND < ½)
0,29 => 15,40 / 0,68 => 20,40 / 0,55 => 19,26 / 0,34 => 15,73 / 0,86 => 21,39
a) b)
Resultados:
a)
b)
( )
𝑓(𝑥) = .𝑥 − 𝐹(𝑋) = 𝑋 = 190 + √𝑅𝑁𝐷. 800
( )
𝑓(𝑥) = .𝑥 + 𝐹(𝑋) = 1 − 𝑋 = 230 − √−𝑅𝑁𝐷. 800 + 800