声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 4602|回复: 15

[FFT] 有谁给我看看我这个信号分解成IMF第一个分量对不对啊?

  [复制链接]
发表于 2011-3-7 12:18 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
未命名.jpg
上面一个是IMF第一个分量,下面是差残,总觉得不对啊,论坛中的大牛帮忙指正一下啊

点评

反对: 1.0
反对: 1
下次发帖请选择正确的帖子分类,这样有更多的人来帮忙~  发表于 2011-4-10 16:34
回复
分享到:

使用道具 举报

发表于 2011-3-7 14:56 | 显示全部楼层
不对,是不是你程序中残差把IMF给覆盖了造成的?
发表于 2011-3-7 21:35 | 显示全部楼层
本帖最后由 Crocodilesky 于 2011-3-7 21:38 编辑

EMD分解是从高到低。从你的残差看,似乎分解不对。
发表于 2011-3-7 21:38 | 显示全部楼层
你的输入信号是什么 ?只有知道了输入信号的成分才能知道你的IMF对不对
 楼主| 发表于 2011-3-8 08:48 | 显示全部楼层
未命名.jpg 输入是这样的,大家给看一下吧,是怎么回事啊,我对HHT刚研究不久!
发表于 2011-3-10 14:44 | 显示全部楼层
回复 5 # flq103 的帖子

你这是什么信号啊?
发表于 2011-3-10 15:23 | 显示全部楼层
你这个是不是分段余弦信号啊,要是的话分解的没错,直接满足IMF条件,不用EMD了~,直接做Hilbert变化就行~
发表于 2011-3-11 08:01 | 显示全部楼层
什么信号?送过来瞧一瞧。
发表于 2011-3-11 15:27 | 显示全部楼层
我觉得你第一个图,是将5个IMF分量画在一起了,这个残差可能是信号和那5个IMF想减之后得到的。
 楼主| 发表于 2011-3-12 15:36 | 显示全部楼层
p=zeros(1,20001);
t=0:0.0001:2;
aa=find(t>0.325&t<=0.556);
    p(aa)=0.9;  
aa=find (t>1.358&t<=1.647);
    p(aa)=1.1;  
aa=find ((t>0.556&t<=1.358)|(t>0&t<=0.325)|(t>1.647&t<=2));
    p(aa)=1;
s=p.*cos(2*pi*50*t).*(t>0&t<=2) +0.1*sin(2*pi*1.5*50*t).*(t>0&t<=0.5)+0.25*sin(2*pi*3*50*t).*(t>=0.5&t<=1)+0.2*sin(2*pi*4*50*t).*(t>=1&t<=2);
plot(t,s)上面是该程序,高手们看看吧!!
发表于 2011-3-14 10:55 | 显示全部楼层
本帖最后由 gczhang 于 2011-3-14 10:55 编辑

我试了一下,结果供参考.


EMD.png
发表于 2011-3-14 10:59 | 显示全部楼层
本帖最后由 gczhang 于 2011-3-14 10:59 编辑

总残差图,非常小。
EMD_残差.png
 楼主| 发表于 2011-3-17 10:36 | 显示全部楼层
本帖最后由 flq103 于 2011-3-17 10:38 编辑

回复 12 # gczhang 的帖子

那我分析出来残差怎么是那样的啊,高手帮忙分析一下啊,还有你分析出来的瞬时幅值和瞬时频率怎么样啊
发表于 2011-3-17 11:35 | 显示全部楼层
本帖最后由 gczhang 于 2011-3-17 13:20 编辑

我感觉你的那个不是残差,更像是最后一阶 IMF8。看我得到的IMF8,和你的所谓残差基本一致。
不知怎么的,图片传不上来。
 楼主| 发表于 2011-3-17 21:36 | 显示全部楼层
回复 13 # flq103 的帖子

是这样啊, 我用的是法国那个人编程序,你的HHT程序边界效应是不是处理过?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-25 08:40 , Processed in 0.062725 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表