%
% function tfldemo 
%
% ejemplo de analisis tiempo-frecuencia
% la seņal se compone de una chirp y una sinusoide
% la seņal chirp tiene una duracion de 1000 muestras
% y la sinusoide tambien 1000 muestras separadas ambas
% 250 muestras.
%
% [x,S]=tfldemo(Lven,Dven,NFFT,ven)
%
% Lven ........  longitud de la ventana de analisis
% Dven ........  desplazamiento en muestras de la ventana de analisis
% NFFT ........  longitud de la fft
% ven  ........  tipo de ventana, 0 -rectangular, 1- hamming
%
% x    ........  seņal temporal
% S    ........  TFL en dB
%

function [x,S]=tfldemo(Lven,Dven,NFFT,ven);

if (nargin<4)
 ven=0;
end

if(NFFT<Lven)
 disp('!!! orden de la FFT inferior a la longitud de la ventana');
end

if (ven==0)
 ventana='rectangular';
else
 ventana='hamming';
end

f1=0.0002;
f0=0.15;
f2=0.2;
%NFFT=256;
%Lven=100;

n1=[0:999];
x=zeros(2750,1);
x(1:1000)=sin(pi*f1*n1.^2+2*pi*f0*n1);
x(1250:2249)=sin(2*pi*f2*n1);
S=tfl(x,Lven,Dven,NFFT,ven);
tam=size(S);
imagesc([0:Dven*tam(2)],[0:0.5/(NFFT-1):0.5],S,[-100,100]);
colormap(flipud(gray));
set(gca,'YDir','normal');
set(gca,'Ygrid','on');
set(gca,'Xgrid','on');
xlabel('muestras');
ylabel('frecuencia');
title(['ventana de longitud: ', num2str(Lven), ' tipo ',ventana]);
colorbar('vert');