% figure with wavelet in TF space
%
% want to print out wavelet on colour printer to see how we can do inverse
%        TF transform; we want bivariate Gaussian with correct TF phase

fc = 20;     %  fcent of wavelet
dt = .25;    % "full support" of wavelet as defined by daugman

g = gabor(fc,dt,0,0,1024);    % gabor wavelet
Y = tf(g,400,16,25);   % redundant TF "Gabor" transform
contour(abs(Y));
title('FREQUENCY')
ylabel('TIME')

mat2cps(Y)

