EPSの東北地方太平洋沖地震特集に投稿していた論文がようやく受理された。ふぅ。
Kalman filterの更新則について、高須さんが提案している形式を早速こちらのプログラムにも導入してみた。実行速度はほぼかわらず。また、結果もほぼかわらず(もともと不安定がなかったのであたりまえ)。ただ、このコードはのちのちいろんなところで使いまわす予定なので、大規模問題に適用する場合に備えてこちらを採用するようにする。
この形式、ほかにやっている例はないかなぁとひとしきり探してみたが、確かにみあたらない。発想としてはきわめて自然なのだが、なんでだろ?みんな情報フィルタ系を使うのかしら。
Kalman filterの更新則について、高須さんが提案している形式を早速こちらのプログラムにも導入してみた。実行速度はほぼかわらず。また、結果もほぼかわらず(もともと不安定がなかったのであたりまえ)。ただ、このコードはのちのちいろんなところで使いまわす予定なので、大規模問題に適用する場合に備えてこちらを採用するようにする。
この形式、ほかにやっている例はないかなぁとひとしきり探してみたが、確かにみあたらない。発想としてはきわめて自然なのだが、なんでだろ?みんな情報フィルタ系を使うのかしら。
個人的には、EKFの初期値依存性の方がより気になります。線形化しているから仕方がないとは言え、もっといい方法がないですかね。