Sabit Kalman filtre tahmincisi nasıl türetilir?


14

Kalman filtreleri ile ilgili bölümünde DSP kitabım, görünüşte maviden çıktı, bir sistem için sabit Kalman filtresinin

{x(t+1)=Ax(t)+w(t)y(t)=Cx(t)+v(t)

öngörücüye sahip

x^(t+1|t)=(AAK¯C)x^(t|t1)+AK¯y(t)

ve durağan durum vektör kovaryansı ve Kalman kazancı

ˉ K = ˉ P CT(C ˉ P CT+R)-1

P¯=AP¯ATAP¯CT(CP¯CT+R)1CP¯AT+Q
K¯=P¯CT(CP¯CT+R)1

burada ve R sırasıyla giriş gürültüsünün w ve ölçüm gürültüsünün v kovaryanslarını belirtir .QRwv

Minimum varyans tahmincisinden buna nasıl ulaşacağımı göremiyorum. Birisi bana açıklayabilir mi yoksa beni ifadeyi türeyen bir kaynağa yönlendirebilir mi? Bu, zaman-varyantı en az varyans filtresi, olabilir türetmek:

x^(t+1|t)=(AK(t)C)x^(t|t1)+K(t)y(t)
P(t+1|t)=A(P(t|t1)P(t|t1)CT(CP(t|t1)CT+R)1CP(t|t1))AT+Q
K(t)=AP(t|t1)CT(CP(t|t1)CT+R)1

Buradan yukarıdaki sabit filtreye nasıl gidileceğinden emin değilim.

P¯=P(t+1|t)=P(t|t1)K(t)=AK¯AKK¯


Hayır, öngörücüyü sistem denklemlerinden "görmek" mümkün değildir. Kalman filtreleriyle ilgili bir ders kitabını sizin için türetmemizi istemek yerine daha iyi olacağını düşünüyorum. Anderson ve Moore tarafından Optimum Filtreleme başlamak için iyi bir yer olabilir. Eğer doğru hatırlıyorsam, bölüm 5'te elde edilmiştir.
Lorem Ipsum

@yoda: Teşekkürler. Benim sorum birisinin beni dersimin önerdiği ders kitabından daha iyi bir kaynağa yönlendirip yönlendiremeyeceğiydi, bu yüzden bu bir cevap.
Andreas

@yoda: Bu arada, belirsiz olduğumda: Durum-uzay sisteminden değil, minimum varyans Kalman filtresinden türetmek istiyorum. Soruyu, zamanla değişmeyen bir Kalman filtresi türetebileceğimi daha net hale getirmek için güncelledim, sabit olanı değil.
Andreas

1
Yukarıdakilerden hangi metni alıyorsunuz? Birisi ona erişebiliyorsa, tüm içeriği görebilmemiz yararlı olabilir.
Jason R

Yanıtlar:


5

Türevleriniz doğru.

P¯=P(t|t1)K(t)=AK¯

Bu senin karışıklığın mı?

  1. t|t1
  2. Derivasyonunuz zamanın değiştiğini gösterdiğinde bu nasıl "durağan" olabilir?

  1. Kitabın bölümünde kötü gösterim seçimi

P¯=AP¯ATAP¯CT(CP¯CT+R)1CP¯AT+QP¯

  1. "Sabit" kelimesinin yanlış anlaşılması.

PKP¯K¯

  • Önceki kendilerinin değerleri
  • ACAC
  • QR

KPy


Sonuç:

Türettiğiniz "zaman değişkeni" denklemler kitaptaki denklemlere eşitti. Ayrıca, notasyonel farklılıklar, neyin değişip neyin değişmediği konusunda küçük bir yanlış anlama vardı.


1
Soruyu sorduğumda yaşadığım sorunun ne olduğunu hatırlamıyorum, ama şimdi mantıklı. Teşekkürler!
Andreas

Bunu tam olarak anlamıyorum. O zaman sabit olmayan bir Kalman filtresi için denklemler nasıl olurdu?
Sandu Ursu
Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.