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") | ![]() | ![]() |