Разделение аудиосигналов на полосы Баркса

Question

Я хочу применить следующие шаги в Matlab

1) Читайте аудио данные

2) Применить оконную функцию к этим данным

3) Применить fft к результату оконной функции

4) Сделать спектр мощности БПФ

5) разделить спектр мощности по полосам Баркса

[y,Fs,bits] = wavread('file1.wav');
Nsamps = length(y);
ham= hamming(Nsamps);
ham_fft = abs(fft(ham));
ham_fft = ham_fft(1:Nsamps/2);      %Discard Half of Points
plot(f, ham_fft,f1)
xlim([0 1000])  
xlabel('Frequency (Hz)')
ylabel('Amplitude')
title('FFT of Hamming Window of  a an audio');

Я не получаю правильную форму волны, хочу нарисовать спектр и затем разделить на полосы лая.