《中衡切比雪夫低通滤波器【设计明细】.pdf》由会员分享,可在线阅读,更多相关《中衡切比雪夫低通滤波器【设计明细】.pdf(9页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、摘要摘要随着信息和数字时代的到来,数字信号处理已成为当今一门极其重要的学科和技术领域。在现代通信系统中,由于信号中经常混有各种复杂成分,因此很多信号的处理都是基于滤波器而进行的。所以,数字滤波器在数字信号处理中起着举足轻重的作用。而数字滤波器的设计都要以模拟滤波器为基础的,这是因为模拟滤波器的理论和设计方方法都已发展的相当成熟,且有典型的模拟滤波器供我们选择。,如巴特沃思滤波器、切比雪夫滤波器等。本次课程设计将运用 MATLAB设计一个基于切比雪夫低通滤波器,并出所设计滤波器的幅度及幅度衰减特性。关键词关键词:模拟低通 滤波 切比雪夫1 1 课题描述课题描述数字滤波器是数字信号处理的重要工具之
2、一,它通过数值运算处理改变输入信号所含频率成分的相对比例或者滤出某些频率成分的数字器件或程序,而数字滤波器处理精度高、体积小、稳定、重量轻、灵活、不存在阻抗匹配问题,可以实现模拟滤波器无法实现的特殊功能。故本课题使用 MATLAB信号处理箱和运用切比雪夫法设计数字低通滤波器。2 2 设计原理设计原理2.12.1 切比雪夫滤波器介绍切比雪夫滤波器介绍在巴特沃兹滤波器中,幅度响应在通带和阻带内都是单调的。因此,若滤波器的技术要求是用最大通带和阻带的逼近误差来给出的话,那么,在靠近通带低频端和阻带截止频率以上的部分都会超出技术指标。一种比较有效的途径是使逼近误差均匀地分布于通带或阻带内,或同时在通带
3、和阻带内都均匀分布,这样往往可以降低所要求的滤波器阶次。通过选择一种具有等波纹特性而不是单调特性的逼近方法可以实现这一点。切比雪夫型滤波器就具有这种性质:其频率响应的幅度既可以在通带中是等波纹的,而在阻带中是单调的(称为I 型切比雪夫滤波器),也可以在通带中是单调的,而在阻带中是等波纹的(称为 II 型切比雪夫滤波器)。I 型切比雪夫滤波器的幅度平方函数是12|HC(j)|212CN(/c)(2.1)=式中为 N 阶切比雪夫多项式,定义为CN(x)cos(N cos1x)(2.2)求CN1(x)从定义切比雪夫多项式可以直接得出由公式。将三角恒等式代入(2.2)式,得CN1(x)CN(x)和CN
4、1(x)的递推=2xCN(x)CN1(x)2CN(x)(2.3)从(2.2)式我们注意到,当0 x1 时,|HC(j)|2CN(x)cos1x是虚数,所以像双曲余弦一样单调地增加。参考(2.1),对/p/p于 01 呈现出在 1 和 1/(1)之间的波动;而对于1 单调2地减小。需要用三个参量来确定该滤波器:,p和 N。在典型的设计中,用容许的通带波纹来确定,而用希望的通带截止频率来确定的阶次 N,以便阻带的技术要求得到满足。定义允许的通带最大衰减p用下式表示:为了求切比雪夫滤波器在椭圆上极点的位置,我们首先要这样确定,在大圆和小圆上以等角度等间隔排列的那些点:这些点对于虚轴呈对称分布,并且没
5、有一个点落在虚轴上;但当 N 为奇数时要有一个点落在实轴上,而当 N 为偶数时,就都不会落在实轴上。切比雪夫滤波器的极点落在椭圆上,起纵坐标由相应的大圆上点的纵坐标来表示,起横坐标由相应的小圆上点的横坐标来表示。c。然后选择合适1.11.1 滤波器的分类滤波器的分类(1)从功能上分;低、带、高、带阻。(2)从实现方法上分:FIR、IIR(3)从设计方法上来分:Chebyshev(切比雪夫),Butterworth(巴特沃斯)(4)从处理信号分:经典滤波器、现代滤波器1.21.2 模拟滤波器的设计指标模拟滤波器的设计指标设ha(j)是一个模拟滤波器的频率响应,则基于平方幅度响应J()Ha(j)的
6、低通滤波器技术指标为:0 Ha(j)21,s(2-2-1)A2112 Ha(j)1,p(2-2-2)2其中为通带波动系数,p和s是通带和阻带边缘频率。A 为阻带衰减系数。这些指标如图所示。从图知必须满足1,p221(2-2-3)Ha(j)1 s2A,其中参数和 A 是数字滤波器指标。1.31.3 切比雪夫切比雪夫 I I 型滤波器型滤波器1.3.11.3.1 切比雪夫低通滤波器的设计原理切比雪夫低通滤波器的设计原理切比雪夫滤波器的幅频特性具有等波纹特性。它有两种形式:振幅特性在通带内是等波纹的,在阻带内是单调递减的切比雪夫I 型滤波器,振幅特性在阻带内是等波纹的,在通带内是单调递减的切比雪夫
7、II 型滤波器,如图所示分别画出了滤波器的幅频特性和衰减函数。以切比雪夫 I 型为例介绍其设计原理幅度平方函数用H(j)2表示H(j)21(2-3-1)12CN()p式中,为小于 1 的正数,表示通带内幅度波动的程度,越大,波动幅度也越大。p称为通带截止频率。令,称为对pp的归一化频率。定义允许的通带内最大衰减p用下式表示p10lgmax Ha(j)2min Ha(j)2,p(2-3-2)式中max ha(j)21(2-3-3)min ha(j)因此21(2-3-4)12p10lg(12)(2-3-5)2100.1p1(2-3-6)这样,可以根据通带内最大衰减p,可求出参数。阶数 N 影响过渡
8、带的宽度,同时也影响通带内波动的疏密,因为 N 等于通带内的最大值和最小值的总个数。设阻带的起点频率为s,则有H(js)211CN(s)p22(2-3-7)令ssp,由s1,有11Hajs2CNs chNarchs可以解出1archN 11(2-3-8)Hajsarchs121(2-3-9)1(2-3-10)1 1s pcharchNHajs23dB 截止频率用c表示,2Hajc1(2-3-11)2按照(2-3-1)式,有1 1c pcharch(2-3-12)N经过一系列推论得归一化系统函数为Gap1i2去归一化的系统函数为Hs GapN1p pi1N(2-3-13)psppN2N1p pi
9、pi1N(2-3-14)1.3.21.3.2 切比雪夫低通滤波器的设计步骤切比雪夫低通滤波器的设计步骤(1)确定低通滤波器的技术指标:边带频率p,通带最大衰减p、阻带最大衰减s、阻带截至频率s,它们满足p 10lg1Ha jp2(2-3-15)(2-3-16)2s10lg(2)求滤波器阶数 N 和参数1Ha jss1sp(2-3-17)100.1s1K1(2-3-18)0.1p101archk11N(2-3-19)archs这样,先由(2-18)式求出K11,代入(2-3-19),求出阶数 N,最后取大于或等于 N 的最小整数。(3)求归一化系统函数GapGap1N(2-3-20)2N1p p
10、ii1(4)将去归一化,得到实际的HasHasGpps(2-3-21)p2.3.3 用 MATLAB设计切比雪夫低通滤波器(1)z,p,Gcheb1apN,Rs该格式用于计算N阶切比雪夫I型归一化模拟低通滤波器系统的零极点和增益因子。返回长度为 N 的列向量 Z 和 P,分别给出 N 个零点和极点的位置。Rs是阻带最小衰减。(2)N,wso cheb1ordwp,ws,Rp,As该格式用于计算切比雪夫 I 型数字滤波器的阶数 N 和阻带截止频率 wso。调用参数分别为数字滤波器的通带频率和阻带边界频率的归一化值。(3)N,wsocheb1ordwp,ws,Rp,As,s该格式用于计算切比雪夫
11、I 型模拟滤波器的阶数 N 和阻带截止频率 wso。wp 和 ws 是实际模拟角频率。(4)B,A cheby(N,Rs,wso,ftype)该格式用于计算N阶切比雪夫I型数字滤波器系统函数的分子和分母多项式系数向量 B 和 A。调用参数 N 和 wso 分别为切比雪夫 I 型数字滤波器的阶数和阻带截止频率的归一化值。(5)B,A cheby(N,Rp,wso,ftype,s)该格式用于计算N阶切比雪夫I型模拟滤波器系统函数的分子和分母多项式系数向量 B 和 A。调用参数 N 和 wso 分别为切比雪夫 I 型模拟滤波器的阶数和阻带截止频率的归一化值。2 2总结总结本次课程设计,让我对设计滤波器的原理和步骤有了更深入的理解。滤波器在数字信号处理中占据着重要的地位,因此掌握模拟低通滤波器的设计对以后的学习有着不菲的收获。首先,设计滤波器前一定要选好设计的方法。这次课程设计采用切比雪夫 I 型来设计一个模拟低通滤波器。其次,确定好参数。最后,用MATLAB软件编程式设计过程大为简化。