FIR数字低通滤波器的设计-课程设计.docVIP

1.数字滤波器相关知识 1.1概念及分类 数字滤波器由数字乘法器、加法器和延时单元组成的一种算法或装置。数字滤波器的功能是对输入离散信号的数字代码进行运算处理，以达到改变信号频谱的目的。数字滤波器是一个按预定的算法，将输入转换为所要求的输出离散时间信号的特定功能装置从性能上来说，IIR滤波器传递函数包括零点和极点两组可调因素，对极点的惟一限制是在单位圆内。因此可用较低的阶数获得高的选择性，所用的存储单元少，计算量小，效率高。但是这个高效率是以相位的非线性为代价的。选择性越好，则相位非线性越严重。FIR滤波器传递函数的极点固定在原点，是不能动的，它只能靠改变零点位置来改变它的性能。所以要达到高的选择性，必须用较高的阶数；对于同样的滤波器设计指标，FIR滤波器所要求的阶数可能比IIR滤波器高5-10倍，结果，成本较高，信号延时也较大；如果按线性相位要求来说，则IIR滤波器就必须加全通网络进行相位校正，同样要大大增加滤波器的阶数和复杂性而FIR滤波器却可以得到严格的线性相位。从结构上看，IIR滤波器必须采用递归结构来配置极点，并保证极点位置在单位圆内。由于有限字长效应，运算过程中将对系数进行舍入处理，引起极点的偏移。这种情况有时会造成稳定性问题，甚至产生寄生振荡。相反，FIR滤波器只要采用非递归结构，不论在理论上还是在实际的有限精度运算中都不存在稳定性问题，因此造成的频率特性误差也较小。此外FIR滤波器可以采用快速傅里叶变换算法，在相同阶数的条件下，运算速度可以快得多。另外，也应看到，IIR滤波器虽然设计简单，但主要是用于设计具有分段常数特性的滤波器，如低通、高通、带通及带阻等，往往脱离不了模拟滤波器的格局。而FIR滤波器则要灵活得多，尤其是他易于适应某些特殊应用，如构成数字微分器或希尔波特变换器等，因而有更大的适应性和广阔的应用领域。从实现方法上分:FIR、IIR低通滤波器容许低频信号通过, 但减弱(或减少)频率高于截止频率的信号的通过。对于不同滤波器而言,每个频率的信号的减弱程度不同图1.4.1 低通滤波器的幅频特性 2. 窗函数设计法 2.1设计思想 先给定理想filter的频响Hd()，设计一个FIR的filter的频响为H()，使H()逼近Hd()。 2.2设计步骤 先用傅氏反变换求出理想filter的单位抽样响应hd(n)，然后加时间窗w(n)对hd(n)截断，以求得FIR filter的单位抽样响应h(n)。 具体来用窗函数法设计FIR滤波器的步骤： （1）根据过渡带宽及阻带衰减要求，选择窗函数的类型并估计窗口长度N（或阶数M=N-1），窗函数类型可根据最小阻带衰减As独立选择，因为窗口长度N对最小阻带衰减As没有影响，在确定窗函数类型以后，可根据过渡带宽小于给定指标确定所拟用的窗函数的窗口长度N，设待求滤波器的过渡带宽为Δw，它与窗口长度N近似成反比，窗函数类型确定后，其计算公式也确定了，不过这些公式是近似的，得出的窗口长度还要在计算中逐步修正，原则是在保证阻带衰减满足要求的情况下，尽量选择较小的N，在N和窗函数类型确定后，即可调用Matlab中的窗函数求出窗函数wd（n）。 （2）根据待求滤波器的理想频率响应求出理想单位脉冲响应hd(n)，如果给出待求滤波器频率应为Hd，则理想的单位脉冲响应可以用下面的傅里叶反变换式求出： 采用离散傅里叶反变换（IDFT）即可求出。 （3）用窗函数wd(n)将截断，并进行加权处理，得到 如果要求线性相位特性， 则h(n)还必须满足： 根据上式中的正、 负号和长度N的奇偶性又将线性相位FIR滤波器分成四类。 要根据所设计的滤波特性正确选择其中一类。 例如， 要设计线性相位低通特性可选择h(n)=h(N-1-n)一类，而不能选h(n)=-h(N-1-n)一类。 （4）验算技术指标是否满足要求，为了计算数字滤波器在频域中的特性，可调用freqz子程序，如果不满足要求，可根据具体情况，调整窗函数类型或长度，直到满足要求为止。 2.3各种窗函数 2.3.1基本概念 （1）窗谱：窗函数的频响的幅度函数亦称作窗谱。 （2）对窗函数要求: a）希望窗谱主瓣尽量窄，以获得较陡的过渡带，这 是因为过渡带等于主瓣宽度。 b）尽量减少窗谱最大旁瓣的相对幅度，这样可使肩峰 和波纹减少。 2.3.2几种窗函数 表2.1 窗函数 主瓣宽度/（） 加窗后过渡带宽 矩形窗 boxcar(n) 2 2 三角窗 triang(n) 4 4 汉明窗 hamming(n)