bj.m

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

[y,fs,bits]=wavread('b.wav');
m=fft(y,512);  %取y的512点fft变换。
p=(0:255)/256*(fs/2);  

wp=0.2;ws=0.4;rp=0.5;rs=40;
[N,Wc]=buttord(wp,ws,rp,rs);
[num,den]=butter(N,Wc);
[h,f]=freqz(num,den,256,1000);
lp=filter(num,den,y);

u=fft(lp,512);  %取y的512点fft变换。
lpp=(0:255)/256*(fs/2); 

plot(y-lp);
%subplot(223);plot(p,abs([m(1:256)]));%信号的频谱图
%subplot(222); plot(lp);%滤波后的信号时域波形
%subplot(224);plot(lpp,abs([u(1:256)]));%滤波后的信号的频谱图