r - 如何将频谱密度转换为幅度

Question

我有从 R“stats”包中的频谱函数调用产生的标准周期图。它在 Y 轴上产生光谱密度。我希望实际检查关键频率信号的幅度。

如何将频谱密度转换为幅度？R中是否有周期图/分析可以自动生成频率与幅度图？感谢任何建议。

Answer 1

也许您使用的术语与我不同。帮助页面说从频谱函数返回的值是一个列表，其前两个元素是：

freq    
vector of frequencies at which the spectral density is estimated. 
(Possibly approximate Fourier frequencies.) The units are the reciprocal 
of cycles per unit time (and not per observation spacing): see ‘Details’

spec    
Vector (for univariate series) or matrix (for multivariate series) of 
estimates of the spectral density at frequencies corresponding to freq.

那么$spec元素就是你所说的“幅度”向量吗？（您还没有说“关键频率信号”是什么，所以我只是从示例中选择了第四个频率?spectrum：

 lh.spec <- spectrum(lh)

 lh.spec$freq[4]
#1] 0.08333333
  lh.spec$spec[4]
#[1] 1.167519