变时数据的傅里叶变换

Question

所以..。我不知道我该在哪里问这个问题，但我希望这里的人至少能给我指明正确的方向。我有一个时间序列，我想做频谱分析，但我找不到任何工具进行FFT，以适应不同的时差之间的数据点(他们都假设dt是常数)。有没有人知道有一个工具能解决这个问题(我特别想找一个周期图或其他方法来确定周期)。

我唯一的想法是在一个特定的时间间隔内做数据点之间的线性插值，给数据一个恒定的dt，但是我担心这样会发现新的光谱分析数据。

这是一小块数据，时间，数据，dt

 time    data        dt
39.630  49662.1     0.170
39.810  49582.5     0.180
40.150  49430.0     0.340
40.320  49413.8     0.170
40.490  49324.0     0.170
40.670  49092.5     0.180
40.830  49025.6     0.160
41.010  49101.5     0.180

有什么建议吗？

Answer 1

FFT后的大多数数学理论都需要固定的采样周期。我建议你做以下几件事：

1. 创建一个时间间隔相等的序列
2. 根据相邻点计算每个瞬间的值(线性或二次插值应该这样做)。为了得到最佳的近似，你可以使用你想要的任意多个点。
3. 根据您希望得到的结果的详细程度，使用参数或非参数方法估计频谱: Burg方法可能有用，LPC/AR模型也可能有用。

查看Mathworks的链接：

http://www.mathworks.com/help/signal/nonparametric-spectral-estimation.html

http://www.mathworks.com/help/signal/parametric-spectral-estimation.html