HHT的问题........
本帖最后由 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变换,得到时间、频率、振幅的三维图。
请教各位大虾,三维图应该怎么画?如何编程?谢谢!!
本帖最后由 wdhd 于 2016-3-11 11:10 编辑
我也想知道啊,下了论坛里的程序,可惜不会用,盼高手指导啊!!
本帖最后由 wdhd 于 2016-3-11 11:10 编辑
没人来帮忙啊,失望!
本帖最后由 wdhd 于 2016-3-11 11:11 编辑
我也画过三维图,效果不好,imagesc函数画出的好看些,用颜色表示第三维,可以一试
本帖最后由 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分别表示时间 瞬时频率矩阵(行向量)和瞬时振幅矩阵(行向量)
绝对没问题 我用过无数遍了 本帖最后由 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编辑过]
本帖最后由 wdhd 于 2016-3-11 11:11 编辑
请教一下上面几位成功者,该程序EMD分解中端点效应的抑制,是什么方法啊?先谢谢了
页:
[1]