スペクトル解析を用いて、時系列データの中で卓越する周期を一瞬で図示 - Qiita

# coding:utf-8 from scipy import fftpack import numpy as np import matplotlib.pyplot as plt class Spectra(object): def __init__(self, t, f, time_unit): """ - t : 時間軸の値 - f : データの値 - time_unit : 時間軸の単位 - Po : パワースペクトル密度 """ assert t.size == f.size # 時間軸の長さとデータの長さが同じであることを確認する assert np.unique(np.diff(t)).size == 1 # 時間間隔が全て一定であることを確認する self.time_unit = time_unit # 時間の単位 T = (t[1] - t[0]) * t.size s