GABPHASEDERIVREAL - DGT phase derivatives

Usage

[phased,c] = gabphasederivreal(dflag,'dgt',f,g,a,M);
phased    = gabphasederivreal(dflag,'cross',f,g,a,M)
phased    = gabphasederivreal(dflag,'phase',cphase,a,M);
phased    = gabphasederivreal(dflag,'phase',cphase,a,M,difforder);
phased    = gabphasederivreal(dflag,'abs',s,g,a,M);
phased    = gabphasederivreal(dflag,'abs',s,a,M,difforder);
[{phased1,phased2,...}] = gabphasederivreal({dflag1,dflag2,...},...);
[{phased1,phased2,...},c] = gabphasederivreal({dflag1,dflag2,...},'dgt',...);

Description

phased=gabphasederivreal(dflag,method,...) computes the time-frequency derivative dflag of the phase of the dgtreal of a signal using algorithm method.

The following strings can be used in place of dflag:

phased is scaled such that (possibly non-integer) distances are measured in samples. Similarly, the frequencies are scaled such that the Nyquist frequency (the highest possible frequency) corresponds to a value of L/2.

The computation of phased is inaccurate when the absolute value of the Gabor coefficients is low. This is due to the fact the the phase of complex numbers close to the machine precision is almost random. Therefore, phased attain very large random values when abs(c) is close to zero.

The phase derivative computation can be done using four different methods depending on the string method:

phased=gabphasederivreal(dflag,'dgt',f,g,a,M) computes the time-frequency derivative using a DGT of the signal f. The DGT is computed using the window g on the lattice specified by the time shift a and the number of channels M. The algorithm used to perform this calculation computes several DGTs, and therefore this routine takes the exact same input parameters as dgtreal.

[phased,c]=gabphasederivreal(dflag,'dgt',f,g,a,M) additionally returns the Gabor coefficients c, as they are always computed as a byproduct of the algorithm.

phased=gabphasederivreal(dflag,'cross',f,g,a,M) does the same as above but this time using algorithm by Nelson which is based on computing several DGTs.

phased=gabphasederivreal(dflag,'phase',cphase,a,M) computes the phase derivative from the phase cphase of a DGT of the signal. The number of channels M is a required input to disambiguate the existence of a Nyquist channel, which is only present if M is even. The original DGT from which the phase is obtained must have been computed using a time-shift of a using the default phase convention ('freqinv') e.g.:

phased=gabphasederivreal(dflag,'phase',angle(dgt(f,g,a,M)),a,M)

phased=gabphasederivreal(dflag,'abs',s,g,a) computes the phase derivative from the absolute values of DGT coefficients s. The number of channels M is a required input to disambiguate the existence of a Nyquist channel, which is only present if M is even. The spectrogram must have been computed using the window g and time-shift a e.g.:

phased=gabphasederivreal(dflag,'abs',abs(dgt(f,g,a,M)),g,a,M)

Currently the 'abs' method only works if the window g is a Gaussian window specified as a string or a cell array.

phased=gabphasederivreal(dflag,'abs',s,g,a,M,difforder) and phased=gabphasederivreal(dflag,'pgase',cphase,a,M,difforder) use a centered finite diffence scheme of order difforder to perform the needed numerical differentiation. Default is to use a 4th order scheme.

Phase conventions

First derivatives in either direction are subject to phase convention. The following additional flags define the phase convention the original phase would have had:

Please see ltfatnote042 for the description of relations between the phase derivatives with different phase conventions. Note that for the 'relative' convention, the following holds:

gabphasederiv('t',...,'relative') == gabphasederiv('t',...,'freqinv')
gabphasederiv('f',...,'relative') == -gabphasederiv('f',...,'timeinv')
gabphasederiv('tt',...,'relative') == gabphasederiv('tt',...)
gabphasederiv('ff',...,'relative') == -gabphasederiv('ff',...)
gabphasederiv('tf',...,'relative') == gabphasederiv('tf',...,'freqinv')

Several derivatives at once

phasedcell=gabphasederiv({dflag1,dflag2,...},...) computes several phase derivatives at once while reusing some temporary computations thus saving computation time. {dflag1,dflag2,...} is a cell array of the derivative flags and cell elements of the returned phasedcell contain the corresponding derivatives i.e.:

[pderiv1,pderiv2,...] = deal(phasedcell{:});

[phasedcell,c]=gabphasederiv({dflag1,dflag2,...},'dgt',...) works the same as above but in addition returns coefficients c which are the byproduct of the 'dgt' method.

Other flags and parameters work as before.