数字信号处理设计IIR滤波器北华大学.doc

北华大学开放实验报告 实验项目名称：设计IIR滤波器 所属课程名称：数字信号处理 实验类型:设计型 实验日期: 班级:通信 学号: 姓名 一．实验目的 加深对脉冲响应不变法以及双线性变换法设计IIR数字滤波器基本方法的了解。 掌握使用模拟滤波器原型进行脉冲响应变换的方法以及掌握用双线性变换法设计数字低通，高通，带通，带阻滤波器的方法。 了解MATLAB有关脉冲响应变换及双线性变换法的子函数。 二．实验步骤 脉冲响应不变法设计步骤 确定数字滤波器性能指标 将数字滤波器频率指标转换成相应的模拟滤波器频率指标 根据指标设计模拟滤波器将展成部分分式形式 把模拟极点转换成数字极点，得到数字滤波器 可见至间的变换关系为 在MATLAB中有两种方法可以实现上述变换 。 方法：利用residue函数和residuez函数实现脉冲响应不变法，这两个函数的使用方法如下： [r,p,k]=residue(b,a) [b,a]=residue(r,p,k) 实现多项式形式 和部分分式形式 之间的转换。 [r,p,k]=residuez(b,a) [b,a]=residuez(r,p,k) 实现多项式形式 和部分分式形式 之间的转换 双线性变换法的设计步骤 ● 输入给定的数字滤波器设计指标 ● 将数字滤波器指标转换成模拟滤波器设计指标 ● 确定模拟滤波器的最小阶数和截止频率 ● 设计模拟低通原型滤波器的系统传递函数 ● 利用模拟域频率变换法，求解实际模拟滤波器的系统传递函数 ● 用双线性变换法将模拟滤波器转换成数字滤波器 三．实验内容 1、采用脉冲响应不变法设计一个巴特沃斯数字低通滤波器，要求：Wp=0.25π，Rp=1dB ;Ws=0.4π，As=15dB ,滤波器采样频率Fs=2000Hz。 编写的程序如下： %数字滤波器指标 wp=0.25*pi; %滤波器的通带截至频率 ws=0.4*pi; %滤波器的阻带截至频率 Rp=1;As=15; %输入滤波器的通阻带衰减指标 ripple=10^(-Rp/20); %计算通带衰减对应的幅度值 Attn=10^(-As/20); %计算阻带衰减对应的幅度值 %转换为模拟滤波器指标 Fs=2000;T=1/Fs; Omgp=wp*Fs;Omgs=ws*Fs; %模拟原型滤波器计算 [n,Omgc]=buttord(Omgp,Omgs,Rp,As,s) %计算阶数n和截止频率 [z0,p0,k0]=buttap(n); %设计归一化的巴特沃斯模拟原型滤波器 ba1=k0*real(poly(z0)); %求原型滤波器系数b aa1=real(poly(p0)); %求原型滤波器系数a [ba,aa]=lp2lp(ba1,aa1,Omgc); %变换为模拟低通滤波器 %用脉冲响应不变法计算数字滤波器系数 [bd,ad]=impinvar(ba,aa,Fs) %[C,B,A]=dir2par(bd,ad) [H,w]=freqz(bd,ad); dbH=20*log10((abs(H)+eps)/max(abs(H))); subplot(2,2,1),plot(w/pi,abs(H)); ylabel(∣H∣);title(幅度响应);axis([0,1,0,1.1]); set(gca,XTickMode,manual,XTick,[0,0.25,0.4,1]); set(gca,YTickMode,manual,XTick,[0,Attn,ripple,1]);grid subplot(2,2,2),plot(w/pi,angle(H)/pi); ylabel(\phi);title(相位响应);axis([0,1,-1,1]); set(gca,XTickMode,manual,XTick,[0,0.25,0.4,1]); set(gca,YTickMode,manual,XTick,[-1,0,1]);grid subplot(2,