用了快速傅里叶变换函数FFT,我是分析时间序列,有458个数据,已转化为时间序列格式
代码如下
 Fs<-458   
 T<- 1/Fs
 y<- fft(xxx)
 L<- length(y)
 f<- (1:L)/2
 y<- y/L *2
 po<- abs(y[1:L/2])
# plot(f,po)
###滤波
 thold<- 10
 y[thold:(L-thold)] <- 0
  X1<-fft(y, inverse = TRUE) / length(y)   #傅里叶逆变换
 plot(f,X1)
Warning message:
In xy.coords(x, y, xlabel, ylabel, log) : 强制改变时丢弃了虚数部分
我看了下X1确实有复数部分
我就是想滤波后,看一下结果,请问除了逆变换还有什么方法吗?