|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
& g: u6 a( w" U* A+ }
b=data(:,8);
0 T4 ~. W* M+ ^x=b(5000:6000);%采样数据
; s0 h" r. {) D+ @ D4 Fimf=emd(x);
4 E: k P6 _: z q4 N3 {1 s[m,n]=size(imf);
$ u; P5 ?1 w; o# _- j$ G- h$ ^, x[A,f,tt]=hhspectrum(imf(1:m-1,: ));%计算HHT时频谱
$ Z) H6 m, o. S8 [6 y[im,tt1]=toimage(A,f,tt,length(tt));% F0 j. @0 r n/ r- ?. }
disp_hhs(im,tt1);%二维图显示HHT时频图
4 D3 t+ q7 n6 c1 z% @for i=1:m-1;- ~2 z: Q+ ~6 Q
fa=f(i,: );' M9 J1 p, l5 r1 l! o0 _
[FA,TT1]=meshgrid(fa,tt1);%三维图显示HHT时频图9 q' F6 l1 Y; C& ^& b
suRF(FA,TT1,im)
* z; F9 q* U* x/ Dhold on
) v3 \: _, w: w# T \/ z3 aend
+ D- _' }1 F8 W. @) _* Chold off
- m4 U3 z* [& \" s& x( W! W请教大神,为什么二维时频图出不来,三维的可以?非常感谢! |
|