【发布时间】：2014-01-12 17:52:34
【问题描述】：
我是 FFT 的新手，我被要求找到一种方法来分析/处理石油钻井平台收集的一组特定数据。 由于钻机移动（例如潮汐和波浪的上下波动），收集的数据中有很多噪音。 我被要求用 FFT=>filtering=>IFFT 清理收集的数据。
我使用 C++ 和 FFTW 3.3.3 库。
一个例子比什么都好：
我有一个数据库，例如泥流（升/分钟）。泥流每 5 秒收集一次，每个度量的 DB 中都有一个时间戳（例如 1387411235）。
所以我的 FFT 的 IN_data 是一对时间戳/泥流（例如 1387456630/3955.94、1387456635/3954.92 等...）
显示这些数据确实看起来像一个嘈杂的声音信号，相关事件可能会被噪音掩盖。
使用在 Internet 上找到的示例，我可以设法执行 FFT，但我缺乏知识和理解是一个大问题，因为我从未从事过信号处理和傅立叶变换。
我真的不知道如何开始这项工作，使用哪个版本的 FFTW 例程（c2c、r2c 等...），如果有任何数据前处理和/或后处理处理要做。 我在互联网上阅读了很多示例和教程，但我是法国人（抱歉我在这里犯了错误），这对我来说并不总是有意义，尤其是对于 OUT_data 单元、OUT_data 类型、输入和输出数据数组大小，窗口化（顺便说一句），简而言之我迷路了......
我想我的问题对于习惯 FFTW 的人来说会很简单，但对我来说现在就非常复杂了。
所以我的问题：
- 在两种方式（FFT 和 IFFT）中使用什么 FFTW 例程（IN_data 和 OUT_data 数组的种类、类型和大小）。
- 如何解释结果数组（FFTW 将返回的单位是什么）。
现在我所做的一个简短示例是：
fftw_plan p;
p = (fftw_plan)fftw_plan_dft_1d(size,in,out,FFTW_FORWARD,FFTW_ESTIMATE);
fftw_execute(p);
fftw_destroy_plan(p);
使用“in”和“out”作为 fftw_complex（我的 In_data 数组的复数元素为每个数据设置为 1，不知道为什么，但教程说要这样做）。
此代码基于互联网上的一个示例，但我缺乏知识/理解是一个很大的拖累，我想知道这里是否有人可以给我解释/工作流程/见解/链接如何拉这个出去。
我正处于新工作的试用期，我真的很想为我的老板实现这个功能，即使这意味着四处寻求帮助，我在这里看到了很多 FFTW 熟练的帖子......
【问题讨论】：
-
不要为自己是法国人而道歉 ;)
-
对于我缺乏技能/知识的朋友，我深表歉意……但还是谢谢……
-
FFTW 对于一个完整的 DSP 菜鸟来说是一个相当复杂的包 - 我建议从更简单的东西开始，例如KissFFT。此外，您需要确定您的数据是均匀采样的，否则基于 FFT 的方法将无法正常工作。
-
感谢 Paul 的回答，测量每 5 秒进行一次，所以我猜采样是统一的（在极少数情况下，未收集数据，但数据库中始终存在相应的时间戳）。
-
好的 - 我的观点是，对于均匀采样的数据，样本之间的间隔必须正好为 5 秒，并且不能丢失数据点，否则您将需要比 FFT 更复杂的方法。