<< Reassigned Time-Frequency Processing Reassigned Time-Frequency Processing tfrrmsc >>

Time Frequency Toolbox >> Time Frequency Toolbox > Reassigned Time-Frequency Processing > tfrrgab

tfrrgab

Reassigned Gabor spectrogram time-frequency distribution.

Calling Sequence

[TFR,RTFR,HAT] = tfrrgab(X)
[TFR,RTFR,HAT] = tfrrgab(X,T)
[TFR,RTFR,HAT] = tfrrgab(X,T,N)
[TFR,RTFR,HAT] = tfrrgab(X,T,N,NH)
[TFR,RTFR,HAT] = tfrrgab(X,T,N,NH,TRACE)
[TFR,RTFR,HAT] = tfrrgab(X,T,N,NH,TRACE,K)
[TFR,RTFR,HAT] = tfrrgab(...,'plot')

Parameters

X :

A Nx elements vector: the signal.

T:

a real Nt vector with elements in [1 Nx] : time instant(s) (default: 1:NX).

N:

a positive integer: the number of frequency bins (default:NX). For faster computation N should be a power of 2.

NH :

an odd positive integer: the length of the gaussian window (default : the first odd number greater than N/4))

TRACE :

A boolean (or a real scalar) if true (or nonzero),the progression of the algorithm is shown (default : %f).

K :

a real sclar: the value at both extremities (default 0.001)

'plot':

if last input parameter value is 'plot', tfrqview is called and the time-frequency representation will be plotted.

TFR :

A real N by Nt array: the time-frequency representation.

RTFR :

A real N by Nt array: the reassigned time-frequency representation.

HAT :

>A complex N by Nt array: the reassignment vectors.

Description

tfrrgab computes the Gabor spectrogram and its reassigned version. This particular window (a Gaussian window) allows a 20 % faster algorithm than the tfrrsp function.

Examples

N=128;
sig=fmlin(N,0.1,0.4);
t=1:128;
Nh=19;
[tfr,rtfr,hat]=tfrrgab(sig,t,N,Nh);
clf;gcf().color_map= jetcolormap(128);
subplot(121)
grayplot(t,linspace(0,0.25,N/2),tfr(1:N/2,:)')
title("TFR")
subplot(122)
grayplot(t,linspace(0,0.25,N/2),tfr(1:N/2,:)')
title("RTFR")

Authors


Report an issue
<< Reassigned Time-Frequency Processing Reassigned Time-Frequency Processing tfrrmsc >>