问题标签 [fftpack]

我曾经scipy.fftpack在 32 位 Python 2.7 上可用，但现在我升级到 64 位 Python 并从这里获得了 SciPy ，我注意到它似乎不包括 FFTPack。

我在哪里可以下载它？

哦，错误是：

假设您使用 dct 函数，然后不处理数据并使用逆变换；反转后的数据与转换前的数据不一样吗？为什么会出现浮点问题？这是报告的问题还是正常行为？

任何人都有解释为什么？当然，一个简单的c*10**-1方法就可以了，但是如果你重复调用函数以在多个维度上使用它，错误会变得更大：

这里是文档的链接。

对于那些使用过 jfftpack 的人，有没有办法根据从 RealDoubleFFT.ft() 命令收到的信息来计算频率？它产生一个可用于绘制 fft 的波表，但我还需要从中获取频率

我正在尝试使用 scipy.fftp 在 python 中平滑和插入一些周期性数据。我已经设法获取数据的 fft，删除高于 wn 的高阶频率（通过执行 myfft[wn:-wn] = 0），然后使用 ifft(myfft) 重建数据的“平滑”版本。ifft 创建的数组与原始数据的点数相同。我如何使用该 fft 创建一个具有更多点的数组。

我还能够以与此处演示的分辨率相同的分辨率手动重新创建 系列使用 FFT 结果重新创建时间序列数据而不使用 ifft

但是当我尝试提高 x 值数组的分辨率时，它也没有给我正确的答案。

提前致谢

尼尔

该scipy.fftpack.rfft函数将 DFT 作为浮点向量返回，在实数部分和复数部分之间交替。这意味着一起乘以 DFT（用于卷积）我将不得不“手动”进行复数乘法，这似乎很棘手。这一定是人们经常做的事情——我想/希望有一个简单的技巧可以有效地做到这一点，但我没有发现？

基本上我想修复这段代码，以便两种方法给出相同的答案：

注意我知道 fftpack 包含一个convolve函数，但我只需要 fft 转换的一半 - 我的过滤器可以提前 fft 一次，然后一遍又一遍地使用。

我是 android 新手，我想将 jfftpack 导入到我在 android studio 中的项目中，我不知道如何导入它。

原始代码是 import ca.uol.aig.realdoublefft ，我不知道将 java 文件放在哪里。我试图把它们放到一个 libs 文件夹中。

这是 jfftpack 源代码：

我在 MATLAB 中编写了一些简单的音频处理算法，将用于 Android 应用程序。我在将 FFT 实现转换为 Java 时遇到了一些麻烦。

我正在尝试使用fftpack的 Java 版本，据报道它很有效，但对我来说也有点神秘，也许是因为我对 FFT 缺乏掌握。

我只是有一个双精度数组，我想从中生成一个对应的复数数组，代表 FFT。我将对 FFT 进行的唯一计算是取其元素的绝对值。

为了清楚起见，这里是我想使用 Java fftpack 复制的基本 MATLAB 代码：

我希望这相当简单，但我无法确定 fftpack 的入口点。

解决了：

LutzL 的以下答案有效，但实际上不需要创建 Complex1D 对象，因为 fftpack 支持双数组作为输入：

我已经阅读了很多关于这个主题的讨论（lomb-scargle 和 fft 之间的比较、在 python 中绘制功率谱、Scipy/Numpy FFT 频率分析等），但仍然无法管理它，所以我需要一些提示。我有一个光子事件列表（检测与时间），数据可在此处获得。列是time、counts、errors和 不同能带中的计数（您可以忽略它们）。我知道来源有一个周期性8.9 days = 1.3*10^-6 Hz。我想绘制在该频率处显示峰值的功率谱密度（可能在对数 x 轴上）。如果我能避免情节的一半（对称），那也很好。这是我到目前为止的代码，到目前为止还没有，但仍然是：

这里产生了（无用的）情节：

我读入了三列数据，所以我有一个 N x 3 数组；前两列是x和y，空间，第三列是值delta。我想得到这个数据的二维快速傅立叶变换，所以，我有 N x 3 数组，从data导入，fftpackscipy

f = fftpack.fft2(data)

我看着f。fN x 3 数组中的哪一列是变换的值，哪一列是 k_x 和 k_y ？

I downloaded numpy1.8 for python 2.7 Installed fine, but when i give:

it gives me a huge error, which ends with

How can i get rid of this error ? I need numpy to perform fft functions.

Any kind of help is appreciated.