Pregunta sobre ponderación de varianza inversa


9

Supongamos que queremos hacer inferencia en una realización no observada de una variable aleatoria , que normalmente se distribuye con media y varianza . Supongamos que hay otra variable aleatoria (cuya realización no observada llamaremos de manera similar ) que normalmente se distribuye con media y varianza . Sea la covarianza de y .xx~μxσx2y~yμyσy2σxyx~y~

Ahora supongamos que observamos una señal en , where y una señal en , donde . Suponga que y son independientes.x

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

¿Cuál es la distribución de condicional en a y b ?xab

Lo que sé hasta ahora: utilizando ponderación de varianza inversa,

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

Dado que x e y se dibujan conjuntamente, b debería llevar cierta información sobre x . Aparte de darme cuenta de esto, estoy atascado. Cualquier ayuda es apreciada!


Esto se ve exactamente como los primeros pasos en la derivación de un filtro de Kalman. Puede mirar la derivación y pensar en la ganancia de Kalman para la actualización de estimación de covarianza del estado. cs.unc.edu/~welch/media/pdf/kalman_intro.pdf
EngrStudent

¡Gracias por la respuesta! Leí el documento en su enlace, pero no veo la conexión con el filtrado de Kalman. ¿Alguna posibilidad de que puedas elaborar? Agradezco la ayuda!
bad_at_math

2
@EngrStudent Si el OP no está familiarizado con el filtro de Kalman, no veo cómo será de mucha ayuda. Tal vez podría explicar cómo abordar el problema sin invocar ninguno de los detalles (o la jerga) involucrados con el KF, aunque quizás haciendo uso de su comprensión para guiar una respuesta sobre los detalles aquí.
Glen_b -Reinstale a Monica el

Publicación cruzada en math.SE here
Glen_b -Reinstate Monica

Respuestas:


2

No estoy seguro de si las fórmulas de ponderación de varianza inversa se aplican aquí. Sin embargo, creo que podría calcular la distribución condicional de dada y suponiendo que , , y siguen una distribución normal multivariada conjunta.xabxyab

En concreto, si se asume (de forma compatible con lo que se especifica en la pregunta) que entonces, dejando y , puedes encontrar que

[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]).
(Tenga en cuenta que en lo anterior se supone implícitamente que y son independientes entre sí y también con e .)uvxy

A partir de esto, puede encontrar la distribución condicional de dada y usando las propiedades estándar de la distribución normal multivariada (vea aquí, por ejemplo: http://en.wikipedia.org/wiki/Multivariate_normal_distribution#Conditional_distributions ).xab

Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.