基于MATLAB的FIR数字滤波器的设计仿真研讨.pdfVIP

—．84．． 第口届京、津、沪、渝反全国城市有线电视技术研讨套论文亲 中国传媒大学信．电工程学院于臻魏景新 摘要：本文研究了切比雪夫等波纹逼近法设计FIR滤波器的特点。在此基础上．利用MATLAB工具箱提供的函数及图形用户 界面工具FDATooI设计分析FIR数宇滤波器．并用SJmullnk实现对FIR数字滤波器可规化建模与动态仿真。 关键词：FIR数宇滤波器切比叠夫等波纹逼近法FD№lSimulink 1 前言 数字滤波器在数字信号处理中起着重要的作用，在信号的滤波、检测和参数的估计等方面，数字滤波器是 使用最为广泛的一种线性系统。数字滤波器按结构可分为[IR数字滤波器和FIR数字滤波器。 IIR数字滤波器是在模拟滤波器的基础上发展起来的，它保留了一些模拟滤波器的优良特性，比如，在实现 上所用的存储单元少，效率高等。但是这些特性的获得是以牺牲线性特性为代价的，而很多实际的电子系统要 求有良好的幅频特性和线性相位特性，在这方面FIR滤渡器有着自身的优点：在设计任意幅频特性滤波器的同 时，能够保证精确、严格的线性相位特性；其单位冲激相应b(t)是有限长的，因而FIR数字滤波器可以用一个因 果稳定的系统来实现，同时它还能够设计出多通带或多阻带滤波器，这些特点使FIR滤波器的应用越来越广 泛。 2 FIR数字滤波器设计方法 HR．数字滤波器的设计方法主要有窗函数法、频率取样法和切比雪夫逼近法等。窗函数法不容易设计预定 给出截止频率的滤波器，也不能解决当滤波器的阶次N给定时，如何设计一个最佳的数字滤波的问题；频率取 样法虽然是一种优化设计方法，但进行优化设计时所使用的变量仅限于过渡带上的几个采样值，因而也不是最 优设计，而切比雪夫逼近法却能对滤波器进行最优设计。因为从FIR数字滤波器的系统函数可以看出，其极点 都是在z平面的原点，零点的分布是任意的。不同的分布将对应不同的频率响应，最优设计实际上就是调节这 小。由于这种设计方法是在一致意义上对凰(矿)作了最佳逼近，因而获得了较好的通带和阻带性能，并能准确 的指定通带和阻带的边缘，是一种更有效的设计方法。 波器，希望指标的精度分布在通带和阻带内，在过渡带内满足一定的要求，同时阶次较低，这就要求逼近函数具 有等波纹特性，而切比雪夫滤波器频率响应满足上述要求，即在通带内是等幅的波纹，在阻带内单调衰减。过渡 迅速。 3滤波器的MATLAB设计 MATLAB的信号处理工具箱包含了各种经典的和现代的数字信号处理技术，是一个非常优秀的算法研究 计FIR滤波器的方法。这里采用含有过渡带的设计方法，它利用等波纹或者最小均方取逼近理想滤波器的频域 第田届京、津、沪、渝反套国城市有线电视技术研讨会论文集 一85一 波器，使实际滤波器的幅频响应最优地接近理想滤波器的幅频响应，使它们之问最大误差最小化，同时在频域 用格式为： [n，fo，aO，w]=remezord(f，a，dev，fs) b=remez(n，fo，aO，w) 在这里用切比雪夫等波纹逼近的方法设计一个最小阶次的低通滤波器，要求通带截止频率为1500Hz，阻 带截止频率为1700Hz，采样频率为8000I-Iz，最小阻带衰减为50dB．通带纹渡小于3dB。 分析所设计出的滤波器的幅频特性和相频特性。 MAⅡAB源程序如下： ％程序remezord．nl ％基于remezord函数的FIR低通滤波器设计 clear all； rp=3； ％通带波纹衰减 m=50； ％阻带波纹衰减 b=8000； ％采样频率 1700]； ％过渡带起止频率 f=[1500 a=[10]； ％理想滤波器的幅度值 ％进行误差计算 dev=[(1叭Ⅱ∥20)一1)／(IO(m／20)+1)lO‘(-rs／20)]； 【