IIR数字滤波器设计与应用.docVIP

　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　IIR数字滤波器设计和应用 一、实验目的 1、熟悉IIR数字滤波器的设计方法。 2、掌握模拟滤波器的matlab实现。 3、熟悉用脉冲响应不变法设计IIR数字滤波器的原理与方法。 4、熟悉用双线性变换法设计IIR数字滤波器的原理与方法。 5、掌握数字滤波器的计算机仿真方法。 6、通过观察对实际心电图信号的滤波作用， 获得数字滤波的感性知识。 二、实验内容和要求 1、用Matlab语言分别设计巴特奥斯低通滤波器和切比雪夫低通 滤波器，其技术指标为： 通带截止频率，通带最大衰减； 阻带起始频率，阻带最小衰减； 要求：求出他们的零点、极点、阶数、增益等，并画出图形作比 较。 （1）巴特奥斯低通滤波器： %Design a Butterworth Analog bandpass filter %Desired performents of the filter wp=5000*2*pi; ws=10000*2*pi; Rp=3; Rs=30; %compute oder and cutoff frequency [N,Wn]=buttord(wp,ws,Rp,Rs,s); N Fc=Wn/(2*pi) %Compute the filter [b,a]=butter(N,Wn,s); %Output w=linspace(1,3000,1000)*2*pi; H=freqs(b,a,w); magH=abs(H); phaH=unwrap(angle(H)); plot(w/(2*pi),20*log10(magH)); xlabel(Frequency(Hz)); ylabel(Magnidute(dB)); title(Butterworth Analog bandpass filter) grid on 运算结果：N =5 Fc = 5.0124e+003 （2）切比雪夫低通滤波器： %Chebyshev ￠?lowpass filter %Desired performents wp=5000*2*pi; ws=10000*2*pi; Rp=3; Rs=30; %Compute Order N ebs=sqrt(10^(Rp/10)-1); A=10^(Rs/20); Wc=wp Wr=ws/wp; g=sqrt(A*A-1)/ebs; N1=log10(g+sqrt(g*g-1))/log10(Wr+sqrt(Wr*Wr-1)); N=ceil(N1) 运算结果：Wc = 3.1416e+004 　　　　N = 4 2、切比雪夫低通滤波器，其技术指标为： 通带截止频率，通带最大衰减； 阻带起始频率，阻带最小衰减； 要求：求出他们的零点、极点、阶数、增益等，并画出图形作比较。 用脉冲响应不变法把模拟滤波器 转换为数字滤波器H(z)，其中采样周期为T=0.1s。 %Impulse invariance method %From analog to digital filter conversion %Model of analog filter b=[1 1]; a=[1 5 6]; Ts=0.1; [R,Ps, K]=residue(b,a); Pz=exp(Ps*Ts); disp(Use direct principle￡o) [bz, az]=residue(R,Pz,K) disp(Use function ?°IMPINVAR?±￡o) [bz1, az1]=impinvar(b,a,1/Ts) 运算结果 Use direct principle： bz = 1.0000 -0.8966 az = 1.0000 -1.5595 0.6065 Use function “IMPINVAR”： bz1 = 0.1000 -0.0897 az1 = 1.0000 -1.5595 0.6065 4、利用双线性变换法设计一个巴特奥斯低通IIR数字滤波器。 (1) 设计的指标参数为：在通带内频率低于0.2π时，最大衰减小于1dB，在阻带内[0.3π，π]频率区间上，最小衰减大于15dB. (2)