|
状態空間モデルを用いて時系列データに潜む振動成分をデータ駆動的に抽出する振動子分解について、パラメータ推定アルゴリズムの改良を行った。従来のカルマンフィルタによる対数尤度の逐次計算に基づく最尤推定アルゴリズムでは、データの時点数が大きいときに計算量が膨大になってしまうという課題があった。この課題に対するいくつかの対処法を検討した結果、Whittle尤度にサブサンプリングを組み合わせたアルゴリズムが計算時間の削減と安定した推定精度の両立という意味で実用的であることがわかった。このアルゴリズムは振動子分解に限らず一般の状態空間モデルのパラメータ推定に適用可能である。また、振動子分解プログラムのMATLABからpythonへの移植を行なった。
スロー地震の解析では0.02-0.05HzのVLFEや2-8Hzの微動をバンドパスフィルタによって取り出すのが一般的である。しかし、最近の研究によって、スロー地震活動はこれらの周波数帯域に限定されないことを示唆する結果が得られてきている。そこで、先行研究に倣って2016年4月10日の紀伊半島沖DONETデータに注目し、振動子分解によるスロー地震のデータ駆動的な抽出を試みた。結果として、0.01Hz付近と0.1Hz付近に振動子のクラスタが存在することが見つかった。
この他に、作用素スケーリング問題に対するSinkhornアルゴリズムのもつ量子情報幾何構造に関する論文が採択された。
|
