冲击响应不变法设计数字Chebyshev低通滤波器.docx
- 文档编号:8967292
- 上传时间:2023-05-16
- 格式:DOCX
- 页数:14
- 大小:142.10KB
冲击响应不变法设计数字Chebyshev低通滤波器.docx
《冲击响应不变法设计数字Chebyshev低通滤波器.docx》由会员分享,可在线阅读,更多相关《冲击响应不变法设计数字Chebyshev低通滤波器.docx(14页珍藏版)》请在冰点文库上搜索。
冲击响应不变法设计数字Chebyshev低通滤波器
课程设计任务书
学生某:
专业班级:
通信工程0603班
指导教师:
许建霞工作单位:
信息工程学院
题目:
冲击响应不变法设计数字Chebyshev低通滤波器
初始条件:
1.计算机一台
2.MATLAB7.0软件
要求完成的主要任务:
用冲击响应不变法设计数字Chebyshev低通滤波器
1、技术指标:
在通带截止频率
处的衰减为不大于1dB。
在阻带截止频率处
衰减不小于15dB。
2、要求:
绘制该数字滤波器的幅频特性曲线、相位特性曲线、单位冲激响应序列、单位阶跃响应序列。
时间安排:
2009年1月5日至2009年1月15日
指导教师签名:
2008年月日
系主任(或责任教师)签名:
年月日
摘要
切比雪夫滤波器的幅度特性具有等纹波特性。
他有两种类型,一种是通带内为等纹波的,在阻带内是单调的成为切比雪夫Ι型滤波器;一种是通带内单调,阻带内等纹波的,称为切比雪夫ΙΙ型。
用冲击相应不变法设计滤波器要求用模拟系统函数Ha(s)求拉普拉斯反变换得到模拟冲击响应然后抽样得到h(n)=ha(nT),再取Z变换得到H(z)。
Abstract
CharacteristicpropertyhastheextentbeingsuretopetewitheachotherintheXuefuwavefilterwaitingforlineswavecharacteristicproperty.Hehastwokindstypes,akindofisthatbelttheinnerisawavewaitforlines,isthatmonotonousbeingissuretobeparedtotheXuefuIOTAtypewavefilterwithinhinderingbelt;Akindofisthatbelttheinnerismonotonous,hinderthesimultaneousinnerfromwaitingforlineswave's,becalledissuretobeparedtotheXuefuIOTAIOTAtype.GetHwithpoundingcorrespondinginvariablefollowthealternationdesigningthatthewavefilterdemandstouseHa(s)simulatingsystemfunctiontoaskLaplaceonthecontraryalternationtobeimitatedpoundrespondtoandthengettingtheh(n)=ha(nT),takingZagainbysampling(z).
Chebyshev低通滤波器
1设计方法
1.1冲击响应不变法原理(Impulseinvariance)
冲击响应不变法的设计原理是使数字滤波器的单位抽样响应序列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)
上式的物理意义为首先将模拟滤波器的系统函数G(s)作周期的延拓,再经过
的映射变换,从而得到数字滤波器的系统函数H(z)。
假设s平面上,s在j
轴上取值,z在Z平面内的单位圆周
上取值,可以得到数字滤波器的频率响应
和模拟滤波器的频率响应
间的关系为
(1-3)
其中
假设模拟滤波器的系统函数G(s)只有单阶极点,且M (1-4) 其拉普拉斯变换脉冲响应g(t)为: g(t)= ,t≥0 0,t<0(1-5) 对g(t)进行等间隔采样,可以得到数字滤波器的单位取样响应函数h(n)为: 当n≥0时,h(n)= ;当n<0时,h(n)=0(1-6) 然后对h(n)进行Z变换,就可以得到数字滤波器的系统传递函数H(z): (1-7) 按照冲击响应不变法的原理,通过模拟滤波器的系统传递函数G(s),可以直接球的数字滤波器的系统函数H(z), 1.2冲击响应不变法设计滤波器的转换步骤归纳 (1)利用ω=ΩT(可由关系式 推导出),将 , 转换为 , ,而 , 不变; (2)求解低通模拟滤波器的传递函数G(s); (3)根据式1-4和1-7将模拟滤波器的传递函数G(s)转换为数字滤波器的传递函数H(z)。 2切比雪夫(chebyshev)滤波器 2.1切比雪夫滤波器简介 切比雪夫滤波器(又译车比雪夫滤波器)是在通带或阻带上频率响应幅度等波纹波动的滤波器。 在通带波动的为“I型切比雪夫滤波器”,在阻带波动的为“II型切比雪夫滤波器”。 切比雪夫滤波器在过渡带比巴特沃斯滤波器的衰减快,但频率响应的幅频特性不如后者平坦。 切比雪夫滤波器和理想滤波器的频率响应曲线之间的误差最小,但是在通频带内存在幅度波动。 这种滤波器来自切比雪夫多项式,因此得名,用以纪念俄罗斯数学家巴夫尼提·列波维其·切比雪夫 2.2切比雪夫滤波器原理 巴特沃兹滤波器在通带内幅度特性是单调下降的,如果阶次一定,则在靠近截止 处,幅度下降很多,或者说,为了使通带内的衰减足够小,需要的阶次N很高,为了克服这一缺点,采用切比雪夫多项式来逼近所希望的 。 切比雪夫滤波器的 在通带X围内是等幅起伏的,所以在同样的通常内衰减要求下,其阶数较巴特沃兹滤波器要小。 切比雪夫滤波器的振幅平方函数为 (2-1) 式中 Ωc—有效通带截止频率 —与通带波纹有关的参量, 大,波纹大0< <1 VN(x)—N阶切比雪夫多项式 (2-2) |x|≤1时,|VN(x)|≤1 |x|>1时,|x|↗,VN(x)↗ 切比雪夫滤波器的振幅平方特性如图所示,通带内,的变化X围为 1(max)→(min) 时,|x|>1,随 ↗, →0(迅速趋于零) 当 =0时, (2-3) N为偶数,cos2( )=1,得到min, ,(2-4) N为奇数,cos2( ,得到max, (2-5) 切比雪夫滤波器的振幅平方特性如图1所示。 图1切比雪夫滤波器的振幅平方特性 2.3Chebyshev有关参数的确定 2.3.1通带截止频率 预先给定 2.3.2ε的确定 与通带波纹有关的参数,通带波纹表示成 (2-6) 所以, 给定通带波纹值 分贝数后,可求得 。 2.3.3阶数N 由阻带的边界条件确定。 、A2为事先给定的边界条件,即在阻带中的频率点处 要求滤波器频响衰减到1/A2以上。 (2-7) (2-8) (2-9) (2-10) 滤波器阶数N对滤波特性有极大的影响,N越大,逼近特性越好,但是相应的结构也越复杂。 一般情况下N等于通带内最大和最小个数的总和。 N的数值可根据阻带衰减来确定。 3切比雪肤低通滤波器的设计 3.1matlab函数说明 1.cheb1ord.m 求切比雪夫1型滤波器的阶次。 2.cheb1ap.m 用来设计原型切比雪夫1型模拟滤波器。 3.cheby1.m 直接设计切比雪夫1型滤波器。 此函数设计N阶切比雪夫1型滤波器,通带波动为RpdB。 在长度为N+1的矢量b和a中返回滤波器系数。 以上3个文件的调用格式和对应的巴特沃兹滤波器的文件类似。 4.impinvar.m 用冲激响应不变法实现 到 及s到z的转换。 3.2Matlab程序及运行 3.2.1matlab程序 %利用模拟切比雪夫滤波器设计数字滤波器 %冲击响应不变法 %数字滤波器指标 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,OmegaS,1,15,'s');%求出满足设计指标的最小阶数n和截止频率Wn [b,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(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所示。 图2matlab程序窗口 3.2.2编程原理说明 MATLAB提供了一个函数[z,p,k]=cheb1ap(N,Rp),来设计一个阶数为N,通带波动为Rp的归一化切比雪夫1型原型滤波器,它在数组z中返回零点,数组p中返回极点,并返回增益k。 我们需要具有任意 的归一化的切比雪夫1型滤波器,这可由归一化滤波器的数组p乘以 得到。 这种滤波器没有零点。 新增益k由旧的增益k乘以非归一化与归一化多项式在s=0出的比值。 上述程序是根据要求所给的性能指标先用cheb1ap设计一个模拟切比雪夫滤波器,然后用函数impinvar将模拟滤波器转换为数字滤波器。 3.3波形记录及分析 在matlab中运行切比雪夫滤波器的设计程序得到它的幅频特性、相频特性、单位冲激响应、单位阶越响应如图3所示 图3切比雪夫滤波器波形 波形分析: 切比雪夫低通滤波器的幅频特性具有等纹波特性,切比雪夫1型如图所示,器通带内为等纹波的,阻带内是单调的。 冲击响应频率坐标的变换是线性的,如果模拟滤波器的频响带限于折叠频率的话,则通过变换后滤波器的频响可不失真的原响应与频率的关系。 4小结 本次课程设计我的任务是用冲击响应不变法设计切比雪夫低通数字滤波器,并应用MATLAB软件进行仿真分析,是数字信号处理的知识与软件应用相结合的一次训练。 所谓冲击相应不变法设计滤波器就是将模拟滤波器转化为数字滤波器,冲击响应不变法的一个重要特点是频域坐标的变换是线性的,例如,对于线性相位的椭圆滤波器,通过冲击响应不变法得到的仍是线性相位的低通滤波器。 如果Ha(s)是稳定的,即其极点在S左半平面,映射到H(z)也是稳定的。 切比雪夫滤波器包括两种,即切比雪夫1型滤波器和2型滤波器,切比雪夫滤波器频率响应的特点是: 在带通内是等幅的纹波,在阻带内单调衰减,过度迅速。 通过本次课程设计,我对平时所学的知识有了更深入的了解,对MATLAB的应用也较以前熟练,知识与应用相结合提高了我的学习热情,希望在一次次的课程设计中可以不断进步,不断提高。 参考文献 [1]X泉,阙大顺编数字信号处理原理与实现,电子工业,2005年6月 [2]苏金明,王永利编Matlab7.0使用指南,电子工业,2004年11月 [3]薛定宇陈某编基于matlab/simulink的系统仿真技术与应用,清华大学2002年 [4]余成波等编,数字信号处理及MATLAB实现(第二版),: 清华大学,2008 [5]陈怀琛编著,数字信号处理教程——MATLAB释义与实现,: 电子工业2004
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 冲击 响应 变法 设计 数字 Chebyshev 滤波器