6

我愿意对时间序列数据应用傅立叶变换以将数据转换为频域。我不确定我用来应用傅立叶变换的方法是否正确？以下是我使用过的数据的链接。

读取数据文件后，我使用绘制了原始数据

t = np.linspace(0,55*24*60*60, 55)
s = df.values
sns.set_style("darkgrid")
plt.ylabel("Amplitude")
plt.xlabel("Time [s]")
plt.plot(t, s)
plt.show()

由于数据是每天的频率，我已经将其转换为秒使用24*60*60并使用 55 天55*24*60*60

该图如下所示：原始时间序列

接下来，我使用以下代码实现了傅立叶变换，并获得了如下图像：

#Applying Fourier Transform
fft = fftpack.fft(s)

#Time taken by one complete cycle of wave (seconds)
T = t[1] - t[0] 

#Calculating sampling frequency
F = 1/T

N = s.size

#Avoid aliasing by multiplying sampling frequency by 1/2 
f = np.linspace(0, 0.5*F, N)

#Convert frequency to mHz
f = f * 1000

#Plotting frequency domain against amplitude
sns.set_style("darkgrid")
plt.ylabel("Amplitude")
plt.xlabel("Frequency [mHz]")
plt.plot(f[:N // 2], np.abs(fft)[:N // 2])  
plt.show()

变身

我有以下问题：

我不确定我的上述方法对于实现傅立叶变换是否正确。

我不确定我用来避免混叠的方法是否正确。

如果，我所做的比如何解释频域图中的三个峰值是正确的。

最后，我将如何仅使用重要的频率进行变换。

4

1 回答 1

2

虽然我不会回答你的前两个问题（对我来说看起来不错，但我喜欢专家的意见），但我可以权衡后两个问题：

如果，我所做的比如何解释频域图中的三个峰值是正确的。

嗯，这意味着你的信号有三个主要分量，频率大约为 0.00025 mHz（可能不是这里最好的单位选择！）、0.00125 mHz 和 0.00275 mHz。

最后，我将如何仅使用重要的频率进行变换。

您可以将低于您决定的截止频率的每个频率归零（例如，绝对值 3 - 这应该涵盖您的峰值）。然后你可以这样做：

below_cutoff = np.abs(fft) < 3
fft[below_cutoff] = 0
cleaner_signal = fftpack.ifft(fft)

那应该这样做，真的！

于 2018-11-05T18:09:07.600 回答