|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
; e( W4 W1 m* O/ _. W; J7 sb=data(:,8);" z) }1 z; ^# ^* l
x=b(5000:6000);%采样数据2 `0 {. n2 @: R+ @/ _+ N
imf=emd(x);, S) ~2 w$ b! x/ ^
[m,n]=size(imf);
1 n0 B6 [2 ]* H# ~' g" O[A,f,tt]=hhspectrum(imf(1:m-1,: ));%计算HHT时频谱
2 D& p" h4 J( }' B- A4 f[im,tt1]=toimage(A,f,tt,length(tt));( M3 R4 O& i M7 k5 V7 c
disp_hhs(im,tt1);%二维图显示HHT时频图
2 Z! M- \3 H# O7 N/ N* Nfor i=1:m-1;
C, B2 B- _- l. o6 z* R" ffa=f(i,: );
% \' i- ?& F3 m. ~0 z% H[FA,TT1]=meshgrid(fa,tt1);%三维图显示HHT时频图
9 p2 U8 f' f5 GsuRF(FA,TT1,im)- W+ K$ }6 j7 S6 }4 x
hold on, I; v- |* O, a0 @. c9 r- e0 A7 {" s
end' u" y- c5 y; Z: c K! u
hold off* p) {6 e, A6 J3 J5 j2 Q
请教大神，为什么二维时频图出不来，三维的可以？非常感谢！ |
|
/1 
发表于 2020-8-18 18:41
收藏
淘帖
支持！
反对！