求助EMD程序如何使用

fengling_335

她用的是小波簇，好的我试试吧

prideheart

原帖由 zhangnan3509 于 2007-6-1 11:32 发表

                               
登录/注册后可看大图

三维，时频，边际

你好，想请教一下，到底怎么来实现你的这个图形？
或者说，是怎么给最后分解出来的变成实际意义的频率？
我看到你在里边说给频率乘上采样频率就可以了，能具体说一下吗？
新手，麻烦了，不好意思，先谢了~~~

zhangnan3509

imagesc(tt1,[0,0.5*fs],E);
set(gca,'YDir','normal')
xray的例子里已经有了啊 如果想用函数解决的话，你可以把disp_hht这个函数改造一下，方法和上面的差不多

maidoudou2008

请问小议EMD程序一文在哪儿啊？没看到

xiaoaq

若是离散数据如X=[12   14     17   29   20     28     27    20.....],程序怎变

baobao1982

http://forum.vibunion.com/forum/viewthread.php?tid=43109

hyhwqy

原帖由 zhlong 于 2007-5-30 17:28 发表

                               
登录/注册后可看大图

t=1:500;
t=t*1/2000;
x=sin(2*pi*20*t);
y=0.4*sin(2*pi*40*t+140);
z=x+y;
plot(t,z)

imf=emd(z);
emd_visu(z,1:length(z),imf)

[A,f,tt]=hhspectrum(imf);
=toimage(A,f);
disp_hhs(im);
colorm ...

您好，这个程序为什么我运行的时候总是出现    ??? Input argument "i" is undefined.
Error in ==> emd_visu at 70
figure(i+1)     运行不出来第三个图形啊

yll

我也画了个这样的图，有点问题请教
??? Input argument "f" is undefined.

Error in ==> mspc at 33
plot(f,ms)
我设置了f，图形就显示不出来了

cboboc

原帖由 zhangnan3509 于 2007-6-1 15:33 发表

                               
登录/注册后可看大图

对应8楼的程序
t=1:500;
t=t*1/400;
x1=sin(2*pi*20*t);
x2=0.4*sin(2*pi*40*t+140);
z=x1+x2
plot(t,z)
y= linspace(1,2,500);%添加趋势项
z=z+y
imf=emd(z);
emd_visu(z,1:length(z),imf)
[A,f,tt]=hhs ...

请问我运行上述程序为什么出现如截图所示的错误呢？ 谢谢

cboboc

hyhwqy你好，我看见你以前问过同样的问题，请问你是怎么解决的呢 ？

maomaoli

请问你的图怎么得来的，程序可以贴上来吗？多谢

chunhuajia

回复#5zhlong:
       你好，我刚学HHT，我用和你一样的方法做希尔伯特谱图
      我有几点不明白的地方：
（1）我的图中有好多亮点，我看了toimage.m和disp_hhs.m，还是不明白，我想问一下这些点是如何在程序中表示的，它们表示的意义是什么，能不能说它们在一定程度上反映了信号的能量呢？
（2）运行结果图中纵坐标（归一化频率）的范围都在0~0.5，为什么要限制在这一范围，是如何计算的？
      谢谢指导！

yanzhang_1

你把i设为1就可以了。会得到3个figure，figure1是emd的结果，figure2，3分别是f2c和c2f

hitdely

请问我运行[A,f,tt]=hhspectrum(imf);个语句的时候怎么总是报错啊 我的imf就是emd的输出用的是法国人的程序

cboboc

这个简单的例子的时频图中怎么会出现波浪的形状呢？照理说说不应该是这样啊？