简单使用
返回离散傅里叶变换的采样频率
import numpy as np
a = np.fft.fftfreq(10, 1)
a
==>array([ 0. , 0.1, 0.2, 0.3, 0.4, -0.5, -0.4, -0.3, -0.2, -0.1])
可以看到后面一半是负频率,一般只用到前面一半
函数原型
def fftfreq(n, d=1.0):
n: 数组长度
d: 输入信号的采样周期,例如1/16000或1/32000
np.fft.fft()函数一般会配合着np.fft.fftfreq()函数使用
注意fftfreq的入参n要与fft返回的数组长度相同
fft = np.fft.fft(np.arange(10))
fft
==>array([45.+0.00000000e+00j, -5.+1.53884177e+01j, -5.+6.88190960e+00j,
-5.+3.63271264e+00j, -5.+1.62459848e+00j, -5.-1.33226763e-15j,
-5.-1.62459848e+00j, -5.-3.63271264e+00j, -5.-6.88190960e+00j,
-5.-1.53884177e+01j])
与fftfreq返回的数组对应,表示相应频率对应的幅值