从信号的 fft 中找到局部最大值

Question

我正在尝试找到信号的 fft 峰值，以用于进一步分析信号。我正在使用数据的 SpanSelect 并进行 fft，表示为频谱。我真的想让绘图具有交互性，并且用户单击一个点进行进一步分析，但我没有找到一种方法可以做到这一点，所以想要一种找到局部频率峰值的方法。频谱可能如下所示：

因此，我想要一种方法来返回峰值为 38 Hz 的频率。有没有办法做到这一点？

Answer 1

用于argrelextrema寻找局部最大值：

import numpy as np
from scipy.signal import argrelextrema
from matplotlib.pyplot import *
np.random.seed()
x = np.random.random(50)
m = argrelextrema(x, np.greater) #array of indexes of the locals maxima
y = [x[m] for i in m]
plot(x)
plot(m, y, 'rs')
show()