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.  ])