切比雪夫低通滤波器设计(共17页).doc

《切比雪夫低通滤波器设计(共17页).doc》由会员分享，可在线阅读，更多相关《切比雪夫低通滤波器设计(共17页).doc（17页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上课程设计任务书20102011学年第一学期专业： 通信工程 学号： 姓名： 课程设计名称： 数字信号处理课程设计 设计题目： 切比雪夫型低通滤波器的软件实现 完成期限：自 年 月 日至 年 月 日共 周一设计目的1巩固所学的理论知识。2提高综合运用所学理论知识独立分析和解决问题的能力。3更好地将理论与实践相结合。4掌握信号分析与处理的基本方法与实现。5熟练使用MATLAB语言进行编程实现。二设计内容用MATLAB编程设计切比雪夫I型低通滤波器，各参数要求如下：fp=100Hz，p=2dB，fs=120Hz，s=60dB。三、设计要求给出所设计滤波器的幅度及幅度衰减特性

2、并分析是否满足设计需要。四设计条件计算机、MATLAB语言环境五、参考资料1数字信号处理（第三版），丁玉美，高西全.西安电子科技大学出版社，2000.2MATLAB及在电子信息课程中的应用，陈怀堔，吴大正，高西全.电子工业出版社，2006.3MATLAB 7.0从入门到精通，求是科技.人民邮电出版社，2006.4数字信号处理（第三版）学习指导，高西全，丁玉美.西安科技大学出版社，2001.指导教师（签字）： 教研室主任（签字）： 批准日期： 年 月 日摘要随着信息和数字时代的到来，数字信号处理已成为当今一门极其重要的学科和技术领域。在现代通信系统中，由于信号中经常混有各种复杂成分，因此很多信号

3、的处理都是基于滤波器而进行的。所以，数字滤波器在数字信号处理中起着举足轻重的作用。而数字滤波器的设计都要以模拟滤波器为基础的，这是因为模拟滤波器的理论和设计方方法都已发展的相当成熟， 且有典型的模拟滤波器供我们选择。，如巴特沃思滤波器、切比雪夫滤波器等。本次课程设计将运用MATLAB设计一个基于切比雪夫低通滤波器，并出所设计滤波器的幅度及幅度衰减特性。关键词：模拟低通 滤波 切比雪夫目录专心-专注-专业1课题描述数字滤波器是数字信号处理的重要工具之一，它通过数值运算处理改变输入信号所含频率成分的相对比例或者滤出某些频率成分的数字器件或程序，而数字滤波器处理精度高、体积小、稳定、重量轻、灵活、不

4、存在阻抗匹配问题，可以实现模拟滤波器无法实现的特殊功能。故本课题使用MATLAB信号处理箱和运用切比雪夫法设计数字低通滤波器。 2设计原理2.1切比雪夫滤波器介绍在巴特沃兹滤波器中，幅度响应在通带和阻带内都是单调的。因此，若滤波器的技术要求是用最大通带和阻带的逼近误差来给出的话，那么，在靠近通带低频端和阻带截止频率以上的部分都会超出技术指标。一种比较有效的途径是使逼近误差均匀地分布于通带或阻带内，或同时在通带和阻带内都均匀分布，这样往往可以降低所要求的滤波器阶次。通过选择一种具有等波纹特性而不是单调特性的逼近方法可以实现这一点。切比雪夫型滤波器就具有这种性质：其频率响应的幅度既可以在通带中是等

5、波纹的，而在阻带中是单调的（称为I型切比雪夫滤波器），也可以在通带中是单调的，而在阻带中是等波纹的（称为II型切比雪夫滤波器）。I型切比雪夫滤波器的幅度平方函数是 = (2.1)式中为N阶切比雪夫多项式，定义为 (2.2)从定义切比雪夫多项式可以直接得出由和求的递推公式。将三角恒等式代入 (2.2)式，得 =2x (2.3)从 (2.2)式我们注意到，当0x1时，是虚数，所以像双曲余弦一样单调地增加。参考(2.1)，对于01呈现出在1和1/（）之间的波动；而对于1单调地减小。需要用三个参量来确定该滤波器：，和N。在典型的设计中，用容许的通带波纹来确定，而用希望的通带截止频率来确定。然后选择合适

6、的阶次N，以便阻带的技术要求得到满足。定义允许的通带最大衰减用下式表示： 为了求切比雪夫滤波器在椭圆上极点的位置，我们首先要这样确定，在大圆和小圆上以等角度 等间隔排列的那些点：这些点对于虚轴呈对称分布，并且没有一个点落在虚轴上；但当N为奇数时要有一个点落在实轴上，而当N为偶数时，就都不会落在实轴上。切比雪夫滤波器的极点落在椭圆上，起纵坐标由相应的大圆上点的纵坐标来表示，起横坐标由相应的小圆上点的横坐标来表示。1.1 滤波器的分类（1）从功能上分；低、带、高、带阻。 （2）从实现方法上分:FIR、IIR （3）从设计方法上来分：Chebyshev(切比雪夫）,Butterworth（巴特沃斯）

7、 （4）从处理信号分：经典滤波器、现代滤波器 1.2 模拟滤波器的设计指标设是一个模拟滤波器的频率响应，则基于平方幅度响应的低通滤波器技术指标为： （2-2-1） （2-2-2）其中为通带波动系数，和是通带和阻带边缘频率。A为阻带衰减系数。这些指标如图所示。从图知必须满足 （2-2-3）其中参数和A是数字滤波器指标。1.3 切比雪夫I型滤波器1.3.1 切比雪夫低通滤波器的设计原理切比雪夫滤波器的幅频特性具有等波纹特性。它有两种形式：振幅特性在通带内是等波纹的，在阻带内是单调递减的切比雪夫I型滤波器，振幅特性在阻带内是等波纹的，在通带内是单调递减的切比雪夫II型滤波器，如图所示分别画出了滤波器

8、的幅频特性和衰减函数。以切比雪夫I型为例介绍其设计原理幅度平方函数用表示 （2-3-1）式中，为小于1的正数，表示通带内幅度波动的程度，越大，波动幅度也越大。称为通带截止频率。令 ，称为对的归一化频率。定义允许的通带内最大衰减用下式表示 （2-3-2） 式中 （2-3-3） （2-3-4） 因此 （2-3-5） （2-3-6） 这样，可以根据通带内最大衰减，可求出参数。 阶数N影响过渡带的宽度，同时也影响通带内波动的疏密，因为N等于通带内的最大值和最小值的总个数。设阻带的起点频率为，则有 （2-3-7） 令，由1，有 （2-3-8） 可以解出 （2-3-9） （2-3-10）3dB截止频率用表

9、示， （2-3-11） 按照（2-3-1）式，有 （2-3-12） 经过一系列推论得归一化系统函数为（2-3-13）去归一化的系统函数为 （2-3-14） 1.3.2 切比雪夫低通滤波器的设计步骤（1）确定低通滤波器的技术指标：边带频率，通带最大衰减、阻带最大衰减、阻带截至频率，它们满足 （2-3-15）（2-3-16）（2）求滤波器阶数N和参数 （2-3-17） （2-3-18） （2-3-19）这样，先由（2-18）式求出，代入（2-3-19），求出阶数N，最后取大于或等于N的最小整数。（3）求归一化系统函数 （2-3-20）（4）将去归一化，得到实际的 （2-3-21）2.3.3用MAT

10、LAB设计切比雪夫低通滤波器 （1）该格式用于计算N阶切比雪夫I型归一化模拟低通滤波器系统的零极点和增益因子。返回长度为N的列向量Z和P，分别给出N个零点和极点的位置。Rs是阻带最小衰减。（2）该格式用于计算切比雪夫I型数字滤波器的阶数N和阻带截止频率wso。调用参数分别为数字滤波器的通带频率和阻带边界频率的归一化值。（3） 该格式用于计算切比雪夫I型模拟滤波器的阶数N和阻带截止频率wso。wp和ws是实际模拟角频率。（4） 该格式用于计算N阶切比雪夫I型数字滤波器系统函数的分子和分母多项式系数向量B和A。调用参数N和wso分别为切比雪夫I型数字滤波器的阶数和阻带截止频率的归一化值。（5） 该

11、格式用于计算N阶切比雪夫I型模拟滤波器系统函数的分子和分母多项式系数向量B和A。调用参数N和wso分别为切比雪夫I型模拟滤波器的阶数和阻带截止频率的归一化值。 2 设计内容2.1 用MATLAB编程实现编程如下：wp=2*pi*100;ws=2*pi*120;Rp=2;As=60;%设置指标参数N1,wp1=cheb1ord(wp,ws,Rp,As,s)%计算切比雪夫I型模拟低通滤波器的阶数和通带边%界频率B1,A1=cheby1(N1,Rp,wp1,s); %计算切比雪夫I型模拟低通滤波器系统函数系数fk=0:120;wk=2*pi*fk;Hk=freqs(B1,A1,wk);%模拟低通滤波

12、器的系统函数subplot(2,2,1);plot(fk/10,abs(Hk);%画出幅频特性xlabel(频率(Hz);ylabel(幅度);subplot(2,2,2);plot(fk/10,20*log10(abs(Hk);grid on;%画出衰减函数xlabel(频率(Hz);ylabel(幅度(dB);axis(0,15,-50,5);图形如下：2.2 设计结果分析（1）求得阶数N1 =13，通带边界频率为wp1=628.3185，表明阶数越大，模拟低通滤波器越理想，且接近矩形。（2）从图上可以看出：切比雪夫I型滤波器是在通带频率响应幅度等波纹的滤波器，通带内具有等波纹起伏特性，阻

13、带内单调下降且衰减更大；幅度在低频处通过，随着频率的增加，逐渐衰减至零；综上所看：本设计符合要求。3 总结本次课程设计，让我对设计滤波器的原理和步骤有了更深入的理解。滤波器在数字信号处理中占据着重要的地位，因此掌握模拟低通滤波器的设计对以后的学习有着不菲的收获。首先，设计滤波器前一定要选好设计的方法。这次课程设计采用切比雪夫I型来设计一个模拟低通滤波器。其次，确定好参数。最后，用MATLAB软件编程式设计过程大为简化。4 参考文献1数字信号处理（第三版），丁玉美，高西全.西安电子科技大学出版社，2000.2MATLAB及在电子信息课程中的应用，陈怀堔，吴大正，高西全.电子工业出版社，2006.

14、3MATLAB 7.0从入门到精通，求是科技.人民邮电出版社，2006.4数字信号处理（第三版）学习指导，高西全，丁玉美.西安科技大学出版社，2001.T=0.1;Fs=1/T;fp=3;fs=4;Rp=1;As=15;wp=fp/Fs*2*pi;ws=fs/Fs*2*pi;Wp=(2/T)*tan(wp/2);Ws=(2/T)*tan(ws/2);%设计切比雪夫低通滤波器原型ep=sqrt(10(Rp/10)-1);A=10(As/20);Wc=Wp;Wr=Ws/Wp;g=sqrt(A*A-1/ep);N=ceil(log10(g+sqrt(g*g-1)/log10(Wr+sqrt(Wr*(Wr-1);%获取零极点参数z,p,k=cheb1ap(N,Rp);p=p*Wc;a=real(poly(p);aNu=a(N+1);k=k*aNu/a(N);b0=k;B=real(poly(z);b=k*B;%双线性变换bz,az=bilinear(b,a,Fs);%绘图H=freqz(bz,az,200,whole);plot(abs(H);