viernes, 5 de octubre de 2012

Método de la transformada inversa

Generación de variables aleatorias

La variabilidad de eventos y actividades se presentan a través de funciones de densidad para fenómenos continuos, y mediante distribuciones de probabilidad para fenómenos de tipo discreto. La simulación de estos eventos o actividades se realiza con la ayuda de la generación de variables aleatorias.
 
 
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 función acumulada F(x) y la generación de números pseudoaleatorios uniformes R = {  ri } ~U (0,1).


El método consiste en:
  1. Definir la función de Densidad f(x) que representa la variable a modelar.
  2. Calcular la función acumulada F(x)
  3.  Igualar a la función acumulada con R, F(x)=R ,y luego despejar la variable aleatoria x para obtener la función acumulada inversa: x =F-1(R).
  4. Generar las variables aleatorias xi, sustituyendo los valores de los números pseudoaleatorios R = { ri } ~U (0,1) en la función acumulada inversa encontrada.
         (observe la portada del libro de texto, Simulación de Raúl Coss Bu)


El método de la transformada inversa también puede emplearse para simular variables aleatorias de tipo discreto, como en las distribuciones de Poisson, de Bernoulli, binomial, geométrica, discreta general, etc. La generación se lleva a cabo a través de la probabilidad acumulada P(x) y la generación de números pseudoaleatorios ri ~U (0,1).

Metodología para generar variables aleatorias continuas.
Metodología para generar variables aleatorias discretas.


Distribución Uniforme

A partir de la función de la densidad de las variables aleatorias uniformes entre a y b.

  • Se obtiene la función acumulada


  • Igualando la función acumulada F(x) con el número pseudoaleatorio ri ~U (0,1), y despejando x se obtiene:


Xi=a + (b - a) F(x)i
Xi=a + (b - a) r


Ejemplo 1:

Los datos del tiempo de servicio en la caja de un banco se comportan de forma exponencial con media de 3 minutos/cliente. Una lista de números pseudoleatorios ri ~U (0,1) y la ecuación generadora exponencial xi = -3In (1 - ri) nos permite simular el comportamiento de la variable aleatoria
Distribución de Bernoulli

A partir de la distribución de probabilidad de las variables aleatorias de Bernoulli con media
p(x) = px (1 – p)1 – x para x=0,1
Se calculan las probabilidades para x=0 y x=1, para obtener


Acumulando los valores de p(x) se obtiene:

Generando números pseudoaleatorios ri ~U (0,1) se aplica la regla:


La tabla siguiente muestra la demanda diaria de cepillos dentales en un supermercado.
Simular el comportamiento de la demanda mediante el método de la transformada inversa.



A partir de la información histórica se calculan las probabilidades puntuales y las acumuladas para x=0,1,2,3



La regla para generar esta variable aleatoria estaría dada por:

Con la lista de números pseudoaleatorios ri ~U (0,1) y la regla anterior es posibles simular la demanda diaria de cepillos dentales, tal como se muestra