%_________________________________________
%  fenetre.m
%  PROGRAMME DE TRANSFORMEE DE FOURIER A FENETRE
%  Valerie Perrier - LMD novembre 1993
%
%  la fenetre est une gaussienne
%  Le signal a analyser doit s'appeler Y
%_________________________________________
%
clf
ii=sqrt(-1);
%
%..signal a analyser........
%
n=length(y)-1;
n2=n/2;
dx=1/n;
t=0:dx:1;
x=-1:dx:1;
%
% calcul de la fenetre (Gaussienne placee en b=i.b0, de taille a)
%
a=input('taille de la fenetre < 1 :');
b0=a/4;
l=fix(4/a);
W=zeros(n2,n);
for i=1:l
	b=i*b0;
	xx=(x-b)./a;
 	g=exp(-xx.^2/2);%Gaussienne
%
% calcul des coefficients de Fourier pour y.g(x-b)
%
	gy(1:n)=y(1:n).*g(n+1:2*n) ;
	fy=fft(gy,n); 
	W(1:n2,i)=abs(fy(1:n2))';
end
%
% trace
%
colormap(gray(60));
subplot(3,1,1),plot(t(1:n),y(1:n)),title('fonction a analyser')
subplot(3,1,2),pcolor(W(1:n2,1:l)) , title('coefficients Fourier Fenetre')
subplot(3,1,3),contour(W(1:n2,1:l)),title('lignes isocontour')