Salta al contenuto
Note per Studenti Esercizio - esponenziale più variabile discreta e probabilità condizionate (esercizi del corso, cap. 2, es. 9)

Esercizio - esponenziale più variabile discreta e probabilità condizionate (esercizi del corso, cap. 2, es. 9)

In questa pagina 3

Testo (raccolta di esercizi del corso, capitolo 2, esercizio 9). Sia yy una variabile esponenziale di parametro λ=2\lambda=2 e xx una variabile aleatoria (indipendente da yy) che assume valore 11 con probabilità p=0,4p=0{,}4 e valore 22 con probabilità 1−p=0,61-p=0{,}6. Sia z=x+yz=x+y.

  1. Scrivere l'espressione di P(z>A∣x=2)P(z>A\mid x=2) e calcolarne il valore per A=1A=1 e A=3A=3.
  2. Scrivere l'espressione di P(z>A∣x=2) px(2)+P(z<A∣x=1) px(1)P(z>A\mid x=2)\,p_x(2)+P(z<A\mid x=1)\,p_x(1) e calcolarne il valore per A=−2A=-2.

Teoria usata: Variabili aleatorie continue - uniforme, gaussiana, esponenziale e funzione QUna v.a. continua ha funzione di distribuzione F_x(a) = P(x ≤ a) e densità f_x(a) = dF_x/da ≥ 0 con integrale 1; P(A ≤ x ≤ B) = F_x(B) − F_x(A) = ∫A^B f_x e P(x = A) = 0. Uniforme U([a,b]): densità 1/(b−a), media (a+b)/2, varianza (b−a)²/12. Esponenziale unilatera di media m: (1/m)e^(−a/m), a ≥ 0. Gaussiana N(m, σ²): P(x ≥ A) = Q((A − m)/σ), con Q(A) coda della gaussiana standard. Somme di gaussiane indipendenti sono gaussiane (medie e varianze si sommano); teorema del limite centrale.Variabili aleatorie continue - uniforme, gaussiana, esponenziale e funzione Q →, Probabilità per le telecomunicazioni - eventi e variabili aleatorie discreteGli eventi hanno probabilità tra 0 e 1 e P(x∪y) = P(x) + P(y) − P(x∩y). Una variabile aleatoria discreta ha un alfabeto A_x e una densità p_x(a) = P(x = a) con somma 1. Si definiscono probabilità condizionata, indipendenza (p{x,y}(a,b) = p_x(a)p_y(b)), valore atteso E(x) = m_x, potenza M_x = E(x²), varianza σ_x² con M_x = m_x² + σ_x²; l'aspettazione è lineare e per variabili indipendenti E(xy) = E(x)E(y).Probabilità per le telecomunicazioni - eventi e variabili aleatorie discrete →, Formula delle probabilità totali e formula di BayesSe (A_i) è una partizione di Ω, P(B) = Σ P(B ∣ A_i) P(A_i) (probabilità totali); la formula di Bayes inverte il condizionamento: P(A_k ∣ B) = P(B ∣ A_k) P(A_k) / P(B).Formula delle probabilità totali e formula di Bayes →.

Convenzione. "Esponenziale di parametro λ\lambda" indica la densità fy(a)=λe−λaf_y(a)=\lambda e^{-\lambda a} per a≥0a\ge0 e 00 per a<0a<0 (media 1λ=0,5\frac1\lambda=0{,}5): è la stessa forma della lezione, con media m=1λm=\frac1\lambda. La coda è P(y>a)={1a<0e−λaa≥0 = e−λmax⁡{a,0}.P(y>a)=\begin{cases}1&a<0\\ e^{-\lambda a}&a\ge0\end{cases}\ =\ e^{-\lambda\max\{a,0\}}. Il massimo nell'esponente serve perché y≥0y\ge0 sempre: per soglie negative la probabilità di superarle è 11.

(1) P(z>A∣x=2)P(z>A\mid x=2)

Condizionando a x=2x=2 la variabile xx non è più aleatoria e z=2+yz=2+y. Inoltre xx e yy sono indipendenti, quindi il condizionamento non cambia la densità di yy. Allora P(z>A∣x=2)=P(2+y>A)=P(y>A−2)=e−λmax⁡{A−2,0}.P(z>A\mid x=2)=P(2+y>A)=P(y>A-2)=e^{-\lambda\max\{A-2,0\}}.

  • Per A=1A=1: A−2=−1<0A-2=-1<0, il massimo vale 00, e P=e0=1P=e^{0}=1. Ha senso: z=2+y≥2>1z=2+y\ge2>1 sempre.
  • Per A=3A=3: A−2=1A-2=1 e P=e−2⋅1=e−2=0,1353P=e^{-2\cdot1}=e^{-2}=0{,}1353.

(2) Una probabilità di decisione corretta

Per P(z<A∣x=1)P(z<A\mid x=1) si procede allo stesso modo con x=1x=1, z=1+yz=1+y: P(z<A∣x=1)=P(1+y<A)=P(y<A−1)=1−e−λmax⁡{A−1,0}.P(z<A\mid x=1)=P(1+y<A)=P(y<A-1)=1-e^{-\lambda\max\{A-1,0\}}. Con px(2)=0,6p_x(2)=0{,}6 e px(1)=0,4p_x(1)=0{,}4 l'espressione richiesta è P(z>A∣x=2) px(2)+P(z<A∣x=1) px(1)=0,6 e−λmax⁡{A−2,0}+0,4(1−e−λmax⁡{A−1,0}).P(z>A\mid x=2)\,p_x(2)+P(z<A\mid x=1)\,p_x(1)=0{,}6\,e^{-\lambda\max\{A-2,0\}}+0{,}4\Big(1-e^{-\lambda\max\{A-1,0\}}\Big). Per A=−2A=-2: max⁡{A−2,0}=max⁡{A−1,0}=0\max\{A-2,0\}=\max\{A-1,0\}=0, quindi 0,6⋅1+0,4⋅(1−1)=0,60{,}6\cdot1+0{,}4\cdot(1-1)=0{,}6.

Significato. La somma è la probabilità di decisione corretta di un ricevitore a soglia che osserva z=x+yz=x+y e decide "x=2x=2" se z>Az>A e "x=1x=1" se z<Az<A: per la formula delle probabilità totali, P(corretta)=P(decido 2, x=2)+P(decido 1, x=1)P(\text{corretta})=P(\text{decido }2,\ x=2)+P(\text{decido }1,\ x=1), cioè P(z>A∣x=2)P(x=2)+P(z<A∣x=1)P(x=1)P(z>A\mid x=2)P(x=2)+P(z<A\mid x=1)P(x=1). Con A=−2A=-2 la soglia è sotto tutti i valori di zz, quindi si decide sempre x=2x=2 e si indovina solo quando x=2x=2, cioè con probabilità 0,60{,}6 (la probabilità a priori). È utile guardare tutta la curva in AA (vedi anche l'esempio sotto).

Grafico interattivo: Probabilità di decisione corretta P(A) = 0,6·e^(−2·max(A−2,0)) + 0,4·(1 − e^(−2·max(A−1,0))): vale 0,6 per A ≤ 1 (si decide sempre x = 2), cresce fino al massimo 0,946 in A = 2 e poi scende verso 0,4 (si decide sempre x = 1)

  • Per A≤1A\le1: z=x+y≥1≥Az=x+y\ge1\ge A, quindi z>Az>A quasi sempre e si decide sempre 22: P=0,6P=0{,}6.
  • Per 1≤A≤21\le A\le2: z=2+y>2≥Az=2+y>2\ge A sempre se x=2x=2 (probabilità 0,60{,}6 intera), mentre se x=1x=1 si decide bene quando 1+y<A1+y<A: P(A)=0,6+0,4(1−e−2(A−1))P(A)=0{,}6+0{,}4(1-e^{-2(A-1)}), crescente in AA.
  • Per A≥2A\ge2 si comincia a perdere anche i casi x=2x=2: P(A)=0,6e−2(A−2)+0,4(1−e−2(A−1))P(A)=0{,}6e^{-2(A-2)}+0{,}4(1-e^{-2(A-1)}), decrescente. Il massimo è in A=2A=2, con P=0,6+0,4(1−e−2)=0,946P=0{,}6+0{,}4(1-e^{-2})=0{,}946.

(Verifica con Python, anche con 10610^6 estrazioni simulate: P(z>1∣x=2)=1P(z>1\mid x=2)=1, P(z>3∣x=2)=0,1353P(z>3\mid x=2)=0{,}1353, P(A=−2)=0,6P(A=-2)=0{,}6, P(A=2)=0,9459P(A=2)=0{,}9459; le simulazioni danno 0,13510{,}1351, 0,59990{,}5999 e 0,94590{,}9459. La soluzione del corso riporta 11, e−2=0,14e^{-2}=0{,}14 e 0,60{,}6.)

Errori comuni

  • Dimenticare il max⁡{⋅,0}\max\{\cdot,0\} e scrivere e−λ(A−2)e^{-\lambda(A-2)} anche per A<2A<2 (verrebbe maggiore di 11).
  • Usare λ\lambda come media: qui la media è 1λ=0,5\frac1\lambda=0{,}5.
  • Nel punto 2, usare P(z>A)P(z>A) senza condizionare, o dimenticare i pesi px(2)p_x(2) e px(1)p_x(1).
  • Scrivere la coda di yy per zz: P(z>A∣x=2)P(z>A\mid x=2) dipende da A−2A-2, non da AA.

Versione ripasso

Testo. yy esponenziale di parametro λ=2\lambda=2 (densità λe−λa\lambda e^{-\lambda a}, media 12\frac12), xx indipendente con P(x=1)=0,4P(x=1)=0{,}4, P(x=2)=0,6P(x=2)=0{,}6, z=x+yz=x+y. (1) P(z>A∣x=2)P(z>A\mid x=2) per A=1,3A=1,3. (2) P(z>A∣x=2)px(2)+P(z<A∣x=1)px(1)P(z>A\mid x=2)p_x(2)+P(z<A\mid x=1)p_x(1) per A=−2A=-2.

(Variabili aleatorie continue - uniforme, gaussiana, esponenziale e funzione QUna v.a. continua ha funzione di distribuzione F_x(a) = P(x ≤ a) e densità f_x(a) = dF_x/da ≥ 0 con integrale 1; P(A ≤ x ≤ B) = F_x(B) − F_x(A) = ∫_A^B f_x e P(x = A) = 0. Uniforme U([a,b]): densità 1/(b−a), media (a+b)/2, varianza (b−a)²/12. Esponenziale unilatera di media m: (1/m)e^(−a/m), a ≥ 0. Gaussiana N(m, σ²): P(x ≥ A) = Q((A − m)/σ), con Q(A) coda della gaussiana standard. Somme di gaussiane indipendenti sono gaussiane (medie e varianze si sommano); teorema del limite centrale.Variabili aleatorie continue - uniforme, gaussiana, esponenziale e funzione Q →, Formula delle probabilità totali e formula di BayesSe (A_i) è una partizione di Ω, P(B) = Σ P(B ∣ A_i) P(A_i) (probabilità totali); la formula di Bayes inverte il condizionamento: P(A_k ∣ B) = P(B ∣ A_k) P(A_k) / P(B).Formula delle probabilità totali e formula di Bayes →)

  • Coda: P(y>a)=e−λmax⁡{a,0}P(y>a)=e^{-\lambda\max\{a,0\}}.
  • (1) Dato x=2x=2, z=2+yz=2+y: P(z>A∣x=2)=P(y>A−2)=e−λmax⁡{A−2,0}P(z>A\mid x=2)=P(y>A-2)=e^{-\lambda\max\{A-2,0\}}; A=1→1A=1\to1, A=3→e−2=0,135A=3\to e^{-2}=0{,}135.
  • (2) P(z<A∣x=1)=1−e−λmax⁡{A−1,0}P(z<A\mid x=1)=1-e^{-\lambda\max\{A-1,0\}}; somma =0,6 e−λmax⁡{A−2,0}+0,4(1−e−λmax⁡{A−1,0})=0{,}6\,e^{-\lambda\max\{A-2,0\}}+0{,}4\big(1-e^{-\lambda\max\{A-1,0\}}\big); A=−2→0,6A=-2\to0{,}6.
  • Significato: probabilità di decisione corretta del ricevitore a soglia (decide 22 se z>Az>A). Per A≤1A\le1 è 0,60{,}6, cresce fino a 0,9460{,}946 in A=2A=2 e scende a 0,40{,}4.
  • Errori: dimenticare il max⁡{⋅,0}\max\{\cdot,0\}; usare λ\lambda come media; non condizionare o non pesare con pxp_x.

Teoria collegata