Python SciPy fftpack.rfftfreq用法及代码示例

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

用法:

scipy.fftpack.rfftfreq(n, d=1.0)#

DFT 采样频率(用于 rfft、irfft)。

返回的浮点数组包含以周期/单位为单位的频率箱(开头为零)，给定窗口长度 n 和样本间距 d：

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

参数 ：：

n： int

窗口长度。

d： 标量，可选

样本间距。默认值为 1。

返回 ：：

out： ndarray

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

例子：

>>> import numpy as np
>>> from scipy import fftpack
>>> sig = np.array([-2, 8, 6, 4, 1, 0, 3, 5], dtype=float)
>>> sig_fft = fftpack.rfft(sig)
>>> n = sig_fft.size
>>> timestep = 0.1
>>> freq = fftpack.rfftfreq(n, d=timestep)
>>> freq
array([ 0.  ,  1.25,  1.25,  2.5 ,  2.5 ,  3.75,  3.75,  5.  ])