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

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