EMD分解并绘制希尔伯特时频谱

对信号进行EMD分解并绘制希尔伯特时频谱的程序

% hilbert谱分析

clear;
clc;

t=1:200;
fs=100;
t=t*1/fs;
% x=sin（2*pi*20*t）;
% y=0.4*sin（2*pi*40*t+140）;
%z=cos（2*pi*5*t）+0.5*cos（2*pi*10*t）+0.7*cos（2*pi*25*t）;
% z=exp（-0.3*t）.*cos（pi*t.*t+2*pi*t）;
%z=0.3*exp（j*8*t）+1*exp（j*16*t）;
[fnamepname]=uigetfile（‘0a.wav‘‘Open Wave File‘）;
file=[pnamefname];
[zfsbits]=wavread（file）;
%N=length（y）;
%z= awgn（y50）;



imf=emd（z）;
 emd_visu（z1:length（z）imf）

[Aftt]=hhspectrum（imf）;
[imtt]=toimage（Af）;
disp_hhs（im[]fs）;  % 绘制hilbert时频谱
      colormap（flipud（gray））;%这条语句加上注释改变图形
title（‘hilbert时频谱‘）;
colorbar;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        618  2009-05-08 16:12  hsa.m

----------- ---------  ---------- -----  ----

                  618                    1