Reassigned Morlet Scalogram time-frequency distribution
[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')
A Nx elements vector or a Nx by 2 array signal.
the time instant(s) with elements in [1 Nx] (default : 1:length(X))
number of frequency bins (default : length(X)). For faster computation N should be a power of 2.
a positive scalar: the time-bandwidth product of the mother wavelet (default : 2.5))
A boolean (or a real scalar) if true (or nonzero),the progression of the algorithm is shown (default : %f).
if one input parameter is 'plot', tfrqview is called and the time-frequency representation will be plotted.
a N by Nt array: the time-frequency representation.
A real N by Nt array: the reassigned time-frequency representation.
a N by Nt complex matrix: the reassignment vectors.
tfrrmsc computes the Morlet scalogram and its reassigned version.
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 | ![]() | ![]() |