function[]=spettrog_musica(C,N,a,b)
%funzione che ricostruisce lo spettrogramma


    atilde=N/a;%quindi N=a*atilde
    btilde=N/b; %quindi N=b*btilde
    M=atilde*btilde;
    
    %C = abs(C(1:floor(btilde/2),:));
%    C = abs(C(1:floor(btilde/2/8),:));
    
%    C=zeros(btilde, atilde);
%    C=reshape(C1,btilde, atilde);
%    Cplot=C(1:round(btilde/2),:);%basta plottare la matrice troncata a metà nella frequanza, infatti per ricostruire 
    %lo spettrogramma basta arrivare alla ferquenza di Nyquist
    minC = min(min(abs(C)));
    maxC = max(max(abs(C)));
    %
    delta_f = 44100/N*b % intervallo in Hz tra una riga e la successiva
    div_b = min([2^(log2(btilde)-8) 16]);
    figure(11), surf(abs(C(1:floor(btilde/2/div_b),:))), caxis([minC maxC/4]), axis([0 size(C,2) 0 size(C,1)/2/div_b minC maxC]) 
                view(2) ,ylabel('frequenza'), xlabel('tempo'),title ('spettrogramma low freq');
    %
    figure(12), surf(abs(C(floor(btilde/2/div_b)+1:2*floor(btilde/2/div_b),:))), caxis([minC maxC/4]), axis([0 size(C,2) 0 size(C,1)/2/div_b minC maxC]) 
                view(2) ,ylabel('frequenza'), xlabel('tempo'),title ('spettrogramma low-middle freq');
    %
    figure(13), surf(abs(C(2*floor(btilde/2/div_b)+1:3*floor(btilde/2/div_b),:))), caxis([minC maxC/4]), axis([0 size(C,2) 0 size(C,1)/2/div_b minC maxC]) 
                view(2) ,ylabel('frequenza'), xlabel('tempo'),title ('spettrogramma middle-low freq');
    %
    figure(14), surf(abs(C(3*floor(btilde/2/div_b)+1:4*floor(btilde/2/div_b),:))), caxis([minC maxC/4]), axis([0 size(C,2) 0 size(C,1)/2/div_b minC maxC]) 
                view(2) ,ylabel('frequenza'), xlabel('tempo'),title ('spettrogramma middle freq');
    %
    