<< tfrrgab Reassigned Time-Frequency Processing tfrrpmh >>

Time Frequency Toolbox >> Time Frequency Toolbox > Reassigned Time-Frequency Processing > tfrrmsc

tfrrmsc

Reassigned Morlet Scalogram time-frequency distribution.

Calling Sequence

[TFR,RTFR,HAT] = tfrrmsc(X)
[TFR,RTFR,HAT] = tfrrmsc(X,T)
[TFR,RTFR,HAT] = tfrrmsc(X,T,N)
[TFR,RTFR,HAT] = tfrrmsc(X,T,N,F0T)
[TFR,RTFR,HAT] = tfrrmsc(X,T,N,F0T,TRACE)
[TFR,RTFR,HAT] = tfrrmsc(...,'plot')

Parameters

X :

A Nx elements vector or a Nx by 2 array signal.

T :

the time instant(s) with elements in [1 Nx] (default : 1:length(X))

N :

number of frequency bins (default : length(X)). For faster computation N should be a power of 2.

F0T :

a positive scalar: the time-bandwidth product of the mother wavelet (default : 2.5))

TRACE :

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

'plot':

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

TFR :

a N by Nt array: the time-frequency representation.

RTFR :

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

HAT :

a N by Nt complex matrix: the reassignment vectors.

Description

tfrrmsc computes the Morlet scalogram and its reassigned version.

Examples

N=128;
sig=fmlin(N,0.1,0.4);
t=1:128;

[tfr,rtfr,hat]=tfrrmsc(sig,t,N,2.1);
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),rtfr(1:N/2,:)')
title("RTFR")

Authors


Report an issue
<< tfrrgab Reassigned Time-Frequency Processing tfrrpmh >>