selezione = 3;
switch selezione
case 1,
  ytot = wavread('intro_KeithJarrett_Country');
  istart = find(ytot(1:85000)>eps); % ci possono essere fino a due secondi di silenzio all'inizio della traccia
  N = 65536;  % numero di campioni del segnale campionato a 44.1 KHz
  a = 128; % passo di traslazione nel tempo
  b = 16; % passo di decimazione in frequenza
case 2,
  ytot = wavread('pianoforte_4'); % note "legate"
  istart = 46000;
  N = 65536*4;  % numero di campioni del segnale campionato a 44.1 KHz
  a = 128*4; % passo di traslazione nel tempo
  b = 32*2; % passo di decimazione in frequenza
case 3,
  ytot = wavread('pianoforte_3'); % note "staccate":  "do" a 262 Hz, "mi" a 330 Hz, "sol" a 396 Hz
  istart = 50000;
  N = 65536*8;  % numero di campioni del segnale campionato a 44.1 KHz
  a = 128*8; % passo di traslazione nel tempo
  b = 32*2; % passo di decimazione in frequenza
case 4,
  ytot = wavread('pianoforte_3'); % note "staccate", solo "do" a 262 Hz
  istart = 50000;
  N = 65536;  % numero di campioni del segnale campionato a 44.1 KHz
  a = 128*2; % passo di traslazione nel tempo
  b = 32; % passo di decimazione in frequenza
otherwise
  error('selezione non presente!');  
end
y = ytot(istart(1)+(0:N-1),:);
%wavplay(y, 44100);
y = y(:,1);
clear istart;
%figure(1); plot((1:N)/44100,y);


sigma = 512;

if(lcm(N,a) ~=1  & lcm(N,b)~=1 )%devo controllare che a e b dividano N
  c=round(N/2); %la gaussiana che costruisco ha centro c nell'intero più vicino a N/2
  g=gaussmf([0:N-1],[sigma c]); % costruzione della gaussiana con una funzione implementata da matlab.
                                % è una gaussiana di centro c e deviazione standard sigma che  è
                                % calcolata  nei punti [0:N-1];
                                % è infatti ipotizzato che la sequenza di dati lunga N abbia gli indici
                                % che corrono da 0 a N-1
  C = stft(y',g,a,b);
  spettrog_musica(C,N,a,b);
else
  print('Errore: i passi di campionamento non dividono N.')
end    





return

