【python3】時系列データの自己相関関数を求める

import statsmodels.api as sm
a = sm.tsa.acf(ts, nlags=n

nlagsは書かないと、すごく少ない数しか結果が出てこないので必須
nlagsの値はデータのサンプル数で置き換えてOK

gwpyでプロットしたいときは、こんな感じ

import statsmodels.api as sm
from gwpy.timeseries import TimeSeries

a = sm.tsa.acf(ts, nlags=int(T*fs))
a = TimeSeries(a, sample_rate=fs)
a.plot()

自己相関関数のプロットは、コレグラムとか言うらしい
へぇ〜

https://qiita.com/hadacchi/items/ea61d6d714c5a7481461

の計算速度の比較を見ると、自作するよりもそれなりに早いらしい

ランキング参加中です

↓クリックしていただけると嬉しいです〜