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

stftb >> stftb > 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 >>