文章/答案/技术大牛

发布

社区首页 >问答首页 >理解时序对象R中的频率(谱分析)

问理解时序对象R中的频率(谱分析)
EN

Stack Overflow用户

提问于 2017-02-17 01:53:37

回答 1查看 856关注 0票数 0

我正在尝试找出我的数据中是否存在6个月的季节性。我想用频谱分析来做这件事，我会在一个特定的频率下查看频谱的Ftest。

如果我基于频率为12的时间序列对象计算x轴，我希望能够帮助我理解频谱分析中x轴的含义(即频率的含义)。

这里生成了一些数据。我的数据是在10年内每月收集的：

data = ts(1:120, frequency = 12, start = c(1959, 1)) # data that is collected each month over 10 years.
library(multitaper)
resSpec <- spec.mtm(data, nFFT = "default", 
                Ftest = TRUE,
                jackknife = FALSE, maxAdaptiveIterations = 100,
                plot = TRUE, na.action = na.fail)

Output from multitaper

如果我正在寻找6个月(每年2次)周期的频谱功率，我是在寻找1/6的频率还是6/12的频率？

time-series

frequency

spectral-density

回答 1

Stack Overflow用户

发布于 2017-02-17 04:31:28

频率输出是每年的循环，而不是每月。您可以这样检查：

curve <- function(i) {
    sin(2*pi/6 * i)
}

data = ts(sapply(1:120, curve), frequency = 12, start = c(1959, 1))

resSpec <- spec.mtm(data, nFFT = "default", 
                Ftest = TRUE,
                jackknife = FALSE, maxAdaptiveIterations = 100,
                plot = TRUE, na.action = na.fail, dtUnits="year", deltat=1/12)

所以你要找的频率是2(周期/年)

票数 0

页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持

原文链接：

https://stackoverflow.com/questions/42281325

复制

相似问题

问理解时序对象R中的频率(谱分析)
EN

回答 1

Stack Overflow用户

社区

活动

圈层

关于

腾讯云开发者

热门产品

热门推荐

更多推荐

问理解时序对象R中的频率(谱分析)EN

回答 1

Stack Overflow用户

社区

活动

圈层

关于

腾讯云开发者

热门产品

热门推荐

更多推荐

问理解时序对象R中的频率(谱分析)
EN