FIR数字滤波器设计课稿.ppt

对格式（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、幅度向量Ao和加权向量W。若设计者事先不能确定要设计的滤波器的阶次，那么，调用remezord后，就可利用这一族参数调用remez, 即 b=remez(N, Fo, Ao, W)，从而设计出所需要滤波器。因此，remez和remezord常结合起来使用。需要说明的是，remezord给出的阶次N有可能偏低，这时适当增加N即可；另外，最好判断一下，若N为奇数，就令其加一，使其变为偶数，这样b的长度为奇数。 13. firls.m 用最小平方法设计线性相位FIR滤波器，可设计任意给定的理想幅频响应； 14. fircls.m用带约束的最小平方法设计线性相位FIR滤波器，可设计任意给定的理想幅频响应； 15. fircls1.m 用带约束的最小平方方法设计线性相位FIR低通和高通滤波器。 16. sgolay.m 用来设计 Savitzky-Golay FIR 平滑滤波器，其原理见9.1.1节 17. firrcos.m 用来设计低通线性相位FIR滤波器，其过渡带为余弦函数形状。 7.5 几种简单形式的滤波器 一、平均滤波器 二、平滑滤波器 三、梳状滤波器 这一类滤波器性能不是很好，但滤波器简单，有时很实用，有的具有一些特殊的用途。 信噪比(SNR)与噪声减少比（NRR） 信噪比： 观察信号 信号 噪声 为了减少噪声，将 通过一个滤波器 噪声减少比（Noise Reduction Ration, NRR): 越小越好！ 可以证明： 一、平均滤波器 点平均器 IIR系统 可以求出： 可见 N 足够大，即可就可以获得足够小的NRR。但是， N 过大会使滤波器具有过大的延迟：