数字信号处理实验十三IIR滤波器设计().docVIP

实验十三 模拟滤波器设计及IIR数字滤波器设计 [实验目的] 1.设计模拟低通滤波器 2.研究并进行模拟滤波器到数字滤波器的转换，得到IIR数字滤波器 [实验原理] IIR设计步骤： 1、根据性能指标确定巴特沃斯滤波器的阶数和截止频率； 2、求出模拟巴特沃斯滤波器在s平面的极点； 3、由s平面左半平面的极点构成传递函数； 4、使用冲激不变法或双线性变换将模拟传递函数转换成数字滤波器的系统函数， 5、计算系统的幅频特性，检验是否符合性能指标。 模拟滤波器的设计 1、根据给定指标得出低通模拟原型滤波器的阶数n和截止频率Wn。 设滤波器的通带截止频率为wp，止带截止频率为ws，通带衰减为Rp，止带衰减为As。 利用函数buttord [N,wc]=buttord(wp,ws,Rp,As,’s’) 2、求出归一化低通模拟原型滤波器，利用函数Buttap [z,p,k]=Buttap(N) 由于所得的结果为零极点型，还必须转换成b/a型，可用函数zp2tf。 [bap,aap]=zp2tf(z,p,k) 3、将归一化模拟低通原型滤波器经频率变换为所要求的低通模拟滤波器可用函数lp2lp (matlab自带） 。 [b,a]=lp2lp(bap,aap,wn)//wn=2*pi*wc 4、求出滤波器的幅频，相频。 直接画出幅频特性和相频特性：freqs(b,a) (matlab自带） 5、使用冲激响应法来实现由模拟滤波器滤波器转变为数字滤波器。 利用Matlab自带的函数： [bz,az]=impinvar(b,a,Fs) [实验内容] 参考教材，完成本章课后习题1的模拟滤波器设计和采用冲激响应不变法将其转换为数字滤波器。编写matlab程序，分别求出模拟滤波器和数字滤波器的幅频特性曲线和相频特性曲线，并加以比较。 解： fp=6000; %fp=6000Hz fs=12000; %fs=12000Hz Fs=1000; wp=2*pi*fp; ws=2*pi*fs; rp=3; rs=25; [N,wc]=buttord(wp,ws,rp,rs,s); [z,p,k]=buttap(N); [bap,aap]=zp2tf(z,p,k); [b,a]=lp2lp(bap,aap,2*pi*wc); freqs(b,a); grid on 解： fp=6000; %fp=6000Hz fs=12000; %fs=12000Hz Fs=3000; wp=2*pi*fp; ws=2*pi*fs; rp=3; rs=25; [N,wc]=buttord(wp,ws,rp,rs,s); [z,p,k]=buttap(N); [bap,aap]=zp2tf(z,p,k); [b,a]=lp2lp(bap,aap,2*pi*wc); [bz,az]=impinvar(b,a,Fs) freqs(bz,az); grid on [实验报告] 1、要求有程序并指出各函数功能，有结果，有图形（坐标要标识清楚）。 2、仔细观察上机操作时出现的各种现象，记录主要情况，作出必要说明和分析。