实验脉冲响应不变法设计IIR数字滤波器.docx

实验5脉冲响应不变法设计IIR数字滤波器一、实验目的1.掌握利用脉冲响应不变法设计IIR数字滤波器的原理及具体方法。2.加深理解数字滤波器和模拟滤波器之间的技术指标转化。3.掌握脉冲响应不变法设计IIR数字滤波器的优缺点及适用范围。二、实验设备及环境计算机、matlab软件环境。三、实验基础理论1.基本原理从时域响应出发，是数字录波器的单位脉冲响应h(n)模仿模拟滤波器的单位冲击响应，h(n)等于的取样值。2变换方法将进行部分分式展开 (2)对进行拉式变换(3)对时域采样得到h(n)(4)对h(n)进行z变换3.设计步骤 (1)确定数字滤波器性能指标，，和。 (2)将数字滤波器频率指标转换成相应的模拟滤波器频率指标 (3)根据指标，，和设计模拟滤波器。 (4)将展成部分分式形式 (5)将模拟极点转换成数字极点，得到数字滤波器可见至H(z)间的变换关系为在MATLAB中有两种方法可以实现上述变换。方法1：利用residue函数和residuez函数实现脉冲响应不变法，这两个函数的使用方法如下[r,p,k]=residue(b,a)[b,a]=reisdue(r,p,k)实现多项式形式和部分分式形式之间的转换。[r,p,k]=residuez(b,a)[b,a]=reisduez(r,p,k)实现多项式形式和部分分式形式之间的转换。方法2：MATLAB提供了impinvar函数采用脉冲响应不变法实现模拟滤波器到数字滤波器的变换，其使用方法如下[bz,az]=impinvar(b,a,fs)采用脉冲响应不变法将模拟滤波器系统函数的系数向量b和a变换成为数字滤波器系统函数的系统向量bz和az，fs为采样频率。[bz,az]=impinvar(b,a)采样频率默认为1的情况下，采用脉冲响应不变法将模拟滤波器变换为数字滤波器。四、实验内容1.设采用率为,采用脉冲响应不变法设计一个三阶巴特沃斯数字低通滤波器，其3dB截止频率为。实验代码及结果： [b, a] = butter(3, 0.5*pi, s)b = 0 0 0 3.8758a = 1.0000 3.1416 4.9348 3.8758 [bz, az] = impinvar(b, a)bz = 0 0.5813 0.2114 0az = 1.0000 -0.3984 0.2475 -0.0432 w = [0:500]*pi/500; [H, w] = freqz(bz, az); subplot(211); plot(w/pi, abs(H)); xlabel(\omega(\pi)); ylabel(|H(e^j^\omega)|); subplot(212); plot(w/pi, angle(H)/pi); xlabel(\omega(\pi)); ylabel(Phase of H(e^j^\omega)(\pi));2.设采用频率为=10kHz，设计数字低通滤波器，满足如下指标：通带截止频率：，通带波动：阻带截止频率：,阻带衰减：要求分别采用巴特沃斯、切比雪夫1型、切比雪夫2型和椭圆模拟原型滤波器及脉冲响应不变法进行设计。结合实验结果，分别讨论采用上述方法设计的数字滤波器是否都能满足给定指标要求，分析脉冲响应不变法设计IIR数字滤波器的优缺点及适用范围。实验代码及结果：Butterworth滤波器： [N, Wn] = buttord(0.2*pi, 0.3*pi, 1, 15, s)N = 6Wn = 0.7087 [b, a] = butter(N, Wn, s)b = 0 0 0 0 0 0 0.1266a = 1.0000 2.7380 3.7484 3.2533 1.8824 0.6905 0.1266 [bz, az] = impinvar(b, a)bz = -0.0000 0.0007 0.0105 0.0167 0.0042 0.0001 0az = 1.0000 -3.3443 5.0183 -4.2190 2.0725 -0.5600 0.0647 w = [0:500]*pi/500;[H, w] = freqz(bz, az);subplot(211);plot(w/pi, abs(H));xlabel(\omega(\pi));ylabel(|H(e