Reassigned Gabor spectrogram time-frequency distribution
[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')
A Nx elements vector: the signal.
a real Nt vector with elements in [1 Nx] : time instant(s) (default: 1:NX).
a positive integer: the number of frequency bins (default:NX). For faster computation N should be a power of 2.
an odd positive integer: the length of the gaussian window (default : the first odd number greater than N/4))
A boolean (or a real scalar) if true (or nonzero),the progression of the algorithm is shown (default : %f).
a real sclar: the value at both extremities (default 0.001)
if last input parameter value is 'plot', tfrqview is called and the time-frequency representation will be plotted.
A real N by Nt array: the time-frequency representation.
A real N by Nt array: the reassigned time-frequency representation.
>A complex N by Nt array: the reassignment vectors.
tfrrgab computes the Gabor spectrogram and its reassigned version. This particular window (a Gaussian window) allows a 20 % faster algorithm than the tfrrsp function.
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 | ![]() | ![]() |