用双线性变换法设计IIR数字低通滤波器.docx

1、用双线性变换法设计IIR数字低通滤波器郑州轻工业学院本科数字信号处理课程设计总结报告设计题目：用双线性变换法设计IIR数字低通滤波器学生姓名：杨彦英、龙芳、何秋慧、赵玉春系 别：计算机与通信工程学院信息与通信工程系专 业：通信工程班 级：08级1班学 号：58号指导教师：杨永双 2010年12月25日郑州轻工业学院课 程 设 计 任 务 书 题目：用双线性变换法设计IIR数字低通滤波器专业、班级 通信工程08-1学号 58姓名 赵玉春主要内容、基本要求、主要参考资料等：1、主要内容根据技术指标求出系统函数Ha(s)；利用s和z之间的转换关系求出函数H(z)，对H(z)进行频谱分析。2、基本要求

2、设计技术指标合适，计算结果正确，结构设计合理，实验数据可靠，程序运行良好。3、主要参考资料1丁玉美，高西全编著数字信号处理西安：西安电子科技大学出版社，20062郑君里等编信号与系统北京：高等教育出版社，1981.3程佩青编著数字信号处理教程北京：清华大学出版社，20074 郑阿奇等编MATLAB实用教程北京：电子工业出版社，20045杨永双等编.数字信号处理实验指导书.郑州：郑州轻工业学院，2007.完 成 期 限： 2010年12月25日 指导教师签名： 杨 永 双 课程负责人签名： 杨 永 双 1、设计题目12、设计内容13、设计思路14、设计原理25、设计原程序26、运行结果37、实验

3、体会58、参考文献6一：设计题目： 用双线性变换法设计IIR数字低通滤波器二：设计内容：1）编制一个用双线性变换法设计IIR数字低通滤波器的源程序；2）根据采样周期、通带和阻带临界频率以及相应的衰减等参数在程序运行时输入，计算阶数N、传递函数Ha(s)；输出Ha(s)分子分母系数；3）根据s和z之间的转换关系，利用系统函数Ha(s)的分子分母系数求出函数H(z)的分子分母系数；4）根据函数H(z)的分子分母系数绘制相应的幅频特性曲线。三：设计思路：1) 确定数字滤波器的性能指标：通带、阻带临界角频率Wp、Ws；通带内的最大衰减A；阻带最小衰减B;采样周期T。2) 计算经过频率预畸的相应参考模拟

4、低通原型的频率，Sp=(2/T)*tan(1/2*Wp)；Ss=(2/T)*tan(1/2*Ws)。3) 计算低通原型阶数N;计算3dB归一化频率Wc，从而求得低通原型的传递函数Ha(s)的分子分母系数.4) 根据s和z之间的转换关系：s=(1-1/z)/(1+1/z)，或者z=(2/T+s)./(2/T-s)，求出数字滤波器传递函数H(z)的分子分母系数。5) 分析滤波器的频谱特性，检查其指标是否满足要求。四：设计原理： 双线性变换法的设计准则是使数字滤波器的频率响应与参考模拟滤波器的频率响应相似。双线性变换法主要依靠双线性变换式：s=(1-1/z)/(1+1/z) z=(1+s)/(1-s

5、) 其中s=a+js, z=r*exp(j*w) 建立起s平面和z平面的单值映射关系，数字频域和模拟频域之间的关系；Ss=tan(1/2*Ws) Ws=arctan(Ss) 由上面的关系可知，当Ws时，Ws终止在折叠频率Ws=pi处，整个js轴单值地对应于单位圆的一周。从上式可以看出双线性变换法不存在频谱混淆的问题。同时还可以知道两者之间的关系不是线性关系。这种非线性关系使得通带截止频率、过渡带的百鸟园频率的相对应位置都发生了非线性畸变。用双线性变换法设计数字滤波器时，一般总是先将数字滤波器各个临界频率经过Ss=tan(1/2*Ws) Ws=arctan(Ss)的频率预畸，求得相应的参考模拟滤

6、波器的各临界频率，然后设计参考模拟滤波器的传递参数，最后通过双线性变换式求得数字滤波器的传递函数。五：设计源程序：Wp=0.2*pi;%数字通带截止频率Ws=0.35*pi;%数字阻带截止频率A=1;%通带最大衰减B=10;%阻带最小衰减T=1;Sp=(2/T)*tan(1/2*Wp);%技术指标数字转换为模拟Ss=(2/T)*tan(1/2*Ws); %技术指标数字转换为模拟Msp=Ss/SpKsp=sqrt(10(B/10)-1)/(10(A/10)-1)N1=log(Ksp)/log(Msp) %求阶数N1N=ceil(log(Ksp)/log(Msp) %对N1向上取整Wc=Wp*(1

7、0(A/10)-1)(-1/2*N) %求3dB截止频率Wcfor K=0:N-1 Sk=Wc*exp(j*pi.*(0.5+(2*(K+1)./(2*N) Zk=(2/T+Sk)./(2/T-Sk) C=real(poly(Sk) D=real(poly(Zk)endnum=(Wc).N; %模拟滤波器的系统函数的分子系数den=C; %模拟滤波器的系统函数的分母的系数den1=D; %数字滤波器的系统函数的分母的系数freqs(num,den,512)freqz(num,den1,512)六：运行结果：Msp =1.8860Ksp =5.8957 N1 = 2.7964N =3Wc =4.

8、7689Sk =-4.1300 + 2.3844iZk = -0.4332 + 0.2205iC =-0.4332 + 0.2205i D =1.0000 0.4332 Sk = -4.1300 - 2.3844iZk =-0.4332 - 0.2205iC =1.0000 4.1300D =1.0000 0.4332 Sk =-0.0000 - 4.7689iZk =-0.7008 - 0.7133iC = 1.0000 0.0000 D =1.0000 0.7008七：实验体会： 在本次课程设计中，经过多次的运用matlab编程，加深了对matlab语言的认识。同时通过通过设计双线性变换法

9、球IIR数字低通滤波器的设计又一次的加深了对以前知识的理解。以前在老师您将这方面的知识时听的不是很懂，经过这次的课程设计，自己在这方面的认识又进入了一个新的境界！同时经过这次的课程设计，让我深刻的认识到在做课程设计时，需要的是团队精神，光靠一个人是很难完成这项任务。因为在写编程时，有时候一个小小的错误一个人是很难发现的。大家一起看时很容易的就发现了！通过我们大家的共同努力，我们很快的完成了本次的课程设计，并且验收合格。总之通过本次的课程设计让我收获非常的大，不仅仅是知识层面的。同时还有合作精神！八、参考文献 1丁玉美，高西全编著数字信号处理西安：西安电子科技大学出版社，20062郑君里等编信号与系统北京：高等教育出版社，1981.3程佩青编著数字信号处理教程北京：清华大学出版社，20074 郑阿奇等编MATLAB实用教程北京：电子工业出版社，20045杨永双等编.数字信号处理实验指导书.郑州：郑州轻工业学院，2007.