www.pudn.com > hilbert.zip > test_e.m, change:2009-07-09,size:403b


%要先修改disp_hhs.m中的采样率,这里为1000 
clc;clear; 
t=0:0.001:1; 
y=cos(300*pi*t)+2*cos(100*pi*t+0.5*sin(20*pi*t));%信号1 
 
x=y'; 
m=x; 
t1=t; 
 
imf=emd_lily01(x); %在这里修改不同的端点处理方法(emd_lily01、emd_lily02、emd_lily04、emd_lily和emd_lily_NB) 
 
emd_visu(x,t,imf); 
[A,fa,tt]=hhspectrum(imf); 
[E,tt1]=toimage(A,fa); 
disp_hhs(E); 
colormap(flipud(gray))            %显示Hilbert谱