python - scipy.fft 挂起某些声音文件

Question

scipy.fft 在运行这个简单的脚本时似乎挂起:

import scipy
from scipy.io import wavfile

sound = 'sounds/silence/iPhone5.wav'
fs, data = wavfile.read(sound)
print scipy.fft(data)

在某些文件上。尝试 this例如文件。

我注意到的一些事情:

1. 从交互式解释器运行单个命令不会挂起。
2. 与其他声音文件一起运行并不总是挂起脚本(但不只是这个文件不起作用)
3. 有时我得到WavFileWarning: chunk not understanded，但它似乎与它发生的时间无关
4. 如果我使用 Ctrl+C 终止脚本，我得到的结果就像它从未卡住一样。
5. 使用 wave 或 audiolab 打开文件会导致相同的结果。

这是错误还是我做错了什么？

Answer 1

查看系统挂起文件的data.shape值。如果您的数据长度恰好是一个质数，或者是几个大质数的乘积，则 FFT 算法无法加快 DFT 的计算速度。如果您用零填充，或将数据修剪为最接近的 2 的幂，一切都应该运行得非常非常快。