% Setup the basic parameters of the Gabor system.L=480;  % Length of signal.a=40;   % Time shift.M=12;   % Number of modulations.% Number of Hermite functions to useR=4;% Create a multiwindow Gabor window using% the R first Hermite functions.g=pherm(L,0:R-1);% Create a Gabor frame object.gf=tfr_create('dgt',a,M,g,'candual');% Alternatively, the following call would have created% the same object, in a more flexible manner.%% gf=tfr_create('dgt',a,M,{'pherm',0:R-1},'candual');% Test that analysis and synthesis works.f=rand(L,1);c=tfr_ar(gf,f);fr=tfr_sr(gf,c);disp('');disp('Norm of reconstruction, should be close to zero.');norm(f-fr)% Get the analysis and synthesis windowsg  = tfr_get_win(gf,'a',L);gd = tfr_get_win(gf,'s',L);% Clear the object again.tfr_clear(gf);% Plot some spetrograms of the Hermite functions.figure(1);subplot(2,2,1);sgram(g(:,1),'nf','tc');zoom(2);subplot(2,2,2);sgram(g(:,2),'nf','tc');zoom(2);subplot(2,2,3);sgram(g(:,3),'nf','tc');zoom(2);subplot(2,2,4);sgram(g(:,4),'nf','tc');zoom(2);% Plot some spetrograms of the dual of the % Hermite functions.figure(2);subplot(2,2,1);sgram(gd(:,1),'nf','tc');zoom(2);subplot(2,2,2);sgram(gd(:,2),'nf','tc');zoom(2);subplot(2,2,3);sgram(gd(:,3),'nf','tc');zoom(2);subplot(2,2,4);sgram(gd(:,4),'nf','tc');zoom(2);