实验三用双线性变换法设计iir数字滤波器

《实验三用双线性变换法设计iir数字滤波器》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、用双线性变换法设计IIR数字滤波器实验三用双线性变换法设计IIR数字滤波器一、实验目的1、熟悉用双线性变换法设计IIR数字滤波器的原理与方法；2、掌握数字滤波器的计算机仿真方法；3、通过观察对实际心电图信号的滤波作用，获得数字滤波的感性知识。二、实验内容及原理1、用双线性变换法设计一个巴特沃斯低通IIR数字滤波器。设计指标参数为：在通带内截止频率低于0.2时，最大衰减小于1dB；在阻带内[0.3,]频率区间上，最小衰减大于15dB。2、以0.02为采样间隔，打印出数字滤波器在频率区间[0，/2]上的幅频响应特性曲线。3、用所设计的滤波器对实际心电

2、图信号采样序列进行仿真滤波处理，并分别打印出滤波前后的心电图信号波形图，观察总结滤波作用与效果。教材例中已求出满足本实验要求的数字滤波系统函数：,式中A=0.09036,三、实验结果心电图信号采样序列-4-用双线性变换法设计IIR数字滤波器一级滤波后的心电图信号：二级滤波后的心电图信号：三级滤波后的心电图信号：-4-用双线性变换法设计IIR数字滤波器滤波器的幅频响应曲线：四、实验代码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

3、,-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];n=0:55;subplot(1,1,1);stem(n,x,'.');axis([055-10050]);xlabel('n');ylabel('x(n)');title('心电图信号采样序列x(n)');N=56;A=[0.090362*0.090360.09036];B=[1-1.26860.7051];B1=[1-1.01060.3583];B2=[1-0.90440.2155];y1=f

4、ilter(A,B,x);n=0:55;figure;subplot(1,1,1);stem(n,y1,'.');xlabel('n');ylabel('y1(n)');title('一级滤波后的心电图信号');y2=filter(A,B1,y1);n=0:55;figure;subplot(1,1,1);-4-用双线性变换法设计IIR数字滤波器stem(n,y2,'.');xlabel('n');ylabel('y2(n)');title('二级滤波后的心电图信号');y3=filter(A,B2,y2);n=0:55;figure;subpl

5、ot(1,1,1);stem(n,y3,'.');xlabel('n');ylabel('y3(n)');title('三级滤波后的心电图信号');A=[0.090362*0.090360.09036];B1=[1-1.26860.7051];B2=[1-1.01060.3583];B3=[1-0.90440.2155];[H1,w]=freqz(A,B1,100);[H2,w]=freqz(A,B2,100);[H3,w]=freqz(A,B3,100);H4=H1.*(H2);H=H4.*(H3);mag=abs(H);db=20*log10

6、((mag+eps)/max(mag));figure;subplot(1,1,1);plot(w/pi,db);axis([00.5-5010]);xlabel('w/pi');ylabel('20lg

7、H(jw)

8、');title('滤波器的幅频响应曲线');五、实验总结双线性变换法的特点：对频率的压缩符合下列公式：-4-用双线性变换法设计IIR数字滤波器这样的变换叫做双线性变换。用双线性变换法来设计数字滤波器，由于从s面映射到s1面具有非线性频率压缩的特点，因此不可能产生频率混叠现象，而且转换成的H(z)是因果稳定的，这是双线性变换法的最大

9、优点。其缺点是w与之间的非线性关系直接影响数字滤波器频香逼真的模仿模拟滤波器的频响。数字滤波器的输入和输出均为数字信号，通过一定的运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分。数字滤波器可以通过模拟其网络传输函数进行实现。如图中所示，滤波器对其高于截止频率的频段产生很高的衰减，所得信号较之原信号剔除了高频的成分。-4-