Bagaimana cara mendapatkan prediktor filter Kalman stasioner?

14

Dalam babnya tentang filter Kalman, buku DSP saya menyatakan, tampaknya tiba-tiba, bahwa filter Kalman stasioner untuk suatu sistem

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

memiliki prediktor

x^(t+1|t)=(SEBUAH-SEBUAHK¯C)x^(t|t-1)+SEBUAHK¯y(t)

dan kovarians vektor keadaan stasioner dan gain Kalman

ˉ 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

di mana Q dan R menunjukkan kovarian dari kebisingan input w dan kebisingan pengukuran v .

Saya tidak bisa melihat bagaimana mencapai ini dari prediktor varians minimum. Bisakah seseorang menjelaskannya kepada saya, atau mengarahkan saya ke sumber daya yang mendapatkan ungkapan? Ini adalah filter varian minimum varian waktu, yang dapat saya peroleh:

P(t+1|t)=A(P(t|t-1)-P(t|

x^(t+1|t)=(SEBUAH-K(t)C)x^(t|t-1)+K(t)y(t)
K ( t ) = A P ( t | t - 1 ) C T ( C P ( t | t -
P(t+1|t)=SEBUAH(P(t|t-1)-P(t|t-1)CT(CP(t|t-1)CT+R)-1CP(t|t-1))SEBUAHT+Q
K(t)=SEBUAHP(t|t-1)CT(CP(t|t-1)CT+R)-1

Saya hanya tidak yakin tentang cara pergi dari sini ke filter stasioner di atas.

Pembaruan: Saya dapat melihat bahwa mengganti dan ke dalam filter varian waktu menghasilkan hasil filter stasioner, tetapi mengapa kalikan dengan ? Apakah ini hanya gejala dari pilihan notasi yang tidak menguntungkan, yang berarti bahwa atau tidak benar-benar menunjukkan perolehan Kalman?K(t)=A ˉ K AK ˉ KP¯=P(t+1|t)=P(t|t-1)K(t)=SEBUAHK¯SEBUAHKK¯

Andreas
sumber
Tidak, tidak mungkin untuk "melihat" prediktor dari persamaan untuk sistem. Saya pikir akan lebih baik jika Anda membaca buku teks pada filter Kalman daripada meminta kami untuk menurunkannya untuk Anda (yang hanya akan memuntahkan sesuatu dari buku teks). Penyaringan optimal oleh Anderson dan Moore mungkin merupakan tempat yang baik untuk memulai. Itu diturunkan di bab 5, jika saya ingat dengan benar.
Lorem Ipsum
@yoda: Terima kasih. Pertanyaan saya adalah apakah seseorang dapat mengarahkan saya ke sumber yang lebih baik daripada buku teks yang saya rekomendasikan, jadi itu jawaban.
Andreas
@Yoda: Omong-omong, kalau-kalau saya tidak jelas: Saya tidak meminta derivasi dari sistem ruang-negara, tetapi dari varians minimum Kalman filter. Saya telah memperbarui pertanyaan untuk membuatnya lebih jelas bahwa saya dapat memperoleh filter Kalman invarian-waktu, hanya saja bukan yang stasioner.
Andreas
1
Dari teks mana Anda mendapatkan hal di atas? Jika seseorang memiliki akses ke sana, mungkin berguna agar kita dapat melihat konteks lengkapnya.
Jason R

Jawaban:

5

Derivasi Anda benar.

K(t)=A ˉ KP¯=P(t|t-1) danK(t)=SEBUAHK¯

Apakah ini kebingungan Anda:

  1. Mengapa mereka tidak memiliki istilah di Kalman Gain dan Covariance Matrix Expressions?t|t-1
  2. Bagaimana ini bisa menjadi "stasioner" ketika derivasi Anda menunjukkan bahwa waktunya bervariasi?

  1. Pilihan yang buruk dari notasi pada buku bagian

Penampilan Mari kita lihat ekspresi: . Fakta bahwa adalah fungsi dari dirinya sendiri menunjukkan hubungan rekursif. Dengan kata lain, ia menggunakan nilai masa lalu. Jadi, BUKAN sama untuk semua instance waktu - Ini berubah di setiap iterasi.P¯=SEBUAHP¯SEBUAHT-SEBUAHP¯CT(CP¯CT+R)-1CP¯SEBUAHT+QP¯

  1. Kesalahpahaman dari kata "stasioner".

Ketika penulis buku mengatakan "stasioner" dia tidak berarti bahwa dan memiliki nilai yang sama setiap saat. Sebagai gantinya, penulis ingin menekankan bahwa ekspresi bagi mereka untuk nilai adalah sama untuk semua realisasi statistik. Stationaritas adalah konsep statistik yang berarti bahwa Statistik sistem sama sepanjang waktu. Lihatlah ekspresi untuk dan lagi. Mereka hanya bergantung pada \PKP¯K¯

  • Nilai-nilai sebelumnya dari diri mereka sendiri
  • Matriks transisi dan yang bersifat deterministik dan dalam kasus Anda invarian waktu ( dan adalah sama setiap saat)SEBUAHCSEBUAHC
  • Q dan yang merupakan matriks kovarians derau. 2 matriks ini menggambarkan statistik dari suara-suara dan sama dalam semua realisasi dan contoh waktu.R

Gain Kalman, , dan kovariansi negara Matrix akan memiliki nilai yang sama untuk semua realisasi proses acak ini. ( Catatan Sisi: Tidak satu pun dari 2 istilah ini bergantung pada pengukuran, . Jadi mereka dapat dihitung sebelumnya. )KPy


Kesimpulan:

Persamaan "varian waktu" yang Anda peroleh sama dengan persamaan dalam buku. Selain itu, perbedaan notasi, ada sedikit kesalahpahaman di pihak Anda tentang apa yang berubah dan apa yang tidak.

ssk08
sumber
1
Saya tidak ingat apa masalah yang saya miliki ketika saya mengajukan pertanyaan, tetapi sekarang masuk akal. Terima kasih!
Andreas
Saya tidak mengerti ini. Seperti apa persamaan untuk filter Kalman yang tidak stasioner?
Sandu Ursu