audioread

%Kamil Janik 226973
clear all;
[y,fs]=audioread('piatek_n.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


%rysowanie sygnalu z pliku

%sygnal przypomina ksztaltem sinusoide z szumem bialym(B)
   figure()
   subplot(211)
    plot(t,y); grid on
    title('Sygnal')
    xlabel('Czas[s]')
    ylabel('Amplituda')

    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]')

    figure()
    subplot(211); %
    %rysowanie hisogramu rozkladu prawdopodobienstwa(A)
nbits = 50; %wykorzystuje bitów na próbkę sygnału audio y
hist(y,nbits);
title('rozklad prawdopodobiensta')
xlabel('wartosc probki');
ylabel('liczba probek w przedziale');

% opoznienie sygnalu i korelacja(E)
  %Sygnal jest sinusoida + szum gausowski
% Czestotliwosc=420 Odczytano z wykresy widma amplitudowego

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()
subplot(211);
plot(r, ry);
xlabel('przesunicie [s]');
ylabel('autokorelacja');