Domanda sulla ponderazione della varianza inversa


9

Supponiamo di voler fare l'inferenza su una realizzazione non osservata di una variabile casuale ~ x , che è distribuita normalmente con media μ x e varianza σ 2 x . Supponiamo che ci sia un'altra variabile casuale ˜ y (la cui realizzazione non osservata chiameremo allo stesso modo y ) che è normalmente distribuita con media μ y e varianza σ 2 y . Lasciate σ x y sia la covarianza di ~ x e ~ y .xx~μxσx2y~yμyσy2σxyx~y~

Supponiamo ora di osservare un segnale su , a = x + ˜ u , dove ˜ uN ( 0 , ϕ 2 x ) e un segnale su y , b = y + ˜ v , dove ˜ vN ( 0 , ϕ 2 anni ) . Supponiamo che ˜ u e ˜ v siano indipendenti.x

a=x+u~,
u~N(0,ϕx2)y
b=y+v~,
v~N(0,ϕy2)u~v~

Qual è la distribuzione di alla condizione che un e B ?xab

Quello che so finora: utilizzando la ponderazione della varianza inversa, e Var(x

E(x|a)=1σx2μx+1ϕx2a1σx2+1ϕx2,
Var(x|a)=11σx2+1ϕx2.

Dal momento che ed y sono disegnati congiuntamente, b dovrebbe portare alcune informazioni su x . Oltre a rendermene conto, sono bloccato. Qualsiasi aiuto è apprezzato!xybx


Questo è esattamente come i primi passi sulla derivazione di un filtro Kalman. Potresti guardare la derivazione e pensare al guadagno di Kalman per l'aggiornamento della stima della covarianza statale. cs.unc.edu/~welch/media/pdf/kalman_intro.pdf
EngrStudent

Grazie per la risposta! Ho letto il documento nel tuo link, ma non vedo la connessione con il filtro Kalman. Qualche possibilità che potresti elaborare? Apprezzo l'aiuto!
bad_at_math il

2
@EngrStudent Se l'OP non ha familiarità con il filtro Kalman, non vedo come sarà di grande aiuto. Forse potresti invece spiegare come affrontare il problema senza invocare nessuno dei dettagli (o il gergo) coinvolti con il KF, anche se forse facendo uso della tua comprensione di esso per guidare una risposta sui dettagli qui.
Glen_b

Trascritto su math.SE qui
Glen_b -Reinstate Monica il

Risposte:


2

xabxyab

[xyuv]N([μxμy00],[σx2σxy00σxyσy20000ϕx20000ϕy2])
a=x+ub=y+v
[xab]N([μxμxμy],[σx2σx2σxyσx2σx2+ϕx2σxyσxyσxyσy2+ϕy2]).
uvxy

xab

Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.