clear;
[x,fs]=wavread('mbi02becz.wav');
dr=6;
y=decimate(x,dr);
fp=fs/dr;
N=length(y);  
t=0:1/fp:(N-1)/fp;
subplot(221);
plot(t,y);
xlabel('czas[s]');
ylabel('sygnal');

Nf=2^14;
N21=Nf/2+1; 
v=fft(y,Nf); 
wy=abs(v); 
f=linspace(0,fp/2,N21);
subplot(222);
plot(f,wy(1:N21));
xlabel('czestotliwosc[Hz]');
ylabel('modul widma');


dr=4;
y=decimate(x,dr);   
fp=fs/dr;
N=length(y);        
t=0:1/fp:(N-1)/fp;      
subplot(223);
plot(t,y);
xlabel('czas[s]');
ylabel('sygnal');

Nf=2^14;                
N21=Nf/2+1;            
v=fft(y,Nf);         
wy=abs(v);              
f=linspace(0,fp/2,N21);   
subplot(224);
plot(f,wy(1:N21));
xlabel('czestotliwosc[Hz]');
ylabel('modul widma');