Name

hhspectrum — compute Hilbert-Huang spectrum

Calling Sequence

   [A,f,tt] = hhspectrum(x)
   [A,f,tt] = hhspectrum([],x)
   [A,f,tt] = hhspectrum(t,x)
   [A,f,tt] = hhspectrum(x,l)
   [A,f,tt] = hhspectrum([],x,l)
   [A,f,tt] = hhspectrum(t,x,l)
   [A,f,tt] = hhspectrum(x,l,aff)
   [A,f,tt] = hhspectrum([],x,aff)
   [A,f,tt] = hhspectrum(t,x,l,aff)
   

Parameters

inputs:

- x :

matrix with one signal per row

- t :

time instants

- l :

estimation parameter for instfreq (integer >=1 (1:default))

- aff :

if 1, displays the computation evolution

outputs:

- A :

instantaneous amplitudes

- f :

instantaneous frequencies

- tt :

truncated time instants

Description

need the Time-Frequency Toolbox (stftb)

Examples


s = rand(1,512,'normal');
imf = emd(s);
[A,f,tt] = hhspectrum(imf(1:$-1,:));
[im,tt]=toimage(A,f);
disp_hhs(im);

s = rand(10,512,'normal');
[A,f,tt] = hhspectrum([],s,2,1);
[im,tt]=toimage(A,f);
disp_hhs(im);

   

See also

emd, toimage, disp_hhs

Authors

H. Nahrstaedt - Aug 2010
G. Rilling, last modification 3.2007
gabriel.rilling@ens-lyon.fr