|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
9 P1 p+ M- m3 O- l1 s/ p
b=data(:,8);6 I% f8 x7 O6 L U4 y0 d; i
x=b(5000:6000);%采样数据) L* _/ w& D+ O' P/ A
imf=emd(x);
/ D9 B* p) F/ ]: E[m,n]=size(imf);
7 X# O+ n& G8 z5 W; R[A,f,tt]=hhspectrum(imf(1:m-1,: ));%计算HHT时频谱
; W$ t) U+ ]3 g0 B[im,tt1]=toimage(A,f,tt,length(tt));# l8 I# r# b" K1 p; W
disp_hhs(im,tt1);%二维图显示HHT时频图" O! z% r% ~1 r; {) h4 j: m2 i" q
for i=1:m-1;
0 q4 W8 H& u; y0 Rfa=f(i,: );: D1 E% X5 K2 w& n3 I: g% ]5 d/ G
[FA,TT1]=meshgrid(fa,tt1);%三维图显示HHT时频图- j/ `2 [' ~' q- N3 E
suRF(FA,TT1,im)
% `0 Y5 N4 n# t3 rhold on) A: ?1 H( P$ n( {8 O( h9 P5 S
end9 }+ [3 ~& m% a7 f
hold off: |8 J1 ^; ]1 b
请教大神，为什么二维时频图出不来，三维的可以？非常感谢！ |
|
/1 
发表于 2020-8-18 18:41
收藏
淘帖
支持！
反对！