Tx Tx Tx Tx Tx Tx Tx Tx Tx Tx Tx Tx Tx Tx Tx Tx Tx Tx Tx实践教学*1* *1* *1* *1* *1* *1* *1* *1* *1* *1* *1* *1* *1* *1* *1* *1* *1* *1* *1*Tx Tx Tx Tx Tx Tx Tx Tx Tx Tx Tx Tx Tx Tx Tx Tx Tx Tx Tx兰州理工大学计算机与通信学院2013 年春季学期信号处理课程设计题目：切比雪夫I型低通滤波器设计专业班级：通信工程三班姓 名：学 号：指导教师：成 绩：摘要本次课程设计将完成一个数字切比雪夫低通 IIR 滤波器的设计，利用双 线性变换和冲激响应不变法完成设计，并利用 MATLAB 进行仿真。已知数字滤波器的性能指标为：通带截止频率为:=0.4k,通带波动为R = 1dB,=0.45k,阻带波动为R = 15dB，要求设计 p p s P满足以上技术指标的切比雪夫I型低通滤波器。绘制出理想冲激响应和实际冲激 响应结果图。并且给出幅度响应结果图。关键字：数字滤波器 切比雪夫 双线性变换 冲激响应不变目录前言 1一数字滤波器 211 数字滤波器的概念 212数字滤波器的分类 21. 3 IIR数字滤波器设计原理3二切比雪夫滤波器 5三. 双线性变换法 8四. 脉冲响应不变法 12五. 切比雪夫低通滤波器的设计 155.1 程序流程图 155.2 设计步骤 15六. 总结 18七. 参考文献 19致谢 20附录 21前言随着信息时代和数字世界的到来，数字信号处理已成为当今一门极其重要的 学科和技术领域。目前数字信号处理在通信、语音、图像、自动控制、雷达、军 事、航空航天、医疗和家用电器等众多领域得到了广泛的应用。在数字信号处理 中起着重要的作用并已获得广泛应用的是数字滤波器(DF, Digital Filter),根据 其单位冲激响应函数的时域特性可分为两类：无限冲激响应IIR(Infinite Impulse Response)滤波器和有限冲激响应FIR(Finite Impulse Response)滤波器。与FIR 滤波器相比，IIR的实现采用的是递归结构，极点须在单位圆内，在相同设计指 标下,实现 IIR 滤波器的阶次较低,即所用的存储单元少,从而经济效率高。 MATLAB是英文 MATrix LABoratory(矩阵实验室)的缩写。它是美国的 MathWorks 公司推出的一套用于科学计算和图形处理可视化、高性能语言与软件 环境。MATLAB的信号处理工具箱是专门应用于信号处理领域的专用工具箱，它 的两个基本组成就是滤波器的设计与实现部分以及谱分析部分。工具箱提供了丰 富而简便的设计，使原来繁琐的程序设计简化成函数的调用。只要以正确的指标 参数调用相应的滤波器设计程序或工具箱函数，便可以得到正确的设计结果，使 用非常方便。一数字滤波器11 数字滤波器的概念滤波器是指用来对输入信号进行滤波的硬件和软件。数字滤波器是对数 字信号实现滤波的线性时不变系统。数字滤波器可以理解为是一个计算程序 或算法，将代表输入信号的数字时间序列转化为代表输出信号的数字时间序 列，并在转化过程中，使信号按预定的形式变化。数字滤波实质上是一种运 算过程，实现对信号的运算处理。数字滤波器和模拟滤波器相比，因为信号 的形式和实现滤波的方法不同，数字滤波器具有比模拟滤波器精度高、稳定、 体积小、重量轻、灵活、不要求阻抗匹配等优点。输入数字信号(数字序列) 通过特定的运算转变为输出的数字序列，因此，数字滤波器本质上是一个完 成特定运算的数字计算过程，也可以理解为是一台计算机。描述离散系统输 出与输入关系的卷积和差分方程只是给数字信号滤波器提供运算规则，使其 按照这个规则完成对输入数据的处理。时域离散系统的频域特性:Y (e j) = X (e j) H (e j)(1-1)其中Y(ej)、X(ej)分别是数字滤波器的输出序列和输入序列的频域特性 (或称为频谱特性), H (e j )是数字滤波器的单位取样响应的频谱，又称为 数字滤波器的频域响应。输入序列的频谱X(ej)经过滤波后X(ej)H(ej), 因此，只要按照输入信号频谱的特点和处理信号的目的，适当选择H(ej)， 使得滤波后的X(ej)H(ej)满足设计的要求，这就是数字滤波器的滤波原 理。12 数字滤波器的分类按照不同的分类方法，数字滤波器有许多种类，但总起来可以分成两大 类：经典滤波器和现代滤波器。经典滤波器的特点是其输入信号中有用的频 率成分和希望滤除的频率成分占有不同的频带，通过一个合适的选频滤波器 滤除干扰，得到纯净信号，达到滤波的目的。但是，如果信号和干扰的频谱 相互重叠，则经典滤波器不能有效地滤除干扰，最大限度地恢复信号，这时 就需要现代滤波器，例如维纳滤波器、卡尔曼滤波器、自适应滤波器等最佳 滤波器。现代滤波器是根据随机信号的一些统计特性，在某种最佳准则下， 最大限度地抑制干扰，同时最大限度地回复信号，从而达到最佳滤波的目的。经典数字滤波器从滤波特性上分类，可以分为：低通滤波器、高通滤波器、带通滤波器、带阻滤波器。|Ha (j糾a低通高通)QW (jQ)|带通c-Q|H(j Q)|al带阻0Q图1 各种理想滤波器的幅频特性数字滤波器根据其冲激响应函数的时域特性，可分为两种，即无限长冲 激响应(IIR)数字滤波器和有限长冲激响应(FIR)数字滤波器。IIR数字滤波 器的特征是，具有无限持续时间冲激响应，需要用递归模型来实现， 其差分方程为：y (n)二迓 a x(n - i) + 迓 by (n - i)(1-2)iii=0i =1系统函数为：艺 b Z - rrH (z) = =0(1-3)1 + a Z - kik=1设计IIR滤波器的任务就是寻求一个物理上可实现的系统函数H(z)，使其 频率响应H(z)满足所希望得到的频域指标，即符合给定的通带截止频率、阻 带截止频率、通带衰减系数和阻带衰减系数。13 IIR 数字滤波器设计原理IIR数字滤波器是一种离散时间系统，其系统函数为kk=o1 一迓 a z - kkk=1n( z)x(z)(14)假设MWN,当MN时，系统函数可以看作一个IIR的子系统和一个(M-N)的FIR子系统的级联。IIR数字滤波器的设计实际上是求解滤波器的系数飞和二，它 是数学上的一种逼近问题，即在规定意义上(通常采用最小均方误差准则) 去逼近系统的特性。如果在S平面上去逼近，就得到模拟滤波器；如果在z平 面上去逼近，就得到数字滤波器。方法冲激响应不变数字原频率模拟高通、带双线性变换型低通交换通和带阻模拟原i 型低通图2.数字滤波器设计的两种方法设计高通、带通、带阻等数字滤波器通常可以归纳为如图所示的两种常 用方法。方法1：首先设计一个模拟原型低通滤波器，然后通过频率变换成所需 要的模拟高通、带通或带阻滤波器，最后再使用冲激不变法或双线性变换成 相应的数字高通、带通或带阻滤波器。方法2 :先设计一个模拟原型低通滤波器，然后采用冲激响应不变法或 双线性变换法将它转换成数字原型低通滤波器，最后通过频率变换把数字原 型低通滤波器变换成所需要的数字高通、带通或带阻滤波器。方法一的缺点是，由于产生混叠是真，因此不能用冲激不变法来变换成 高通或阻带滤波器，故一般采用第二种方法进行设计。二切比雪夫滤波器目的：构造一个模拟低通滤波器。为了从模拟滤波器出发设计IIR数字滤波器，必须先设计一个满足技术 指标的模拟滤波器，亦即要把数字滤波器的指标转换成模拟滤波器的指标， 因此必须先设计对应的模拟原型滤波器。模拟滤波器的理论和设计方法己发展得相当成熟，且有一些典型的模拟 滤波器供我们选择，如巴特沃斯(Butterworth)滤波器、切比雪夫(Chebyshev) 滤波器、椭圆(Cauer)滤波器、贝塞尔(Bessel)滤波器等，这些典型的滤波器 各有特点。这里介绍切比雪夫滤波器。切比雪夫滤波器特点：误差值在规定的频段上等波纹变化。巴特沃兹滤波器在通带内幅度特性是单调下降的，如果阶次一定，则在 靠近截止处，幅度下降很多，或者说，为了使通带内的衰减足够小，需 要的阶次N很高，为了克服这一缺点，采用切比雪夫多项式来逼近所希望的 阳価才。切比雪夫滤波器的|甩(血才在通带范围内是等幅起伏的，所以在 同样的通常内衰减要求下，其阶数较巴特沃兹滤波器要小。切比雪夫滤波器的振幅平方函数为A(n 2)= |h(jo)|2=1 + 2V2()(2-1)N Q式中Qc为有效通带截止频率，E表示与通带波纹有关的参量，E值越大通带不动愈大。VN(x)是N阶切比雪夫多项式，定义为cos( Nar cos x)Ixl 1aq=o 1 + 2切比雪夫滤波器的振幅平方特性如图所示:N为偶数，cos2-f )=1,得到 min，H (jO)|2aQ=0N为奇数，cos2(，得到max,H (jO)|2= 1(23)(24)aQ=0他)(a)礦奇数(b)N为偶数图3切比雪夫滤波器的振幅平方特性有关参数的确定:a、通带截止频率：Of预先给定;b、e与通带波纹有关的参数，通带波纹表示成6 二 10 lg|H (j0)|2amax|H (jI2amin二 20lg|H (j1aH (jamaxmin(2-5)所以,8 二 101g( + e2),e 2 = 100.18 1给定通带波纹值占分贝数后，可求得e 2。C、阶数N：由阻带的边界条件确定。0、A 2为事先给定的边界条件，即在s阻带中的频率点0处，要求滤波器频响衰减到1 /A2以上。s0 = 0 时,s(2-6)11 1 时,得NV (x)= cosh Nar cosh x)N()ar cosh A2 一 1 ear cos2-9)因此，要求阻带边界频率处衰减越大，要求N也越大，参数N, 给定后，查阅有关模拟滤波器手册，就可求得系统函数Ha(s)。三双线性变换法为了克服冲激响应法可能产生的频率响应的混叠失真，这是因为从S平 面到Z平面是多值的映射关系所造成的。为了克服这一缺点，可以采用非线 性频率压缩方法，将整个频率轴上的频率范围压缩到-n/Tn/T之间，再 用z=esT转换到Z平面