UNIDAD II
NUMEROS PSEUDOALEATORIOS
En todos los experimentos de Simulación,
existe la necesidad de generar valores de
variables aleatorias que representan a una
cierta distribución de probabilidad.
Este proceso se hace a partir de la generación
de números aleatorios en el intervalo (0, 1),
llamados rectangulares, que representan el
valor de una variable aleatoria distribuida
uniformemente en (0,1) a los cuales se hará
referencia como ri
Algunas de las formas para obtener a los
números aleatorios son:
Provisión externa.
Generación interna a partir de un proceso
físico.
Generación interna por medio de una relación
de recurrencia.
Características que deben de poseer
los números aleatorios:
1. Uniformemente distribuidos.
2. Estadísticamente independientes
3. Reproducibles
4. Periodos largos (sin repetición dentro
de una longitud n determinada de la
sucesión)
5. Generados a través de un método
rápido.
6. Generados a través de un método
que no requiera mucha capacidad de
almacenamiento en la computadora.
Técnicamente, un NA (ri), se define como una
muestra aleatoria independiente tomada de una
Distribución continua Uniforme cuya función de
densidad de probabilidad esta dada por:
1, 0 ≤ r ≤1
f(r) = 0, en cualquier otro caso.
Generar un conjunto ri de números
aleatorios es una tarea relativamente
sencilla.
Sin embargo lo difícil es diseñar un
algoritmo que genere un conjunto de ri con
periodo de vida suficientemente grande, y
que además pase sin problema las pruebas
de uniformidad e independencia, es decir,
se deben evitar los problemas siguientes.:
• Que los números del conjunto ri no
estén uniformemente distribuidos, es
decir que existan demasiados ri en un
subíntralo y en otro muy pocos o
ninguno.
• Que la media del conjunto sea muy
alta o muy baja, es decir que este por
arriba o por debajo de ½.
• Que la varianza del conjunto sea muy
alta o muy baja, que se localice por
arriba o por debajo de 1/12.
Los algoritmos determinísticos para
generar los ri, se clasifican en :
Congruenciales
No congruenciales
Algoritmo de cuadrados medios o central
Fue uno de los primeros métodos para
generar números aleatorios, creado por J. Von
Neumann en 1946. Consiste en el siguiente
procedimiento:
1.- Elegir un NA de “n” dígitos, llamado
semilla.
2.- Elevarlo al cuadrado y añadir ceros
a la izquierda en donde se necesite
para producir “2n” dígitos.
3.- El nuevo número aleatorio se elige
de la parte media (central) de la etapa
anterior, con “n” dígitos.
4.- Repetir los pasos 2 y 3 para obtener
los siguientes NA de la Serie.
Entre las desventajas de este método
está que tienen períodos breves (pocos
números aleatorios), no satisfacen
completamente las pruebas estadísticas
de aleatoriedad.
Ver ejemplos cuadrado_1 y cuadrado_2.
Resolver cuadrado_3.
Basados en una relación fundamental de
Congruencia:
“Dos enteros a y b son congruentes módulo
“m”, si su diferencia es múltiplo entero de
m”. Esto es:
a ≡ b (mód m)
1.- (a – b) es divisible entre m.
2.- a y b dan el mismo residuo al ser
divididos entre |m|.
Para este generador, la relación congruencial
de recurrencia se expresa de la siguiente
forma:
ni+1 = (ani + c) mód m.
Donde:
n : es la semilla o valor inicial (n0 > 0).
0
a : factor multiplicativo (a > 0).
c : constante aditiva (c > 0).
m : el módulo (m>n , m>a y m>c).
0
ni+1: residuo de dividir (ani+c) /m.
ri+1: el número uniforme obtenido de: (ni+1)
/m.
Para que el algoritmo pueda lograr el máximo
periodo de vida, es preciso que los
parámetros cumplan ciertas condiciones
como:
m = 2g, donde g debe ser entero
a = 1 + 4k, donde k debe sr entero.
c relativamente primo a m.
Periodo logrado n=m = 2g
Ejemplos:
Utilizar el generador:
ni+1 = (9ni+13) mód 32, con n0 = 7 para
obtener {ni} y {ri}.
Ver mixto_1
b)Generar 4 números entre 0 y 1 con los
siguientes parámetros:
n0 = 37, a = 19, c= 33 y m = 100.
Ver mixto_2.
Resolver los ejercicio mixto_3 y mixto_4
Produce una sucesión de enteros {ni}, entre 0
≤ni ≤ m y una sucesión de números
uniformes {Ri}, entre 0 ≤ ri ≤ 1.
La ecuación de Congruencia recurrente
empleada es:
ni+1 = (ani) mód m.
donde:
no: es la semilla o valor inicial (n0 > 0).
a : Factor multiplicativo (a > 0).
m : Módulo
ni+1: residuo de dividir (ani) entre el módulo m.
Ri+1: Número uniforme, se obtiene del cociente
de (ni+1) entre m.
Para que el algoritmo pueda lograr el
máximo periodo de vida, es preciso que los
parámetros cumplan ciertas condiciones
como:
m = 2g, donde g debe ser entero
a = 3 + 8k, o a = 5 + 8k , donde k>0.
n0 debe ser un entero impar.
relativamente primo a m.
Periodo logrado n=m/4 = 2g-2
Ejemplos:
a) Generar suficientes números entre 0 y 1 con
los siguientes parámetros:
n0 = 17, k = 2, g = 5,
Hasta encontrar el periodo de vida máximo.
Obtenga :
un generador con los siguientes valores de sus
parámetros:
a=3, m=102 y n0=17. La ecuación de
recurrencia se expresará como:
ni+1 = ( 3ni ) mód 102.
Ver multiplicativo_1.
Resolver multiplicativo_2
