数字信号处理课程方案实验报告(打印).docVIP

北京工商大学 数字信号处理课程设计实验报告 班 级： 信息081 学 号： 姓 名： 同组同学姓名： 成 绩： 2011年7月 设计任务 1. 设计目的： （1）熟悉和巩固模拟滤波器的设计方法和原理 （2）掌握Butterworth/Chebyshev滤波器设计方法 （3）实现滤波器设计的有关经典算法 （4）熟练掌握使用高级语言程序设计各种要求的数字滤波器 （5）熟练掌握双线性变换方法 2.设计技术指标： （1）按要求设计Butterworth型数字低通滤波器， （2）性能指标如下： ① 通带截止频率； ② 通带最大衰减αp=3dB； ③ 阻带起始频率； ④ 阻带最小衰减αs=20dB； 3.设计要求： （1） 根据模拟滤波器的性能指标，确定数字滤波器指标； （2） 程序应具有通用性； （3）采用双线性变换法，设计满足上述性能指标要求的Butterworth型数字低通滤波器； （4）由软件直接给出数字滤波器的级联型结构； （5）确定数字滤波器的频率响应（幅值响应和相位响应）及群延迟等； 4.软件要求： 基本要求：软件能根据滤波器性能技术指标要求，能自动设计Butterworth型、Chebyshev型或其他滤波器； 5.其它要求： （1）学习并巩固模拟滤波器（Butterworth,Chebyshev型滤波器）设计的基本概念和基本理论； （2）熟悉Butterworth,Chebyshev型滤波器设计的有关公式； （3）掌握归一化设计方法中各种表格的应用方法； 二、实验原理与计算方法 １、双线性变换法设计IIR低通数字滤波器的基本原理和算法 双线性变换法设计数字滤波器，采用了二次映射的方法，就是先将整个s平面压缩到s1平面的一个的横形条带范围内，然后再将这个条带映射到z平面上，就能建立s平面到z平面的一一对应关系。对于低通数字滤波器，映射关系为 　　 　 (1) 其中T为抽样周期。 用双线性变换法设计低通IIR数字滤波器的基本步骤，首先根据设计要求确定相应的模拟滤波器的传递函数，再应用(1)式得数字滤波器的传递函数 　　　　　　　　　 　　　　　　　　　　　　　　 (2)通常可以给定的参数为：低通数字滤波器通带边界频率、阻带边界频率和对应的通带衰减函数、阻带衰减函数。s1平面中的模拟角频率与数字角频率的关系为线性关系，在计算模拟滤波器的阶数N、极点si和传递函数之前，应作预畸变处理 　　　 　　　　 (3) 　　模拟滤波器的阶数N、极点si和传递函数的计算方法与冲激响应不变法相同，可以采用Butterworth逼近或Chebyshev逼近。 ２、双线性变换法设计IIR高通带通数字滤波器的基本原理和算法由于双线性变换法获得的数字滤波器频率响应特性中不会出现混叠现象，因此可以适用于高通、带通和带阻滤波器的设计。IIR数字滤波器的设计通常要借助于模拟低通滤波器的设计，给定的对数字滤波器(DF)的技术要求转换为一个低通模拟滤波器(AF)的技术要求，根据这种要求用某种逼近设计出原型的低通模拟滤波器(LP AF)，计算出模拟滤波器的阶数N、极点si和传递函数，再按照双线性变换的变换关系，将模拟滤波器的传递函数转换为数字滤波器的传递函数。表1中列出将给定的对数字滤波器(DF)的技术要求直接转换为对一个低通模拟滤波器(AF)的技术要求的频率预校正关系表1 双线性变换和频率预校正的计算公式 变换类型 变换关系 频率预校正 备　　　　注 低通变换 其中，T为抽样周期，f为模拟频率 其中，分别为数字带通滤波器通带的上下边界频率数字带滤波器的上下边界频率。 高通变换 带通变换 例：数字高通滤波器的设计 首先将给定的数字高通滤波器的技术指标根据公式转换为模拟低通滤波器的技术指标，利用cheb1ord(Wp,Ws,ap,as,s)函数求出chebyshev模拟低通滤波器的阶数N，再利用cheb1ap(N,ap)函数求出模拟低通滤波器系统函数的零极点，zp2tf(z,p,k)函数将零极点转换为系统函数系数；然后利用lp2hp由模拟低通滤波器的系统函数得到模拟带通滤波器的系统函数，bilinear函数则用于实现双线性变换法由模拟带通滤波器系统函数计算数字数字带通滤波器系统函数的系数。 （2）数字频带变换法 首先将给定的对数字滤波器(DF)的技术要求转换为一个低通模