时频Gabor变换

function [tfr,dgr,gam]=tfrgabor(sig,N,q,h,trace) %TFRGABOR Gabor representation of a signal. % [TFR,DGR,GAM]=TFRGABOR(SIG,N,Q,H,TRACE) computes the Gabor % representation of signal X, for a given synthesis window H, on a % rectangular grid of size (N,M) in the time-frequency plane. M and N % must be such that % N1 = M * N / Q % where N1=length(X) and Q is an integer corresponding to the % degree of oversampling. % % SIG : signal to be analyzed (length(SIG)=N1). % N : number of Gabor coefficients in time (N1 must be a multiple % of N) (default : divider(N1)). % Q : degree of oversampling ; must be a divider of N % (default : Q=divider(N)). % H : synthesis window, which was originally chosen as a Gaussian % window by Gabor. Length(H) should be as closed as possible % from N, and must be >=N (default : Gauss(N+1)). % H must be of unit energy, and CENTERED. % TRACE : if nonzero, the progression of the algorithm is shown % (default : 0). % TFR : Square modulus of the Gabor coefficients. When % called without output arguments, TFRGABOR runs TFRQVIEW. % DGR : Gabor coefficients (complex values). % GAM : biorthogonal (dual frame) window associated to H.