%  show gabors of long duration (with fourier trans)

clg
axis('normal') % anisotropic
s = gabor(20,.25/8,0,0,100);
[s,t] = pad(s,150);
subplot(221)
axis([-.5 .5 -.7 1.2])
plot(t,s)
grid
%%%plot(t,s)  % need to redo to get border and axes
[s,f] = pad(s,1024);
S=fftshift(abs(fft(s)));
subplot(223)
axis([-.5 .5 0 4.99])
plot(f,S)  % plot previous fft in third place (under prev. time plot)
grid

hold on
subplot(211)
title('an expansion onto short Gabor functions gives good physical resolution but ...')
hold off

s = gabor(2,.25/8,0,0,100);
[s,t] = pad(s,150);
subplot(222)
axis([-.5 .5 -.7 1.2])
plot(t,s)
grid

[s,f] = pad(s,1024);
S=fftshift(abs(fft(s)));
subplot(224)
axis([-.5 .5 0 4.99])
plot(f,S)
grid

subplot(212)
title('poor Fourier resolution')

hold off

