2次元平面内を、時間とともにふらふらと動く点の軌跡、すなわち 2次元軌道 を考えます。たとえば、静止立位時の足圧中心（center of pressure: CoP）のゆらぎ（下図）、視覚固視中に生じる眼球の微小な揺らぎ、あるいは画面上でのポインタやカーソルの動きな…

型のパワースペクトルと一口に言っても、指数 の値によって、その統計的・物理的性質は大きく異なります。教科書や論文で典型的に議論されるのは、概ね 程度の範囲であり、この範囲で、弱定常過程、ホワイトノイズ・非整数ガウスノイズ（fractional Gaussian…

ブラウン運動と非整数ブラウン運動（fractional Brownian motion）について、しっかり理解したい人向けに、基礎の基礎から説明してみたいと思います。 1. 生物学者ブラウンの発見から、数学的ブラウン運動へ 2. 実世界で観測した時系列は「離散的」 3. 理想…

近年、統計解析や時系列解析は、多くの分野で「道具」として広く使われるようになっています。その一方で、そのような手法がどのような構造を仮定し、何を計算しているのかを十分に理解しないまま使っている学生や実務者が増えているとも感じます。統計や時…

複数のセンサで観測した信号が、実は複数の独立した「音源」の混合であった、という状況は、信号処理の分野ではよく見られます。そのような問題は ブラインド音源分離（Blind Source Separation, BSS） と呼ばれます。 BSS の代表的な手法としては ICA（Inde…

瞬時周波数や振幅を同時に推定できる TKEO（Teager–Kaiser Energy Operator） は、振幅変調・周波数変調（AM–FM）信号解析の古典的な手法です。しかし、実データへの適用を考えると、従来の TKEO 法には明確な問題があります。ごくわずかな加法ノイズが含ま…

振動する時系列信号を解析するとき、私たちはしばしば「この信号は、いまどれくらいの強さで振動しているのか」「振動の速さは時間とともにどう変わっているのか」といった問いを立てます。Teager–Kaiser energy operator（TKEO）は、まさにそのような問いに…

以前説明した「自己回帰 (AR)パワースペクトルの分解：波素分析」（詳細は以下のリンク参照） 自己回帰 (AR)パワースペクトルの分解：波素分析 - ケィオスの時系列解析メモランダム をRで実際にやってみようという記事です。ARスペクトル分解（波素分析）は…

今回はC言語プログラムを使って、DMA（Detrending Moving Average Analysis）の高速アルゴリズムの処理速度が如何に速いかを見てみます。素朴に Savitzky-Golay フィルタの畳み込みを計算をする DMA アルゴリズムと、我々のグループが提案している高速DMAア…

時系列の長時間相関（スケーリング）を調べるとき、トレンド除去を含む detrended moving-average analysis (DMA) は、DFA の次に良く使われている手法です。実は、DMAには、DFAよりも優れた部分が複数あり、上位互換と言える方法です。しかし、知名度（主に…