MATLAB在数字滤波器教学中应用杨小凤.docVIP

MATLAB在数字滤波器教学中应用杨小凤 　　摘要：针对数字信号处理课程的重要教学内容――数字滤波器，具有概念抽象、理论性强、设计方法灵活等特点，以培养应用型人才为出发点，对课程教学进行了改革与探索，利用MATLAB演示数字滤波器的设计过程，提高了教学质量，取得了较好的教学效果。 　　关键词：教学改革；数字信号处理；应用型人才培养 　　中图分类号：G642 文献标志码：A 文章编号：1674-9324（2018）18-0207-03 　　一、引言 　　数字信号处理（Digital Signal Processing，DSP）是应用最广泛、发展最迅速的学科之一，是通信工程、电子信息工程等专业的一门重要的专业基础课[1]。数字滤波器的设计是这门课程的重要教学内容，具有概念抽象、理论性强、设计方法灵活等特点，对于应用型本科高校的学生来说，由于数学基础比较薄弱，不能很好地掌握这部分内容的基本原理、物理意义和设计方法，更难将其应用于今后的学习和工作中。 　　为了达到培养应用型人才的教学目标，笔者对DSP课程中关于数字滤波器设计的教学进行了若干教学改革，根据因材施教的原则和应用性的导向进行教学，授课时淡化公式的数学推导，仅在必要时推导一些重要的公式，把概念和结论讲清楚，其余的数学内容可以让有兴趣的学生课后自学，同时着重强调概念的物理意义、设计方法和应用领域，做到数学概念、物理意义、设计应用三者并重。教?W中利用MATLAB软件强大地数值分析、信号处理和图形显示的功能，将其与数字滤波器设计的教学相结合，不仅可以使学生从繁重的数学运算与推导中解脱出来，还可以使设计达到最优化，从而快速有效地实现各种数字滤波器的设计、仿真与性能比较。 　　二、IIR数字滤波器的教学 　　无限长冲激响应（Infinite Impulse Response，IIR）滤波器主要有如下几种类型：巴特沃斯滤波器（Butterworth）、切比雪夫I型滤波器（Chebyshev I）、切比雪夫II型滤波器（Chebyshev II）和椭圆滤波器（Ellipse）[2]。上课时应说明这几种滤波器在性能上的区别，以便学生设计滤波器时根据性能要求选择不同类型的滤波器。对相同的指标要求，巴特沃斯滤波器阶数最高，切比雪夫次之，椭圆滤波器阶数最低；当阶数相同时，巴特沃斯滤波器过渡带最宽，椭圆滤波器过渡带最窄。 　　IIR数字滤波器的设计结果是滤波器的系统函数H（z），设计步骤为：首先确定滤波器的阶数N和3dB截止频率wc，然后得到滤波器系统函数，运用MATLAB进行设计也相应地需要调用两个函数。 　　设计巴特沃斯数字滤波器： 　　[N，wc]=butterord（wp，ws，Rp，As）；[B，A]=butter 　　（N，wc，’ftype’）； 　　设计切比雪夫（I型）数字滤波器： 　　[N，wc]=cheb1ord（wp，ws，Rp，As）；[B，A]=cheby1 　　（N，Rp，wp，’ftype’）； 　　设计椭圆数字滤波器： 　　[N，wc]=ellipord（wp，ws，Rp，As）；[B，A]=ellip（N， 　　Rp，As，wc，’ftype’）； 　　其中，wp为通带截止频率；ws为阻带截止频率；Rp为通带最大衰减；As为阻带最小衰减；’ftype’为滤波器类型，ftype=high时，设计高通滤波器，ftype=stop时，设计带阻滤波器，省略ftype为低通或带通滤波器；A、B分别为滤波器系统函数的分母多项式系数和分子多项式系数。讲解各个参数的含义时可以引导学生查阅MATLAB帮助文件，培养学生触类旁通的能力。下面通过两个案例说明如何运用MATLAB进行IIR数字滤波器设计的教学。 　　例1：设计一个低通数字滤波器，要求在通带[0，0.2π]内衰减不大于3dB，在阻带[0.6π，π]内衰减不小于40dB，且幅频特性单调下降。 　　解：因为题目要求滤波器的幅频特性单调下降，因此选择巴特沃斯滤波器。MATLAB设计程序如下，滤波器幅频响应如图1所示，从图1可以看出边界频率处的衰减均符合题目指标要求。 　　wp=0.2；ws=0.6；rp=3；rs=40； 　　[N，wc]=buttord（wp，ws，rp，rs）； 　　[B，A]=butter（N，rp，wc）； 　　hfvt=fvtool（B，A）；%画幅频响应图 　　例2：设计一个数字低通滤波器，要求滤波器阶数最低。等效的模拟低通滤波指标是fp=100Hz，ap=1dB，fs=250Hz，as=60dB，采样周期T=1ms。 　　解：因为题目要求滤波器的幅频特性单调阶数最低，因此选择椭圆滤波器。MATLAB设计程序如下，滤波器幅频响应如图2所示，从图2