matlab设计数字滤波器VIP专享VIP免费

如何能得知FIR的系数呢？假如是一些没有特殊要求的平滑（低通）滤波的场合，可以简单地采用平均值滤波的方法，所谓平均值滤波，即各个系数相等。设平均值滤波器的阶数是N，则可以取bn=1/N。实际使用时，只要每个系数相等即可，不必一定等于1/N。假如面对一些对通频带、截止频率、阶数有特殊要求的场合，有两种方法，一种是根据系数计算步骤（可参考 Fundamentals of Digital Signal Processing, Joyce Van de Vegte），二是用Matlab的Fdatool工具进行滤波器设计。选择第二种方法明显的有优越性。关于fdatool的使用，Matlab有提供详细的指导，在 Matlab命令行输入 help fdatool，再点入 doc fdatool即可查看。 >> help fdatool FDATOOL Filter Design & Analysis Tool. FDATOOL launches the Filter Design & Analysis Tool (FDATool). FDATool is a Graphical User Interface (GUI) that allows you to design or import, and analyze digital FIR and IIR filters. If the Filter Design Toolbox is installed, FDATool seamlessly integrates advanced filter design methods and the ability to quantize filters. See also FVTOOL, SPTOOL. >> doc fdatool Toolbox下的 Filter Design 图1 Fdatool设计滤波器 图2 冲激响应 图3 滤波器系数 这些系数是按下Design Filter后自动生成的。对于 FIR，它只有Numerator，即分子。这些系数表示的就是这个滤波器的系统函数的分子中各项的系数，也就是差分方程里面的b系数。这些系数可以导出成C语言的头文件，这需要按下菜单栏的Targets，选择 Generate C header，然后弹出如下窗口，作出相应的选择后按下Generate即可，生成的头文件中的系数可以以常系数数组的方式烧入 DSP处理器的Flash中。 图 4 生成C 头文件 数字滤波器的Matlab实现 数字滤波器在音频处理领域也显示出相当的作用，最直观的一个例子是数字音频效果器。比如要用数字的方式产生延时效果，其原理框图如下： 图5 延时系统 它的表达式是： y[n]=G2*x[n]+G1*c[n] c[n]=x[n-k]+G3*c[n-k] 用对G1、G2、G3分别赋值 0.8、0.8和 0.8后，用Matlab的绘图功能画出其幅频和相频响应如下： 图6延时系统额频率响应 其中蓝色部分是幅频响应，绿色部分是相频响应，可见梳状滤波器的特点是幅频响应是一个个尖峰。在时域里，梳状滤波器对...