全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
3956 1
2018-09-14
用了快速傅里叶变换函数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确实有复数部分
我就是想滤波后,看一下结果,请问除了逆变换还有什么方法吗?

二维码

扫码加我 拉你入群

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

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

全部回复
2018-9-14 11:26:35
前面有考虑复数是误差,但是复数和实数没有相差很大啊
大概是2.052763e-04+6.77946e-05i
          1.987376e-04+6.56543e-05i
说明结果是不可以用的吗?
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群