Función de Densidad: Generar Con La Serie U (0,1) : 0,19 0,41 0,83 0,16 0,50

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 4

Método de la transformada inversa

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.

Los pasos a seguir para efectuar el método son:


1. Definir la función f(x) que represente la variable a modelar.
2. Calcular la función acumulada F(X).
3. Igualar la función acumulada F(X) a un número pseudo aleatorio ri U(0,1), y luego despejar la variable aleatoria
x para obtener la función acumulada inversa F(X)-1.
4. Generar las variables aleatorias x, a través de números pseudo aleatorios ri U(0,1) en la función acumulada
inversa.

Ejemplo:

1. Definición: distribución uniforme en el intervalo A-B.

función de densidad 1
1 𝑠𝑖 𝑎 ≤ 𝑥 ≤ 𝑏
𝑓(𝑥) = 𝑓(𝑥) = 𝑏 − 𝑎
𝑏−𝑎
0 𝑠𝑖 𝑎 > 𝑥 > 𝑏

2. Cálculo de la función acumulada.

1 1 1 𝑥 𝑥−𝑎
𝐹(𝑥) = 𝑑𝑥 = 𝑑𝑥 = [𝑥] 𝐹(𝑥) =
𝑏−𝑎 𝑏−𝑎 𝑏−𝑎 𝑎 𝑏−𝑎

3. Igualar a ri y despejar la variable aleatoria x.

𝑥−𝑎
𝐹(𝑥) = 𝑟 𝑟 = 𝑟 . (𝑏 − 𝑎) = 𝑥 − 𝑎 𝑥 = 𝑎 + 𝑟 . (𝑏 − 𝑎)
𝑏−𝑎

4. Generar variables aleatorias

Ejercicios:

1. 𝑓(𝑥) = 𝜆𝑒 para 𝑥≥0

⎧ ; 𝑝𝑎𝑟𝑎 0 ≤ 𝑥 < 1

2. 𝑓(𝑥) = Generar con la serie U(0,1): 0,19 0,41 0,83 0,16 0,50

⎪ ; 𝑝𝑎𝑟𝑎 1 ≤ 𝑥 < 6

(18 − 𝑥); 𝑝𝑎𝑟𝑎 14 ≤ 𝑥 < 18


3. 𝑓(𝑥) = Generar con la serie U(0,1): 0,29 0,68 0,55 0,34 0,86
(𝑥 − 18); 𝑝𝑎𝑟𝑎 18 ≤ 𝑥 < 22
4. Generar 5 números aleatorios de cada uno de los siguientes gráficos de probabilidad:

a) b)

Algunas soluciones:

1. 𝑓(𝑥) = 𝜆𝑒 para 𝑥≥0

Calculo de la función acumulada:

. . 𝑥 . . 𝝀.𝒙
𝐹(𝑋) = 𝜆. ∫ 𝑒 𝑑𝑥 = −𝑒 = −𝑒 +𝑒 = −𝒆 +𝟏
0

Integral indefinida:

.
𝑑𝑢 𝜆 𝒙
𝜆. 𝑒 𝑑𝑥 = 𝜆. 𝑒 = 𝑒 𝑑𝑢 = −1. 𝑒 = −𝒆
−𝜆 −𝜆
𝑢 = −𝜆. 𝑥
𝑑𝑢 = −𝜆. 𝑑𝑥
𝑑𝑢
𝑑𝑥 =
−𝜆

Igualación a ri (lo denominaremos RND) y despeje de x:

𝑅𝑁𝐷 = .
−𝑒 +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

Se calcula la función acumulada para f(x)1:

𝑥 √𝑥 𝑥 𝑥 𝑥
𝐹(𝑋) = 𝑑𝑥 = 𝑑𝑥 = 𝑥 𝑑𝑥 = = . 𝑥 = 𝑥 =
16 4 0 0 6

Igualación a RND y despeje de x:

𝑅𝑁𝐷 = 𝑥
6

𝑅𝑁𝐷. 6 = 𝑥

𝑋 = (𝑅𝑁𝐷. 6) (para 0 ≤ RND < 1/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

Igualación a RND y despeje de x:

𝑅𝑁𝐷 = 1
.𝑥
6

𝑅𝑁𝐷. 6 = 𝑥

𝑋 = 𝑅𝑁𝐷. 6 (para 1/6 ≤ RND < 1)

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).

0,19 𝑋 = 𝑅𝑁𝐷. 6 = 0,19 . 6 = 1,14


0,41 𝑋 = 𝑅𝑁𝐷. 6 = 0,41 . 6 = 2,46
0,83 𝑋 = 𝑅𝑁𝐷. 6 = 0,83 . 6 = 4,98
0,16 𝑋 = (𝑅𝑁𝐷. 6) = 0,97
0,50 𝑋 = 𝑅𝑁𝐷. 6 = 0,50 . 6 = 3,00
(18 − 𝑥); 𝑝𝑎𝑟𝑎 14 ≤ 𝑥 < 18
3. 𝑓(𝑥) = Generar con la serie U(0,1): 0,29 0,68 0,55 0,34 0,86
(𝑥 − 18); 𝑝𝑎𝑟𝑎 18 ≤ 𝑥 < 22

Resultados:

( )
𝐹(𝑋) = − Valor de corte: 𝑋 = 18 − 32. (para 0 ≤ RND < ½)

𝐹(𝑋) = + (𝑥 − 18) 𝑋 = 18 + 32. 𝑅𝑁𝐷 − (para ½ ≤ RND < 1)

0,29 => 15,40 / 0,68 => 20,40 / 0,55 => 19,26 / 0,34 => 15,73 / 0,86 => 21,39

4. Generar 5 números aleatorios de cada uno de los siguientes gráficos de probabilidad:

a) b)

Resultados:

a)

𝑓(𝑥) = .𝑥 − 𝐹(𝑋) = . (𝑥 − 170) 𝑋 = √𝑅𝑁𝐷. 1600 + 170

b)

( )
𝑓(𝑥) = .𝑥 − 𝐹(𝑋) = 𝑋 = 190 + √𝑅𝑁𝐷. 800

( )
𝑓(𝑥) = .𝑥 + 𝐹(𝑋) = 1 − 𝑋 = 230 − √−𝑅𝑁𝐷. 800 + 800

También podría gustarte