|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
/ u9 ]) v- I3 E
b=data(:,8);
. A1 o/ I- m2 d, }8 d- Cx=b(5000:6000);%采样数据
$ [5 z" S8 s5 n) L# X& O* S$ Cimf=emd(x);* ^' S. ^* o7 c' U
[m,n]=size(imf);
7 v, b2 ^. p8 F- [5 b6 d[A,f,tt]=hhspectrum(imf(1:m-1,: ));%计算HHT时频谱/ G% h q/ [0 R! i5 I3 p
[im,tt1]=toimage(A,f,tt,length(tt));
" V& u# g- b9 S% Ldisp_hhs(im,tt1);%二维图显示HHT时频图
9 r) ?7 r5 B" I7 V- F% wfor i=1:m-1;7 c4 k7 |% M- P1 k; m" J" Q
fa=f(i,: );
4 ~! k2 B8 u' j7 @5 z& n[FA,TT1]=meshgrid(fa,tt1);%三维图显示HHT时频图
) Q* t8 z s. k: A7 ysuRF(FA,TT1,im)
+ Q! V! Z- o& ohold on
) }" ], B6 ~0 u, X. z' r% t- b( A/ t' c5 Vend
: d/ ~* V: _" U6 `* shold off
$ J. ~# m F+ Q! {& {3 y请教大神,为什么二维时频图出不来,三维的可以?非常感谢! |
|