lps.m

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

[y,fs,bits]=wavread('b.wav');
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); 
plot(lp);