假设我有一个时间序列t,其中有100次测量,每个条目代表每天的测量值。我假设信号中存在一些周期性-可能每天,每周或每月重复一次。
将时间序列转换为傅立叶域可能有助于找到这样的周期性?
如何使用numpy的fft模块查找时间序列中最可能的时间段?
来自scipy.signal的函数argrelextrema不会检测到平坦的极值.例:
import numpy as np
from scipy.signal import argrelextrema
data = np.array([ 0, 1, 2, 1, 0, 1, 3, 3, 1, 0 ])
argrelextrema(data, np.greater)
(array([2]),)
Run Code Online (Sandbox Code Playgroud)
检测到第一个max(2),未检测到第二个max(3,3).
这种行为的解决方法是什么?谢谢.