声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2129|回复: 6

[HHT] 求教这个HHT三维谱是怎么画出来的?

[复制链接]
发表于 2014-10-14 19:51 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
Cena马给了这个图谱,但是他始终没有说明这个图是怎么画的,求大神指教
1.jpg
2.jpg
回复
分享到:

使用道具 举报

 楼主| 发表于 2014-10-28 09:31 | 显示全部楼层
算了,没人回答。cena马也太小气了。我自己已经解决了。
发表于 2016-5-26 21:09 | 显示全部楼层
你好 请问你能分享一下你图谱的程序吗?
发表于 2016-5-27 08:51 | 显示全部楼层
是三维图吗?matlab自带命令能实现吗?
发表于 2016-5-31 18:02 | 显示全部楼层
求指教,怎么出来的
发表于 2016-6-1 10:41 | 显示全部楼层
楼主为什么不回复
发表于 2016-6-1 10:44 | 显示全部楼层
  1. function plot_hht_3d(imf,numfreq,fs,ANGLE)
  2. if nargin<3
  3.     fs=1;
  4.     ANGLE=[-37.5,30];
  5. end
  6. if nargin<4
  7.     if size(fs,2)>1
  8.         ANGLE=fs;
  9.         fs=1;
  10.     else
  11.         ANGLE=[-37.5,30];
  12.     end   
  13. end
  14. N=size(imf,2);
  15. [A,f,tt]=hhspectrum(imf);
  16. [m,n]=size(f);
  17. MaxFreq=max(max(f));
  18. MaxFreq=ceil(MaxFreq/0.5)*0.5;
  19. if nargin<2
  20.     numfreq=512;
  21. end
  22. df=linspace(0,MaxFreq,numfreq);
  23. Spectrum=zeros(numfreq,n);
  24. Temp=f;
  25. Temp=min(round((Temp/MaxFreq)*numfreq)+1,numfreq);
  26. for k=1:m
  27.     for u=1:n
  28.         Spectrum(Temp(k,u),u)=Spectrum(Temp(k,u),u)+A(k,u);
  29.     end
  30. end
  31. df=df*fs;
  32. figure
  33. clf
  34. mesh(tt,df,Spectrum);
  35. set(gca,'XLim',[0,N]);
  36. xlabel('采样点数/n');
  37. if fs==1
  38.     ylabel('归一化频率');
  39. else
  40.     ylabel('频率/Hz');
  41. end
  42. zlabel('幅值');
  43. title('三维联合时频图');
  44. colormap jet;
  45. shading interp;
  46. view(ANGLE(1),ANGLE(2));
  47. set(gca,'YLim',[0,fs/2]);
  48. end
复制代码

要要切克闹!!!!
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-5 13:51 , Processed in 0.072632 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表