冲击响应不变法设计数字Cheshev低通滤波器

1、课程设计任务书课程设计任务书学生姓名：学生姓名： 专业班级：专业班级： 通信工程通信工程 06030603 班班 指导教师：指导教师： 许建霞许建霞 工作单位：工作单位： 信息工程学院信息工程学院 题题 目目: : 冲击响应不变法设计数字 chebyshev 低通滤波器 初始条件：初始条件：1. 计算机一台2. matlab 7.0 软件要求完成的主要任务要求完成的主要任务: : 用冲击响应不变法设计数字 chebyshev 低通滤波器1、技术指标：在通带截止频率处的衰减为不大于 1db。2 . 0p在阻带截止频率处衰减不小于 15db。3 . 0t2、要求：绘制该数字滤波器的幅频特性曲线、相

2、位特性曲线 、单位冲激响应序列、单位阶跃响应序列。时间安排：时间安排：2009 年 1 月 5 日至 2009 年 1 月 15 日指导教师签名：指导教师签名： 20082008 年年 月月 日日系主任（或责任教师）签名：系主任（或责任教师）签名： 年年 月月 日日目录目录摘要 .1abstract .21 设计方法.31.1 冲击响应不变法原理（impulse invariance） .31.2 冲击响应不变法设计滤波器的转换步骤归纳.42 切比雪夫（chebyshev）滤波器.42.1 切比雪夫滤波器简介 .42.2 切比雪夫滤波器原理.52.3 chebyshev有关参数的确定 .72.

3、3.1 通带截止频率.72.3.2 的确定 .72.3.3 阶数 n .73 切比雪肤低通滤波器的设计.83.1 matlab函数说明.83.2 matlab程序及运行.83.2.1 matlab 程序 .83.2.2 编程原理说明.103.3 波形记录及分析.114 小结 .12参考文献 .13摘要切比雪夫滤波器的幅度特性具有等纹波特性。他有两种类型，一种是通带内为等纹波的，在阻带内是单调的成为切比雪夫 型滤波器；一种是通带内单调，阻带内等纹波的，称为切比雪夫 型。用冲击相应不变法设计滤波器要求用模拟系统函数 ha（s）求拉普拉斯反变换得到模拟冲击响应然后抽样得到 h(n)=ha(nt),再

4、取 z 变换得到 h（z） 。abstractcharacteristic property has the extent being sure to compete with each other in the xuefu wave filter waiting for lines wave characteristic property. he has two kinds types , a kind of is that belt the inner is a wave wait for lines , is that monotonous becoming is sure to be

5、compared to the xuefu iota type wave filter within hindering belt; a kind of is that belt the inner is monotonous , hinder the simultaneous inner from waiting for lines waves, be called is sure to be compared to the xuefu iotaiota type. get h with pounding corresponding invariable follow the alterna

6、tion designing that the wave filter demands to use ha (s) simulating system function to ask laplace on the contrary alternation to be imitated pound respond to and then getting the h (n) = ha (nt) , taking z again by sampling (z).chebyshevchebyshev 低通滤波器低通滤波器1 设计方法1.1 冲击响应不变法原理（impulse invariance）冲击

7、响应不变法的设计原理是使数字滤波器的单位抽样响应序列 h（n） ，模仿模拟滤波器的冲击响应 g（t） 。设系统传递函数为 g（s）的模拟滤波器的单位冲击响应 g（t） ，并将冲击响应 g（t）进行等间隔采样，使得数字滤波器的单位抽样响应 h（t）刚好等于g（t）的采样值，即： （1-1）其中 ts 为采样周期。因为 g(s)是模拟滤波器的系统传递函数，故他是该系统冲击响应函数g（t）的拉普拉斯变换；又设 h（z）是数字滤波器的系统传递函数，从而可的它是数字滤波器的单位抽样响应函数 h（n）的 z 变换。模拟信号的拉普拉斯变换与其采样序列 z 变换的关系为： （1-2） 上式的物理意义为首先将模

8、拟滤波器的系统函数 g（s）作周期的延拓，再经过的映射变换，从而得到数字滤波器的系统函数 h（z） 。假设 s 平面esstz 上，s 在 j轴上取值，z 在 z 平面内的单位圆周上取值，可以得到数字滤ej波器的频率响应和模拟滤波器的频率响应间的关系为)(ejh)(jg （1-3）其中st假设模拟滤波器的系统函数 g（s）只有单阶极点，且 mn，系统传递函数可以用部分分式表示：)(1)(ksjjkjgtshe)(1)(kssezjsgtzhst0)()()()()(nssstnthntttgeztgnh （1-4） 其拉普拉斯变换脉冲响应 g（t）为：g（t）= ，t0stnkkea10 ，

9、t0 （1-5） 对 g（t）进行等间隔采样，可以得到数字滤波器的单位取样响应函数h（n）为：当 n0 时，h(n)=；当 n0 时，h(n)=0 （1-sntnkkea16）然后对 h(n)进行 z 变换，就可以得到数字滤波器的系统传递函数 h（z）： （1-7） 按照冲击响应不变法的原理，通过模拟滤波器的系统传递函数 g（s） ，可以直接球的数字滤波器的系统函数 h（z） ，1.2 冲击响应不变法设计滤波器的转换步骤归纳（1） 利用 =t（可由关系式推导出） ，将，转换为，stez psps而，不变；ps（2） 求解低通模拟滤波器的传递函数 g（s） ；（3）根据式 1-4 和 1-7 将

10、模拟滤波器的传递函数 g（s）转换为数字滤波器的传递函数 h（z） 。nkkkssasg1)(nktsknntsnnnkknnzeazeaznhzhkk11011)()(2 切比雪夫（chebyshev）滤波器 2.1 切比雪夫滤波器简介切比雪夫滤波器（又译车比雪夫滤波器）是在通带或阻带上频率响应幅度等波纹波动的滤波器。在通带波动的为 “i 型切比雪夫滤波器 ”，在阻带波动的为“ii 型切比雪夫滤波器 ”。切比雪夫滤波器在过渡带比巴特沃斯滤波器的衰减快，但频率响应的幅频特性不如后者平坦。切比雪夫滤波器和理想滤波器的频率响应曲线之间的误差最小，但是在通频带内存在幅度波动。这种滤波器来自切比雪夫多

11、项式，因此得名，用以纪念俄罗斯数学家巴夫尼提列波维其切比雪夫2.2 切比雪夫滤波器原理 巴特沃兹滤波器在通带内幅度特性是单调下降的，如果阶次一定，则在靠近截止 处，幅度下降很多，或者说，为了使通带内的衰减足够小，需要的阶次 n 很高，为了克服这一缺点，采用切比雪夫多项式来逼近所希望的 。切比雪夫滤波器的 在通带范围内是等幅起伏的，所以在同样的通常内衰减要求下，其阶数较巴特沃兹滤波器要小。 切比雪夫滤波器的振幅平方函数为 （2-1） 式中 c有效通带截止频率 与通带波纹有关的参量， 大,波纹大 0 1 时, |x|, vn(x) 切比雪夫滤波器的振幅平方特性如图所示,通带内， 的变化范围为1(m

12、ax) (min)时，|x|1,随 ， 0 (迅速趋于零)当 =0 时， （2-3） n 为偶数，cos2( )=1，得到 min， （2-4）n 为奇数，cos2( ，得到 max， （2-5） 切比雪夫滤波器的振幅平方特性如图 1 所示。图 1 切比雪夫滤波器的振幅平方特性2.3 chebyshev 有关参数的确定 2.3.1 通带截止频率预先给定 2.3.2 的确定与通带波纹有关的参数,通带波纹表示成 （2-6）所以, , 给定通带波纹值 分贝数后，可求得 。 2.3.3 阶数 n由阻带的边界条件确定。 、a2为事先给定的边界条件,即在阻带中的频率点处 ,要求滤波器频响衰减到 1/a2以

13、上。 （2-7）（2-8）（2-9）（2-10）滤波器阶数 n 对滤波特性有极大的影响，n 越大，逼近特性越好，但是相应的结构也越复杂。一般情况下 n 等于通带内最大和最小个数的总和。n 的数值可根据阻带衰减来确定。3 切比雪肤低通滤波器的设计3.1 matlab 函数说明1. cheb1ord.m求切比雪夫 1 型滤波器的阶次。2cheb1ap.m用来设计原型切比雪夫 1 型模拟滤波器。3cheby1.m 直接设计切比雪夫 1 型滤波器。此函数设计 n 阶切比雪夫 1 型滤波器，通带波动为 rpdb。在长度为 n+1 的矢量 b 和 a 中返回滤波器系数。以上 3 个文件的调用格式和对应的巴

14、特沃兹滤波器的文件类似。4impinvar.m用冲激响应不变法实现到及 s 到 z 的转换。3.2 matlab 程序及运行3.2.1 matlab 程序%利用模拟切比雪夫滤波器设计数字滤波器%冲击响应不变法%数字滤波器指标wp=0.2*pi;%通带截止频率ws=0.3*pi;%阻带截止频率rp=1;%通带最大衰减as=15;%阻带最大衰减t=1; %设定周期为 1s%将数字滤波器指标反转换为模拟滤波器的参数%性能指标rip=10(-rp/20);atn=10(-as/20);omgp=wp*t;omgs=ws*t;%切比雪夫 1 型模拟滤波器的设计n,wn=cheb1ord(omegap,o

15、megas,1,15,s); %求出满足设计指标的最小阶数n 和截止频率 wnb,a=cheby1(n,1,wn,low,s); %chebyshev 模拟低通滤波器的生成freqs(b,a); %设计模拟的频率响应bz,az=impinvar(b,a,t); %用脉冲响应不变法映射为数字的 h,w=freqz(bz,az,512,t); %设计数字的频率响应ma=20*log10(abs(h),pha=20*log10(unwrap(angle(h),hi=impz(bz,az);ni=step(bz,az);%幅度，相位，冲击响应，阶跃相应%图形的生成subplot(2,2,1),plot

16、(w,ma);title(幅频特性(db);xlabel(w(/pi);ylabel(db);subplot(2,2,2),plot(w,pha);title(相频特性);xlabel(w(/pi);ylabel(pha(/pi);subplot(2,2,3),plot(hi);title(单位冲击响应);xlabel(n);ylabel(h(n);subplot(2,2,4),plot(ni); title(单位阶越响应);xlabel(n);ylabel(h(n);bz,az; %加坐标网格;程序运行窗口如图 2 所示。 图 2 matlab 程序窗口3.2.2 编程原理说明matlab

17、提供了一个函数z，p，k=cheb1ap（n，rp） ，来设计一个阶数为n，通带波动为 rp 的归一化切比雪夫 1 型原型滤波器，它在数组 z 中返回零点，数组 p 中返回极点，并返回增益 k。我们需要具有任意的归一化的切比雪夫c1 型滤波器，这可由归一化滤波器的数组 p 乘以得到。这种滤波器没有零点。c新增益 k 由旧的增益 k 乘以非归一化与归一化多项式在 s=0 出的比值。上述程序是根据要求所给的性能指标先用 cheb1ap 设计一个模拟切比雪夫滤波器，然后用函数 impinvar 将模拟滤波器转换为数字滤波器。3.3 波形记录及分析在 matlab 中运行切比雪夫滤波器的设计程序得到它

18、的幅频特性、相频特性、单位冲激响应、单位阶越响应如图 3 所示图 3 切比雪夫滤波器波形波形分析波形分析：切比雪夫低通滤波器的幅频特性具有等纹波特性，切比雪夫 1型如图所示，器通带内为等纹波的，阻带内是单调的。冲击响应频率坐标的变换是线性的，如果模拟滤波器的频响带限于折叠频率的话，则通过变换后滤波器的频响可不失真的原响应与频率的关系。4 4 小结本次课程设计我的任务是用冲击响应不变法设计切比雪夫低通数字滤波器，并应用 matlab 软件进行仿真分析，是数字信号处理的知识与软件应用相结合的一次训练。所谓冲击相应不变法设计滤波器就是将模拟滤波器转化为数字滤波器，冲击响应不变法的一个重要特点是频域坐标的变换是线性的，例如，对于线性相位的椭圆滤波器，通过冲击响应不变法得到的仍是线性相位的低通滤波器。如果 ha（s）是稳定的，即其极点在 s 左半平面，映射到 h（z）也是稳定的。切比雪夫滤波器包括两种，即切比雪夫 1 型滤波器和 2 型滤波器，切比雪夫滤波器频率响应的特点是：在带通内是等幅的纹波，在阻带内单调衰减，过度迅速。通过本次课程设计，我对平时所学的知识有了更深入的