analiza wav

clear all
%subplot(111);
[y,fs,bits]=wavread('mbi04popis.wav'); % odczytanie pliku wav

disp(['Czestotliwosc probkowania fs=',num2str(fs),'[Hz]']);
N=length(y);

    dt=1/fs; % okres probkowania
    t=0:dt:(N-1)*dt; % wektor czasu

    %wykres sygnalu
    subplot(211)
    plot(t,y); grid on
    title('Sygnal')
    xlabel('Czas[s]')

    df=fs/N % rozdzielczosc widma w Hz
    widmo=abs(fft(y)); % widmo amplitudowe (modul(
    f=0:df:fs-df;
    subplot(212)
    plot(f,widmo); grid on
    title ('Widmo amplitudowe')
    xlabel('Czestotliwosc[Hz]')