频带部分重叠的多分量线性调频信号用HHT分解不了

zhangnan3509

我做出来的结果和楼主一样，但是分不开，是不是版本的问题，我用的是2005版，而不是新版

[ 本帖最后由 zhangnan3509 于 2007-11-22 22:16 编辑 ]

zhlong

imf=emd(x,'STOP',[0.0001,0.0005,0.0005],'MAXITERATIONS',1000000);
试试这个。

zhangnan3509

imf=emd(x,'STOP',[0.05,0.5,0.05],'MAXITERATIONS',2000);
比起这个 分解细了和运算量大了

kevin19821

相当有价值

xueer534

在输入信号形式已知，而具体参数未知的情况下，如何设定缺省值是个问题。

form

原帖由 xueer534 于 2007-11-16 21:11 发表
比如说5~25Hz和15~35Hz的的线性调频信号在同一时间段叠加，就无法完全分解呀
是不是说EMD不适用于此类信号

现在这样来，一个频率升的，一个频率降的两个频带有重叠的信号；还以前面的信号为例，大家看看效果如何
t=0:0.001:0.5;
x=sin(2*pi*(5+40*t).*t);
y=sin(2*pi*(40+15*t).*t);
z=x+y;
运行了很长时间，好像分的不是很好

zhangnan3509

效果不如zhlong的明显，尤其是前两阶变化趋势上，再加上阶数要比zhlong的多出3阶，虚假分量有点多啊

form

顺带问下你的附件图片怎么保存为gif格式的，只有5K多

form

是啊，第二个分量还差不多，第一个频率没变似的，就尾部变了些

zhangnan3509

呵呵 这是zhlong介绍的软件，名叫snagit，我把保存下来的bmp文件，转化成gif文件，大小可以节省很多，上传之后打开很方便

zhangnan3509

照理说不应该啊，我把参数变了之后，效果和上面的一样。你是怎么操作的？直接处理的？还是改了默认值？改成什么数值？

[ 本帖最后由 zhangnan3509 于 2007-11-26 11:02 编辑 ]

form

:@D 那个软件好呀，节省空间而且很清楚，这么一看我那个图好糙~~
和哪个一样？我按照zhlong 的设置的细致些imf=emd(x,'STOP',[0.0001,0.0005,0.0005],'MAXITERATIONS',1000000);就会得到那个结果。

[ 本帖最后由 form 于 2007-11-26 11:17 编辑 ]

zhlong

首先你的两个信号都是频率线性增大的，没有频率减小的。
第二，用短时傅立叶等时频分析方法分析你的信号，发现里面怎么有一个是频率不变的分量？（其实是频率变化本来就很小，即第二分量y，频率从40Hz增加到47.5Hz，视觉上看起来就很小）


[ 本帖最后由 zhlong 于 2007-11-26 11:29 编辑 ]

zhangnan3509

form分解出来的第一阶imf频率也是几乎没有变化

zhlong

频率本来变化就很小，40～47.5Hz，而采样频率为1000Hz，所以看起来几乎没变化很正常。