syp.m

我的数字信号处理课程设计

[y,fs,bits]=wavread('b.wav');
m=fft(y,512); 
p=(0:255)/256*(fs/2);  
plot(p,abs([m(1:256)]));