hht程序，帮忙翻译一下，一定要具体详细啊。新手菜鸟

q296844597

1. clear;
   
2. N=2048;
   
3. t=1:N;
   
4. fs=1000;
   
5. z=sin(2*pi*40*t/fs)+sin(2*pi*15*t/fs)+sin(2*pi*5*t/fs);
   
6. plot(t,z);                                     %做频谱图
   
7. xlabel('时间');
   
8. ylabel('幅值');
   
9. title('原始信号');
   
10. imf=emd(z);                           %EDM分解
    
11. cemd_visu(z,1:length(z),imf)
    
12. [A, fa, tt] = hhspectrum(imf);                  % 对IMF分量求取瞬时频率与振幅
    
13. [E, tt1] = toimage(A,fa,tt,length(tt));    %绘图toimage是把一维的变成二维的显示出来
    
14. disp_hhs(E);                                % disp_hhs(E,[],fs);是非归一化频率
    
15. colormap(flipud(gray))
    
16. for k=1:size(E,1);
    
17.     bjp(k) = sum(E(k,:))*1/fs;  
    
18. end;
    
19. f = (0:N-3)/N*(fs/2);         %进行对应的频率转换，定义横坐标频率范围的
    
20. figure(4)                 %创建画图窗口
    
21. plot(f,bjp);            %做频谱图
    
22. xlabel('频率');
    
23. ylabel('幅值');
    
24. title('边际谱');

复制代码

q296844597

包括小参数，也说明一下，感激不尽。

竹叶清风123

具体的函数可以上网查一下用法，理解的会更深刻。

江海客

这不是有工具箱的吗

cufflink

程序不完整，有些参数不好写。