我在数学课上了解了傅立叶变换,并认为我已经理解了它们。现在,我正在尝试使用 R(统计语言)并在实践中解释离散 FFT 的结果。这就是我所做的:
x = seq(0,1,by=0.1)
y = sin(2*pi*(x))
calcenergy <- function(x) Im(x) * Im(x) + Re(x) * Re(x)
fy <- fft(y)
plot(x, calcenergy(fy))
并得到这个情节:
如果我理解正确,这代表能量密度谱的“一半”。由于变换是对称的,我可以将所有值镜像到 x 的负值以获得完整的光谱。
但是,我不明白的是,为什么我会得到两个尖峰?这里只有一个正弦频率。这是混叠效应吗?
另外,我不知道如何从这个图中得到频率。让我们假设正弦函数的单位是秒,那么密度谱中 1.0 处的峰值是 1Hz 吗?
再说一遍:我理解 FFT 背后的理论;实际应用是问题:)。
谢谢你的帮助!