Python SciPy fft.rfftfreq用法及代碼示例

本文簡要介紹 python 語言中 scipy.fft.rfftfreq 的用法。

用法:

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

返回離散傅裏葉變換采樣頻率(用於 rfft、irfft)。

返回的浮點數組 f 包含頻率 bin 中心，以每單位樣本間隔的周期為單位(開頭為零)。例如，如果樣本間隔以秒為單位，則頻率單位為周期/秒。

給定窗口長度 n 和樣本間距 d：

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

與 fftfreq 不同(但與 scipy.fftpack.rfftfreq 類似)，奈奎斯特頻率分量被認為是正的。

參數 ：：

n： int

窗口長度。

d： 標量，可選

采樣間隔(采樣率的倒數)。默認為 1。

xp： array_namespace，可選

返回數組的命名空間。默認值為“無”，其中使用NumPy。

device： 設備，可選

返回數組的設備。僅當xp.fft.rfftfreq實現設備參數時有效。

返回 ：：

f： ndarray

包含采樣頻率的長度為 n//2 + 1 的數組。

例子：

>>> signal = np.array([-2, 8, 6, 4, 1, 0, 3, 5, -3, 4], dtype=float)
>>> fourier = scipy.fft.rfft(signal)
>>> n = signal.size
>>> sample_rate = 100
>>> freq = scipy.fft.fftfreq(n, d=1./sample_rate)
>>> freq
array([  0.,  10.,  20., ..., -30., -20., -10.])
>>> freq = scipy.fft.rfftfreq(n, d=1./sample_rate)
>>> freq
array([  0.,  10.,  20.,  30.,  40.,  50.])