限制scipy.signal.spectrogram仅计算特定频率
限制scipy.signal.spectrogram仅计算特定频率
提问于 2018-01-04 22:18:07
我正在遵循scipy.signal.spectrogram的教程示例。有没有可能限制我们想要计算频谱的频率，比如说在1000-2000 in的范围内。根据FFT，我们将在一半的采样频率下获得它。如果不是，那么我是否可以切片输出以获得t,f,Sxx的3D矩阵，其中我将频率限制为1000-2000 Hz？在我的例子中的用法是，我的信号只包含一定范围的频率，并且保持所有频率对于我的分析来说将是昂贵的。
回答 1
Stack Overflow用户
回答已采纳
发布于 2018-01-05 10:41:44
除非你要分析一小部分频谱，否则用FFT计算整个频谱(并丢弃你不想要的部分)通常仍然是获得部分频谱的一种相当有效的方法。
然后，您可以使用numpy.where在f中找到感兴趣的频率(遵循this answer by askewchan)。下一步将是简单地对输出进行切片，在1D阵列f中提取所需的频率值，并在2D阵列Sxx中提取相关的功率谱值。这可以按如下方式完成：
...
f, t, Sxx = signal.spectrogram(x, fs)
fmin = 1000 # Hz
fmax = 2000 # Hz
freq_slice = np.where((f >= fmin) & (f <= fmax))
# keep only frequencies of interest
f = f[freq_slice]
Sxx = Sxx[freq_slice,:][0]页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/48097164
复制相关文章
相似问题
腾讯云开发者
