求教这个HHT三维谱是怎么画出来的？

airazor

Cena马给了这个图谱，但是他始终没有说明这个图是怎么画的，求大神指教

airazor

算了，没人回答。cena马也太小气了。我自己已经解决了。

zou1017

你好 请问你能分享一下你图谱的程序吗？

Catsayer

是三维图吗？matlab自带命令能实现吗？

归零

求指教，怎么出来的

jiangwj

楼主为什么不回复

sovereign

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

复制代码

要要切克闹！！！！