matlab - 将音频信号划分为 Barks Band

Question

我想在matlab中应用以下步骤

1) 读取音频数据

2）对该数据应用窗口函数

3）对窗口函数的结果应用fft

4）制作fft的功率谱

5) 将功率谱划分为 Barks 频带

[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');

我没有得到正确的波形，想绘制频谱，然后分成树皮带。