实验六IIR数字滤波器设计及应用.docxVIP

YUNNAN NORMAL UNIVERSITY 本科学生实验报告 学号 k!^ k￡^ k￡^ rTw rr^ rr^ rr^ rr^ rr rr rj rTw rfw T ri ri i rT ri ri 姓名 学院物电学院 专业、班级 k! k! k￡^ k￡^ k￡^ kP k{ k{ k{ k{ k! k! kJ k￡ 2^ 2^ 叫?叫?叫?叫?叫?叫?叫? 实验课程名称 数字信号分析与处理 教师及职称 vt k{^ k!^ k!^ k!^ k￡^ k!^ k!^ ^9^ k1^ J rj rj rj rj rj rj i i T i i i rT 开课学期2015 至2016学年上学期 填报时间 2016年 5 月 12 H 云南师范大学教务处编印 一、验设计方案 实验序号实验时间实验六买验名称IIR数字滤波器设计及应用2016/5/12 实验序号 实验时间 实验六 买验名称 IIR数字滤波器设计及应用 2016/5/12 实验室 同析楼三栋313实验室 1.实验目的 加深理解IIR数字滤波器的特性，掌握IIR数字滤波器的设计原理与设计方法，以及IIR 数字滤波器的应用。 2?实验原理、实验流程或装置示意图 【例231】设计一个数字滤波器，要求在0~0.2龙(rad)内衰耗不大于3dB, 在0.6〃■?71 (rad)内衰耗不小于60dB。 [解] 若模拟低通滤波器原型采用巴特沃思滤波器： [N,Wc]=buttord(0.2,0.6,3,60); [b,a]=butter(N,Wc); reqz(b,a);axis([0,l，-120,0]); 设计结果如图2.3.1所示。可见IIR是数字滤波器相位为非线性。 -50-1000 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 -50 -100 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Normalized Frequency (xti rad/sample) 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1Normalized Frequency (xk 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Normalized Frequency (xk rad/sample) ?100 -200 -300 -400 ■500 图2.3」 Buttertworth低通滤波器 【例2.3.2】设计一个数字带通滤波器，要求在100?200Hz通带内纹波不大 于3dB,通带两边各50Hz外是阻带，衰耗不小于40dB。抽样频率为lOOOHzo ［解］ Wp二[100 200]/500;Ws=[ 100-50 200+50]/500; Rp=3;Rs=40; |N,Wn|=cheb 1 ord(Wp,Ws,Rp,Rs); [b,a]=cheby l(N,Rp,Wn); freqz(b,a,512,1000); title(fchebyshev Type I Bandpass Filter1); axis([0,500,-80,0]); 设计结果如图2.3.2所示。 Chebyshev Type I Bandpass Filter 00 50 100 150 200 250 300 350 400 450 500Frequency (Hz)20■40-60-(8P)① P3U6E Chebyshev Type I Bandpass Filter 0 0 50 100 150 200 250 300 350 400 450 500 Frequency (Hz) 20 ■ 40 - 60 - (8P)① P3U6E 乏 -80 500 0 -500 -1000 0 50 100 150 200 250 300 350 400 450 500 Frequency (Hz) 图2.3.2 Chebyshev I型滤波器 3?实验设备及材料 计算机，MATLAB软件 4.实验方法步骤及注意事项 注意事项： 在使用MATLAB时应注意中英输入法的切换，在中文输入法输入程序时得到的程序 是错误的； MATLAB中两个信号相乘表示为x.*u,中间有个同样两个信号相除也是如此； 使用MATLAB编写程序吋，应新建一个m文件，而不是直接在Comandante窗口下编 写程序； 在使用MAT LAB编程时，应该养成良好的编写习惯。 5?实验数据处理方法 图像法、比较法 6.参考文献 《信号分析与处理》 《MATLAB数值计算与方法》 二、报告 1.实验现象与结果 实验内容 第一题(1) 实验程序代码 k=0:100; x=cos(pi*k/4)+cos(2*pi*k