说明：  Hilbert变换只是单纯地求信号的瞬时振幅，频率和相位，有可能出现没有意义的负频率；HHT变换先将信号进行EMD分解，得到的是各个不同尺度的分量，对每一个分量进行Hilbert变换后得到的是有实际意义的瞬时频率。
(imagesc(t,fliplr((1:size(im,1))/(2*size(im,1))),im,[inf,0]) set(gca, YDir , normal ) xlabel([ time ]) ylabel([ normalized frequency ]) title( Hilbert-Huang spectrum ) function [f,z]=fftfenxi(t,y) L=length(t) N=2^nextpow2(L) )