GABREASSIGNREAL - Reassign time-frequency distribution for real signals

Usage

sr = gabreassignreal(s,tgrad,fgrad,a,M);

Description

gabreassignreal(s,tgrad,fgrad,a,M) reassigns the values of the positive time-frequency distribution s using the phase gradient given by fgrad and tgrad. The lattice is determined by the time shift a and the number of channels M.

fgrad and tgrad can be obtained by the routine gabphasegrad.

Examples:

The following example demonstrates how to manually create a reassigned spectrogram. An easier way is to just call resgram:

% Create reassigned vector field of the bat signal.
a=4; M=100;
[phased,c] = gabphasederivreal({'t','f'},'dgt',bat,'gauss',a,M,'relative');
[tgrad, fgrad] = deal(phased{:});

% Perform the actual reassignment
sr = gabreassignreal(abs(c).^2,tgrad,fgrad,a,M);

% Display it using plotdgt
plotdgt(sr,a,143000,50);