实验七 (实验课) 幻灯片.pptVIP

Matlab及其应用 实验九  用窗口法设计FIR滤波器 一、实验目的 二、实验原理与方法 二、实验原理与方法 二、实验原理与方法 实验用MATLAB语言工具函数简介 实验用MATLAB语言工具函数简介 例9.1 运行结果： 三、实验内容及要求 四、思考题 五、实验报告要求 要 求： 六、相关程序： 六、相关程序： （1）了解一个实际滤波器设计过程， （2）掌握用窗口法设计FIR滤波器的原理和窗函数 对滤波器性能的影响。 Hd ( e jω ) 过渡带Δω 阻带衰减δ2 2. 求理想的单位取样响应 hd ( n )。 3. 由阻带衰减δ2 确定窗的类型 w( n )。 查表7.2.2 4. 由过渡带Δω及窗的类型确定 N 值 1. 确定所给理想频率响应函数及其技术指标。 h( n ) = hd ( n ) w( n ) 5. 求实际可设计的单位取样响应 h( n )。 6. 检验技术指标是否满足要求。 H ( e jω ) → Hd ( e jω ) ？？？ ( 若不满足，改变N 及窗的类型 ) 窗函数设计步骤： ( 1 ) 形成一过渡带 过渡带的宽度 ≈ 窗口频谱的主瓣宽度 ( 2 ) 在过渡带两旁产生肩峰和余振，它们取决于窗口频谱的旁瓣； 旁瓣越多，余振也越多；旁瓣相对值越大，肩峰则越强。 对频率响应会产生以下三点影响： h( n ) = hd ( n ) w( n ) 旁瓣与主瓣的相对关系只决定于窗口函数形状。 增加N，只能相对应减小过渡带宽，而不能改变肩峰值。 ( 3 ) 增加截取长度N， 只能缩小窗口频谱的主瓣宽度 而不能改变旁瓣的相对值； 肩峰值的大小直接决定通带内的平稳和阻带的衰减，对滤波器性能有很大关系。 例如： 矩形窗的情况下，肩峰达8.95%，致使阻带最小衰减只有21分贝，这在工程上往往是不够的。 怎样才能改善阻带的衰减特性呢？ 只能改变窗口的形状，所以希望的窗口频谱中应该尽量减少旁瓣，使能量集中在主瓣，这样可以减少肩峰和余振，提高阻带的衰减。 而且要求主瓣宽度尽量窄，以获得较陡的过渡带。 然而这两个要求总不能兼得，往往需要用增加主瓣宽度带换取决对主瓣的抑制，于是提出了海明窗、凯塞窗、切比雪夫窗等窗口函数。 在MATLAB的数字信号处理工具箱中提供了函数fir1。 fir1是采用经典窗函数法设计线性相位FIR数字滤波器，且具有标准 低通、带通、高通和带阻等类型。 语法格式： B=fir1(n,Wn) B=fir1(n,Wn,’ftype’) B=fir1(n,Wn,window) B=fir1(n,Wn,’ftype’,window) 其中，n为FIR滤波器的阶数; ‘ftype’为滤波器类型，缺省时为低通或带通滤波器， 为’high’时是高通滤波器, 为’stop’时是带阻滤波器。 window为窗函数，列向量，其长度为n+1， 缺省时，自动取hamming窗。 输出参数B为FIR滤波器系数向量，长度为n+1。 语法格式： B=fir1(n,Wn) B=fir1(n,Wn,’ftype’) B=fir1(n,Wn,window) B=fir1(n,Wn,’ftype’,win