p4.m

来自「MATLAB语音信号处理功能」· M 代码 · 共 8 行

S=fft(s,512);
SF=fft(sf,512);
w=(0:255)/256*(Fs/2);
plot(w,abs([S(1:256)' SF(1:256)']));
xlabel('Frequency (Hz)');
ylabel('Mag.of Fourier transform');
grid;
legend({'before','after'});