请选择 进入手机版 | 继续访问电脑版
楼主: Melody810
2849 1

[问答] 【求助】R语言 用傅里叶变换滤波后怎么画图 [推广有奖]

  • 0关注
  • 0粉丝

初中生

14%

还不是VIP/贵宾

-

威望
0
论坛币
25 个
通用积分
0
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
49 点
帖子
4
精华
0
在线时间
29 小时
注册时间
2015-10-29
最后登录
2018-12-13

Melody810 发表于 2018-9-14 10:37:23 |显示全部楼层

+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确实有复数部分
我就是想滤波后,看一下结果,请问除了逆变换还有什么方法吗?

二维码

扫码加我 拉你入群

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

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

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

Melody810 发表于 2018-9-14 11:26:35 |显示全部楼层
前面有考虑复数是误差,但是复数和实数没有相差很大啊
大概是2.052763e-04+6.77946e-05i
          1.987376e-04+6.56543e-05i
说明结果是不可以用的吗?

使用道具

您需要登录后才可以回帖 登录 | 我要注册

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2022-6-25 12:40