数字信号处理实验06.docVIP

实验六 2019年12月13日 一、实验目的 1.加深对双线性变换法设计IIR数字滤波器基本方法的了解。 2.掌握用双线性变换法设计数字低通、高通、带通、带阻滤波器的方法。. 3.了解MATLAB有关双线性变换法的子函数。 4.学习MATLAB直接设计IIR模拟滤波器和数字滤波器的方法。 5.了解MATLAB有关直接设计IIR模拟和数字滤波器的子函数，明确设计模拟滤波器和数字滤波器的区别。 6.初步了解采样频率的选择与数字滤波器实现的关系。 二、用到的MATLAB函数 1.bilinear 双线性变换——将s域映射到z域的标准方法，将模拟滤波器变换为离散等效滤波器。 2.butter 巴特沃兹模拟或数字滤波器设计。 3.cheby1 切比雪夫I型滤波器设计。 4.cheby2 切比雪夫II型滤波器设计。 5.ellip 椭圆滤波器设计。 三、实验原理 1.双线性变换法 双线性变换法是将整个s平面映射到整个z平面，其映射关系为： s 四、实验内容 1.用双线性变换法设计切比雪夫Ⅱ型数字滤波器，列出传递函数并描绘模拟和数字滤波器的幅频响应曲线。 设计一个数字低通，要求：通带wp=0.2π，Rp 传递函数系数如下： 分子：bd = 0.1624 0.0056 0.0056 0.1624 分母：ad = 1.0000 -1.4073 0.9056 -0.1623 设计一个数字高通，要求：通带wp=0.35π，Rp 传递函数系数如下： 分子：bd = 0.5221 -1.3531 1.3531 -0.5221 分母：ad = 1.0000 -1.4510 1.0283 -0.2711 2.采用双线性变换法设计一个切比雪夫Ⅱ型数字带通滤波器，要求：通带fp1=200Hz，fp2 传递函数系数： bd = 0.1397 0.0000 -0.0911 -0.0000 0.0911 0.0000 -0.1397 ad = 1.0000 -0.0000 1.1454 -0.0000 0.7275 -0.0000 0.1205 3.采用双线性变换法设计一个椭圆数字带阻滤波器，要求：下通带wp1=0.35π，上通带W 传递函数系数： bd = 0.4452 -0.0000 1.1373 -0.0000 1.1373 0.0000 0.4452 ad = 1.0000 0.0000 1.1776 0.0000 0.8975 -0.0000 0.0900 4.用MATLAB直接法设计切比雪夫Ⅱ型数字低通滤波器，要求： 通带wp=0.2π 系统传递函数系数： b = 0.1160 -0.0591 0.1630 -0.0591 0.1160 a = 1.0000 -1.8076 1.5891 -0.6201 0.1153 5.用MATLAB直接法设计巴特沃斯型数字带通滤波器，要求： 通带fp1=3.5kHz， 系统传递函数系数： b = 0.0612 0 -0.1837 0 0.1837 0 -0.0612 a = 1.0000 0.0000 0.9995 0 0.6006 -0.0000 0.1111 6.人体心电图信号在测量过程中往往受到工业高频干扰，所以必须经过低通滤波处理后，才能作为判断心脏功能的有用信息。下面给出一实际心电图信号采样序列样本x(n)，其中存在高频干扰，编写MATLAB代码，以x(n x -60,-84,-90,-66,-32,-4,-2,-4,8,12,12,10,6,6,6,4,0,0,0,0,0,-2,-4,0,0,0, -2,-2,0,0,-2,-2,-2,-2,0} 五、MATLAB代码 1. clear; close all; % 双线性变换法设计切比雪夫II型数字低通滤波器 wp = 0.2 * pi; ws = 0.35 * pi; Rp = 1; As = 15; ripple = 10^(-Rp/20); Attn = 10^(-As/20); % 转换为模拟滤波器指标 Fs = 10; T = 1/Fs; Omgp = (2/T) * tan(wp/2); Omgs = (2/T) * tan(ws/2); % 模拟原型滤波器计算 [n, Omgc] = cheb2ord(Omgp, Omgs, Rp, As, s); [z0, p0, k0] =