设计数字滤波器.docVIP

14题：设计一个工作频率80kHz的巴特沃斯低通数字滤波器，要求通带边界频率为4kHz，通带最大衰减为0.5dB，阻带边界频率为20kHz，阻带最小衰减为45dB，调用MATLAB工具箱函数buttord和butter设计，并显示数字滤波器系统函数H(z)的系数，绘制损耗函数和相频特性曲线。

程序：

fs=80000;

wp=2__pi__4000/fs;ws=2__pi__20000/fs;rp=0.5;rs=45;

[N,wc]=buttord(wp/pi,ws/pi,rp,rs);

[B,A]=butter(N,wc);

[h,w]=freqz(B,A);

f=w/pi__40;

subplot(121);

plot(f,20__log10(abs(h)));

grid;

title(损耗函数曲线);__label(频率/KHz);ylabel(幅度/dB)

subplot(122);

plot(f,angle(h));

grid;

title(相频特性曲线);__label(频率/KHz);ylabel(相位/\pi)

程序运行结果：N=4

17题：设计一个工作频率5MHz的椭圆带通数字滤波器，要求通带边界频率为560kHz和780kHz，通带最大衰减为0.5dB，阻带边界频率为375kHz和1000kHz，阻带最小衰减为50dB，调用MAT