FIR数字滤波器原理及创新设计.ppt

7.6.1.2 窗函数法设计FIR数字滤波器 1. b=fir1(N,Wn) fir1函数用来设计 FIR 滤波器，其中N为滤波器的阶数，因此h(n)的长度为N+1；Wn是截止频率，其取值在0～1之间，这是以抽样频率为基准频率的标称值，故1对应抽样频率；b对应设计好的滤波器的系数h(n)。若Wn是一标量，则用来设计低通滤波器；如果Wn是1×2的向量，则可以用来设计带通滤波器；如果Wn是1×L的向量，则可以用来设计L通带滤波器，注意这时的调用方式为： b=fir1(N,Wn,’DC-1’) 或 b=fir1(N,Wn,’DC-0’) 前者保证第一个带为通带，后者保证第一个带为阻带。 fir1 函数还有以下多种调用形式： b = fir1(N,Wn,ftype) 当ftype 中的输入参数为 ‘high’ 字串，即可用来设计高通滤波器；当ftype 中的输入参数为 ‘stop’ 字串，即可用来设计带阻滤波器，此时Wn是1×2的向量。 b = fir1(N,Wn,window) b = fir1(N,Wn,ftype,window) 参量 ‘window’ 表示设计 FIR 滤波器所采用的窗函数类型，以列向量形式表示。向量window 的长度必须为 N+1。若 window 缺省，则 fir1 默认使用哈明窗。 2． b = fir2(N,f,m) 该函数采用窗函数法设计具有任意频率响应的FIR数字滤波器。其中 f是频率向量，其值在0～1（标称值）之间，1 对应抽样频率，其第一个点必须是 0，最后一个点必须是 1，而且频率点必须是递增的。m 是对应于频率点 f 处的期望的幅频响应， f和m的长度必须相等。 如同fir1函数，b = fir2(n,f,m,window)可以根据windows的值来选取不同的窗函数，缺省时自动选用哈明窗。 fir1和fir2函数可以用来设计低通、高通、带通、带阻和通用多带 FIR 滤波器。 例7.3 利用矩形窗和哈明窗设计一个FIR低通滤波器，已知ωc =0.25π，N=10。运行结果如图7.16所示。 图7.16 分别加矩形窗和哈明窗后的幅频响应曲线 例7.4 利用窗函数法完成数字带通滤波器的设计，并画出所设计的滤波器的幅频响应图。滤波器的性能指标如下： 低端阻带边界频率：ωs1=0.2π, 高端阻带边界频率：ωs2=0.8π； 阻带最小衰减：As=60dB。 低端通带边界频率：ωp1=0.35π, 高端通带边界频率：ωp2=0.65π； 通带最大衰减：Rp=1dB。 图7.17 例7.4中带通滤波器的幅频响应曲线 解：因为 As=60dB，所以应该选用Blackman窗。 程序运行结果如下： M = 75（滤波器阶数） Rp =0.0028（实际通带波动） As =75（最小阻带衰减） 从运行结果来看，75 阶 Blackman 窗的 FIR 数字滤波器的实际阻带衰减为 75dB，通带波动为 0.0028dB，显然满足上面所提的技术要求，其幅频响应曲线如上图所示。 例7.5 利用Kaiser窗设计一个满足下列指标的数字高通滤波器，并画出幅频响应曲线图。 通带边界频率ωp=0.6π, 阻带边界频率ωs=0.45π， 阻带波纹δ=0.03。 图7.18 例7.5的高通滤波器的幅频响应曲线 7.6.2 用Matlab进行等波纹FIR滤波器的设计 1、b = remez(N,f,a,w,ftype) 它的几种调用形式为： （1） b = remez(N,f,a) 用来设计一个N阶（h(n)的长度为 N+1）的FIR 数字滤波器。f是频率向量，其单位为π，范围为 0≤ f ≤1，这些频率是顺序递增的；数组 a 对应f中为各指定频率上理想的幅频响应，f 和 a 的长度必须相等，且为偶数；每个频率带中所用的权函数等于 1，这说明在每个频率带中的容限（δi）是相同的。数组 b 返回滤波器的系数（即冲激响应）。 （2） b = remez(N,f,a,w) 数组 w 是每个频率带的加权向量，其它参数与（1）中类似。 （3） b = remez(N,f,a,ftype) 当 ftype 是字符串 ‘hilbert’ 或 ‘differentiator’ 时，这个函数可以用来设计数字希尔伯特变换器或数字微分器。 2、[N,fo,ao,w] = remezord(f,a,dev) 该函数用来确定当用切比雪夫最佳一致逼近设计FIR 滤波器时所需要的阶数、归一化边界频率、频带幅度以及加权系数。输入参量 f 是频率带边界频率，范围为 0 到 fs/2（抽样频率的一半）；a 是所期望的幅度值，f 的长度是 a 的长度的 2 倍减