ZADANIE GOTOWE DO ODDANIA

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

disp(['Czestotliwosc probkowania fs=',num2str(fs),'[Hz]']);
N=length(y); %ilość próbek


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

    figure(2)
    subplot(211); %
    %rozkład prawdopodobieństwa
nbits = 100; %wykorzystuje bitów na próbkę sygnału audio y
hist(y,nbits);% histogram
title('rozklad prawdopodobiensta')
xlabel('wartosc probki');
ylabel('liczba probek w przedziale');
    %wykres sygnalu
   figure(1)
    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
    f=0:df:fs-df;
    subplot(212)
    plot(f,widmo); grid on
    title ('Widmo amplitudowe')
    xlabel('Czestotliwosc[Hz]')

    x = randn(1,N);% rozkład normalny liczb losowych

max = 1000;
ry = xcorr(y,x,max); % opóźnienie sygnału w sekundach - korelacja wzajemna z sygnałem
r = -max/fs:1/fs:max/fs;
figure(3)
subplot(211);
plot(r, ry);
xlabel('przesunicie [s]');
ylabel('autokorelacja');