Salta al contenuto
Note per Studenti Esercizio 22 · somma di un numero geometrico di Bernoulli

Esercizio 22somma di un numero geometrico di Bernoulli

In questa pagina 5

Testo (II parziale, 15 gennaio 2026, esercizio 1, Tema A; punti 8 = 2+3+3). Siano (Xk)k≥1(X_k)_{k \ge 1} v.a. indipendenti e identicamente distribuite, con distribuzione di Bernoulli di parametro p∈(0,1)p \in (0, 1), e sia T∼Geo(q)T \sim \text{Geo}(q), una variabile geometrica di parametro q∈(0,1)q \in (0, 1), indipendente dalle (Xk)k≥1(X_k)_{k \ge 1}. Poniamo Y:=∑k=1TXkY := \sum_{k=1}^T X_k.

(a) Determinare la densità discreta di YY condizionata all'evento {T=n}\{T = n\}, per n∈Nn \in \mathbb N. È una distribuzione nota?

(b) Calcolare il valor medio di YY.

(c) Calcolare la funzione caratteristica di YY, φY(θ):=E[eiθY]\varphi_Y(\theta) := E[e^{i\theta Y}].

Sugg.: usare la proprietà della torre e nello svolgimento porre a:=φBe(p)(θ)a := \varphi_{\text{Be}(p)}(\theta).


Che cosa rappresenta YY

YY è una somma con un numero aleatorio di addendi. Un esempio concreto: un trasmettitore invia pacchetti finché il canale non si chiude; il numero di pacchetti inviati è T∼Geo(q)T \sim \text{Geo}(q) (il canale si chiude con probabilità qq dopo ogni pacchetto); ogni pacchetto arriva corretto con probabilità pp, indipendentemente. YY = numero di pacchetti arrivati corretti.

La difficoltà: TT è aleatorio, quindi YY non è una binomiale. L'idea è condizionare a TT: fissato il numero di addendi, tutto diventa noto.

(a) Legge di YY dato T=nT = n

Sull'evento {T=n}\{T = n\} si ha Y=X1+⋯+XnY = X_1 + \dots + X_n. Quindi

P(Y=j∣T=n)=P(X1+⋯+Xn=j∣T=n)=P(X1+⋯+Xn=j).P(Y = j \mid T = n) = P(X_1 + \dots + X_n = j \mid T = n) = P(X_1 + \dots + X_n = j).

L'ultimo passaggio usa l'indipendenza di TT dalle XkX_k: sapere che T=nT = n non cambia la legge di X1,…,XnX_1, \dots, X_n, quindi si può togliere il condizionamento.

X1+⋯+XnX_1 + \dots + X_n è la somma di nn Bernoulli indipendenti di parametro pp, cioè il numero di successi in nn prove indipendenti: è una binomialeBernoulli Be(p): un solo tentativo, vale 1 con probabilità p; binomiale Bin(n, p): numero di successi in n prove indipendenti, P(X = k) = C(n,k) p^k (1−p)^(n−k), media np e varianza np(1−p); uniforme discreta: n valori equiprobabili.Bernoulli, binomiale e uniforme discreta → Bin(n,p)\text{Bin}(n, p).

P(Y=j∣T=n)=(nj)pj(1−p)n−j,j=0,1,…,n;Y∣T=n ∼ Bin(n,p)\boxed{P(Y = j \mid T = n) = \binom nj p^j(1 - p)^{n - j}, \quad j = 0, 1, \dots, n; \qquad Y \mid T = n \ \sim\ \text{Bin}(n, p)}

Sì, è una distribuzione nota: la binomiale (Densità e media condizionataLa densità condizionata di X dato Y = y è p(x, y)/p_Y(y) nel discreto e f(x, y)/f_Y(y) nel continuo; la sua media è E[X | Y = y], e E[X | Y] è la v.a. che si ottiene sostituendo y con Y. Proprietà della torre: E[E[X | Y]] = E[X]. Per un vettore gaussiano, X dato Y = y è gaussiana con media m₁ + (Cov/Var Y)(y − m₂) e varianza σ₁²(1 − ρ²).Densità e media condizionata →).

(b) Valor medio di YY

Proprietà della torre (Densità e media condizionataLa densità condizionata di X dato Y = y è p(x, y)/p_Y(y) nel discreto e f(x, y)/f_Y(y) nel continuo; la sua media è E[X | Y = y], e E[X | Y] è la v.a. che si ottiene sostituendo y con Y. Proprietà della torre: E[E[X | Y]] = E[X]. Per un vettore gaussiano, X dato Y = y è gaussiana con media m₁ + (Cov/Var Y)(y − m₂) e varianza σ₁²(1 − ρ²).Densità e media condizionata →): E[Y]=E[E[Y∣T]]E[Y] = E\big[E[Y \mid T]\big].

Passo 1. Dato T=nT = n, Y∼Bin(n,p)Y \sim \text{Bin}(n, p), che ha media npnp: E[Y∣T=n]=np⟹E[Y∣T]=pT.E[Y \mid T = n] = np \quad \Longrightarrow \quad E[Y \mid T] = pT.

Passo 2. Media di pTpT, con E[T]=1qE[T] = \frac1q (Distribuzione geometricaGeo(p) è il numero della prova in cui arriva il primo successo in prove indipendenti: P(X = n) = (1−p)^(n−1) p per n ≥ 1, P(X > n) = (1−p)^n (lunga attesa), media 1/p, varianza (1−p)/p², ed è senza memoria.Distribuzione geometrica →): E[Y]=E[pT]=p E[T]=pq.E[Y] = E[pT] = p\,E[T] = \frac pq.

E[Y]=pq\boxed{E[Y] = \frac pq}

Interpretazione: si fanno in media 1q\frac1q prove, ognuna con successo di probabilità pp: in media pq\frac pq successi. È la formula di WaldPer una somma di un numero aleatorio T di v.a. i.i.d. indipendenti da T: E[somma] = E[T] per E[X_1]. E[Y]=E[T] E[X1]E[Y] = E[T]\,E[X_1].

Attenzione: il risultato sembra ovvio, ma la linearità della media da sola non basta, perché il numero di addendi è aleatorio: serve il condizionamento.

(c) Funzione caratteristica di YY

Funzione caratteristica di Bernoulli e binomiale

Per X∼Be(p)X \sim \text{Be}(p) (Funzione caratteristicaLa funzione caratteristica φ_X(θ) = E[e^(iθX)] = E[cos θX] + iE[sin θX] esiste per ogni v.a. e ogni θ, vale 1 in 0, ha modulo ≤ 1, determina la legge (è la trasformata di Fourier della densità), trasforma le somme di v.a. indipendenti in prodotti e dà i momenti: E[X^k] = φ^(k)(0) / i^k. Per N(μ, σ²) vale e^(iμθ − σ²θ²/2).Funzione caratteristica →): a:=φBe(p)(θ)=E[eiθX]=(1−p)eiθ⋅0+p eiθ⋅1=peiθ+1−p.a := \varphi_{\text{Be}(p)}(\theta) = E[e^{i\theta X}] = (1 - p)e^{i\theta\cdot 0} + p\,e^{i\theta\cdot 1} = pe^{i\theta} + 1 - p.

Per Bin(n,p)\text{Bin}(n, p), somma di nn Bernoulli indipendenti, la funzione caratteristica è il prodotto: φBin(n,p)(θ)=an\varphi_{\text{Bin}(n,p)}(\theta) = a^n.

Torre

Con la versione generale della proprietà della torre, E[g(Y)]=E[E[g(Y)∣T]]E[g(Y)] = E\big[E[g(Y) \mid T]\big], e g(y)=eiθyg(y) = e^{i\theta y}:

φY(θ)=E[E[eiθY∣T]].\varphi_Y(\theta) = E\big[E[e^{i\theta Y} \mid T]\big].

Dato T=nT = n, Y∼Bin(n,p)Y \sim \text{Bin}(n, p), quindi E[eiθY∣T=n]=anE[e^{i\theta Y} \mid T = n] = a^n, cioè E[eiθY∣T]=aTE[e^{i\theta Y} \mid T] = a^T. Resta da calcolare E[aT]E[a^T]:

φY(θ)=E[aT]=∑k=1∞ak P(T=k)=∑k=1∞ak(1−q)k−1q.\varphi_Y(\theta) = E[a^T] = \sum_{k=1}^\infty a^k\,P(T = k) = \sum_{k=1}^\infty a^k(1 - q)^{k - 1}q.

La serie geometrica

Si raccoglie qaqa e si pone j=k−1j = k - 1:

∑k=1∞ak(1−q)k−1q=qa∑k=1∞((1−q)a)k−1=qa∑j=0∞((1−q)a)j=qa1−(1−q)a.\sum_{k=1}^\infty a^k(1 - q)^{k-1}q = qa\sum_{k=1}^\infty \big((1 - q)a\big)^{k-1} = qa\sum_{j=0}^\infty\big((1 - q)a\big)^j = \frac{qa}{1 - (1 - q)a}.

Convergenza. La serie geometrica (anche complessa) converge se la ragione ha modulo <1< 1. Qui ∣a∣=∣φBe(p)(θ)∣≤1|a| = |\varphi_{\text{Be}(p)}(\theta)| \le 1 (ogni funzione caratteristica ha modulo al più 11) e 1−q<11 - q < 1, quindi ∣(1−q)a∣≤1−q<1|(1 - q)a| \le 1 - q < 1 ✓.

φY(θ)=qa1−(1−q)a,a=peiθ+1−p\boxed{\varphi_Y(\theta) = \frac{qa}{1 - (1 - q)a}, \qquad a = pe^{i\theta} + 1 - p}

È la funzione E[sT]=qs1−(1−q)sE[s^T] = \frac{qs}{1 - (1 - q)s} della Distribuzione geometricaGeo(p) è il numero della prova in cui arriva il primo successo in prove indipendenti: P(X = n) = (1−p)^(n−1) p per n ≥ 1, P(X > n) = (1−p)^n (lunga attesa), media 1/p, varianza (1−p)/p², ed è senza memoria.Distribuzione geometrica → calcolata in s=as = a.

Controlli.

  • θ=0\theta = 0: a=1a = 1, φY(0)=q1−(1−q)=1\varphi_Y(0) = \frac{q}{1 - (1 - q)} = 1 ✓.
  • Media: E[Y]=−iφY′(0)E[Y] = -i\varphi_Y'(0). Derivando φY=qa1−(1−q)a\varphi_Y = \frac{qa}{1 - (1 - q)a} rispetto ad aa si ottiene q(1−(1−q)a)2\frac{q}{(1 - (1 - q)a)^2}, che in a=1a = 1 vale qq2=1q\frac{q}{q^2} = \frac1q; e dadθ=ipeiθ\frac{da}{d\theta} = ipe^{i\theta}, che in 00 vale ipip. Quindi φY′(0)=ipq\varphi_Y'(0) = \frac{ip}{q} e E[Y]=−i⋅ipq=pqE[Y] = -i\cdot\frac{ip}{q} = \frac pq ✓, come in (b).

Un'osservazione in più (non richiesta). Posto r=p+q−pqr = p + q - pq, sostituendo aa e semplificando si verifica che φY(θ)=π+(1−π) βeiθ1−(1−β)eiθ,π=q(1−p)r,β=qr.\varphi_Y(\theta) = \pi + (1 - \pi)\,\frac{\beta e^{i\theta}}{1 - (1 - \beta)e^{i\theta}}, \qquad \pi = \frac{q(1 - p)}{r}, \quad \beta = \frac qr. Per l'unicità della funzione caratteristica: YY vale 00 con probabilità π\pi (nessun successo prima che il canale si chiuda) e, altrimenti, è una Geo(β)\text{Geo}(\beta). Riconoscere la forma della funzione caratteristica basta a identificare la legge.

Errori comuni

  • Scrivere Y∼Bin(T,p)Y \sim \text{Bin}(T, p) come se fosse una legge: lo è solo condizionatamente a TT; la legge di YY è un'altra.
  • Calcolare E[Y]=E[T]⋅pE[Y] = E[T]\cdot p "per linearità" senza giustificazione: il risultato è giusto, ma va motivato con la torre (o con Wald).
  • Partire la serie da k=0k = 0: la geometrica del corso parte da 11.
  • Dimenticare di verificare che la serie converge.

Teoria collegata