python和idl_IDL和Python之间的FFT差异

我将一些简单的IDL代码传递给

Python.但是,从SciPy / NumPy包中返回的FFT值与IDL包不同,我无法找到原因.

将所有内容简化为8个元素的简单示例我发现SciPy / NumPy例程返回的值比IDL大8(2 ^ 3)倍(我认为是规范化问题).

以下是两种语言的示例代码(从here复制)：

IDL

signal = ([-2., 8., -6., 4., 1., 0., 3., 5.])

fourier = fft(signal)

print, fourier

回报

( 1.62500, 0.00000) ( 0.420495, 0.506282) ( 0.250000, 0.125000) ( -1.17050, -1.74372) ( -2.62500, -0.00000) ( -1.17050, 1.74372) ( 0.250000, -0.125000) ( 0.420495, -0.506282)

Python

from scipy.fftpack import fft

import numpy as N

…

signal = N.array([-2., 8., -6., 4., 1., 0., 3., 5.])

fourier = fft(signal)

print fourier

回报

[ 13. +0.j , 3.36396103 +4.05025253j, 2. +1.j , -9.36396103-13.94974747j, -21. +0.j , -9.36396103+13.94974747j, 2. -1.j , 3.36396103 -4.05025253j]

我用NumPy包做到了,得到了相同的结果.我试图打印fft(信号,8)以防万一,但它返回相同,如预期的那样.

然而,并非所有,回到我真正的256个元素阵列,我发现差异不再是8或256,而是256 * 8！这只是疯了.

虽然我解决了这个问题,但我需要知道为什么会有这种差异.

解决了：这只是规范化,在某些时候我将IDL 256阵列分成8倍,我忘了删除.在Dougal的回答中,有我错过的文档.