syp.m

来自「我的数字信号处理课程设计」· M 代码 · 共 4 行

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