实验四IIR数字滤波器设计实验报告.docVIP

实验四IIR数字滤波器设计实验报告 HUNAN UNIVERSITY OF TECHNOLOGY 数字信号处理 实验报告 实验四 IIR 数字滤波器的设计 学 生 姓 名 张 志 翔 班 级 电子信息工程 1203 班 学 号 12401720522 指 导 教 师 xx429 实验四 IIR 数字滤波器的设计 一、 实验目的： 1. 掌握双线性变换法及脉冲响应不变法设计 IIR 数字滤波器的具体 设计方法及其原理，熟悉用双线性变换法及脉冲响应不变法设计低通、 高通和带通 IIR数字滤波器的 MATLA 编程。 2. 观察双线性变换及脉冲响应不变法设计的滤波器的频域特性，了 解双线性变换法及脉冲响应不变法的特点。 3. 熟悉 Butterworth 滤波器、切比雪夫滤波器和椭圆滤波器的频率 特性。 二、 实验原理 ： 1 .脉冲响应不变法 用数字滤波器的单位脉冲响应序列 模仿模拟滤波器的冲激响应 ， 让 正好等于 的采样值，即，其中 为采样间隔，如果以 及 分别表示 的 拉式变换及的 Z 变换，则 1 2 H(z ) ze sT 〒 H a (s j—m) I m I 2 .双线性变换法 S 平面与 z 平面之间满足以下映射关系: 1 Is s 平面的虚轴单值地映射于 z 平面的单位圆上， s 平面的左半平面完 全映射到 z 平面的单位圆内 双线性变换不存在混叠问题。可通过预畸而得到校正 三、实验内容及步骤： 1 T s (s j ； z re j ) 双线性变换是一种非线性变换 ，这种非线性引起的幅频特性畸变 实验中有关变量的定义 ： fc 通带边界频率； fr 阻带边界频率； 3 通带波动； At 最小阻 带衰减； fs 采样频率； T 采样周期 (1) =0.3KHz, S =0.8Db, =0.2KHz, At =20Db,T=1ms; 设计一个切比雪夫高通滤波器，观察其通带损耗和阻带衰减是否满足 要求。 MATLA 源程序： wp=2*1000*ta n(2*pi*300/(2*1000)); ws=2*1000*ta n(2*pi*200/(2*1000)); [N,wn]=cheb1ord(wp,ws,0.8,20,s); % 给定通带( wp )和阻带 (ws) 边界角频率，通带波动波动 0.8 ，阻带最小衰减 20dB, 求出最低阶数 和通带滤波器的通带边界频率 Wn [B,A]=cheby1(N,0.5,wn,high,s);% 给定通带( wp )和阻带 (ws) 边界角频率，通带波动 [nu m,de n]=bili near(BA 1000); [h,w]=freqz (nu m,de n); f=w/(2*pi)*1000; plot(f,20*log10(abs(h))); axis([0,500,-80,10]); grid;xlabel( 频率 );ylabel( 幅度 /dB) 程序结果 分析：由图可知，切比雪夫滤波器幅频响应是通带波纹，阻带单调衰 num = 0.0304 -0.1218 0.1827 -0.1218 0.0304 den = 1 1.3834 1.4721 0.8012 0.2286 系统函数： =0.0304 -0.1218Z 1 0.1827 Z 2 -0.1218 Z 3 0.0304 Z 4 =1.0000+1.3834 Z 1 +1.4721 Z 2 + 0.8012 Z 3 +0.2286 Z 4 幅频响应图: Figure 1 File Edit View Insert Tools Desktop Window Help □曰?空| 脅?熄貳 TD| 口图|?口 m蚩 减的。 3 =0.8 ， fr=0.2kHz,At=30Db, 满足设计要求 (2) fc=0.2kHz, S =1dB,fr=0.3kHz,At=25dB,T=1ms; 分别用脉冲响 应不变法及双线性变换法设计一 Butterworth 数字低通滤波器，观察 所设计数字滤波器的幅频特性曲线， 记录带宽和衰减量，检查是否满 足要求。比较这两种方法的优缺点。 MATLA 源 程序： T = 0.001;fs = 1000;fc = 200;fr = 300; wp1 = 2*pi*fc;wr1 = 2*pi*fr; [N1,w n1] = buttord(wp1,wr1,1,25,s) [B1,A1] = butter(N1,w n1,s); [nu m1,de n1] = imp in var(B1,A1,fs);% 脉冲响应不变法 [h1,w] = freqz (nu m1,de