Python SciPy fft.fftfreq用法及代码示例

本文简要介绍 python 语言中 scipy.fft.fftfreq 的用法。

用法:

scipy.fft.fftfreq(n, d=1.0, *, xp=None, device=None)#

返回离散傅里叶变换采样频率。

返回的浮点数组 f 包含频率 bin 中心，以每单位样本间隔的周期为单位(开头为零)。例如，如果样本间隔以秒为单位，则频率单位为周期/秒。

给定窗口长度 n 和样本间距 d：

f = [0, 1, ...,   n/2-1,     -n/2, ..., -1] / (d*n)   if n is even
f = [0, 1, ..., (n-1)/2, -(n-1)/2, ..., -1] / (d*n)   if n is odd

参数 ：：

n： int

窗口长度。

d： 标量，可选

采样间隔(采样率的倒数)。默认为 1。

xp： array_namespace，可选

返回数组的命名空间。默认值为“无”，其中使用NumPy。

device： 设备，可选

返回数组的设备。仅当xp.fft.fftfreq实现设备参数时有效。

返回 ：：

f： ndarray

包含样本频率的长度为 n 的数组。

例子：

>>> signal = np.array([-2, 8, 6, 4, 1, 0, 3, 5], dtype=float)
>>> fourier = scipy.fft.fft(signal)
>>> n = signal.size
>>> timestep = 0.1
>>> freq = scipy.fft.fftfreq(n, d=timestep)
>>> freq
array([ 0.  ,  1.25,  2.5 , ..., -3.75, -2.5 , -1.25])