我最近不得不处理一个非常大的文件来获得功率谱。该文件有大约 60m 的测量值(可以说它不是那么大,但它已经到了那个地步)。虽然最后(最初我尝试使用 Labview)我设法找到了一些工具(Matlab,R),可以让我在我的 PC 上执行我需要的计算,但我不禁想知道什么是数学上正确的解决该问题的方法,以便我可以分段执行。
我想,我的问题归结为:
我可以对非常大的信号片段进行 FFT,然后将结果组合起来以获得总信号的 FFT 吗?
奖金问题:
周期图会发生什么?
我解释周期图的方式是它对一小部分数据的 fft。那么,如果我获取信号的子集(重叠)并执行 fft 并创建周期图,我会得到截然不同的结果吗?
更新信号 No:2 (500kHz) 的信息:
提前为更改采样率道歉。我刚刚获得了更高采样率的更多结果,因为我怀疑噪声频率太高而无法捕获。
- 信号说明:
该信号来自风洞中的测量值。
当风洞速度设置为特定值时,记录的值是皮托管的电信号。
之所以进行测量,是因为噪声(可能是电的)、信号中正在蔓延,并且需要研究有贡献的频率,以便应用适当的滤波器(不影响风湍流测量的滤波器) .
- 元数据:
数据采集率:500kHz(我刚拿到另一组500kHz的数据)
数据记录率:500kHz
记录时间:大约 1 分钟(这是典型的长度)。
- 预处理:
原始信号在没有预处理的情况下被记录。
对于 FFT,唯一的预处理涉及从信号中去除平均值(以去除直流分量)。(不应用窗口)。
- 数据样本
以下是两个摘录(来自 500kHz 数据的 1000 个数据点)。
第一个是疑似噪声源开启的数据(实际上是为风洞供电的逆变器)。虽然逆变器已打开,但风速设置为零。
第二组是逆变器关闭的摘录(显然风速设置为零),在 y 轴上的比例相同。
下图与在 y 轴上使用自动缩放的第二张图完全相同。
- 较长区间摘录 0.06[s](30000 点)
这是为了更好地了解信号。我切换到点(而不是线)以提供更清晰的画面。以下两张图同样是电噪声 a) 开启和 b) 关闭