python - 科学 : fourier transform of a few selected frequencies

Question

我在信号上使用 scipy.fft，使用移动窗口绘制随时间变化的频率幅度(这里是 an example，时间在 X 上，频率在 Y 上，并且振幅是颜色)。

但是，只有少数频率让我感兴趣(仅约 3、4 个频率)。使用 FFT 似乎我不能只选择我想要的频率(显然频率范围是由算法决定的)，所以我计算了很多无用的东西，我的程序甚至因 MemoryError< 而崩溃 如果信号太长。

我该怎么办？我是否必须使用自定义傅立叶变换 - 在这种情况下，欢迎提供良好实现的链接 - 或者是否有 scipy 方式？

---

编辑

@jfaller 回答后，我决定(尝试)实现 Goertzel 算法。我想出了这个:https://gist.github.com/4128537但它不起作用(频率 440 没有出现，没关系峰值，我懒得应用适当的窗口)。任何帮助!？我不擅长 DSP。

Answer 1

您确实希望使用 Goertzel 算法:http://en.wikipedia.org/wiki/Goertzel_algorithm .基本上，它是一个单点的 FFT，如果您只需要信号中有限数量的频率，它是高效的。如果您无法从维基百科中提取算法，请返回，我会帮助您。另外，如果你用谷歌搜索一些资源，就会发现有用 python 编写的 DTMF 解码器(按键式电话解码器)。你可以看看他们是怎么做到的。