Salta al contenuto
Note per Studenti Esercizio 20 · vettore gaussiano, combinazione lineare e correlazione

Esercizio 20vettore gaussiano, combinazione lineare e correlazione

In questa pagina 5

Testo (III appello, 24 giugno 2026, esercizio 5, Tema A; punti 6 = 2+2+2). Sia (X,Y)(X, Y) un vettore aleatorio Gaussiano con media m=(1,−1)m = (1, -1) e matrice di covarianza

Σ=(4222).\Sigma = \begin{pmatrix} 4 & 2 \\ 2 & 2 \end{pmatrix}.

(a) Scrivere la densità congiunta di (X,Y)(X, Y).

(b) Determinare la legge di U:=X−2YU := X - 2Y e calcolare P(U>5)P(U > 5).

(c) Calcolare Cov(X,U)\text{Cov}(X, U) e il coefficiente di correlazione ρ(X,U)\rho(X, U).


Lettura dei dati

  • E[X]=1E[X] = 1, E[Y]=−1E[Y] = -1;
  • Var(X)=4\text{Var}(X) = 4, Var(Y)=2\text{Var}(Y) = 2 (diagonale);
  • Cov(X,Y)=2\text{Cov}(X, Y) = 2 (fuori diagonale). Qui XX e YY non sono indipendenti: ρ(X,Y)=24⋅2=12\rho(X, Y) = \frac{2}{\sqrt{4\cdot 2}} = \frac{1}{\sqrt2}.

(a) Densità congiunta

Formula del vettore gaussianoX = (X₁, ..., Xₙ) è un vettore gaussiano N(m, Σ) se ogni combinazione lineare a·X è gaussiana (equivalentemente X = m + AZ con Z gaussiane standard indipendenti); se Σ è invertibile ha densità exp(−½(x−m)ᵀΣ⁻¹(x−m)) / √((2π)ⁿ det Σ). Proprietà chiave: AX + b ~ N(Am + b, AΣAᵀ), le marginali sono gaussiane e componenti non correlate sono indipendenti.Vettori gaussiani → in dimensione 22:

f(x,y)=12πdet⁡Σexp⁡(−12(x−m1, y−m2) Σ−1(x−m1y−m2)).f(x, y) = \frac{1}{2\pi\sqrt{\det\Sigma}}\exp\left(-\frac12(x - m_1,\ y - m_2)\,\Sigma^{-1}\begin{pmatrix} x - m_1 \\ y - m_2 \end{pmatrix}\right).

Passo 1: determinante. det⁡Σ=4⋅2−2⋅2=4>0\det\Sigma = 4\cdot 2 - 2\cdot 2 = 4 > 0: la matrice è invertibile e la densità esiste.

Passo 2: inversa. Per una 2×22 \times 2 (abbd)\begin{pmatrix} a & b \\ b & d \end{pmatrix}: si scambiano gli elementi diagonali, si cambia segno a quelli fuori diagonale e si divide per il determinante:

Σ−1=14(2−2−24)=(12−12−121).\Sigma^{-1} = \frac14\begin{pmatrix} 2 & -2 \\ -2 & 4 \end{pmatrix} = \begin{pmatrix} \frac12 & -\frac12 \\ -\frac12 & 1 \end{pmatrix}.

Controllo: ΣΣ−1=(4⋅12−2⋅12−4⋅12+22⋅12−2⋅12−2⋅12+2)=(1001)\Sigma\Sigma^{-1} = \begin{pmatrix} 4\cdot\frac12 - 2\cdot\frac12 & -4\cdot\frac12 + 2 \\ 2\cdot\frac12 - 2\cdot\frac12 & -2\cdot\frac12 + 2 \end{pmatrix} = \begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix} ✓.

Passo 3: forma quadratica. Con u=x−1u = x - 1 e v=y+1v = y + 1 (scarti dalle medie):

(u,v)(12−12−121)(uv)=12u2−12uv−12uv+v2=12u2−uv+v2.(u, v)\begin{pmatrix} \frac12 & -\frac12 \\ -\frac12 & 1 \end{pmatrix}\begin{pmatrix} u \\ v \end{pmatrix} = \frac12 u^2 - \frac12 uv - \frac12 uv + v^2 = \frac12 u^2 - uv + v^2.

(I due termini fuori diagonale danno ciascuno −12uv-\frac12 uv.)

Passo 4: costante. 2πdet⁡Σ=2π⋅2=4π2\pi\sqrt{\det\Sigma} = 2\pi\cdot 2 = 4\pi.

f(x,y)=14πexp⁡(−12(12u2−uv+v2)),u=x−1, v=y+1\boxed{f(x, y) = \frac{1}{4\pi}\exp\left(-\frac12\left(\frac12 u^2 - uv + v^2\right)\right), \qquad u = x - 1,\ v = y + 1}

Per esteso: f(x,y)=14πexp⁡(−(x−1)24+(x−1)(y+1)2−(y+1)22)f(x, y) = \frac{1}{4\pi}\exp\left(-\frac{(x - 1)^2}{4} + \frac{(x - 1)(y + 1)}{2} - \frac{(y + 1)^2}{2}\right).

Il termine misto uvuv c'è perché la covarianza è diversa da zero: la densità non si fattorizza, coerentemente con la non indipendenza.

(b) Legge di U=X−2YU = X - 2Y e P(U>5)P(U > 5)

La legge

Ogni combinazione lineare delle componenti di un vettore gaussiano è gaussiana (è la definizione di vettore gaussianoX = (X₁, ..., Xₙ) è un vettore gaussiano N(m, Σ) se ogni combinazione lineare a·X è gaussiana (equivalentemente X = m + AZ con Z gaussiane standard indipendenti); se Σ è invertibile ha densità exp(−½(x−m)ᵀΣ⁻¹(x−m)) / √((2π)ⁿ det Σ). Proprietà chiave: AX + b ~ N(Am + b, AΣAᵀ), le marginali sono gaussiane e componenti non correlate sono indipendenti.Vettori gaussiani →). Quindi UU è gaussiana e bastano media e varianza.

Media (linearità): E[U]=E[X]−2E[Y]=1−2⋅(−1)=3.E[U] = E[X] - 2E[Y] = 1 - 2\cdot(-1) = 3.

Varianza (varianza di una combinazione lineareCov(X, Y) = E[(X − E X)(Y − E Y)] = E[XY] − E[X]E[Y] misura quanto X e Y variano insieme; è bilineare, Cov(X, X) = Var(X), Var(X + Y) = Var X + Var Y + 2Cov(X, Y); ρ = Cov / (σ_X σ_Y) sta in [−1, 1] e vale ±1 solo per legami lineari. Indipendenti ⇒ non correlate, ma non viceversa (tranne per i vettori gaussiani).Covarianza e coefficiente di correlazione →, con a=1a = 1, b=−2b = -2): Var(U)=a2Var(X)+b2Var(Y)+2ab Cov(X,Y)=4+4⋅2+2⋅1⋅(−2)⋅2=4+8−8=4.\text{Var}(U) = a^2\text{Var}(X) + b^2\text{Var}(Y) + 2ab\,\text{Cov}(X, Y) = 4 + 4\cdot 2 + 2\cdot 1\cdot(-2)\cdot 2 = 4 + 8 - 8 = 4.

U∼N(3,4)\boxed{U \sim N(3, 4)}

Il termine con la covarianza è negativo (perché aa e bb hanno segni opposti e la covarianza è positiva) e compensa in parte le varianze: XX e YY tendono a muoversi insieme, e nella differenza X−2YX - 2Y le loro oscillazioni si cancellano parzialmente. Senza covarianza si sarebbe trovato 1212.

La probabilità

σU=4=2\sigma_U = \sqrt4 = 2. Standardizzando con Z=U−32∼N(0,1)Z = \frac{U - 3}{2} \sim N(0, 1) (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(U>5)=P(Z>5−32)=P(Z>1)=1−Φ(1).P(U > 5) = P\left(Z > \frac{5 - 3}{2}\right) = P(Z > 1) = 1 - \Phi(1).

P(U>5)=1−Φ(1)≈0,159\boxed{P(U > 5) = 1 - \Phi(1) \approx 0{,}159}

(c) Covarianza e correlazione tra XX e UU

Covarianza, per bilinearità:

Cov(X,U)=Cov(X,X−2Y)=Cov(X,X)−2 Cov(X,Y)=Var(X)−2 Cov(X,Y)=4−2⋅2=0.\text{Cov}(X, U) = \text{Cov}(X, X - 2Y) = \text{Cov}(X, X) - 2\,\text{Cov}(X, Y) = \text{Var}(X) - 2\,\text{Cov}(X, Y) = 4 - 2\cdot 2 = 0.

Coefficiente di correlazione, con Var(X)=4\text{Var}(X) = 4 e Var(U)=4\text{Var}(U) = 4:

ρ(X,U)=Cov(X,U)Var(X) Var(U)=016=0.\rho(X, U) = \frac{\text{Cov}(X, U)}{\sqrt{\text{Var}(X)\,\text{Var}(U)}} = \frac{0}{\sqrt{16}} = 0.

Cov(X,U)=0,ρ(X,U)=0\boxed{\text{Cov}(X, U) = 0, \qquad \rho(X, U) = 0}

XX e UU sono non correlate. Di più: (X,U)(X, U) è un vettore gaussiano (trasformazione lineare di (X,Y)(X, Y)), quindi XX e UU sono anche indipendenti (Vettori gaussianiX = (X₁, ..., Xₙ) è un vettore gaussiano N(m, Σ) se ogni combinazione lineare a·X è gaussiana (equivalentemente X = m + AZ con Z gaussiane standard indipendenti); se Σ è invertibile ha densità exp(−½(x−m)ᵀΣ⁻¹(x−m)) / √((2π)ⁿ det Σ). Proprietà chiave: AX + b ~ N(Am + b, AΣAᵀ), le marginali sono gaussiane e componenti non correlate sono indipendenti.Vettori gaussiani →).

Perché proprio −2-2. Il coefficiente di YY in UU non è casuale: U=X−cYU = X - cY è non correlata con XX quando Var(X)−c Cov(X,Y)=0\text{Var}(X) - c\,\text{Cov}(X, Y) = 0, cioè c=Var(X)Cov(X,Y)=42=2c = \frac{\text{Var}(X)}{\text{Cov}(X, Y)} = \frac42 = 2. È lo stesso trucco con cui si dimostra la formula del vettore gaussiano condizionatoLa 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 →: si scompone una v.a. in una parte "spiegata" dall'altra e una parte indipendente da essa.

Errori comuni

  • Inversa 2×22 \times 2 sbagliata: dimenticare di cambiare segno agli elementi fuori diagonale, o di dividere per il determinante.
  • Forma quadratica con −12uv-\frac12 uv invece di −uv-uv: i termini fuori diagonale sono due.
  • Dimenticare 2ab Cov(X,Y)2ab\,\text{Cov}(X, Y) nella varianza di UU (si otterrebbe 1212 invece di 44).
  • Standardizzare dividendo per 44: si divide per σU=2\sigma_U = 2.

Teoria collegata