% S_chirprate.m
%
% P. Flandrin

clf

N = 256;
f0 = .0;
ff = .5;
sig = fmlin(N,f0,ff);

K = 12;

Nbin = 128;

Lmin = 3;
Lmax = N-1;
Lh = odd(logspace(log10(Lmin),log10(Lmax),K));

pause

for k = 1:K
	h = window(Lh(k),'Kaiser');
	[s,T,F] = tfrsp(sig,1:N,2*Nbin,h);
	imagesc(flipud(s(1:Nbin,N/4:3*N/4)));
	colormap(jet)
	set(gca,'YTick',[]);set(gca,'XTick',[])
	title(['spectrogramme, Lh = ',int2str(Lh(k))])
	pause(.1)
end
