设计数字滤波器
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
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('损耗函数曲线');xlabel('频率/KHz');ylabel('幅度/dB') subplot(122); plot(f,angle(h)); grid;
title('相频特性曲线');xlabel('频率/KHz');ylabel('相位/\pi')
程序运行结果：N=4
-250
-200
-150
-100
-50
50
损耗函数曲线
频率/KHz
幅度/d B
-4
-3
-2
-1
01
2
3
4相频特性曲线
频率/KHz
相位/
17题：设计一个工作频率5MHz 的椭圆带通数字滤波器，要求通带边界频率为
560kHz 和780kHz ，通带最大衰减为0.5dB ，阻带边界频率为375kHz 和1000kHz ，阻带最小衰减为50dB ，调用MATLAB 工具箱函数ellipord 和ellip 设计，并显示数字滤波器系统函数H(z)的系数，绘制损耗函数和相频特性曲线。
程序：
fpl=560000;fpu=780000;fsl=375000;fsu=1000000;fs=5000000; wp=[2*fpl/fs,2*fpu/fs]; ws=[2*fsl/fs,2*fsu/fs];rp=0.5; rs=50; [N,wp0]=ellipord(wp, ws, rp, rs); [B,A]=ellip(N,rp,rs,wp0); [h,w]=freqz(B,A); f=w/pi*2500; subplot(121);
plot(f,20*log10(abs(h))); grid;
title('损耗函数曲线');xlabel('频率/KHz');ylabel('幅度/dB') subplot(122); plot(f,angle(h)); grid;
title('相频特性曲线');xlabel('频率/KHz');ylabel('相位/\pi')
程序运行结果：N=4
-100
-90-80-70
-60-50-40
-30-20
-100损耗函数曲线
频率/KHz
幅度/d
B
-4
-3-2
-1012
3
4相频特性曲线
频率/KHz
相位/
18题：设计一个工作频率5KHz 的椭圆带通数字滤波器，要求通带边界频率为
500kHz 和2125kHz ，通带最大衰减为1dB ，阻带边界频率为1050kHz 和1400kHz ，阻带最小衰减为40dB ，调用MATLAB 工具箱函数ellipord 和ellip 设计，并显示数字滤波器系统函数H(z)的系数，绘制损耗函数和相频特性曲线。
程序：
fpl=500;fpu=2125;fsl=1050;fsu=1400;fs=5000;
wp=[2*fpl/fs,2*fpu/fs]; ws=[2*fsl/fs,2*fsu/fs];rp=1; rs=40; [N,wp0]=ellipord(wp, ws, rp, rs); [B,A]=ellip(N,rp,rs,wp0,'stop'); [h,w]=freqz(B,A); f=w/pi*2500; subplot(121);
plot(f,20*log10(abs(h))); grid;
title('损耗函数曲线');xlabel('频率/Hz');ylabel('幅度/dB') subplot(122); plot(f,angle(h)); grid;
title('相频特性曲线');xlabel('频率/Hz');ylabel('相位/\pi')
程序运行结果： N=3
-90
-80-70-60-50-40-30
-20-100损耗函数曲线
频率/Hz
幅度/d
B
-4
-3-2-1012
34相频特性曲线
频率/Hz
相位/