kokos

clear all;

figure(1)

[x, fpx] = wavread('aut_60_wr.wav');
Nx = length(x);
tx = 0:1/fpx:(Nx-1)/fpx;

subplot(331);
plot(tx,x);
xlabel('czas [s]');
title('wczytany sygnal');
ylabel('amplituda')

subplot(332);
nbits = 51;
hist(x,nbits);
title('rozklad prawdopodobiensta')
xlabel('wartosc probki');
ylabel('liczba prbek w przedziale');

Nfx = 2^15;
N21x = Nfx/2+1;
v = fft (x,Nfx);                   % widmo zespolone
wx = abs(v);                       % widmo amplitudowe
fx = linspace (0, fpx/2, N21x);        % generowanie osi czestotliwosci
subplot(333);
plot(fx,wx(1:N21x));
title('widmo aplitudowe sygnalu')
xlabel('czestotliwosc[Hz]');
ylabel('modul amplitudy');

maxfx = 0;
maxwx = 0;
for i = (1:N21x)
  if maxwx < wx(i)
    maxwx = wx(i);
    maxfx = fx(i);
  end
end

display('Czestotliwosc podstawowa: [Hz]')
display(maxfx)

xn = randn(1,Nx);

kmax = 1000;
ry = xcorr(x,xn,kmax);
tr = -kmax/fpx:1/fpx:kmax/fpx;
subplot(312);
plot(tr, ry);
xlabel('przesunicie [s]');
ylabel('autokorelacja');