离散信号处理——应用与实践 教学课件 作者 张延华 黎玉玲 编著 第8章 FIR 滤波器的设计.ppt

8.3.3 带阻滤波器 理想带阻滤波器的频率响应可由下式定义： （8-45） 则其冲激响应序列是式（8-45）的IDTFT，可求出为： （8-46） 8.4 频率采样滤波器设计 根据频率采样法设计FIR滤波器时，则是从要求的滤波器频率响应 出发， 用内插方法和DFT获得（滤波器）冲激响应。 由此可看出， 频率采样法可以用于对任意频率响应波形滤波器的设计， 具有多样性 的特点，因此有很大的应用价值。 假设理想FIR滤波器频率响应 的频率采样（在0到 区间等间隔采样 点） 用滤波器冲激响应 的N点DFT近似， 则有 （8-47） 对 作IDFT，可求出N点单位冲激响应序列 为 （8-48） 滤波器传递函数 是 的z变换 （8-49） 交换求和，令 ，并用（有限）几何级数的 闭式，可得出 （8-50） 上式对应的频率响应为 （8-51） 考虑到 则有 （8-52） 式中 （8-53） 是内插多项式， 它以间隔1/N采样序列 并从中重构 当 时， ，其它情况下为0。 所以，虽然 在其它频率处变化较大，但在采样点等于要求的 这就是频率采样的含义。 8.6 基于MATLAB环境的滤波器设计技术 8.6.1 MATLAB滤波器设计函数 firrcos 升余弦方法 cremez 任意响应滤波器设计 fircls, fircls1 约束最小二乘法 firls, remez, remezord 带过渡带的多带方法 fir1, fir2, kaiserord 窗函数方法 MATLAB设计函数 方法 表8-3 MATLAB FIR滤波器设计函数 1 标准设计函数 标准通带设计函数fir1 fir1函数实现了加窗线性相位FIR数字滤波器设计的经典 方法， 主要用于标准的通带滤波器设计， 如低通、高通、带通和 带阻数字滤波器设计，调用句法是： b = fir1(n,Wn) b = fir1(n,Wn,ftype) b = fir1(n,Wn,window) b = fir1(n,Wn,ftype,window) b = fir1(...,normalization) 多带设计函数fir2 fir2函数也实现加窗的FIR数字滤波器设计， 但它针对具有任意 形状的分段线性频率响应。 这个特点在fir1中是受到限制的， 因为fir1 函数只能设计低通、高通、带通和带阻线性相位FIR滤波器。 fir2函数 的调用句法是： b = fir2(n,f,m) b = fir2(n,f,m,window) b = fir2(n,f,m,npt) b = fir2(n,f,m,npt,window) b = fir2(n,f,m,npt,lap) b = fir2(n,f,m,npt,lap,window) 约束最小二乘FIR滤波器设计函数fircls和fircls1 fircls函数用约束最小二乘法实现多带线性相位FIR数字滤波器的 设计。fircls函数的调用句法是： b = fircls(n,f,amp,up,lo) fircls(n,f,amp,up,lo,design_flag) 最小二乘线性相位FIR滤波器设计函数firls firls函数在一组期望的频带上用最小二乘法设计线性相位FIR 数字滤波器， 并且使理想分段线性函数 和滤波器的幅度响应之间的 加权积分均方误差为最小。 firls函数的调用句法是： b = firls(n,f,a) b = firls(n,f,a,w) b = firls(n,f,a,ftype) b = firls(n,f,a,w,ftype) Parks-McClellan 最优FIR滤波器设计函数remez remez函数的主要调用句法是： b = remez(n,f,a) b = remez(n,f,a,w) b = remez(n,f,a,ftype) b = remez(n,f,a,w,ftype) b = remez(...,{lgrid}) 8.2.2 基于窗函数的设计 基于窗函数的设计思想就是从冲激响应序列 着手，选择一 个合适的窗函数和理想滤 波器。 一个带宽为 的理想低通滤波器， 在通带上具有单位 增益和线性相位， 在阻带上具有零响应，可以用下式定义： （8-20） 式中 称为截止频率，系数 称为采样延迟。 由于 的IDTFT是其冲激响应序列 ，它具有无限 持续时间且可由下式给出： （8-21） 注意 关于采样延迟 对称，且是非因果和无限长序列。 为了从 得到一个FIR滤波器，对 进行对称（双边） 截断得到一个长度 为N的因果线性相 位FIR滤波器的冲激响应序 列 ，因此必须有： （8-22） 上述操作通常称为“加窗”