R语言学习 - 由I/Q数据得到频谱

idata和qdata是从excel表中导出的数据

plotFFT <- function(idata, qdata)
{
  sample_rate <- 1.6e9;
  idata <- data.matrix(idata);
  qdata <- data.matrix(qdata);
  num_samples <- dim(idata)[1];
  x <- 1:num_samples;
  z <- idata + qdata * 1i;
  f <- x * (sample_rate/num_samples);

  plot(f[1:num_samples/2], abs(fft(z))[1:num_samples/2]);
}