fft1.m

来自「一些数字信号处理课程的小程序」· M 代码 · 共 5 行

S=fft(s,512);
SF=fft(sf,512);
w=(0:255)/256*(fs/2);
plot(w,abs(S(1:256)),w,abs(SF(1:256)),'r');
grid on;