Trabajo Metodo de Aceptacion y Rechazo
Trabajo Metodo de Aceptacion y Rechazo
Trabajo Metodo de Aceptacion y Rechazo
Presentado a:
LUIS CARLOS FORERO
PÁ
G.
Introducción...........................................................................................................................................4
1. Variable aleatoria...............................................................................................................................5
1.1. Variable aleatoria discreta............................................................................................................5
1.2. Variable aleatoria continua...........................................................................................................5
2. Generacion de variables aleatorias....................................................................................................5
2.1. Métodos más empleados para la generación de variables aleatorias........................................6
2.2. Tipos de generadores de variables aleatorias...............................................................................7
3. Método de aceptación y rechazo......................................................................................................10
3.1. Demostración...................................................................................................................................10
3.2. Algoritmo de aceptación y rechazo...............................................................................................10
4. Simulación del método de aceptación y rechazo............................................................................14
conclusiones.............................................................................................................................................17
bibliografía................................................................................................................................................18
INTRODUCCIÓN
Una variable aleatoria es una función que asocia un número a cada elemento del
espacio muestral, en donde el recorrido o rango de la variable es el conjunto de valores
que toma.
Se llama variable aleatoria a toda función que asocia a cada elemento del espacio
muestral E un número real.
Una variable aleatoria discreta es aquella que sólo puede tomar valores enteros.
Una variable aleatoria continua es aquella que puede tomar todos los valores posibles
dentro de un cierto intervalo de la recta real.
Ejemplo: La altura de los alumnos de una clase, las horas de duración de una pila 1 .
1
(VITUTOR, Variables Aleatorias (En línea))
2.1. MÉTODOS MÁS EMPLEADOS PARA LA GENERACIÓN DE VARIABLES
ALEATORIAS
2.1.1. Método de la transformada inversa.
Cuando f(x) es una función acotada y x tiene un rango finito, como a x b, se utiliza
este método para encontrar los valores de las variables aleatorias. El método consiste
en normalizar el rango de f mediante un factor de escala c, luego definir a x como una
función lineal de r, después se generan parejas de números aleatorios r1 , r2 y por
último si el número encontrado se elige al azar dentro del rango (a,b) y r b, se utiliza
este método para encontrar los valores de las variables aleatorias. El método consiste
en normalizar el rango de f mediante un factor de escala c, luego definir a x como una
función lineal de r, después se generan parejas de números aleatorios r1 , r2 y por
último si el número encontrado se elige al azar dentro del rango (a,b) y r cf (x) se
acepta, en caso contrario se rechaza. El problema de este método es la cantidad de
intentos que se realizan antes de encontrar una pareja exitosa.
Ejemplo. Desarrollar cinco iteraciones del generador X n+1 = 3Xn mod 100, con X0=51.
Xn RND
51 0.51
53 0.53
59 0.59
77 0.77
31 0.31
2.2.2.1. Selección de los parámetros del generador congruencial multiplicativo.
A continuación damos los parámetros para dos casos especiales que ya han sido
probados para este generador (Varela, Kelton respectivamente).
a = es la constante multiplicativa.
c = es la constante aditiva.
m = es la magnitud del módulo.
X0 = es la semilla.
Aquí mod representa a la operación aritmética módulo entre los enteros a y b tal que el
resultado de (a mod b) es el residuo entero de la división a entre b.
3.1. DEMOSTRACIÓN
4
(Poli Libros Mexico, (En línea))
Se genera y un valor y de la variable Y (de forma
independiente).
Si u>g(y), se va al paso 1.
Como el número de iteraciones del método hasta aceptar un valor sigue una
distribución geométrica de parámetro (cuya esperanza es C), entonces si queremos
optimizar el método habrá que intentar que C sea próxima a uno y que h sea sencilla de
generar.
Además,
Ejemplo:
Si y sólo si:
Si y sólo si:
Si y sólo si
Ejemplo:
Tomándose:
Ejemplo:
Ahora,
Como , entonces,
Consideremos X una variable aleatoria cuyo soporte es un intervalo finito [a, b] sobre el
que la función de densidad es acotada y no nula. Sea c = máx. x ∈ [a , b] f(x). Los
métodos generales de generación de variables aleatorias procedimiento es el siguiente:
La probabilidad de rechazo será el cociente entre el área de esta zona y el área del
rectángulo. La región de aceptación tiene área 1 por ser f(x) una función de densidad.
0,5 0,20
0,15
0
0 0,2 0,4 0,6 0,8 1 0,10
0,00 [0 - 0.1] (0.1 - 0.2] (0.2 - 0.3] (0.3 - 0.4] (0.4 - 0.5] (0.5 - 0.6] (0.6 - 0.7] (0.7 - 0.8] (0.8 - 0.9] (0.9 - 1.0]
Frecuencias Probabilidades
0,5 0,20
0,15
0
0 0,2 0,4 0,6 0,8 1 0,10
0,00 [0 - 0.1] (0.1 - 0.2] (0.2 - 0.3] (0.3 - 0.4] (0.4 - 0.5] (0.5 - 0.6] (0.6 - 0.7] (0.7 - 0.8] (0.8 - 0.9] (0.9 - 1.0]
Frecuencias Probabilidades
Resultados de la simulación
U Y f(Y)/t(Y) Clases Simulaciones válidas (frecuencia)
relativas (%) teóricas
0,27836 0,64659 0,97695 [0 - 0.1] 9,55% 100,00% 0,00% 0,0013
0,70519 0,54876 0,97362 (0.1 - 0.2] 10,10% 96,04% 0,83% 0,0157
0,89257 0,05578 0,00448 (0.2 - 0.3] 9,95% 76,38% 4,87% 0,0535
0,93204 0,74735 0,77097 (0.3 - 0.4] 10,25% 50,73% 10,47% 0,1087
0,66478 0,18476 0,12128 (0.4 - 0.5] 10,10% 18,32% 17,10% 0,1646
0,17960 0,87281 0,31123 (0.5 - 0.6] 9,65% 4,66% 19,07% 0,2006
0,99542 0,50102 0,90605 (0.6 - 0.7] 11,45% 2,62% 23,11% 0,2000
0,89294 0,60799 0,99933 (0.7 - 0.8] 9,20% 22,83% 14,72% 0,1568
0,30003 0,64265 0,98070 (0.8 - 0.9] 9,65% 58,03% 8,39% 0,0830
0,15524 0,13447 0,05270 (0.9 - 1.0] 10,10% 93,07% 1,45% 0,0158
0,42751 0,49771 0,90005 Total 100,00% 51,75% 100,00% 1,0000
0,87500 0,14364 0,06289
0,70247 0,36457 0,56612
0,72045 0,82932 0,48080
0,88721 0,26554 0,29225
0,33205 0,42225 0,72713
0,48143 0,61018 0,99891
0,33882 0,31709 0,43024
0,06519 0,88285 0,27325
0,63818 0,77975 0,66545
0,74429 0,10516 0,02695
0,50432 0,50155 0,90701
0,61945 0,66463 0,95547
0,62156 0,64887 0,97462
0,98747 0,58387 0,99732
0,21515 0,38726 0,63094
0,93043 0,29755 0,37613
0,25327 0,43088 0,74973
0,87994 0,27218 0,30907
0,93493 0,43744 0,76653
0,38620 0,81643 0,53063
.
.
.
CONCLUSIONES
UNAL. URL:
www.virtual.unal.edu.co/cursos/sedes/manizales/4060015/Lecciones/Capitulo
%20VI/valeatorias.htm
ECCI. URL:
http://www.kramirez.net/ProbaEstad/Material/Presentaciones/GeneracionVariable
sAleatorias.pdf