“X”和“Y”的长度与FFT不同R中

Question

您好我接过例子来自：stat.ethz.ch

## Now let's look at some artificial data: 
x <- seq(100000)/1000 # pretend we're sampling at 1 kHz 

## We'll put in two frequency components, plus a dc offset 
f1 <- 5 # Hz 
f2 <- 2 # Hz 
y <- 0.1*sin(2*pi*f1*x) + sin(2*pi*f2*x) + 50 
fft.y <- fft(y) 
delta <- x[2] - x[1] 
f.Nyquist <- 1/2/delta 
f <- f.Nyquist*c(seq(length(x)/2), -rev(seq(length(x)/2)))/(length(x)/2) 

par(mfrow=c(2,2)) 
plot(x,y, type='l', xlim=c(0,20)) 
plot(f, Mod(fft.y), type='l', log='y') 

## Now let's zoom in and mark the points were I expect to see peaks: 
plot(f, Mod(fft.y), type='l', log='y', xlim=c(-10,10)) 

现在我有一个数据帧分析 - DF。 df的第一列（V1）是日期，第二列是值（V2）。我设定了分数，但我总是得到一个更多的fft值（与奈奎斯特做的是什么）。因此得到：“'x'和'y'长度不同”。不知道问题出在哪里！

y <- df$V2 
fft.y <- fft(y) 
delta <- 10 # I know that there are 10sec between values 
f.Nyquist <- 1/2/delta 
f <- f.Nyquist*c(seq(length(df$V1)/2), -rev(seq(length(df$V1)/2)))/(length(df$V1)/2) 

DF的样子：

07032012-185821;20.0 
07032012-185831;12.0 
07032012-185841;14.0 

，并有大约20000个obersvations

感谢您的帮助！

Answer 1

大胆猜测，我认为当观察的长度很奇怪时会出现问题。基本上，

length(c(seq(length(df$V1)/2), -rev(seq(length(df$V1)/2)))) 

不等于

length(df$V2) 

这是由于整数前的舍入：seq(length(1:9)/2)是相同seq(length(1:8)/2)虽然原始向量的长度不同。

编辑： 它可以通过固定：

if (length(df$V2) %% 2 == 1) 
{ 
    f <- f.Nyquist*c(seq(length(df$V1)/2), 0, -rev(seq(length(df$V1)/2)))/(length(df$V1)/2) 
} else 
{ 
    f <- f.Nyquist*c(seq(length(df$V1)/2), -rev(seq(length(df$V1)/2)))/(length(df$V1)/2) 
}