Introduction - If you have any usage issues, please Google them yourself
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)