Analiza widma sygnału parc.wav przed i po decymacji (dr=4) - wykresy i FFT

Witam. Zrobiłem trochę nieporęcznie wykresy i widma przed i po decymacji sygnału parc.wav.

clear;
oneplot;
hold off;

% sygnal parc

[x,fp] = auload('mbi02parc.wav');
N = length(x);
fp
t = linspace(0,(N-1)/fp,N);
plot(2,2,1);
clg;
title('sygnał parc.wav');
xlabel('czas [s]');
ylabel('amplituda');
plot(t,x);
'po pierwszym wykresie'
pause

% sygnal po decymacji

dr = 4;
Nd = floor(N/dr);
fd = fp/dr;
for i = 1:Nd
xd(i)=x(i*dr);
end;
td = linspace(0,(Nd-1)/fd,Nd);
plot(2,2,3);
clg;
xlabel('czas[s]');
ylabel('amplituda');
title('sygnał po decymacji dr = 4');
plot(td,xd);
'po drugim wykresie'
pause

%widmo wygnalu

M=2^14;
v=fft(x,M);
wx=abs(v);
f=linspace(0,fp,M);
plot(2,2,2);
clg;
title('widmo parc.wav');
xlabel('czestotliwosc [Hz]');
ylabel('| fft(x) |');
plot(f,wx);
'po trzecim wykresie'
pause

% widmo sygnału po decymacji

vd=fft(xd,M);
wxd=abs(vd);
plot(2,2,4);
clg;
xlabel('czestotliwosc [Hz]');
ylabel('| fft(x) |');
title('widmo sygnału po decymacji dr = 4');
plot(f,wxd);
'po czwartym wykresie'


Problem polega na tym że czwarty wykres powinien mieć 12khz a ja nie mam pojęcia jak to zrobić. W octavie oczywiście