Representasi ruang-negara yang berbeda untuk Regresi Otomatis dan filter Kalman


8

Saya melihat bahwa ada berbagai cara untuk menulis model AR ke dalam representasi state-space, sehingga kita dapat menerapkan filter Kalman untuk memperkirakan sinyal. Lihat Contoh 1, 2 dan 3 di sini .

Saya bertanya-tanya apa perbedaan antara representasi ruang-negara yang berbeda pada estimasi oleh filter Kalman?

Terima kasih!


Ini adalah tempat yang tepat untuk itu, bukan Ilmu Komputasi . Jika Anda belum mendapatkan jawaban, coba perbarui posting yang menunjukkan upaya Anda selama seminggu terakhir - sudahkah Anda mencoba meneliti sendiri? Pilihan lain adalah menambahkan hadiah ...
Lorem Ipsum

Diskusi di sana tampaknya lebih teoretis daripada di sini. Filter Kalman adalah metode estimasi optimal untuk sistem dinamis stokastik. Jadi itu sangat cocok dengan ilmu komputasi. Saya belum menemukan sesuatu yang bermanfaat.
Tim

Sudahkah Anda mencoba menempatkan hadiah? Anda hanya perlu mendapatkan lebih banyak perhatian pada pertanyaan Anda, dan ada beberapa cara untuk melakukannya ...
Lorem Ipsum

Jawaban:


8

Sayangnya saya tidak tahu banyak tentang filter Kalman, tapi saya pikir saya bisa membantu Anda dengan hal-hal luar angkasa.

Pada Contoh 1, model AR adalah definisi output rekursif DSP lama yang baik:

yt=α+ϕ1yt1+ϕ2yt2+ηt

Dalam hal ini kita menuliskan model state-space dengan korespondensi langsung dengan persamaan di atas:

(ytyt1)=(ϕ1ϕ210)(yt1yt2)+(α0)+(10)ηt

Perhatikan bahwa dalam kasus ini, status sistem adalah nilai keluaran saat ini dan sebelumnya.

Dalam contoh kedua, Anda memisahkan negara Anda cdari nilai output Anda. Ini berarti bahwa status sekarang dapat berupa apa saja, meskipun mereka masih langsung memetakan ke nilai output. Dengan cara ini kita dapatkan

yt=μ+ct

ct=ϕ1ct1+ϕ2ct2+ηt

Dan oleh karena itu

(ctct1)=(ϕ1ϕ210)(ct1ct2)+(10)ηt

Anda juga harus mengenali ini sebagai representasi ruang-negara standar dari sistem linier, karena Anda persamaan untuk evolusi negara dan output yang bergantung pada negara adalah dua persamaan yang berbeda . Pemisahan ini sepele dalam hal model AR, tetapi notasi terakhir ini adalah bagaimana kita memikirkan semua model ruang-ruang linear secara umum.

Contoh ketiga adalah yang aneh. Jika Anda mengalikan semua koefisien Anda akan menyadari bahwa itu sebenarnya setara dengan contoh pertama dan kedua. Jadi mengapa melakukannya? Saya ternyata contoh 2 (menjadi representasi ruang-negara yang tepat dari sistem) disebut Bentuk Canonical yang Dapat Dikontrol dari sistem ini. Jika Anda membaca atau menganalisis sistem dengan hati-hati, Anda akan menyadari bahwa kami dapat menempatkan sistem ini ke keadaan apa pun yang kami sukai dengan memberikan nilai-nilai yang berlakuϕ1 dan ϕ2 dengan input tunggal α. Karena itu kami menyebut sistem seperti itu dapat dikendalikan , dan sangat mudah untuk melihat dari bentuk persamaan state-space ini.

Anda harus memperhatikan bahwa dua sistem linier dapat identik hingga perubahan basis. Ini berarti bahwa kita dapat memilih basis yang berbeda untuk mewakili sistem linier yang sama. Anda dapat meyakinkan diri sendiri bahwa itulah yang kami lakukan untuk beralih dari contoh kedua ke ketiga. Khususnya, kami menyukai transformasi linier ini untuk memindahkan matriks transisi keadaan, sehingga kami akan mendapatkan keadaan yang tidak diketahuis

yt=(10)αt

αt=(stst1)=(ϕ1ϕ210)(st1st2)+(α0)+(10)ηt

Sekarang kita bisa menggunakan perubahan basis untuk mencari tahu seperti apa keadaan ini s harus sehubungan dengan negara y. Dan kita bisa menghitungnya

(stst1)=(ytϕ2yt1)

Bentuk ini (transpos dari Controllability Canonical Form) disebut sebagai Observability Canonical Form karena jika kita dapat menempatkan sistem dalam formulir ini, kita dapat dengan mudah menyimpulkan keadaan sistem yang dapat diamati dengan hanya melihat output. Untuk beberapa deskripsi bentuk kanonik, Anda dapat membaca dokumen ini , dan tentu saja melihat-lihat di web. Perhatikan bahwa dalam dokumen keadaan terbalik, yang tidak mengubah apa pun tentang representasi sistem, hanya menyusun ulang baris / kolom matriks.


2

Singkatnya, itu semua tergantung pada apa yang Anda coba perkirakan, yaitu apa yang Anda ketahui tentang sinyal dan apa yang tidak. Filter Kalman akan mencoba memperkirakan negara berdasarkan definisi Anda tentang apa itu negara. Masalah konvensionalnya adalah ketika kita mencoba memperkirakan koefisien AR.

Mari kita ambil contoh AR(2) model tanpa suku konstan μ.

yk=a1yk1+a2yk2+ηk

Untuk memperkirakan sistem di atas, yang perlu Anda lakukan adalah memperkirakan koefisien AR, a1 dan a2.

Pengaturan Ruang Umum Filter Negara Kalman:

xk=Fk1xk1+wk
yk=Hkxk+vk
wk=WGN(0,Qs) dan vk=WGN(0,Qo)

Dalam hal ini, kita perlu memperkirakan a1 dan a2. Jadi wajar untuk mengatur negara sebagai koefisien ini.xk=[a1,a2]T Untuk contoh ini, koefisien ini konstan (Fk=Fk1=I) dan juga tidak ada noise dalam koefisien ini -> wk=0Qs=0.

Karena semua yang kita amati adalah yk, mereka menjadi ukuran untuk sistem kami. Karena kita telah mendefinisikan apa vektor keadaan, untuk persamaan pengukuran kita harus sama dengan model AR yang diberikan, kita mengganti noise pengukuran kitavk dengan ηk dan Hk=[yk1,yk2].

xk=xk1=[a1a2]
yk=Hkxk+ηk=[yk1yk2][a1a2]+ηk

Sekarang, Anda dapat menggunakan filter Kalman untuk memperkirakan kondisi Anda dan akibatnya sinyal Anda.

Catatan: Satu-satunya hal yang aneh di sini adalah matriks Anda Hk tergantung pada pengukuran Anda yk. Beberapa orang memiliki kesalahpahaman bahwa Keuntungan Kalman dan Negara Covariance Matriks selalu pengukuran independen dan bahwa mereka dapat dihitung sebelumnya. Kasus ini jelas menunjukkan bahwa ini bukan masalahnya. Baik Kalman Gain dan State Covariance Matrix diperkirakan dengan fungsi Hk, yang dalam hal ini tergantung pada pengukuran.


Saya tidak setuju. Saya pikir Anda kompromi pengamatan negara dengan memasukkan pengukuran dalam matriks
Dengan menggunakan situs kami, Anda mengakui telah membaca dan memahami Kebijakan Cookie dan Kebijakan Privasi kami.
Licensed under cc by-sa 3.0 with attribution required.