实验四IIR数字滤波器设计及软件实现实验报告.docx

如文档对您有帮助，欢迎下载支持，谢谢！ 如文档对您有帮助，欢迎下载支持，谢谢！ PAGE PAGE 10 实验四 IIR 数字滤波器设计及软件实现 实验报告 一、实验目的 （ 1）熟悉用双线性变换法设计 IIR 数字滤波器的原理与方法； （ 2）学会调用 MATLAB信号处理工具箱中滤波器设计函数（或滤波器设计分析工具fdatool ）设计各种 IIR 数字滤波器，学会根据滤波需求确定滤波器指标参数。 （ 3）掌握 IIR 数字滤波器的 MATLAB实现方法。 （ 3）通过观察滤波器输入输出信号的时域波形及其频谱，建立数字滤波的概念。 二、实验原理 设计 IIR 数字滤波器一般采用间接法（脉冲响应不变法和双线性变换法） ，应用最广泛的是双线性变换法。 基本设计过程是： ①先将给定的数字滤波器的指标转换成过渡模拟滤波 器的指标； ②设计过渡模拟滤波器；③将过渡模拟滤波器系统函数转换成数字滤波器的系 统函数。 MATLAB信号处理工具箱中的各种 IIR 数字滤波器设计函数都是采用双线性变换法。第六章介绍的滤波器设计函数 butter 、cheby1 、cheby2 和 ellip 可以分别被调用来直接设计巴特沃斯、切比雪夫 1、切比雪夫 2 和椭圆模拟和数字滤波器。本实验要求读者调用如上函数直接设计 IIR 数字滤波器。 本实验的数字滤波器的 MATLAB实现是指调用 MATLAB信号处理工具箱函数 filter 对给定的输入信号 x(n) 进行滤波，得到滤波后的输出信号 y(n ） 。 三、实验内容及步骤 （1）调用信号产生函数 mstg 产生由三路抑制载波调幅信号相加构成的复合信号 st ，该函数还会自动绘图显示 st 的时域波形和幅频特性曲线， 如图 10.4.1 所示。 由图可见， 三路信号时域混叠无法在时域分离。 但频域是分离的， 所以可以通过滤波的方法在频域分离， 这就是本实验的目的。 图 10.4.1 三路调幅信号 st 的时域波形和幅频特性曲线 （ 2）要求将 st 中三路调幅信号分离，通过观察 st 的幅频特性曲线，分别确定可以分 离 st 中三路抑制载波单频调幅信号的三个滤波器 （低通滤波器、 带通滤波器、 高通滤波器） 的通带截止频率和阻带截止频率。要求滤波器的通带最大衰减为 0.1dB, 阻带最小衰减为60dB。 提示：抑制载波单频调幅信号的数学表示式为 s(t) cos(2 f0t)cos(2 fct) 1 [cos(2 ( 2 fc f0 )t ) cos(2 ( fc f0 )t )] 其中， cos(2 fct ) 称为载波， f c 为载波频率， cos(2 f0 t) 称为单频调制信号， f 0 为调制正 弦波信号频率，且满足 f c f0 。由上式可见，所谓抑制载波单频调幅信号，就是 2 个正弦 信号相乘，它有 2 个频率成分：和频 fc f 0 和差频 f c f0 ，这 2 个频率成分关于载波频率 f c 对称。所以， 1 路抑制载波单频调幅信号的频谱图是关于载波频率 f c 对称的 2 根谱线，其中没有载频成分，故取名为抑制载波单频调幅信号。容易看出，图 10.4.1 中三路调幅信号的载波频率分别为 250Hz、500Hz、1000Hz。如果调制信号 m(t) 具有带限连续频谱，无直流 成分，则 s(t ) m(t ) cos(2 fct) 就是一般的抑制载波调幅信号。其频谱图是关于载波频率 f c 对称的 2 个边带（上下边带） ，在专业课通信原理中称为双边带抑制载波 (DSB-SC) 调幅 信号 , 简称双边带 (DSB) 信号。如果调制信号 m(t) 有直流成分， 则 s(t) m(t) cos(2 fct ) 就 是一般的双边带调幅信号。其频谱图是关于载波频率 f c 对称的 2 个边带（上下边带） ，并包含载频成分。 （ 3）编程序调用 MATLAB滤波器设计函数 ellipord 和 ellip 分别设计这三个椭圆滤波器，并绘图显示其幅频响应特性曲线。 （ 4）调用滤波器实现函数 filter ，用三个滤波器分别对信号产生函数 mstg 产生的信 号 st 进行滤波，分离出 st 中的三路不同载波频率的调幅信号 y 1(n) 、y 2(n) 和 y 3(n) ， 并绘图显示 y1(n) 、y2(n) 和 y3(n) 的时域波形，观察分离效果。 四、信号产生函数 mstg 清单 function st=mstg %产生信号序列向量 st, 并显示 st 的时域波形和频谱 %st=mstg 返回三路调幅信号相加形成的混合信号，长度 N=1600 N=1600 %N 为信号