是否有一个简单的解决方案来获得特定信号的傅里叶变换?我已经用GeneCycle尝试了它,但结果似乎是错误的。
信号:https://www.sendspace.com/file/bmwjvr
我尝试过的代码:
trajectory <- testmatrix[,2]
library(GeneCycle)
f.data <- GeneCycle::periodogram(trajectory)
harmonics <- 2:10
plot(f.data$freq[harmonics]*length(trajectory),
f.data$spec[harmonics]/sum(f.data$spec),
xlab="Harmonics", ylab="Amplitute Density", type="h")