时频Gabor变换

时频Gabor变换
function[tfr,dgr,gam]=tfrgabor(sig,N,q,h,trace)%TFRGABORGaborrepresentationofasignal.% [TFR,DGR,GAM]=TFRGABOR(SIG,N,Q,H,TRACE)computestheGabor% representationofsignalX,foragivensynthesiswindowH,ona% rectangulargridofsize(N,M)inthetime-frequencyplane.MandN% mustbesuchthat% N1=M*N/Q% whereN1=length(X)andQisanintegercorrespondingtothe% degreeofoversampling.%% SIG:signaltobeanalyzed(length(SIG)=N1).% N:numberofGaborcoefficientsintime(N1mustbeamultiple% ofN) (default:divider(N1)).% Q:degreeofoversampling;mustbeadividerofN% (default:Q=divider(N)).% H:synthesiswindow,whichwasoriginallychosenasaGaussian% windowbyGabor.Length(H)shouldbeasclosedaspossible% fromN,andmustbe＞=N(default:Gauss(N+1)).% Hmustbeofunitenergy,andCENTERED.% TRACE:ifnonzero,theprogressionofthealgorithmisshown% (default:0).% TFR:SquaremodulusoftheGaborcoefficients.When% calledwithoutoutputarguments,TFRGABORrunsTFRQVIEW.% DGR:Gaborcoefficients(complexvalues).% GAM:biorthogonal(dualframe)windowassociatedtoH.