2.1 Valor Esperado E[x] a partir de pmf – Ejercicio

Referencia: LoRaWan – Descriptor estadístico de un punto en Girni,

El indicador de fuerza de la señal recibida (RSSI) para un dispositivo a una distancia desde el gateway es un valor entero que se encuentra alrededor el valor esperado o media. Las condiciones del ambiente generan un componente aleatorio que se muestra en la figura.

Las mediciones se pueden obtener para el enlace de subida (dispositivo a gateway) o enlace de bajada (gateway a dispositivo). A partir de los valores se pueden obtener la pmf del enlace de subida o de bajada

Ejercicio

A partir de los datos de los datos obtenidos para un punto con ubicación  etiquetada LOS04, calcular el valor esperado para los enlaces de subida (UP) y de bajada (Down) a partir de los siguientes datos, presentados en forma de un diccionario de Python:

{"pmf":
 {"pmf_down":
  {"intervalo":[-70,-69,-68,-67,-66,-65,-64],
   "freq_relativa":[0.036101083,0.1371841155,0.1732851986,0.1624548736,0.2707581227,0.202166065,0.0180505415]
   },
  "pmf_up":
  {"intervalo":[-83,-82,-81,-80,-79,-78,-77,-76,-75,-74,-73,-72,-71,-70,-69,-68,-67,-66,-65,-64,-63,-62,-61,-60,-59,-58,-57],
   "freq_relativa":[0.036101083,0.0,0.0252707581,0.0,0.0072202166,0.0,0.1010830325,0.0036101083,0.1083032491,0.0,0.1227436823,0.0,0.0433212996,0.0324909747,0.0180505415,0.0252707581,0.0685920578,0.1335740072,0.0288808664,0.0,0.1046931408,0.1046931408,0.0180505415,0.0072202166,0.0072202166,0.0,0.0036101083]
   }
  },
 "dispositivo":{"pmf_down":"cc24","pmf_up":"cc24"}
 }

PAM a PSK – valor esperado

Y(t) = a \cos \big( 2\pi t + \frac{\pi}{2}X(t) \big)

La pmf de x(t) es 0.5 para cada valor de [-1,1]

c) Encuentre la media y autocorrelacion de Y(t)

E[Y(t)] = E\big[ a \cos \big( 2\pi t + \frac{\pi}{2}X \big)\big]

Referencia:  León-García 3.3.1 p. 107. Valor esperado de funciones de variable aleatoria

Si z =g(x)

E[g(x)] = \sum_k g(x_k)p_x(X_k)

se tiene entonces que:

= \big[a \cos \big( 2\pi t + \frac{\pi}{2}(-1) \big)\big]\frac{1}{2} + \big[a \cos \big( 2\pi t + \frac{\pi}{2}(1) \big)\big] \frac{1}{2} = \frac{a}{2}\cos \big( 2\pi t - \frac{\pi}{2} \big) + \frac{a}{2}\cos \big( 2\pi t + \frac{\pi}{2} \big) = \frac{a}{2} \sin \big( 2\pi t \big) - \frac{a}{2} \sin \big( 2\pi t \big) = 0 E[Y(t)] = 0

Media, Varianza correlacion

Procesos aleatorios Discretos en tiempo

Media:
m_X(n) = E[X(n)]
Varianza:
VAR[X(n)] = E[(X(n) - m_X(n))^2]

Autocorrelación
R_X(n_1,n_2) = E[X(n_1,n_2)]

Autocovarianza:
C_X(n_1,n_2) = E[\{X(n_1) - m_X(n_1) \} \{ X(n_2) - m_X(n_2)\}
C_X(n_1,n_2) = R_X(n_1,n_2) - m_X(n_1) m_X(n_2)

1. Valor Esperado de una función de variable aleatoria

Referencia: Gubner 2.4 p83 , Ross 2.4.3 p42, León-García 3.3.1 p.107

Dada una variable aleatoria X, se puede definir una nueva variable aleatoria Z = g(X), donde g(x) es una función de valor real de la variable real x.

Para calcular E[Z] se puede proceder como:

E[g(X)] = \sum_i g(x_i) p_X(x_i) E[g(x)] = \int_{-\infty}^{\infty} g(x) f(x) dx

dado que la fórmula es mas fácil de usar que encontrar la pmf de Z, la formula se la conoce como la «ley del estadístico inconsciente» o LOTUS (Law Of The Unconscious Statistician).

Una aplicación simple es :

E[aX] = \sum_i ax_i p_X(x_i) = = a \sum_i x_i p_X(x_i) = a E[X]

Ejemplo

Referencia: León- García 3.17 p107

Sea X el ruido en el voltaje que está uniformemente distribuido en SX = {-3,-1,+1,+3} con pX (k) =  1/4 para k en SX. Encuentre E[Z] donde Z=X2.

Solución: Se busca primero encontrar la pmf (probability mass function) de Z, el SZ ={9,1,1,9} = {1,9}, por lo que:

pZ(9) = P[X ∈ {-3,+3}] 
      = pX(-3) + pX(3) 
      = 1/4 + 1/4 = 1/2
pZ(1) = pX(-1) + pX(1) = 
      = 1/4 + 1/4 = 1/2
entonces:
E[Z] = 1(1/2) + 9(1/2) = 5 

usando la fórmula para E[Z]:

E[Z] = E[g(X)] = \sum_i g(x_i)p_X(x_i) = \sum_i i^2 p_X(x_i) = \frac{1}{4} [(-3)^2 + (-1)^2+1^2+2^2] = = \frac{20}{4} = 5

con lo que se obtuvo el mismo resultado.


Ross Corolario 2.2. Siendo a y b constantes, entonces:

E[aX + b] = aE[X] +b

3. Valor Esperado de variables aleatorias contínuas

Referencia: Ross 2.4.2 p39, Gubner 4.2 p149, León-García 4.3 p 155, p16

Si X es una variablea aleatoria contínua que tiene una función densidad de probabilidad f(x), el valor esperado de X se define como:

E[X]= \int_{-\infty}^{\infty}xf(x) \delta x

Ejemplo

Referencia: Ross 4.7.

Sea X una variable aleatoria contínua uniforme [a,b], encuentre el valor esperado:

Solución:

E[X]= \int_{-\infty}^{\infty} xf(x) \delta x = \int_{a}^{b} x\frac{1}{b-a}\delta x = \left. = \frac{1}{(b-a)} \frac{x^2}{2} \right|_{a}^{b} = \frac{1}{(b-a)} \frac{b^2-a^2}{2} = = \frac{(b+a)(b-a)}{2(b-a)} = \frac{a+b}{2}

que es el promedio simple entre a y b cuando la función es uniforme.


Ejemplo

Ross 4.8/Leon-García 4.20. Un convertidor analógico-digital o cuantizador con resolución de paso Δ voltios redondea en la entrada el valor más cercano al múltiplo de Δ voltios como se muestra en la figura.

La entrada es un voltaje Vin de tipo aleatorio y la salida del convertidor es A/D es Vout, y su desempeño se mide por el error cuadratico medio:

E[|Vin – Vout|2]

Se supone que el error Vin – Vout se puede aproximar a una función aleatoria uniforme [-Δ/2,Δ/2] dado que siempre el valor siempre cae en el intervalo dado. Determine el valor esperado para la señal de entrada, y también el valor del error cuadrático medio.

Solución: Para el intervalo centrado en el origen

f(x)= \frac{1}{b-a} E[X] = \int_{-\Delta /2}^{\Delta /2} xf(x) \delta x = =\int_{-\Delta /2}^{\Delta /2} x\frac{1}{\Delta /2 - (-\Delta /2)} \delta x = \int_{-\Delta /2}^{\Delta /2} \frac{x}{\Delta} \delta x = \left. \frac{x^2}{2 \Delta} \right|_{-\Delta /2}^{\Delta /2} = \frac{1}{2\Delta} \big[ {\big( \frac{-\Delta}{2}\big)}^2 - {\big( \frac{\Delta}{2}\big)}^2 \big] = 0

Para el error cuadrático medio en [a,b], a=-Δ/2, b=Δ/2

E[|V_{in}-V_{out}|^2] \approx E[X^2] = \int_{-\infty}^{\infty} x^2 f(x) \delta x = \int_{a}^{b} x^2 \frac{1}{b - a} \delta x = \left. \frac{1}{b-a} \frac{x^3}{3} \right|_{a}^{b} = \frac{1}{b-a} \frac{b^3-a^3}{3} = = \frac{(b-a)(b^2+ba+a^2)}{3(b-a)} = \frac{b^2+ba+a^2}{3} = \frac{({\frac{\Delta}{2})}^2 + {(\frac{\Delta}{2})}{(\frac{-\Delta}{2})} + {(\frac{-\Delta}{2})}^2}{3} = \frac{1}{3} \frac{\Delta ^2 -\Delta^2 + \Delta ^2}{4} = \frac{\Delta ^2}{12}

Teoría de Colas – Cajeras María y Alicia

Referencia: Ross problema 8.3 p.568

Un administración de un supermercado puede contratar a María o Alicia.

María atiende con una tasa de servicio exponencial de 20 clientes por hora y puede ser contratada por $3/hora.

Alicia atiende con una tasa exponencial de 30 clientes por hora, y puede ser contradada por $C/hora.

El administrador estima que en promedio cada cliente tiene un costo de $1 por hora y debería ser contabilizado para el modelo.

Suponga que los clientes llegan a una tasa Poisson de 10 clientes por hora.

a) ¿Cuál es el costo promedio por hora si María es contratada? o ¿Si Alicia es contratada?

b) Encuentre C si el promedio de costo por hora es el mismo para María y Alicia.


Solución:

a) ¿Cuál es el costo promedio por hora si María es contrada? o ¿Si Alicia es contratada?

λ = 10 clientes/hora

María: 

μMaría = 20 clientes/hora

LMaría = λ/(μ -λ) = 10/(20-10) = 1 

Costo/hora María =  $3 + ($1 * (número promedio de clientes en la cola cuando trabaja María))
      =  $3 + ($1 * LMaría)
      =  $3 + ($1 * 1) = $ 4 

Alicia:
μAlicia = 30 clientes/hora
LAlicia = λ/(μ -λ) = 10/(30-10) = 1/2

Costo/hora Alicia =  $C + ($1 * (número promedio de clientes en la cola cuando trabaja Alicia))
      =  $C + ($1 * LAlicia)
      =  $C + ($1 * 1/2) = $ (C+ 0.5)
      
b) Encuentre C si el promedio de costo por hora es el mismo para María y Alicia.
Costo/hora María = Costo/hora Alicia
 4 = C + 0.5
 C = 3.5 $/hora

 que es lo máximo que el administrador podría pagar a Alice sin incurrir en costos mayores que María, dadas las condiciones del problema.

Teoria de colas – Llegadas Esperadas

Referencia: Ross problema 8.1 p.568

Para una cola M/M/1, calcule

a) el numero esperado de llegadas durante un periodo de servicio
b) la probabilidad que no lleguen clientes durante un periodo de servicio
pista: «condición»


Solución:
a) el numero esperado de llegadas durante un periodo de servicio

E[numero de llegadas] = E[λ S] = λ[1/μ] = λ/μ

b) la probabilidad que no lleguen clientes durante un periodo de servicio

P{0 llegadas} = E[P{0 llegadas | periodo de servicio S}] =
   = E[P{N(S)=0}] = E[e-λS]
E[e^{-\lambda S}] = \int_{0}^{x} e^{-\lambda S}\mu e^{-\mu S}dS = \mu \int_{0}^{\infty} e^{-\lambda S -\mu S}dS = \mu \int_{0}^{\infty} e^{-S(\lambda+\mu)}dS = \left.\frac{-\mu}{\lambda+\mu} e^{-S(\lambda+\mu)}\right|_0^{\infty} = \frac{-\mu}{\lambda+\mu}\left[ e^{-\infty(\lambda+\mu)}-e^{-0(\lambda+\mu)}\right] = \frac{-\mu}{\lambda+\mu} [0-1] = \frac{\mu}{\lambda+\mu}

Teoría de colas – Dos servidores en paralelo

Referencia: Ross 8.17 p.571

Los clientes llegan a una sucursal de dos servidores según a un proceso de Poisson con tasa de dos por hora.

Si al llegar el cliente el servidor 1 esta libre, se atiende con éste.
Si al llegar el cliente el servidor 1 esta ocupado y el servidor 2 esta libre, inicia con el servidor 2.
Clientes que al llegar encuentran los dos servidores ocupados se retiran.

Un cliente que fué atendido por el servidor 1, pasa a ser atendido por el servidor 2 siempre que esté libre, caso contrario sale de la sucursal.
Un cliente al completar el servicio con el servidor 2 sale de la sucursal.
Los tiempos de servicio en los servidores 1 y 2 son variables aleatorias exponenciales con tasas de cuatro y seis por hora.

a) ¿Qué fracción de clientes no entran a la sucursal?

b) ¿Cuál es el valor promedio de tiempo que un cliente que entra, permanece en el sistema?

c) ¿Cuál fracción de los clientes que entraron son atendidos por el servidor 1?


Solución:
Se definen los estados como (servidor1,servidor2), ocupado=0, ocupado=1:

Estado (0,0)
Estado (0,1)
Estado (1,0)
Estado (1,0)

con lo que se contruye el diagrama de estados y transiciones:

Usando los valores del enunciado se ubican los valores para las tasas de llegada y atención

El siguiente paso consiste en escribir son las ecuaciones de balanceo:

λ P00 = μ2 P012 + λ) P01 = μ1 P10 + μ1 P111 + λ) P10 = λ P00 + μ2 P112 + μ1) P11 = λ P01 + λ P10
1 = P00 + P01 + P10 + P11

usando los valores para λ, μ1 y μ2:

2 P00 = 6 P01
8 P01 = 4 P10 + 4 P11
6 P10 = 2 P00 + 6 P11
10 P11 = 2 P01 + 2 P10
1 = P00 + P01 + P10 + P11

de la ecuación(1)

P01 = (2/6) P00 = (1/3) P00
P01 = (1/3) P00

usando (2) (1/4) y reordenando
2 P01 = P10 + P11
P11 = 2 P01 - P10

sumando con (4)
  P11 = 2 P01 - P10
5 P11 =   P01 + P10

6 P11 = 3 P01
2 P11 = P01
usando (1)
2 P11 = (1/3) P00
  P11 = (1/6) P00

usando (4)
5 P11 =   P01 + P10
5 (1/6) P00 = (1/3) P00 + P10
  P10 = [(5/6) - (1/3)]P00 = [(5-2)/6] P00 = (3/6) P00 = (1/2) P00
  P10 = (1/2) P00 

usando (5) con los resultados dependientes de P00:

1 = P00 + (1/3) P00 +(1/2) P00 + (1/6) P00
1 = [(6+2+3+1)/6] P00 = [12/6] P00 = 2 P00
P00 = 1/2

quedando:
P00 = 1/2
P01 = 1/6
P10 = 1/4
P11 = 1/12

usando numpy de python, se reorganiza las ecuaciones y se crean las matrices:

2 P00 - 6 P01                 = 0
        8 P01 - 4 P10 - 4 P11 = 0
 2 P00        - 6 P10 + 6 P11 = 0
        2 P01 + 2 P10 - 10 P11 = 0
   P00 +  P01 +   P10 +    P11 = 1
import numpy as np
A=np.array([
    [2,-6,0,0],
    [0,8,-4,-4],
    [2,0,-6,6],
    [1,1,1,1]])
B=np.array([0,0,0,1])
P=np.linalg.solve(A,B)
print(P)
[ 0.5         0.16666667  0.25        0.08333333]

que son los resultados anteriores encontrados:

P00 = 0.5
P01 = 0.16666667
P10 = 0.25
P11 = 0.08333333

a) ¿Qué fracción de clientes no entran a la suscursal?

solo ocurre cuando ambos servidores estan ocupados (1,1)

P11 = 1/12 = 0.08333333

b) ¿Cuál es el valor promedio de tiempo que un cliente que entra, permanece en el sistema?

W = L/λ

  = (valor esperado)/(proporcion de los clientes que si entran)
  
  = [0*P00 + 1*P01 + 1*P10 + 2*P11]/ [λ(1-P11)]
  
  = [0*(1/2) + 1*(1/6) + 1*(1/4) + 2*(1/12)]/ [2(1-(1/12)]
  
  = [(0+3+2+2)/12]/[22/12] = (7/12)/(22/12) =  7/22 

c) ¿Cuál fracción de los clientes que entraron son atendidos por el servidor 1?

(atiende el servidor 1 cuando esta libre) / (los que entraron)

[P00 + P01] / [1 - P11] = [1/2 + 1/6]/[1 - 1/12]
  
  = [(3+1)/6] / [11/12] = (4/6)/(11/12) =   8/11 

Teoría de colas – Dos Servidores en secuencia

Referencia: Ross problema 8.15 p.570

Considere un sistema de servicio secuencial de dos servidores A y B.
Los clientes que llegan entran al sistema solo si el servidor A esta libre.
Si un cliente entra, entonces es atendido inmediatamente por el servidor A.
Cuando la atención del servidor A se completa, el cliente pasa a ser atendido por el servidor B siempre que esté libre, si B esta ocupado, el cliente sale del sistema.
Una vez que el cliente sale del servidor B, el cliente se parte del sistema.

Suponga que las llegadas de cientes son tipo Poisson con tasa de llegada de dos clientes por hora, y que los servidores A y B atienden a tasas exponenciales de cuatro y dos clientes por hora.

a) ¿Cuál es la proporción de clientes que entran al sistema?

b) ¿Cuál es la proporción de clientes que entraron al sistema son atendidos por el servidor B?

c) ¿Cuál es el número promedio de clientes en el sistema?

d) ¿Cuál es el monto promedio de tiempo que un cliente que entró se mantiene dentro del sistema?


Solución

los estados propuestos son:

00 P00
10 P10
01 P01
11 P11

El diagrama de transición entre estados propuesto, considerando de un solo evento por unidad de tiempo:

con lo que se puede plantear las ecuaciones de balanceo («sale=entra»):

2 P00 = 2 P01
4 P10 = 2 P00 + 2 P11
(2+2) P01 = 4 P10 + 4 P11
(4+2) P11 = 2 P01
P00 + P10 + P01 + P11 = 1

Resolviendo:

P00 = P01              ecuación(1)
2 P10 = P00 + P11      (2)
P01 = P10 + P11        (3)
3 P11 = P01            (4)
P00 + P10 + P01 + P11 = 1      (5)

de la ecuación (2) restando la ecuación (3)

2 P10 - P01 = P00 + P11 - P10 - P11
2 P10 - P01 = P00 - P10
3 P10 - P01 = P00

usando la ecuación (1)

3 P10 - P00 = P00
3 P10 = 2 P00
P10 = 2/3 P00

usando la ecuación (4) y (1)

P11 = (1/3) P01
P11 = (1/3) P00

usando la ecuación (5)

P00 + P10 + P01 + P11 = 1
P00 + P00 + (2/3) P00 + (1/3)P00 = 1
3 P00 = 1
P00 = 1/3

y se obtienen los otros valores para p

P01 = 1/3
P10 = (2/3)(1/3 )= 2/9
P11 = (1/3)(1/3) = 1/9

Con lo que se puede responder:

a) ¿Cuál es la proporción de clientes que entran al sistema?

P00 + P01 = 1/3 + 1/3 
         = 2/3

c) ¿Cuál es el número promedio de clientes en el sistema?

0*P00 +  1*P10 + 1*P01 + 2*P11 = 0 + 2/9 + 1/3 + 2*1/9 = 
     = (0+2+3+2)/9 = 7/9

d) ¿Cuál es el monto promedio de tiempo que un cliente que entró se mantiene dentro del sistema?

W = L/λA
W = (promedio de clientes en el sistema)/(tasa de los que entran)
W= (7/9)/[2(2/3)] = 7/12

Usando Numpy

resolviendo con numpy, se convierten las ecuaciones a matrices, sustituyendo una ecuación con la última (suma de probabilidades =1), se obtiene:

2 P00 = 2 P01
4 P10 = 2 P00 + 2 P11
(2+2) P01 = 4 P10 + 4 P11
(4+2) P11 = 2 P01
P00 + P10 + P01 + P11 = 1

reorganizando:

2 P00 -2 P01                = 0
2 P00        - 4 P10 + 2 P11 = 0
      -4 P01 + 4 P10 + 4 P11 = 0
       2 P01        - 6 P11 = 0
  P00  +  P01 +   P10 +  P11 = 1

se obtienen las matrices A y B para usar con la instrucción solver de numpy:

import numpy as np
A=np.array([
    [2,-2,0,0],
    [2,0,-4,2],
    [0,-4,4,4],
    [1,1,1,1]])
B=np.array([0,0,0,1])

p=np.linalg.solve(A,B)

print(p)
[ 0.33333333  0.33333333  0.22222222  0.11111111]

con lo que los valores buscados son:

P00=0.33333333
P01=0.33333333
P10=0.22222222
P00=0.11111111

2. Cadena Markov – Ejemplo: Clima con dos días

Referencia: Ross p.193, 197

(4.4) Pronostico del clima con dos días

Suponga que si llueve o no hoy depende de las condiciones del clima de los dos últimos días.

Si ha llovido en los dos últimos días, entonces que llueva mañana tiene probabilidad de 0.7;
si llovió hoy pero no ayer, entonces que llueva mañana tiene probabilidad de 0.5;
si llovio ayer pero no hoy, entonces que llueva mañana tiene probabilidad de 0.4;
si no ha llovido ni ayer ni hoy, entonces que llueva mañana tiene probabilidad de 0.2

Solución propuesta

Se puede hacer un modelo de cadena de Markov haciendo que el estado del clima mañana sea determinado por las condiciones de hoy y ayer:

Estado 0: si ha llovido hoy y ayer
Estado 1: si ha llovido hoy pero no ayer
Estado 2: Si ha llovido ayer pero no hoy
Estado 3: si no ha llovido ayer y tampoco hoy

Para considerar los posibles estados añadiendo los valores cuando está soleado, dado que el problema enuncia solo cuando lloverá.

con lo que se puede construir el diagrama de transición de estados:

y la matriz de transición es entonces:

p= \begin{pmatrix} 0.7 & 0 & 0.3 & 0 \\ 0.5 & 0 & 0.5 & 0 \\ 0 & 0.4 & 0 & 0.6 \\ 0 & 0.2 & 0 & 0.8 \end{pmatrix}

(4.9) Pronóstico para pasado mañana

Siguiendo con el ejemplo anterior.
Dado que llovió el Lunes y Martes, ¿Cuál es la probabilidad que lloverá el Jueves?

Solución propuesta

La matriz de transición permite hacer el pronóstico hasta el miércoles, para el día jueves es necesario hacer una matriz de dos pasos dada por p(2)

p^{(2)}= \begin{pmatrix} 0.7 & 0 & 0.3 & 0 \\ 0.5 & 0 & 0.5 & 0 \\ 0 & 0.4 & 0 & 0.6 \\ 0 & 0.2 & 0 & 0.8 \end{pmatrix} \text{.} \begin{pmatrix} 0.7 & 0 & 0.3 & 0 \\ 0.5 & 0 & 0.5 & 0 \\ 0 & 0.4 & 0 & 0.6 \\ 0 & 0.2 & 0 & 0.8 \end{pmatrix} = \begin{pmatrix} 0.49 & 0.12 & 0.21 & 0.18 \\ 0.35 & 0.20 & 0.15 & 0.30 \\ 0.20 & 0.12 & 0.20 & 0.48 \\ 0.10 & 0.16 & 0.10 & 0.64 \end{pmatrix}

Que llueva el jueves es equivalente no importa si llovió o nó el miércoles, por lo que para la probabilidad del jueves se usa la probabilidad de la fila para estado 0, sumando si llovió o no el miércoles.

p(lloverá jueves|llovió lunes y martes)= p(2)00+p(2)01= = 0.49 + 0.12 = 0.61.


Usando python para determinar p(2)

# 4.9 pronostico para pasado mañana
import numpy as np

# INGRESO
p = np.array([[0.7,   0, 0.3, 0  ],
              [0.5,   0, 0.5, 0  ],
              [  0, 0.4,   0, 0.6],
              [  0, 0.2,   0, 0.8]])

#PROCEDIMIENTO
p2 = np.linalg.matrix_power(p,2)

# SALIDA
print(p2)
print('Para el jueves:', p2[0,0]+p2[0,1])
[[ 0.49  0.12  0.21  0.18]
 [ 0.35  0.2   0.15  0.3 ]
 [ 0.2   0.12  0.2   0.48]
 [ 0.1   0.16  0.1   0.64]]
Para el jueves: 0.61