数字信号滤波器的设计.docVIP

三峡大学 课 程 设 计 报 告 专业班级 09电信 课 程 数字信号处理 学 号 2 学生姓名 杰 指导教师 2012年 5月 平时成绩（20%） 报告成绩（40%） 答辩成绩（40%） 总成绩 一、设计要求： 1、课程设计题目 题目见附录，每人至少两道题。实验一2小题中必选一道，实验二和实验三4小题中必选一题，全做更好（答辩时都要参加回答问题）。 2、设计安排 1） 课程设计的时间为一周 2）第1天明确设计要求，选定设计题目 3）第1——4天学生进行设计 4）第4，5天教师验收，然后学生撰写和打印课程设计说明书 3、指导与答疑 每天都有教师现场答疑，学生有疑难问题可找教师答疑，设计时应充分发挥主观能动性，不应过分依赖教师。 4、课程设计说明书应包括封面（设计题目，学生的姓名，班级和学号），设计原理，设计过程，程序（matlab），运行结果，心得体会等几个部分，要求用A4纸打印，左侧装订。 二、设计步骤： 实验一：用双线性变换法和脉冲响应不变法设计IIR数字滤波器 采用双线性变换法设计一个巴特沃斯数字低通滤波器，要求：通带截止频率100Hz， 阻带截止频率200Hz，通带衰减指标Rp小于2dB，阻带衰减Rs大于15dB，滤波器采样频率Fs=500Hz。绘制频率响应曲线。 解：1）IIR滤波器的设计原理：双线性不变法是使数字滤波器的频率响应与模拟滤波器的频率响应相似的一种变换方法。首先把整个s平面压缩变换到某一中介平面s1的一条横带里（宽度为），再通过标准变换关系将此横带变换到整个z平面上去。 2)IIR低通滤波器的设计步骤： （1）按一定规则将给出的数字的滤波器的技术指标转换为模拟滤波器的技术指标； （2）根据转换后的技术指标设计模拟低通滤波器H（s）； （3）再按一定规则将H（s）转换成H（z）； 3）由题知，，，，采样频率 ，下面求阶数N 模拟巴特沃斯低通滤波器的幅度平方函数为 将已知条件代入上式可得 由此联立求解，得出 因而把代入上式得所需滤波器的阶数N为 所以取大于此数的N=2。 4）求系统函数 因为N=2，故查表得出二阶模拟低通滤波器的系统函数为 又因为在低频处，模拟滤波器与数字滤波器有较确切的对应关系，即，所以得到，故得到数字滤波器的系统函数为 令，即得该数字滤波器的频率响应为 其中 通过MATLB编程绘出该巴特沃斯低通滤波器的频率响应曲线如下 wp=100;ws=200;Fs=500; rp=2;rs=15; [n,Wn]=buttord(wp/(Fs/2),ws/(Fs/2),rp,rs,z); [num,den]=butter(n,Wn); [H,W]=freqz(num,den); plot(W*Fs/(2*pi),abs(H));grid; xlabel(频率/Hz); ylabel(幅值); title(巴特沃斯低通滤波器的频率响应); 实验二：用窗函数法设计FIR数字滤波器 1、选择合适的窗函数设计一个FIR数字低通滤波器，要求：带通截止频率为Wp=0.2π，阻带截止频率Ws=0.3π，阻带衰减不小于40dB，通带衰减不大于3dB。描绘滤波器的幅频响应曲线。 解：1）设计原理 设所要求的理想数字滤波器的频率响应为是与其对应的单位脉冲响应，因此先由的傅里叶反变换导出，即 由于是矩形频率特性，故一定是无限长的序列，且是非因果的，而我们要设计的是FIR滤波器，其必然是有限长的，所以要用有限长的来逼近无限长的。最有效的方法是截断，即用有限长的来截取，即 这即是窗函数的设计方法。 2）设计步骤 （1）求各对应的数字频率 由题知通带截止频率，阻带截止频率，通带衰减指数，阻带衰减指数。 （2）求 设为理想线性相位滤波器 已知通带截止频率，阻带截止频率，故其对应的数字频率为 由此可得 其中，为线性相位所必需的移位，我们已经知道应满足。 （3）求窗函数 由阻带衰减确定窗形状，由过渡带宽确定 由于，查表知可选海明窗，其阻带最小衰减dB满足要求 所要求的过渡带宽（数字频域） 由于海明窗过渡带宽满足 所以，。 （4）求 由海明窗表达式确定FIR滤波器的 海明窗 所以 3）MATLAB编程求解如下： （1）函数ideal_lp源程序 function hd=ideal_lp(wc,N) int N; alpha = (N-1)/2; n