Salta al contenuto
Note per Studenti Esercizio 23 · teste di due monete e approssimazione normale

Esercizio 23teste di due monete e approssimazione normale

In questa pagina 3

Testo (I appello, 15 gennaio 2026, esercizio 3, Tema A; punti 5 = 2+3). Supponiamo di disporre di due monete eque, indicate con A e B, e di lanciarle ripetutamente per 50 volte. Indichiamo con XAX_A il numero di teste ottenute dai lanci della moneta A, e con XBX_B il numero di teste ottenute dai lanci della moneta B.

(a) Determinare la distribuzione di XAX_A, XBX_B e di Y:=XA+XBY := X_A + X_B.

(b) Stimare la probabilità dell'evento {Y≤47}\{Y \le 47\} tramite approssimazione normale.


(a) Le distribuzioni

XAX_A e XBX_B

Ogni lancio della moneta A è una prova con successo "testa" di probabilità 12\frac12; i 50 lanci sono indipendenti. XAX_A conta i successi in 50 prove: per il modello binomialen prove indipendenti, ciascuna con probabilità di successo p: una sequenza con k successi ha probabilità p^k (1−p)^(n−k), e la probabilità di esattamente k successi è (n su k) p^k (1−p)^(n−k) (modello binomiale); il primo successo alla prova k ha probabilità (1−p)^(k−1) p.Prove ripetute e modello binomiale →,

XA∼Bin(50,12),P(XA=k)=(50k)1250,k=0,…,50.X_A \sim \text{Bin}\left(50, \tfrac12\right), \qquad P(X_A = k) = \binom{50}{k}\frac{1}{2^{50}}, \quad k = 0, \dots, 50.

(Con p=12p = \frac12 i due fattori pk(1−p)50−kp^k(1 - p)^{50 - k} si combinano in 1250\frac{1}{2^{50}}.) Lo stesso vale per XBX_B.

XAX_A e XBX_B sono indipendenti: riguardano lanci di monete diverse, che non si influenzano.

Y=XA+XBY = X_A + X_B

Somma di due binomiali indipendenti con lo stesso pp: è una binomiale con i numeri di prove sommati (Somma di variabili aleatorie indipendentiSe X e Y sono indipendenti, la legge di Z = X + Y è la convoluzione: p_Z(n) = Σ_k p_X(k) p_Y(n − k) nel discreto, f_Z(z) = ∫ f_X(z − y) f_Y(y) dy nel continuo. Casi notevoli: Bin(n,p) + Bin(m,p) = Bin(n+m,p), Poi(λ) + Poi(μ) = Poi(λ+μ), Geo + Geo con densità (n−1)p²(1−p)^(n−2), Exp(λ) + Exp(λ) = Γ(2,λ), gaussiane indipendenti sommano medie e varianze.Somma di variabili aleatorie indipendenti →). Il motivo, senza conti: YY conta le teste in 50+50=10050 + 50 = 100 lanci indipendenti di monete eque.

XA,XB∼Bin(50,12),Y∼Bin(100,12),P(Y=n)=(100n)12100, n=0,…,100\boxed{X_A, X_B \sim \text{Bin}\left(50, \tfrac12\right), \qquad Y \sim \text{Bin}\left(100, \tfrac12\right), \quad P(Y = n) = \binom{100}{n}\frac{1}{2^{100}}, \ n = 0, \dots, 100}

(b) P(Y≤47)P(Y \le 47) con l'approssimazione normale

Perché si può approssimare

Il calcolo esatto richiede ∑n=047(100n)12100\sum_{n=0}^{47}\binom{100}{n}\frac1{2^{100}}: 48 termini con coefficienti enormi. Ma YY è somma di 100 Bernoulli indipendenti, e per il teorema del limite centraleSe X₁, X₂, ... sono i.i.d. con media μ e varianza σ² ∈ (0, ∞), la somma standardizzata (Sₙ − nμ)/(σ√n) converge in distribuzione a N(0, 1): per n grande P(Sₙ ≤ x) ≈ Φ((x − nμ)/(σ√n)). Caso binomiale (De Moivre-Laplace): Bin(n, p) ≈ N(np, np(1 − p)), con correzione di continuità ±0,5. Si dimostra con le funzioni caratteristiche e il teorema di Lévy.Teorema del limite centrale e approssimazione normale → (teorema di De Moivre-Laplace) ha legge approssimativamente gaussiana con la stessa media e varianza. Le condizioni pratiche np=50≥5np = 50 \ge 5 e n(1−p)=50≥5n(1 - p) = 50 \ge 5 sono ampiamente soddisfatte.

Media e varianza

E[Y]=np=100⋅12=50,Var(Y)=np(1−p)=100⋅12⋅12=25,σY=5.E[Y] = np = 100\cdot\tfrac12 = 50, \qquad \text{Var}(Y) = np(1 - p) = 100\cdot\tfrac12\cdot\tfrac12 = 25, \qquad \sigma_Y = 5.

Standardizzazione

Con Z∼N(0,1)Z \sim N(0, 1):

P(Y≤47)=P(Y−505≤47−505)≈P(Z≤−0,6)=Φ(−0,6)=1−Φ(0,6).P(Y \le 47) = P\left(\frac{Y - 50}{5} \le \frac{47 - 50}{5}\right) \approx P(Z \le -0{,}6) = \Phi(-0{,}6) = 1 - \Phi(0{,}6).

L'ultimo passaggio usa la simmetria Φ(−z)=1−Φ(z)\Phi(-z) = 1 - \Phi(z) (Distribuzione gaussiana (normale)N(μ, σ²) ha densità e^(−(x−μ)²/(2σ²)) / √(2πσ²), a campana centrata in μ con larghezza σ; media μ, varianza σ²; si standardizza con Z = (X − μ)/σ ~ N(0, 1) e si calcola P(X ≤ x) = Φ((x − μ)/σ), con Φ(−z) = 1 − Φ(z); aX + b è ancora gaussiana, N(aμ + b, a²σ²).Distribuzione gaussiana (normale) →).

P(Y≤47)≈1−Φ(0,6)≈1−0,7257=0,274\boxed{P(Y \le 47) \approx 1 - \Phi(0{,}6) \approx 1 - 0{,}7257 = 0{,}274}

Con la correzione di continuità

YY è a valori interi: {Y≤47}={Y<48}\{Y \le 47\} = \{Y < 48\}. Senza correzione, la gaussiana darebbe a questi due eventi probabilità diverse (Φ(−0,6)\Phi(-0{,}6) contro Φ(−0,4)\Phi(-0{,}4)). La correzione di continuità taglia a metà strada, in 47,547{,}5 (ogni intero jj "occupa" l'intervallo [j−0,5,j+0,5][j - 0{,}5, j + 0{,}5]):

P(Y≤47)=P(Y≤47,5)≈P(Z≤47,5−505)=Φ(−0,5)=1−Φ(0,5).P(Y \le 47) = P(Y \le 47{,}5) \approx P\left(Z \le \frac{47{,}5 - 50}{5}\right) = \Phi(-0{,}5) = 1 - \Phi(0{,}5).

P(Y≤47)≈1−Φ(0,5)≈1−0,6915=0,309\boxed{P(Y \le 47) \approx 1 - \Phi(0{,}5) \approx 1 - 0{,}6915 = 0{,}309}

Confronto con il valore esatto

metodo valore
normale senza correzione 0,27430{,}2743
normale con correzione 0,30850{,}3085
binomiale esatta (calcolata al computer) 0,30860{,}3086

Con la correzione di continuità l'errore è di un decimillesimo; senza, di circa 0,030{,}03. La soluzione ufficiale accetta entrambe le risposte, ma la seconda è molto migliore: conviene usarla sempre con v.a. a valori interi.

Plausibilità. 4747 è sotto la media 5050 di poco più di mezza deviazione standard: la probabilità deve essere un po' meno di 12\frac12. Infatti la probabilità di "al più 50" sarebbe circa 12\frac12 (un po' di più, per il peso del valore 5050 stesso, P(Y=50)≈0,08P(Y = 50) \approx 0{,}08).

Errori comuni

  • Y∼Bin(50,1)Y \sim \text{Bin}(50, 1) o Bin(50,12)⋅2\text{Bin}(50, \frac12)\cdot 2: la somma di due binomiali indipendenti con lo stesso pp somma i numeri di prove, non le probabilità.
  • Standardizzare dividendo per la varianza 2525 invece che per σ=5\sigma = 5.
  • Correzione di continuità nel verso sbagliato: per P(Y≤47)P(Y \le 47) si usa 47,547{,}5, non 46,546{,}5.
  • Scrivere Φ(−0,6)=−Φ(0,6)\Phi(-0{,}6) = -\Phi(0{,}6).

Teoria collegata