基于MATLAB的数字均衡器设计.docVIP

基于MATLAB的数字滤波器设计 摘要：通过Matlab 强大的信号处理功能，分析数字均衡器的设计要求，对各种数字音频信号进行滤波处理，设计出一种比较合理的数字滤波器，并在此基础上设计一种均衡器，最后对该数字滤波器和均衡器进行综合测试并改进，使其达到要求的指标。 关键词： MATLAB；数字滤波器；均衡器 一、概述 随着数字化技术的快速、深入发展，人们对数字化电子产品所产生的图像、图形以及声音等质量的要求越来越高。在实时数字处理过程中，与D/A 和A/D 转换相关的模拟信号重构过程是决定数字系统输出质量的关键。在声音的拾取过程及通过音响设备的传送过程中，由于设备或器件的原因，其幅度对频率的响应往往不一致，这样就达不到原来的效果，往往需要对目标信号进行滤波处理，以满足用户对信号的要求。MATLAB 是一个数据分析和处理功能十分强大的工程实用软件，它的滤波器设计工具箱为实现声音信号的数字滤波提供了十分方便的函数和命令。本文将介绍基于MATLAB设计出的一种实用的数字滤波器，并对其功能进行扩展，设计出一种均衡器。 滤波器的种类很多，按所通过信号的频段分为低通、高通、带通和带阻滤波器四种。低通滤波器：它允许信号中的低频或直流分量通过，抑制高频分量或干扰和噪声。高通滤波器：它允许信号中的高频分量通过，抑制低频或直流分量。带通滤波器：它允许一定频段的信号通过，抑制低于或高于该频段的信号、干扰和噪声。带阻滤波器：它抑制一定频段内的信号，允许该频段以外的信号通过。上述每种滤波器又可以分为模拟滤波器和数字滤波器。如果滤波器的输入输出都是数字信号，则这样的滤波器称之为数字滤波器，它通常通过一定的运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分来实现滤波。根据数字滤波器冲激响应的时域特性，可将数字滤波器分为两种，即无限长冲激响应（IIR）滤波器和有限长冲激响应（FIR）滤波器。根据数字信号处理的一般理论，IIR滤波器的特征是具有无限持续时间的冲激响应，而FIR滤波器使冲激响应只能持续一定的时间。 从设计方法上讲，数字滤波器又可分为： 巴特沃斯响应（最平坦响应） 巴特沃斯响应能够最大化滤波器的通带平坦度。该响应非常平坦，非常接近DC信号，然后慢慢衰减至截止频率点为-3dB，最终逼近-20ndB/decade的衰减率，其中n为滤波器的阶数。巴特沃斯滤波器特别适用于低频应用，其对于维护增益的平坦性来说非常重要。 贝塞尔响应 除了会改变依赖于频率的输入信号的幅度外，滤波器还会为其引入了一个延迟。延迟使得基于频率的相移产生非正弦信号失真。就像巴特沃斯响应利用通带最大化了幅度的平坦度一样，贝塞尔响应最小化了通带的相位非线性。 切贝雪夫响应 在一些应用当中，最为重要的因素是滤波器截断不必要信号的速度。如果你可以接受通带具有一些纹波，就可以得到比巴特沃斯滤波器更快速的衰减。附录A包含了设计多达8阶的具巴特沃斯、贝塞尔和切贝雪夫响应滤波器所需参数的表格。其中两个表格用于切贝雪夫响应∶一个用于0.1dB最大通带纹波；另一个用于1dB最大通带纹波。 二、设计原理 根据数字滤波器冲激响应的时域特性的比较，FIR数字滤波器优点是对有限字长效应不敏感，严格线性相位实现同样指标，所需计算量远大于IIR滤波器 FIR滤波器 根据过渡带宽及阻带衰减要求，选择窗函数的类型并估计窗口长度N（或阶数M=N-1），窗函数类型可根据最小阻带衰减As独立选择，因为窗口长度N对最小阻带衰减As没有影响，在确定窗函数类型以后，可根据过渡带宽小于给定指标确定所拟用的窗函数的窗口长度N，设待求滤波器的过渡带宽为Δw，它与窗口长度N近似成反比，窗函数类型确定后，其计算公式也确定了，不过这些公式是近似的，得出的窗口长度还要在计算中逐步修正，原则是在保证阻带衰减满足要求的情况下，尽量选择较小的N，在N和窗函数类型确定后，即可调用MATLAB中的窗函数求出窗函数wd（n）。 根据待求滤波器的理想频率响应求出理想单位脉冲响应hd(n)，如果给出待求滤波器频率应为Hd，则理想的单位脉冲响应可以用下面的傅里叶反变换式求出： 在一般情况下，hd(n)是不能用封闭公式表示的，需要采用数值方法表示；从w=0到w=2π采样N点，采用离散傅里叶反变换（IDFT）即可求出。 用窗函数wd(n)将hd(n)截断，并进行加权处理，得到 如果要求线性相位特性， 则h(n)还必须满足： 根据上式中的正、 负号和长度N的奇偶性又将线性相位FIR滤波器分成四类。 要根据所设计的滤波特性正确选择其中一类。 例如， 要设计线性相位低通特性可选择h(n)=h(N-1-n)一类，而不能选h(n)=-h(N-1-n)一类。 验算技术指标是否满