HHT的问题........

6543124 · 发表于 2006-5-19 13:03

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

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

x

本帖最后由 wdhd 于 2016-3-11 11:10 编辑

　　我想做的就是t=0:0.001:0.5;

　　x=(1+0.2*sin(2*pi*7.5*t)).*(cos(2*pi*30*t+0.5*sin(2*pi*15*t)))+sin(2*pi*120*t);

　　这个函数进行emd分解，得到imf图，然后进行hht变换，得到时间、频率、振幅的三维图。

　　请教各位大虾，三维图应该怎么画?如何编程?谢谢!!

loryin · 发表于 2006-5-19 18:23

本帖最后由 wdhd 于 2016-3-11 11:10 编辑

　　我也想知道啊，下了论坛里的程序，可惜不会用，盼高手指导啊!!

loryin · 发表于 2006-5-27 00:07

本帖最后由 wdhd 于 2016-3-11 11:10 编辑

　　没人来帮忙啊，失望!

hyl2323 · 发表于 2006-5-27 10:13

本帖最后由 wdhd 于 2016-3-11 11:11 编辑

　　我也画过三维图，效果不好，imagesc函数画出的好看些，用颜色表示第三维，可以一试

lueac · 发表于 2006-6-3 10:28

本帖最后由 wdhd 于 2016-3-11 11:11 编辑

　　我认为三维图最好用散点图画，就是数据量大的时候有点慢但是这个和HHT那篇原著里的图最像，背景和点的颜色可以改。

　　for i=1:IMFN(IMF个数)

　　XX((i-1)*length(t)+1:i*length(t))=t;

　　YY((i-1)*length(W)+1:i*length(W))=W(i,:);

　　ZZ((i-1)*length(A1)+1:i*length(A1))=A1(i,:);

　　end

　　figure(1);

　　scatter(XX,YY,3,ZZ,'filled');

　　其中t、W、A1分别表示时间瞬时频率矩阵(行向量)和瞬时振幅矩阵(行向量)

lueac · 发表于 2006-6-3 10:29

绝对没问题我用过无数遍了

hdwang · 发表于 2006-6-4 09:26

本帖最后由 wdhd 于 2016-3-11 11:11 编辑

　　lueac你好，我用了你这段程序，显示错误信息

　　??? Undefined function or variable 'XX'.

　　Error in ==> C:\Documents and Settings\wang\桌面\singal.m

　　On line 20 ==> scatter(XX,YY,3,ZZ,'filled');

　　可否详细给出楼主问题答案，

　　谢谢阿

　　

[此贴子已经被作者于2006-6-4 9:29:31编辑过]

loryin · 发表于 2006-6-8 01:45

本帖最后由 wdhd 于 2016-3-11 11:11 编辑

　　请教一下上面几位成功者，该程序EMD分解中端点效应的抑制，是什么方法啊?先谢谢了

账号		自动登录	找回密码
密码			我要加入

[HHT] HHT的问题........

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