贵州大学-MATLAB数字信号处理课程设计

xlabel('频率单位（HZ）') axis([0,800,0,1500])

3.2滤波器的设计

低通滤波器的设计

wp=2*pi*100;ws=2*pi*150;Rp=0.1;As=50; [N,wc]=buttord(wp,ws,Rp,As,'s'); [B,A]=butter(N,wc,'s'); k=0:511;

fk=0:150/512:150;wk=2*pi*fk; HK=freqs(B,A,wk);

plot(fk,20*log(abs(HK)));grid on title('低通滤波器'); ylabel('幅度/db'); xlabel('频率单位（HZ）'); axis([0,200,-1,1]);

通过上述程序的运行得到下述波形：

低通滤波器

10.80.60.40.2

幅度/db

0-0.2-0.4-0.6-0.8-1

频率单位（HZ）