数字信号处理实验二报告材料.docx

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