Python SciPy fftpack.irfft用法及代碼示例

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

用法:

scipy.fftpack.irfft(x, n=None, axis=-1, overwrite_x=False)#

返回實序列 x 的離散傅裏葉逆變換。

的內容x被解釋為的輸出scipy.fftpack.rfft函數。

參數 ：：

x： array_like

要反轉的轉換數據。

n： 整數，可選

傅裏葉逆變換的長度。如果 n < x.shape[axis]，x 被截斷。如果 n > x.shape[axis]，x 是零填充的。默認結果為 n = x.shape[axis]。

axis： 整數，可選

計算 ifft 的軸；默認值在最後一個軸上(即，axis=-1)。

overwrite_x： 布爾型，可選

如果為 True，則 x 的內容可以被銷毀；默認為假。

返回 ：：

irfft： 浮點數數組

離散傅裏葉逆變換。

注意：

返回的實數數組包含：

[y(0),y(1),...,y(n-1)]

其中 n 是偶數：

y(j) = 1/n (sum[k=1..n/2-1] (x[2*k-1]+sqrt(-1)*x[2*k])
                             * exp(sqrt(-1)*j*k* 2*pi/n)
            + c.c. + x[0] + (-1)**(j) x[n-1])

對於 n 是奇數：

y(j) = 1/n (sum[k=1..(n-1)/2] (x[2*k-1]+sqrt(-1)*x[2*k])
                             * exp(sqrt(-1)*j*k* 2*pi/n)
            + c.c. + x[0])

抄送表示前麵表達式的複共軛。

有關輸入參數的詳細信息，請參閱 rfft 。

要處理具有複雜數據類型的 (conjugate-symmetric) frequency-domain 數據，請考慮使用較新的函數 scipy.fft.irfft 。

例子：

>>> from scipy.fftpack import rfft, irfft
>>> a = [1.0, 2.0, 3.0, 4.0, 5.0]
>>> irfft(a)
array([ 2.6       , -3.16405192,  1.24398433, -1.14955713,  1.46962473])
>>> irfft(rfft(a))
array([1., 2., 3., 4., 5.])