% Wczytanie sygnalu w formacie WAV, FUNKCJA DECIMATE
clear;
% wczytanie sygnału w formacie WAV
[x,fpx] = wavread('mbi04czap.wav');
Nx=length(x)
tx=0:1/fpx:(Nx-1)/fpx;

% wykres sygnału przed decymacją

subplot(221)
plot(tx,x)
xlabel('czas[s]')
ylabel('sygnal')
set(gcf,'Position',[10,100,1000,550]) % ustalenie wielkości okna
set(gcf,'PaperPositionMode','auto') % ustalenie wielkości okna

% widmo sygnału przed decymacją
Nfx = 2^15; % Nf powinno być większe niż N
N21x = Nfx/2+1;
v = fft(x,Nfx); % widmo zepsolone
wx = abs(v); % widmo amplitudowa
fx = linspace(0,fpx/2,N21x); % generowanie osi czestotliwosci
subplot(222)
plot(fx,wx(1:N21x))
xlabel('czestotliwosc[Hz]')
ylabel('modul widma')

% decymacja wykres

dr = 8 % rzad decymacji
y = decimate(x,dr);
Ny = length(y);
fpy = fpx/dr;
ty = 0:1/fpy:(Ny-1)/fpy;
subplot(223)
plot(ty,y)
xlabel('czas[s]')
ylabel('sygnal po decymacji')

% widmo sygnału decymacji

Nfy = 2^13
N21y = Nfy/2+1;
v = fft(y,Nfy);
wy = abs(v);
fy = linspace(0,fpy/2,N21y);
subplot(224)
plot(fy,wy(1:N21y));
xlabel('czestotliwosc[Hz]')
ylabel('modul widma sygn. po decymacji')