数字信号处理综合性实验报告.doc

重庆交通大学 综合性设计性实验报告 班 级：xxxxxxxxxxxxxxxxxxxxxxxxx 学 号： xxxxxxxxxxxxxxx 姓 名： xxxxxxxxxx 实验项目名称： 虑波器设计综合实验应用 实验项目性质： 综合性 实验所属课程： 数字信号处理 实验室(中心)： 现代电子实验中心 指 导 教 师 ： xxxxxxxxx 实验完成时间： 2014 年 12 月 25 日 一、实验目的 1、学生自己运用 MATLAB 设计IIR 数字低通滤波器，方法不限。（基本要求） 2、实现信号的滤波。（要求扩展） 二、实验主要内容及过程 任务1、设计一模拟IIR 模拟低通滤波器并转换为数字IIR 低通滤波器。(考核基本要求) 1）模拟滤波器设计采用巴特沃斯或者切比雪夫一型滤波器作为原型。 2）模拟到数字转换采用冲激响应或者双线性映射法。 任务2、利用实现的滤波器对信号进行滤波。（考核扩展功能部分） 完成此步骤需要与任务一综合考虑。分析有用和噪音信号的频率，并参考噪音信号的频率利用采样定理等知识选取合理的滤波器截止频率进行任务一的设计。 三、设计方案 （一）题目 用冲激响应不变法设计 Butterworth 低通数字滤波器，要求通带频率为，通带波纹小于1dB，阻带在内，幅度衰减大于15dB，采样周期T=0.01s。 (1)用冲激响应不变法设计该数字滤波器 H(z)。 (2)使用MATLAB 软件对滤波器性能进行分析。 (3)假设一个信号,其中f1=5Hz,f2=30Hz。试将原信号与经过该滤波器的输出信号进行比较。 （二）设计的主要思路 IIR 数字滤波器经典设计法的一般步骤为： (1)根据给定的性能指标和方法不同，首先对设计性能指标中的频率指标，如数字边界频率进行变换，转换后的模拟频率指标作为模拟滤波器原型设计的性能指标。 (2)估计模拟滤波器最小阶数和截止频率，利用MATLAB 工具函数buttord、cheb1ord、cheb2ord、ellipord 等。 (3)设计模拟低通滤波器原型。利用MATLAB 工具函数buttap、cheb1ap、cheb2ap、ellipap 等。 (4)由模拟原型低通滤波器经频率变换获得模拟滤波器（低通、高通、带通、带阻等）利用MATLAB 工具函数lp2lp、lp2hp、lp2bp、lp2bs。 (5)将模拟滤波器离散化获得IIR 数字滤波器，利用MATLAB 工具函数bilinear 或impinvar。设计IIR 滤波器时，给出的性能指标通常分数字指标和模拟指标两种。数字 性能指标给出通带截止频率，阻带起始频率，通带波纹，阻带衰减等。数字频率和的取值范围为 ，单位弧度。而MATLAB 工具函数常采用归一化频率，和的取值范围为0~1，对应于，此时需进行转换。 模拟性能指标给出通带截止频率，阻带起始频率 ，通带波纹Rp，阻带衰减Rs 等。模拟频率和单位为弧度/秒（rad/s）。 MATLAB 信号处理工具箱中，设计性能指标的转换应根据不同设计方法进行不同处理。 （三）主要功能 该滤波器是一个低通滤波器，能够对信号进行低通滤波，滤除信号的高频分量；它允许一定频段的信号通过，抑制低于或高于该频段的信号、干扰和噪声。 四、实验结果及说明（一）参数设计分析： 根据题目可知，所滤信号为,其中f1=5Hz,f2=30Hz。低通滤波器需要滤除f2=30Hz的分量。 对于这样一个信号来说，信号最高频率为30Hz。由采样定理知，采样频率，即采样频率应大于信号最高频率的两倍即60Hz，因此将滤波器的采样率设置为100Hz能够满足要求，此时滤波器最高频率对应此采样率下信号最高频率50Hz。可算出f1对应，f2对应。 最终将数字滤波器性能设计为：通带频率为，通带波纹小于1dB，阻带在内，幅度衰减大于15dB，采样周期T=0.01s。应当能够满足对信号 进行滤波的要求。 （二）数字频率到模拟频率的转换： 根据理论可知，数字频率和模拟频率之间的关系为： 由此，可以计算出模拟原型滤波器的性能参数要求为：滤波器在rad/s处的幅频响应衰减不多于-1dB，在rad/s处的幅频特性响应衰减应大于-15dB。 即： （三）模拟低通原型滤波器设计过程： 将参数代入公式： 计算可得： 将代入公式： 计算可得： 根据，从巴特沃斯滤波器多项式表中找到归一化（即时）的巴特沃斯低通原型滤波器的系统函数为