《IIR数字滤波器设计和应用》由会员分享,可在线阅读,更多相关《IIR数字滤波器设计和应用(9页珍藏版)》请在金锄头文库上搜索。
1、IIR数字滤波器设计和应用一、实验目的1、熟悉 IIR数字滤波器的设计方法。2、掌握模拟滤波器的 matlab实现。3、熟悉用脉冲响应不变法设计 IIR数字滤波器的原理与方法。4、熟悉用双线性变换法设计 IIR数字滤波器的原理与方法。5、掌握数字滤波器的计算机仿真方法。6、通过观察对实际心电图信号的滤波作用, 获得数字滤波的感性知识。二、实验内容和要求1、用 Matlab语言分别设计巴特奥斯低通滤波器和切比雪夫低通滤波器,其技术指标为:通带截止频率 ,通带最大衰减 ;)5(2KHzpdBP3阻带起始频率 ,阻带最小衰减 ;10s s0要求:求出他们的零点、极点、阶数、增益等,并画出图形作比较。
2、(1)巴特奥斯低通滤波器:%Design a Butterworth Analog bandpass filter%Desired performents of the filterwp=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);%Outputw=linspace(1,3000,1000)*2*pi;H=freqs(b,a,w
3、);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 =5Fc = 5.0124e+003(2)切比雪夫低通滤波器:%Chebyshev lowpass filter %Desired performentswp=5000*2*pi;ws=10000*2*pi;Rp=3;Rs=30;%Compute Order Nebs=sqrt
4、(10(Rp/10)-1);A=10(Rs/20);Wc=wpWr=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+004N = 42、切比雪夫低通滤波器,其技术指标为:通带截止频率 ,通带最大衰减 ;MHzfp3dBP1.0阻带起始频率 ,阻带最小衰减 ;s12s6要求:求出他们的零点、极点、阶数、增益等,并画出图形作比较。3、用脉冲响应不变法把模拟滤波器 转换为数字滤651)(2ssHa波器 H(z),其中采样周期为 T=0.1s。%Impul
5、se 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)bz, az=residue(R,Pz,K) disp(Use function IMPINVAR)bz1, az1=impinvar(b,a,1/Ts)运算结果Use direct principle:bz = 1.0000 -0.8966az = 1.0000
6、 -1.5595 0.6065Use function “IMPINVAR”:bz1 = 0.1000 -0.0897az1 = 1.0000 -1.5595 0.60654、利用双线性变换法设计一个巴特奥斯低通 IIR数字滤波器。(1) 设计的指标参数为:在通带内频率低于 0.2 时,最大衰减小于 1dB,在阻带内0.3,频率区间上,最小衰减大于15dB.(2) 以 0.02 为采样间隔, 打印出数字滤波器在频率区间0, /2上的幅频响应特性曲线。(3) 用所设计的滤波器对实际心电图信号采样序列(在本实验后面给出)进行仿真滤波处理,并分别打印出滤波前后的心电图信号波形图, 观察总结滤波作用与
7、效果。例如假设最后求出的数字滤波器系统函数 H(z)为: )215.094.1()358.016.)(705.268.1( (7) 22 61 zzzzzzzH )1.).(31 式(zk 式 )( 2.3,21,)2)(21kzCBAHkk式中A=0.090 36 B1=1.2686,C1=-0.7051 B2=1.0106,C2=-0.3583 B3=0.9044,C3=-0.2155 由(1 式)和(2 式)可见, 滤波器 H(z)由三个二阶滤波器 H1(z), H2(z)和 H3(z)级联组成, 如图 1 所示。图 1 滤波器 H(z)的组成x(n) H1(z) H2(z) H3(z)
8、y1(n) y2(n) y3(n)y(n)H(z)(2) 编写滤波器仿真程序, 计算 H(z)对心电图信号采样序列 x(n)的响应序列 y(n)。设 yk(n)为第 k级二阶滤波器 Hk(z)的输出序列, y k-1(n)为输入序列, 如图 1 所示。 由(2 式)可得到差分方程:yk(n)=Ayk-1(n)+2Ayk-1(n-1)+Ayk-1(n-2)+Bkyk(n-1)+Ckyk(n-2) (3式) (3) 在通用计算机上运行仿真滤波程序, 并调用通用绘图子程序, 完成实验内容(2)和(3)。人体心电图信号在测量过程中往往受到工业高频干扰,所以必须经过低通滤波处理后,才能作为判断心脏功能的
9、有用信息。下面给出一实际心电图信号采样序列样本 x(n),其中存在高频干扰。在实验中,以 x(n)作为输入序列,滤除其中的干扰成分。程序:x=-4,-2,0,-4,-6,-4,-2,-4,-6,-6,-4,-4,-6,-6,-2,6,12,8,0,-16,-38,-60,-84,-90,-66,-32,-4,-2,-4,8,12,12,10,6,6,6,4,0,0,0,0,0,-2,-4,0,0,0,-2,-2,0,0,-2,-2,-2,-2,0;k=1;close all;figure(1)subplot(2,2,1);n=0:55;stem(n, x, .);axis(0 56 -100
10、50);hold on;n=0:60;m=zeros(61);plot(n, m);xlabel(n);ylabel(x(n);title(心电图信号采样序列X (n));B =0.09036 2*0.09036 0.09036;A=1.2686 -0.7051;A1= 1.0106 -0.3583;A2=0.9044 -0.2155;while(k=3)y=filter(B, A, x); %The function is to filte() the singal xx=y;if k=2;A=A1;endif k=3A=A2;endk=k+1;endsubplot(2,2,3)n=0:55
11、;stem(n, y, .);axis(0 56 -15 5);hold on;n=0:60;m =zeros(61);plot(n, m);xlabel(n);ylabel(y(n);title(三级滤波后的心电图信号 );%求数字滤波器的幅频特性A= 0.09036, 0.18072, 0.09036; B1=1, -1.2686, 0.7051; B2=1, -1.0106, 0.3583; B3=1, -0.9044, 0.2155; H1, w=freqz(A, B1, 100); H2, w=freqz(A, B2, 100); H3, w=freqz(A, B3 ,100); H
12、4=H1.*(H2); H=H4.*(H3); mag=abs(H); db=20*log10( mag+ eps )/max( mag ); subplot(2,2,2) plot(w/pi, db); axis(0,0.5,-50,10); title(滤波器的幅频响应曲线) ;三、思考题1. 答:二者一致。2.答:(1)优点:一定是递归结构、对频率分量的选择性好(零极点可同时起作用) 、相同性能下阶次较低、运算误差大,有可能出现极限环振荡、设计有封闭形式的公式,一次完成、对计算手段的要求较低;(2)缺点:一般是非线性的、系统不一定稳定、不能用FFT作快速卷积、有噪声反馈,噪声大。3.答:脉冲响应不变法具有较好的线性特性,但会产生频率混叠。双线性变换法靠频率的严重非线性关系得到 s平面与 z平面的单值,固不会产生频率混叠,但其线性关系较差。4.答:T 的值对结果有影响,因为整个 s平面频率轴上的频率范围压缩的/T 之间。