質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

新規登録して質問してみよう
ただいま回答率
85.50%
アルゴリズム

アルゴリズムとは、定められた目的を達成するために、プログラムの理論的な動作を定義するものです。

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

Q&A

1回答

2310閲覧

カルマンフィルターのシステムノイズと観測ノイズの同時推定について

shimotani1028

総合スコア5

アルゴリズム

アルゴリズムとは、定められた目的を達成するために、プログラムの理論的な動作を定義するものです。

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

0グッド

0クリップ

投稿2021/09/11 01:25

今カルマンフィルターを使おうとしており、ノイズを状態変数と同時に推定する方法を勉強しています。
最尤推定法で推定する方法はわかったのですが、リアルタイムで同時に推定したいと考えており、以下のURLで自己組織化型状態空間モデルの存在を知りました。これをカルマンフィルターに適用しようとしているのですが、計算方法がわかりません。
https://www.jstage.jst.go.jp/article/jjsai/16/2/16_300/_pdf
https://staff.aist.go.jp/naoyuki.ichimura/research/nng_ssm/papers/sossm_cvim01july.pdf
添付の図のように各種条件を設定したのですが、このままだと手で計算してもxtに設定した2つのノイズ(システムノイズの分散と観測ノイズの分散)成分を求めることができません。具体的には更新しても当たり前ですが、値が変わらないです。GやFの設定が間違っている?、予測フィルタステップ以外にノイズを更新する手続きが必要?な気もしているのですが、よくわかっておりません。
もし具体的な手順をご存知であれば教えていただけないでしょうか。
よろしくお願いいたします。
イメージ説明

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

ozwk

2021/09/13 04:00 編集

フィルタ: 側で RとQに添え字tがついていないですけど時間変化すべきものではないでしょうか。
shimotani1028

2021/09/13 12:56

そうですね。記載漏れです。申し訳ございません。
guest

回答1

0

あまり深く検証してないんですがとりあえず、

カルマンフィルタでxtのノイズの分散成分が変化しないのは
H=(1 0 0)ですのでカルマンゲインのノイズの分散成分が0になるからですね。

「自己組織化型状態空間モデルを用いた運動軌跡のフィルタリング」
の方を読んでいくと

ノイズベクトルが状態に含めた超パラメータに依存するため,2.3 の状態空間モデルは非線形となる.

とあるので線形なモデルを前提とするカルマンフィルタでは計算できないんでしょう(適当ですみません)

投稿2021/09/14 02:35

ozwk

総合スコア13512

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

shimotani1028

2021/09/14 21:38

ご回答ありがとうございます。 非線形であればヤコビアンを求める必要があると思いますが、FやHに変数があるわけではないので特に求めることができません。そもそもこの係数行列が間違っているんですかね。 でも上記の論文のFとHは同様の考え方なんですよね。他に特別な処理が必要な気しているんですが、それがわからないんですよね。
fana

2021/09/15 01:25 編集

(カルマンフィルタは10年くらい前に一度触ったことがあるだけで,すっかり忘れてしまったので) どの式が何なのかというところから「?」な状態の者です. 故に以下は全く見当違いな可能性が濃厚と思われ,その場合には無視するか罵倒するかしてください. --- カルマンフィルタのノイズのパラメータは定数的に与えるようなものだったと記憶しています(そしてその設定の正解がさっぱりわからんという). それに対して,本件ではそれを変数として扱おう,ということなのかな,と. であれば(素人考えでは),このノイズのパラメータも他の推定対象パラメータ(xか?)と同様の扱いにすることになるのかな? とか. すなわち x_t が x_t-1 と「xに対するノイズ」から更新されるモデルであるように, 「ノイズのパラメータ」もまた, 時刻tでのそのパラメータの値が,時刻t-1での値と「そのパラメータに対するノイズ」から更新されるような式 みたいなのを何かしら用意して突っ込むことになるのでは. (ノイズのパラメータに関して,時間方向の変化を速度とか加速度とかいう話に考えることはできそうもない(?)ので,「そのパラメータに対するノイズ」だけがノイズのパラメータを変化させる源になる…のかな? でもこれだとどこまで行っても定数的に設定せねばならないハイパパラメータが減らない気がするが…) で,そうなると(FなのかHなのかわからんのですが)行列のサイズも相応にでかくなることになるのではないでしょうか?(パラメータベクトルの次元を増やすのだから)
shimotani1028

2021/09/17 07:43

アドバイスまことにありがとうございます。 おっしゃる通りノイズをパラメータとして扱っているので、常に更新する式が必要だと思っています。 ”時刻tでのそのパラメータの値が,時刻t-1での値と「そのパラメータに対するノイズ」から更新されるような式”というのが、図の中のFの2行目(システムノイズ)、3行目(観測ノイズ)とGに当てはまっていると考えております。 ハイパーパラメータのハイパーパラメータについては別途最尤推定法とかで決めてしまえば、いいのかなと思っておりました。
fana

2021/09/17 08:13 編集

よくある {位置x,速度v,加速度a, ...} なモデル: x += v * dt v += a * dt で,加加速度(?)とか加加加速度(??)とかまでは考えない(考えたくない)から, どこかで打ち切って「a = a + ノイズ」的にしちゃうのと同じような扱いにならんのかな? とか思ったのですが,話が全然違うのかも. (あるいは,そういうふうにしたような…という記憶自体が間違いかも) 誠に失礼.罵倒願います.
ozwk

2021/09/17 08:16

@fanaさん > そうなると(FなのかHなのかわからんのですが)行列のサイズも相応にでかくなることになるのではないでしょうか?(パラメータベクトルの次元を増やすのだから) 質問文ですでに入ってますよ。 1要素の状態ベクトルに観測ノイズとシステムノイズの2要素を追加して3要素で扱ってます。
fana

2021/09/17 08:20

(頭の中に何故か3次元空間があって,勝手に x がノイズ含まず状態で3次元と見てました)
shimotani1028

2021/09/17 08:23

僕の書き方が悪かったですね。ごめんなさい。 3次元ではなく1次元の簡単な時系列モデルです。
fana

2021/09/17 08:30

いや,どう考えてもそちらに落ち度は全くないと思います.
guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

まだベストアンサーが選ばれていません

会員登録して回答してみよう

アカウントをお持ちの方は

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.50%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問