%STFT de una señal no estacionaria con contenido de 10,25,50 y 100 Hz y varios spikes
%Se emplea una ventana del tipo gaussiana con diferentes valores de 
%amplitud a
%
%Subrutinas empleadas: 
%          - WindowFT
%
%Nota: el parámetro a se introduce en la subrutina Makewindow, 
%que es parte de la subrutina WindowFT

close all; 
clear all;

spike1(1:307)=0;
spike2(1:308)=0;
spike3(1:205)=0;
npuntos=1024;
w=2*pi;
t=0:1/npuntos:1;
t1=t(1:npuntos*0.3);
t2=t(npuntos*0.3:npuntos*0.6);
t3=t(npuntos*0.6:npuntos*0.8);
t4=t(npuntos*0.8:npuntos);
tiempo=1000*t;
spike1(201)= 2;
spike1(251)=1.5;
spike2(15)=-0.8;
spike2(292)=1;
spike3(105)=1.5;
x1=sin(w*100*t1) + spike1;
x2=sin(w*50*t2) + spike2;
x3=sin(w*25*t3);
x4=sin(w*10*t4)+spike3;
x=[x1 x2 x3 x4];
specgm = WindowFT(x,npuntos/4,50,'Gaussian');
z=abs(specgm);
mesh(z);


