【求助】R语言 用傅里叶变换滤波后怎么画图

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

用了快速傅里叶变换函数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确实有复数部分
我就是想滤波后，看一下结果，请问除了逆变换还有什么方法吗？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：傅里叶变换 傅里叶 R语言 时间序列 FFT

前面有考虑复数是误差，但是复数和实数没有相差很大啊
大概是2.052763e-04+6.77946e-05i
          1.987376e-04+6.56543e-05i
说明结果是不可以用的吗？