快速傅里叶变换和汉明窗

Question

我需要帮助在汉明窗口中实现FFT以在图形/图表上搜索峰值（确定光谱相位函数）。此刻我知道如何标记图表上的所有峰值，以及如何在一个图表上加载文件夹中的所有文件。现在最重要的是通过FFT汉明窗口找到一个峰值。对于我来说，对我来说最重要的事情就是将其实现到我的脚本中或者告诉我如何执行此操作。当我理解这一点时，我可以修改它来搜索我想要的峰值。

我的代码：

cormat <- round(cor(test, use = "pairwise.complete.obs"),2)

melted_cormat <- melt(cormat)

ggplot(data = melted_cormat, aes(x=Var1, y=Var2, fill=value)) + geom_tile(color='white')

Answer 1

您的步骤：

计算阵列大小的汉明窗口权重

通过汉明权重乘以数据

制作FFT

在通用语言中，我以这种方式使用HW（N是数组大小）：

  Re[i] := Re[i] * (0.54-0.46*cos(2*Pi*i/N));

但似乎在Matlab中你准备好了function