|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
' L* p: P, x7 Z7 P$ C! hb=data(:,8);/ k! u$ [" }1 G$ i4 V* [5 F
x=b(5000:6000);%采样数据; z4 g4 y {) S% R. J
imf=emd(x);
+ H/ ?% j8 E7 T3 t( s: C[m,n]=size(imf);
4 w% K$ i* u9 L8 F3 E[A,f,tt]=hhspectrum(imf(1:m-1,: ));%计算HHT时频谱/ @/ [* {* _+ @ t' l
[im,tt1]=toimage(A,f,tt,length(tt));# X3 _. C9 Q3 W" U
disp_hhs(im,tt1);%二维图显示HHT时频图
6 o( a% u9 `- g. Y w8 V$ p% S, vfor i=1:m-1;
/ @5 O# B- x9 \3 \, [fa=f(i,: );
* F- s" D/ F4 @+ S[FA,TT1]=meshgrid(fa,tt1);%三维图显示HHT时频图
, b# H! L# j3 JsuRF(FA,TT1,im)2 V& g: o$ C" N, N
hold on
* U- T3 }& E4 ^4 v: h2 p! Z: nend" {# n) r( D6 L D/ g* _6 i( X; D5 s
hold off
& V1 H# Z3 }8 Q# U- F请教大神，为什么二维时频图出不来，三维的可以？非常感谢！ |
|
/1 
发表于 2020-8-18 18:41
收藏
淘帖
支持！
反对！