整理IIR数字滤波器设计及软件实现.pptxVIP

实验四：IIR 数字滤波器设计及软件实现 一、实验原理与方法 1、设计 IIR 数字滤波器一般采用间接法（脉冲响应不变法和双线性变换法）， 应用最广泛的是双线性变换法，其基本设计过程是： 将给定的数字滤波器的指标转换成过渡模拟滤波器的指标； 设计过渡模拟滤波器； 将过渡模拟滤波器系统函数转换成数字滤波器的系统函数。 本实验的数字滤波器的MATLAB 实现是指调用MATLAB 信号处理工具箱函数 filter 对给定的输入信号x(n)进行滤波，得到滤波后的输出信号 y(n）。 二、实验内容 1、调用信号产生函数 mstg 产生由三路抑制载波调幅信号相加构成的复合 信号 st，该函数还会自动绘图显示 st 的时域波形和幅频特性曲线，如图 4.1 所 示。由图可见，三路信号时域混叠无法在时域分离。但频域是分离的，所以可以 通过滤波的方法在频域分离，这就是本实验的目的。 ;这 2 个频率成分关于载波频率 fc 对称。所以，1 路抑制载波单频调幅信号的频谱 图是关于载波频率 fc 对称的 2 根谱线，其中没有载频成分，故取名为抑制载波 单频调幅信号。容易看出，图 4.1 中三路调幅信号的载波频率分别为 250Hz、 500Hz 、1000Hz 。如果调制信号 m(t) 具有带限连续频谱， 无直流成分， 则 s(t) ? m(t) cos(2? fct) 就是一般的抑制载波调幅信号。其频谱图是关于载波频率 fc 对称的 2 个边带（上下边带），在专业课通信原理中称为双边带抑制载波 (DSB-SC) 调幅信号,简称双边带 (DSB) 信号。如果调制信号 m(t)有直流成分，则就是一般 的双边带调幅信号。其频谱图是关于载波频率 fc 对称的 2 个边带（上下边带）， 并包含载频成分。 3、编程序调用MATLAB 滤波器设计函数ellipord 和ellip 分别设计这三个椭 圆滤波器，并绘图显示其幅频响应特性曲线。实验结果如图 4.2、4.3、4.4，程序 见附录 4.1、4.2、4.3。 4、调用滤波器实现函数filter，用三个滤波器分别对信号产生函数mstg 产 生的信号 st 进行滤波，分离出 st 中的三路不同载波频率的调幅信号y1(n)、y2(n) 和y3(n)， 并绘图显示 y1(n)、y2(n)和 y3(n) 的时域波形，观察分离效果。实验结 果如图 4.2、4.3、4.4，程序见附录 4.1、4.2、4.3。 注：信号产生函数mstg 清单 function st=mstg %产生信号序列向量 st,并显示 st 的时域波形和频谱 %st=mstg 返回三路调幅信号相加形成的混合信号，长度 N=800 N=800 %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;;三、实验结果和分析、讨论及结论 1、滤波器参数选取 观察图 4.1 可知，三路调幅信号的载波频率分别为 250Hz、500Hz、1000Hz。 带宽（也可以由信号产生函数 mstg 清单看出）分别为 50Hz、100Hz、200Hz。所 以，分离混合信号 st 中三路抑制载波单频调幅信号的三个滤波器（低通滤波器、 带通滤波器、高通滤波器）的指标参数选取如下： 对载波频率为 250Hz 的条幅信号，可以用低通滤波器分离，其指标为： 通带截止频率 fp ? 280Hz ， 通带最大衰减 ap ? 0.1dB ； 阻带截止频率 fs ? 450Hz ，阻带最小衰减as ? 60dB 。 对载波频率为 500Hz 的条幅信号，可以用带通滤波器分离，其指标为： 带截止频率 fpl ? 440Hz ， fpu ? 560Hz ，通带最大衰减ap ? 0.1dB ；阻带截止 频率 fsl ? 275Hz ， fsu ? 900Hz ，阻带最小衰减as ? 60dB 。 对载波频率为 1000Hz 的条幅信号，可以用高通滤波器分离，其指标为： 带截止频率 fp ? 890Hz ，通带最大衰减ap ? 0.1dB ；阻带截止频率 fs ? 550Hz ， 阻带最小衰减as ? 60dB 。 说明：（1）为了使滤波器阶数尽可能低，每个滤波器的边界频率选择原则是 尽量使滤波器过渡带宽尽可能宽。 与信号产生函数 mstg 相同，采样频率 Fs=10kHz。 为了滤波器阶数最低，选用椭圆滤波器。 2、实验结果 由图 4.2、4.3、4.4 可见，三个分离滤波器指标参数选取正确，损耗函数曲线 达到所给指标。分离出的三路信号 y1(n)、y2(n) 和 y3(n) 的波形是抑制载波的单频 调幅波。 ;;四、思