利用matlab设计iir滤波器.doc

实验目的使用脉冲响应不变法设计IIR滤波器。使用双线性变化法设计IIR滤波器。比较两种IIR滤波器的设计方法有什么不同。使用频率脉冲采样法设计FIR滤波器,并与窗函数设计的FIR滤波器比较有什么不同。二、实验条件PC机,,通带截止频率=,阻带下限频率=,通带最大衰减为3dB,阻带最小衰减为20dB,给定Ts=。程序如下:Ts=;Ap=3;As=20;OmegaP=*pi/Ts;OmegaS=*pi/Ts;%模拟通带、阻带截止频率[n,Wn]=buttord(OmegaP,OmegaS,Ap,As,'s');%确定最小阶数n和反归一化截止频率Wn[b,a]=butter(n,Wn,'s');%b、a分别为模拟滤波器的分子、分母按降幂排列的多项式系数[bz,az]=impinvar(b,a,1/Ts);%脉冲响应不变法得到数字滤波器的分子分母系数omega=[0::pi];%确定坐标轴范围h=freqz(bz,az,omega);%得到模拟滤波器的单位冲激响应系数Ampli=20*log10(abs(h)/abs(h(1)));%求衰减的分贝subplot(2,1,1);plot(omega/pi,Ampli,'k');%显示滤波器的幅度响应xlabel('数字频率/\pi');ylabel('幅度/dB');grid;subplot(2,1,2);theta=phasez(bz,az,omega);%滤波器的相位响应及坐标值plot(omega/pi,theta*360/(2*pi),'k');%显示滤波器的相位响应xlabel('数字频率/\pi');ylabel('相位/度');grid;程序所得图像如下:利用双线性变换法设计一巴特沃斯低通数字滤波器,通带截止频率=,阻带下限频率=,通带最大衰减为3dB,阻带最小衰减为20dB,给定Ts=。程序如下:Ap=3;As=20;OmegaP=*pi;%数字通带截止频率OmegaS=*pi;%数字阻带截止频率[n,Wn]=buttord(OmegaP/pi,OmegaS/pi,Ap,As);%确定最小阶数n和反归一化截止频率Wn[bz,az]=butter(n,Wn);%bz、az分别为数字滤波器的分子、分母按降幂排列的多项式系数omega=[0::pi];%确定坐标轴范围h=freqz(bz,az,omega);%得到滤波器的单位冲激响应系数Ampli=20*log10(abs(h));%求衰减的分贝subplot(2,1,1);plot(omega/pi,Ampli,'k');%显示滤波器的幅度响应xlabel('数字频率/\pi');ylabel('幅度/dB');grid;subplot(2,1,2);theta=phasez(bz,az,omega);%滤波器的相位响应及坐标值plot(omega/pi,theta*360/(2*pi),'k');%显示滤波器的相位响应xlabel('数字频率/\pi');ylabel('相位/度');grid;程序所得图像:比较脉冲响应不变法与双线性变换法的区别:将两种方法的幅度响应做比