1Eva_IT2010_T1 Urnas de Ehrenfest

Tema 1. (Urnas de Ehrenfest) Se tienen dos urnas, con 5 bolas repartidas dentro de ellas, y en cada etapa se escoge una bola al azar y se cambia de urna.

La cadena Xn representa el número de bolas de una de las urnas tras n etapas.

Escriba la matriz de transición de un paso de esta cadena de Markov.

Referencia: FCNM/ICM01420

1Eva_IIT2009_T5 Manejo de inventarios

Tema 5. Considérese la cadena de Markov Xt que representa el número de unidades que se tienen al final del período t y considérese la política (2,2), es decir si el nivel del inventario al final del período es menor a 2 se ordenan 2 unidades, de otra manera no se ordena.  La demanda tiene distribución Poisson con λ=1.

Suponga que el costo de ordenar es 10 +25 z dólares, si no se ordena el costo es 0 y por cada unidad de demanda insatisfecha se tiene un costo de 50 dólares por unidad.

Encontrar las probabilidades de estado estable y el costo promedio esperado a la larga por unidad de tiempo.

Referencia: FCNM/ICM01420

1Eva_IIT2009_T4 Nacimientos en país poco poblado

Tema 4. Los niños nacen en un país poco poblado, con una frecuencia de un nacimiento cada 12 minutos. El tiempo entre nacimientos sigue una distribución exponencial.

Determinar:

a) La cantidad promedio de nacimientos por año.

b) La probabilidad de que no haya nacimientos en cualquier día.

c) La probabilidad de emitir 50 certificados de nacimientos en 3 horas, cuando se emitieron 40 certificados durante las primeras 2 horas del período de 3 horas

Referencia: FCNM/ICM01420

1Eva_IIT2009_T3 Falla de maquina embotelladora

Tema 3. Suponga que el tiempo entre descomposturas de una máquina es exponencial, con promedio de 6 horas. 

Si la máquina ha trabajado sin fallar durante las últimas tres horas,

a. ¿cuál es la probabilidad de que continúe sin fallar durante la próxima hora?

b. ¿Que se descomponga durante la siguiente 1/2 hora?

1Eva_IIT2009_T2 Bebidas gaseosas

Tema 2. Suponga que toda la industria de refresco produce solo dos bebidas gaseosas:  ROJA y AZUL. 
Cuando una persona ha comprado ROJA hay una probabilidad de 90% de que siga comprándola la vez siguiente.
Si una persona compró AZUL, hay 80% de que repita la vez siguiente.

Se pide:

a) Si una persona actualmente es comprador de AZUL. ¿Cuál es la probabilidad de que compre ROJA pasadas dos compras a partir de hoy?

b) Si en la actualidad una persona es comprador de ROJA. ¿Cuál es la probabilidad de que compre ROJA pasadas tres compras a partir de ahora?

c) Supongamos que el 60% de toda la gente toma hoy ROJA y el 40% AZUL. A tres compras a partir de ahora, ¿Qué fracción de los compradores estará tomando ROJA?. (utilice probabilidad total)

d) Determinar las probabilidades de estado estable.

Referencia: FCNM/ICM01420

1Eva_IIT2009_T1 Preguntas Teóricas

Tema 1.

  1. ¿Qué es un proceso estocástico?
  2. ¿Cuáles son las propiedades de una cadena de Markov?
  3. ¿Qué son tiempos de primera pasada de ir del estado i al j?
  4. ¿Qué es un estado absorbente?, ¿Qué es un estado recurrente?
  5.  ¿Qué propiedades debe cumplir una clase?

Referencia: FCNM/ICM01420

Proceso pasos aleatorios-caminata aleatoria

Ejemplo: Caminata de pasos aleatorios o Random Step Process

Referencia: Leon-García. E9.14 p.500

Un contador «up-down» o «sube-baja» genera pulsos +1 ó -1. La entrada del contador está dada por Dn = 2 In-1, donde In es un proceso aleatorio tipo Bernoulli.

D_n = \begin{cases} +1 & \quad \text{, } I_n = 1\\ -1 & \quad \text{, } I_n = 0 \end{cases}

Por ejemplo Dn representaría el cambio en la posición de una partícula que se mueve a lo largo de una linea recta, y cambia entre ±1 cada unidad de tiempo. Ej: n=20, p=0.5

Caminata aleatoria

La media de Dn es:

mD(n) = E[Dn] 
    = E[ 2In -1] = 2E[In]-1
    = 2p - 1

La varianza de Dn se encuentra como:

VAR[Dn] = VAR[2In - 1]
    = 22 VAR[In]
    = 4p(1 - p)

Las probabilidades de los eventos de Dn se calculan como en el ejemplo del tema de la Binomial.


Instrucciones en Python

# Proceso caminata de pasos aleatorios
# Leon-Garcia E 9.14 p.500
# propuesta: edelros@espol.edu.ec
import numpy as np
import matplotlib.pyplot as plt
import scipy.stats as stats

# INGRESO
# n=int(input('cuantos aleatorios: '))
# p=float(input('probabilidad p: '))
n = 20
p = 0.5

# PROCEDIMIENTO

# inicializa vectores
pasos  = np.zeros(n, dtype=int)
camina = np.zeros(n, dtype=int)
ejex   = np.zeros(n, dtype=int)

s = 0 
for i in range(0,n):
    # genera aleatorio con binomial
    pasos[i] = 2*(stats.binom.rvs(1,p))-1
    s = s+pasos[i]
    camina[i] = s
    ejex[i]   = i

# SALIDA
# grafica pasos 
plt.subplot(211)
plt.stem(ejex,pasos)
plt.ylabel('Pasos Dn')
plt.margins(0.05)

# grafica caminata aleatoria
plt.subplot(212)
plt.stem(ejex,camina)
plt.ylabel('Caminata Sn')
plt.margins(0.05)
plt.xlabel('n')

plt.show()

Referencia: León-García 9.14 p.500-501

Procesos Discretos en tiempo

La forma mas simple de un proceso aleatorio – independiente, con secuencias identicamente distribuidas (iid).

Procesos aleatorios iid

La secuencia Xn de un proceso aleatorio iid consiste en una secuencia de variables independientes, identicamente distribuidas con cdf Fx(x), media m y varianza σ2 .

mX(n) = E[Xn] = m     para todo n

lo que indica que su media es constante.

La función de autocovarianza , si n1 ≠ n2 :

Cx(n1, n2) = E[(Xn1 - m)(Xn2 - m)] =
= E[(Xn1 - m)] E[(Xn2 - m)] = 0       

dado que Xn1 y Xn2 son variables independientes y si n1 = n2 = n:

Cx(n1, n2) = E[(Xn - m)2] = σ2 

lo que permite expresarla como:

Cx(n1, n2) = σ2 δn1n2

donde δn1n2 = 1 si n1 = n2 y 0 para cualquier otro caso, que también expresa que su autocovarianza es cero en cualquier lugar exceptuando n1 = n2.

La función de autocorrelación de un proceso iid es:

Rx(n1, n2) = Cx(n1, n2) + m2

Ejemplo: Proceso aleatorio Bernoulli

Leon-García E9.13 pdf/p.499

La secuencia In es una secuencia de variable aleatoria independiente tipo Bernoulli.
Una ejecución se muestra en a figura. Se puede interpretar como el evento que un foco falle y se reemplace en un dia n.

Proceso aleatorio bernoulli

tiene media y varianza:

m1 = p     VAR[In] = p(1-p)

permite realizar cálculos de forma muy sencilla, por ejemplo, la probabilidad que los primeros cuatro bits tengan la secuencia 1001

p[I1=1,I2=0,I3=0,I4=1] = 
= p[I1=1]p[I2=0]P[I3=0]P[I4=1]=
= p2(1-p)2

o que la probabilidad de que el segundo bit sea 0 y el septimo sea 1

P[I2=0, I7=1] = P[I2=0]P[I7=1] = p(1-p)

El código python para generar el proceso bernoulli usado es:

% matplotlib inline
# Generar numeros aleatorios con distribución bernoulli
# Leon-Garcia E9.13 Proceso aleatorioBernoulli 
# propuesta: edelros@espol.edu.ec
import numpy as np
import matplotlib.pyplot as plt
import scipy.stats as stats

# INGRESO
# n=int(input('cuantos aleatorios: '))
# p=float(input('probabilidad p: '))
n=20
p=0.5

# Procedimiento
pasos=np.zeros(n, dtype=int)
camina=np.zeros(n, dtype=int)

# procesa los datos
s=0
for i in range(0,n):
    pasos[i]=stats.binom.rvs(1,p)
    s=s+pasos[i]
    camina[i]=s

# SALIDA
plt.subplot(211)
plt.stem(pasos)
plt.ylabel('bernoulli In')
plt.margins(0.05)

plt.subplot(212)
plt.stem(camina)
plt.ylabel('Sn')
plt.margins(0.05)
plt.show()

Combinatorias

Referencia: León-García p.42, Gubner p.35, Ross p.10

Combinatorias es el estudio de metodos sistemáticos de conteo, las cuatro clases de problemas principales son:

  1. Muestreo ordenado con reemplazo
  2. Muestreo ordenado sin reemplazo
  3. Muestreo no ordenado sin reemplazo
  4. Muestreo no ordenado con reemplazo

Muestreo ordenado con reemplazo

Se llaman variaciones con repetición de m elementos tomados de n en n, a los diferentes grupos que pueden formarse con los m elementos dados, tomados de n en n.
Dos grupos pueden ser distintos entre si, si tienen diferentes elementos en diferente orden.

El número de posibles k-tuplas en distinto orden (x1x2x3… xk) con elementos xi de un grupo de ni elementos diferentes es:

número de k-tuplas en orden distinto = n1n2n3… nk

Ejemplo : Posibles rutas para un paquete de internet

León-García E2.15 p.42.
Una urna tiene cinco pelotas numeradas de 1 al 5. Suponga que se sacan dos pelotas de la urna con reemplazo.

– ¿Cuántos pares diferentes de se pueden obtener?

pares diferentes = 5 * 5 = 52 = 25 pares diferentes

– ¿Cuál es la probabilidad que en se repitan las pelotas?

las formas de repetir son (1,1),(2,2),(3,3),(4,4) y (5,5), en total 5 de un total de 25 pares diferentes, por lo que la probabilidad será 5/25 = 1/5 = 0.20


Muestreo ordenado sin reemplazo

Se seleccionan k elementos en sucesi{on sin reemplazo de una población A de n elementos diferentes. Con k≤n, la primera vez se pueden escoger n1=n elementos diferentes, la segunda n2=n-1, la tercera n3=n-2, … en la última nk=n-(k-1)

número de k-tuplas en orden distinto = n(n-1)(n-2)…(n-k+1)

n(n-1)(n-2)...(n-k+1) = \frac{n!}{(n-k)!}

Ejemplo :

León-García E2.17 p.43.
Una urna tiene cinco pelotas numeradas de 1 al 5. Suponga que se sacan tres pelotas sin reemplazo .

– ¿Cuántos tripletas diferentes de se pueden obtener?

tripletas diferentes = 5 * 4 * 3 = 60 formas diferentes

(5!)/(5-3)! = 5!/2! = 5*4*3*(2!)/2! =5*4*3 = 60 formas diferentes

import scipy.special as sts

sts.perm(5,3)
60.0

Muestreo no ordenado sin reemplazo

Se sacan k elementos de un grupo A de n objetos diferentes sin reemplazo, y que se escriben los resultados sin importar e orden. Sería como colocarlos en otro conjunto B, el orden deja de importar.

En el nuevo conjunto B, existen k! formas ordenadas de seleccionar los objetos, y Ckn será el valor buscado de las combinaciones de tamaño k del conjunto A de n elementos.

C_k^n k! = n(n-1)(n-2)...(n-k+1) = = \frac{n!}{(n-k)!}

que simplificando se convierte el «coeficiente binomial» y se lee «de n toma k elementos»:

C_k^n = \frac{1}{n!} \frac{k!}{(n-k)!} = \left( \begin{array}{c} n \\ k \end{array} \right)

da lo mismo escoger k elementos del conjunto A, que dejar n-k elementos en el conjunto, por lo que:

\left( \begin{array}{c} n \\ k \end{array} \right) =\left( \begin{array}{c} n \\ n-k \end{array} \right)

Ejemplo :

Gubner E 1.38 p.38.

Se requiere conformar un jurado de 12 personas seleccionados de un total de 20 jueces. ¿Cuántas formas posibles existen para conformar el jurado?. No importa el orden.

\left( \begin{array}{c} 20 \\ 12 \end{array} \right) = \frac{20!}{12!8!} = 125970
import scipy.special as sts

sts.comb(20,12,repetition=False)
125970.0

Muestreo no ordenado con reemplazo

Se toman k objetos de un grupo de n objetos diferentes con reemplazo, se escribe el resultado sin importar el orden.

\left( \begin{array}{c} n-1+k \\ k \end{array} \right) = \left( \begin{array}{c} n-1+k \\ n-1 \end{array} \right)

Ejemplo : Dispensadora de frutas

(Gubner Ej 1.38 )

En una maquina dispensadora automática se entregan manzanas, bananas y peras. Por un precio fijo, se puden obtener cinco frutas seleccionadas por el cliente.

El proceso se maneja electrónicamente con una secuencia de 7 bits que los ceros (0) representan las manzanas, bananas y peras en orden y se separan por un bit uno (1) como en el ejemplo:

0100100 son un manzana, dos bananas y dos peras.

El primer grupo de 0’s es manzanas, el segundo grupo de 0’s son bananas y el grupo final de 0’s son peras.
¿Cuántas opciones tienen los clientes?

Solución: Equivale a preguntar cuántas secuencias de 7 bits hay con cinco ceros y dos unos.

\left( \begin{array}{c} 7 \\ 5,2 \end{array} \right) = \left( \begin{array}{c} 7 \\ 5 \end{array} \right) =\left( \begin{array}{c} 7 \\ 2 \end{array} \right)
import scipy.special as sts

c1 = sts.comb(7,2)
c2 = sts.comb(7,5)

print(c1)
print(c2)
21.0
21.0

3. Eventos Independientes

Referencia: León-García p.53, Gubner p.30, Parsen p.13/pdf.34, Ross p.10

Eventos Independientes

Conocer que al ocurrir un evento B no se altera la probabilidad de otro evento A, se dice que A es independiente de B. Para lo cual se cumple que:

P[ A \cap B ] = P[A] P[B]

que implica que:

P[A|B] = P[A]

y que:

P[B|A] = P[B]

Ejemplo: Transmisión de paquetes por Routers

Gubner E1.23 p.31

Un mensaje se transmite en forma de paquete de dato desde la ciudad de Guayaquil a Daule usando los “Router1”, un enlace de fibra óptica y el “router2” mostrados en la figura. Cada router puede descartar un paquete con una probabilidad p=0.01.

¿Cuál sería la probabilidad de transmitir con éxito un paquete entre el origen y destino?

Routers Guayaquil-Daule

Un paquete se transmite con éxito si y solo si ninguno de los routers descarta el paquete.

En lenguaje de eventos se dice que: descartar un paquete por el router 1 es D1 y para el router 2 es D2.

Sea el evento A cuando el paquete se transmite con éxito, ocurre solo cuando el paquete no se descarta en ningún router.

A = D_1^c \cap D_2^c

El problema indica que D1 y D2 son eventos independientes, por lo que D1c y D2c también son independientes, entonces:

P[A] = P[D_1^c \cap D_2^c] = P[D_1^c] P[D_2^c] = [1-P[D_1][1-P[D_2]] = (1-p)^2

dado que p=0.01, entonces P[A] = (1-0.01)2 = 0.9801

otra forma de ver el problema:

D1c D1
D2c (1-p)(1-p) (1-p)p
D2 (1-p)p p2

con números es:

D1c D1 PMarginal
D2c 0.9801 0,0099 0,99
D2 0,0099 0,0001 0,01
Pmarginal 0,99 0,01 1

Ejemplo: urnas con pelotas numeradas y de color

León-García E2.31 p.54

De una urna que contiene dos pelotas negras numeradas 1 y 2, y dos blancas numeradas 3 y 4, se obtienen una pelota.

Sean los eventos en los que se obtiene :

A = {(1,negra), (2,negra)}, «una pelota negra»
B = {(2,negra), (4,blanca)}, «una pelota numerada par»
C = {(3,blanca), (4,blanca)}, «una pelota con número mayor que 2»

a) Los eventos A y B ¿son independientes?
b) Los eventos A y C ¿son independientes?

Desarrollo :
a) las probabilidades de A y B son:

P[A] = P[B] = \frac{2}{4} = \frac{1}{2} P[A \cap B] = P[\{(2,negra)\}] = \frac{1}{4}

entonces:

P[A \cap B] = \frac{1}{4} = P[A]P[B]

por lo que A y B son independientes. Otra forma de escribirlo es:

P[A|B] = \frac{P[A \cap B]}{P[B]} = \frac{P[\{(2,negra)\}]}{P[\{(2,negra), (4,blanca)\}]} = = \frac{1/4}{1/2} = \frac{1}{2}

y para P[A]:

P[A] = \frac{P[A]}{P[S]} = \frac{P[\{(1,negra), (2,negra)\}]}{P[\{(1,negra), (2,negra),(3,blanca), (4,blanca) \}]} = = \frac{1/2}{1}

las ecuaciones implican que P[A]=P[B] debido que la proporción de las salidas en S tienen como resultado que A ocurre el mismo número de veces que B. Por lo que al conocer cuando ocurre B, no se altera la probavilidad de que ocurra A.

b) Los Eventos A y C no son independientes dado que:
P[A \cap C] = P[\emptyset] = 0
P[A|C] = 0 \neq P[A] = \frac{1}{2}

A y C son mutuamente excluyentes dado que A∩C=∅, por lo que ocurra C implica que A no ha ocurrido de forma definitiva.


En general, si dos eventos tienen probabilidad diferente de cero y son mutuamente excluyentes, no pueden ser independientes.

Si fueran independientes y mutiamente excluyentes:

0=P[A∩B] =P[A]P[B]

lo que implica que al menos uno de los eventos deben tener probabilidad cero.