DSP7有限冲激响应数字滤波器设计要点.ppt

* * 7.8 滤波器设计小结 IIR 滤波器的优点： 1. 好的通带与阻带衰减;准确的通带与阻带边缘频率; 2. 滤波时需要的计算量较少 缺点: 不具有线性相位，有可能存在稳定性问题。 FIR 滤波器的优点： 1. 可取得线性相位; 2. 无稳定性问题; 缺点： 滤波时需要的计算量较大 * FIR 窗函数法 频率抽样法 一致逼近法 简单平均 简单平滑 设计方法简单，性能不够好 性能非常好 简单，实用，性能不够好 IIR 梳状滤波器 极零抵消滤波器 特殊用途，周期性 简单实用，速度快 * 产生窗函数的文件有八个： bartlett（三角窗）; 2. blackman（布莱克曼窗） ; 3. boxcar（矩形窗）; 4. hamming（哈明窗）; 5. hanning（汉宁窗）; 6. triang（三角窗）; 7. chebwin（切比雪夫窗）; 8 .kaiser（凯赛窗）; 两端为零 两端不为零 调用方式都非常简单请见help文件 稍为复杂 * 9．fir1.m 用“窗函数法”设计FIR DF。 调用格式： （1）b = fir1(N,Wn); (2) b = fir1(N,Wn,‘high’); (3) b = fir1(N,Wn, ‘stop’); N：阶次，滤波器长度为N＋1； Wn：通带截止频率，其值在0～1之间，1对应 Fs/2 b: 滤波器系数。 * 对格式（1），若Wn为标量，则设计低通滤波器，若 Wn是1×2的向量，则用来设计带通滤波器，若Wn是 1×L的向量，则可用来设计L带滤波器。这时，格式 （1）要改为: b = fir1(N,Wn, DC-1), 或 b = fir1(N,Wn, DC-0) 前者保证第一个带为通带，后者保证第一个带为阻带。 格式（2）用来设计高通滤波器， 格式（3）用来设计带阻滤波器。 在上述所有格式中，若不指定窗函数的类型，fir1自动选择Hamming窗。 * 10．fir2.m 本文件采用“窗函数法”设计具有任意幅 频响应的FIR 数字滤波器。其调用格式是： b = fir1(N, F, M); F是频率向量，其值在0～1之间，M是和F相对应的所希望的幅频相应。如同fir1, 缺省时自动选用 Hamming窗。 例 ：设计一多带滤波器，要求频率在0.2~0.3, 0.6~0.8 之间为1，其余处为零。 设计结果如下： * N=30,90 时幅频响应响应及理想幅频响应； N=30 N=90 * 11. remez.m 设计Chebyshev最佳一致逼近FIR滤波器、Hilbert变换器和差分器。调用格式： (1) b=remez(N, F, A); (2) b=remez(N, F, A, W); （3）b=remez(N,F,A,W,‘Hilbert’); (4) b=remez(N, F, A,W, ‘differentiator) N是给定的滤波器的阶次，b是设计的滤波器的系数，其长度为N＋1；F是频率向量，A是对应F的各频段上的理想幅频响应，W是各频段上的加权向量。 * F、A及W的指定方式和例7.4.1和7.4.2所讨论过 的一样，唯一的差别是F的范围为0～1，而非0～0.5, 1对应抽样频率的一半。需要指出的是，若b的长度为偶数，设计高通和带阻滤波器时有可能出现错误，因此，最好保证b的长度为奇数，也即N应为偶数。 * 例1： 设计低通 FIR DF: b=remez(N, F, A, W) F = (0, 0.6, 0.7, 1) A = (1, 0) W = (1, 10) * 12．remezord.m 本文件用来确定在用Chebyshev最佳一致逼近设计FIR滤波器时所需要的滤波器阶次。其调用格式是： [N, Fo, Ao, W] = remezord(F, A, DEV, Fs)。 F、A的含意同文件remez，DEV是通带和阻带上的偏差；输出的是适合要求的滤波器阶次N、频率向量Fo、幅度向量