fftfreq函数的用法

fftfreq函数的用法：

from scipy.fftpack import fftfreq
from numpy.fft import fftfreq #numpy和scipy都有fftfreq函数
signal = np.array([-2, 8, 6, 4, 1, 0, 3, 5], dtype=float)
fourier = np.fft.fft(signal)
n = signal.size
timestep = 0.1
freq = np.fft.fftfreq(n, d=timestep)
freq
    array([ 0.  ,  1.25,  2.5 , ..., -3.75, -2.5 , -1.25])

我们在画频谱图的时候需要对信号做傅里叶变换，x轴对应的是频率范围，
使用fftfreq的好处就是可以自动生成频率范围，而不用去考虑信号长度是奇数还是偶数的问题。
处理实际信号的时候这样使用：fftfreq(len(signal),1/samplerate),即fftfreq(信号长度，1/采样率)。注意这里的信号长度是信号做傅里叶变换之后的原始长度（画频谱图的时候则取一半长度），得到频率范围后发现频率有负值，这时同样对频率取一半即可。

频谱图如下：