function stft(a1)

% ********** Funci髇 de Matlab para hallar la STFT **********************
% Halla la transformada de Forurier de Corto Tiempo
% Entrada: a1 es el coeficiente de dilatacion de la ventana exploradora
% Funcion a transformar
n=10; % exponente del numero total de muestras: 2^n 1024(casi las 1280)
f1=0.5; f2=40; % frecuencias de la se馻l entr 0,5 y 40Hz
for i=1:2^n,
if i<2^(n-1), V(i)=sin(2*pi*f1*(i-1)/2^(n)); else V(i)=sin(2*pi*f2*(i-1)/2^(n));end
end
Nw=32; % ancho de la ventana
i=1:Nw;w(i)=exp(-a1/2*((i-1)-Nw/2).^2); % vector ventana de Nw elementos
for i=0:2^n-1,
VV(i+1)=V(i+1)*w(mod(i,Nw)+1);
end
% el vector VV tiene los productos de la ventana por cada tramo del vector
% datos
i=0; while i<2^n,
j=0; while j<Nw
A(floor(i/Nw)+1, j+1)=VV(i+1);i=i+1;j=j+1;
end
end
for i=1:2^n/Nw,
for j=1:Nw, B(j)=A(i,j); end
C=abs(fft(B))/sqrt(Nw);
for j=1:Nw, F(i,j)=C(j); end
end
surf(F)