filip32Print

%część 2:

%[x,fpx]=wavread('/home/c/mbi04czap.wav');
[x,fpx]=audioread('/mbi04czap.wav');
Nx=length(x);
tx=0:1/fpx:(Nx-1)/fpx;
subplot(2,2,1);
plot(tx,x);
xlabel('czas[s]');
ylabel('sygnal');






Nfx=2^15;
N21x=Nfx/2+1;
v=fft(x,Nfx);
wx=abs(v);
fx=linspace(0,fpx/2,N21x);
subplot(2,2,2);
plot(fx,wx(1:N21x));
xlabel('czestotliwosc[Hz]');
ylabel('modul widma');






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







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