实验五实验报告分析解析.docxVIP

综合性、设计性实验报告 姓名 陶建敏 学号 201308002107 专业 通信工程 班级 通信 1301 实验课程名称 数字信号处理 指导教师及职称 李玲香 讲师 开课学期 2015 至＿ 2016 学年 上 学期 上课时间 2015 年 12 月 10、17 日 湖南科技学院教务处编印 设计题目  实验五  IIR  滤波器的实现与应用 小组合作  是（  ） 否（  √） 小组成员 实验目的： （ 1）熟悉用双线性变换法设计 IIR 数字滤波器的原理与方法； （ 2）学会调用 MATLAB 信号处理工具箱中滤波器设计函数设计各种 IIR 数字滤波器， 学会根据滤波需求确定滤波器指标参数。 3）掌握 IIR 数字滤波器的 MATLAB 实现方法。 4）通过观察滤波器输入输出信号的时域波形及其频谱，建立数字滤波的概念。 实验环境： 1） Windous XP 或、 Windous 7 之类的操作系统版本； 2） MATLAB7.0 实验相关知识 无限长单位冲击响应（ IIR ）数字滤波器的设计方法有多种，对于一般条件下使用 的数字滤波器，其常用的设计方法是基于模拟滤波器变换原理的经典设计：首先，将数 字滤波器的技术指标转换成相应的模拟滤波器技术指标；设计模拟低通滤波器；然后， 将设计好的模拟滤波器转换成满足给定技术指标的数字滤波器。 MATLAB 的数字信号 处理工具箱中，提供的有关设计函数如下： 1）模拟低通滤波器设计 模拟低通滤波器的逼近有巴特沃思型、切比雪夫型、和考尔型，分别用如下的函数 实现。 [Z,P,K]=buttap(n) ； 返回一个 n 阶、巴特沃思型归一化的模拟低通滤波器的零极点 增益模型 [Z,P,K]=cheb1ap(n,Rp) [Z,P,K]=cheb2ap(n,Rs) [Z,P,K]=elliap(n,Rp,Rs)  ；n 阶、通带内的最大衰减； n 阶、阻带内的最小衰减； n 阶、通带内的最大衰减  Rp、切比雪夫Ⅰ型 Rs、切比雪夫Ⅱ型 Rp、阻带内的最小衰减  Rs、考 尔型。 （ 2）模拟低通滤波器阶数  n 的选择函数 滤波器阶数的选择在整个滤波器的设计中占有十分重要的地位和作用。 根据需要选择合适的滤波器阶数， MATLAB 工具箱中提供了对应于各类模拟低通滤波器的阶数选 择函数，如巴特沃思型的 buttord、切比雪夫型的 cheb1ord、cheb2ord 和考尔型的 ellipord 。 这些函数的调用格式大同小异。 [n,Wn]=buttord(Wp,Ws,Rp,Rs, ’s’) 输入参数： Wp 通带截止频率， Ws 阻带截止频率， Rp 通带最大衰减， Rs 阻带最小衰减。输出参数： n 为符合要求的滤波器最小阶数， Wn 为巴特沃思型模拟低通滤波器 3dB 截止频率。 ’s’：表示模拟域。 （ 3）零极点增益模型到传递函数模型的转换 [num,den]=zp2tf(Z,P,K) 输入参数： Z， P， K 分别表示系统的零极点增益模型的零点、极点和增益； 输出参数： num，den 分别为同一系统传递函数模型的分子和分母多项式系数。 4）模拟域的频率变换 将归一化的模拟低通滤波器转换成所需要类型（低通、高通、带通和带阻）的模拟 滤波器，可分别用如下命令实现： [b,a]=lp2lp(Bap,Aap,Wn) ；把传递函数形式的归一化模拟低通滤波器原型转换成 3dB 截止频率为 Wn 的同型低通滤波器。 [b,a]=lphp(Bap,Aap,Wn) ；转换成高通 [b,a]=lp2bp(Bap,Aap,W0,Bw) ；转换成带通， W0：中心频率， Bw ：带宽 [b,a]=lp2bs(Bap,Aap,W0,Bw) ；转换成带阻 5．模拟滤波器数字化 [bz,az]=bilinear(b,a,Fs) ：采用双线性变换法的映射关系。 其中，Fs 是采样频率。 [bz,az]=impinvar(b,a,Fs) ：采用冲击响应不变法的映射关系。 4．实验内容 设信号 f (t) 2sin(300 t ) sin(600 t) ，试将它的两个频率分量分离，并绘制它们的时域 波形及频谱图。具体要求： 1）分析实验内容，设计实验方案； 2）绘制 f (t) 时域和频域，观察输出波形，选择合理的滤波器参数； 3）绘制数字滤波器的频率特性曲线； 4）编 程 完成 f (t ) 信号 的分离 ，并 输出 分 离信 号 。 实验设计方案 （ 1）: 分析 f(t) 可知原信号有两个频率分量组成：一个频率分量是 150Hz，另一个频率 为 300Hz； 2）: 由 f(t) 由两个频率分量组成， 而要实现将这两个频率分量信号进行分离， 则必须要设计两个滤波器实现对两个