《数字信号处理教学课件》第七章 FIR滤波器的设计.ppt

第七章 FIR滤波器的设计 线性相位FIR滤波器的特点 窗函数设计法 频率抽样设计法 IIR和FIR数字滤波器的比较 5、线性相位FIR滤波器的结构 N为奇数时 h(n)偶对称，取“+” N为偶数时 第七章作业 3 5 20（实验） 其单位抽样响应： 理想高通的频响： 5、线性相位FIR高通滤波器的设计 其单位抽样响应： 理想带通的频响： 6、线性相位FIR带通滤波器的设计 其单位抽样响应： 理想带阻的频响： 7、线性相位FIR带阻滤波器的设计 对理想频率响应等间隔抽样 作为实际FIR数字滤波器的频率特性的抽样值 第三节 频率抽样设计法 1、设计方法 知道H(k)后, 由IDFT定义, 可以用这N个采样值H(k)来惟一确定有限长序列h(n)，即 h(n)为待设计的滤波器的单位脉冲响应。其系统函数H(z)为 以上就是频率采样法设计滤波器的基本原理。此外，由频域内插公式知道，利用这N个频域采样值H(k)同样可求得FIR滤波器的系统函数H(z) 2、线性相位的约束 如果我们设计的是线性相位的FIR滤波器，则其采样值H(k)的幅度和相位一定要满足前面所讨论的二类线性相位滤波器的约束条件。 （1）对于第一类线性相位滤波器， h(n)偶对称，长度N为奇数时， 式中： 第一类线性相位滤波器幅度函数H(ω)关于ω=0, π, 2π为偶对称，即 如果采样值H(k)=H(ej2πk/N)也用幅值Hk（纯标量）与相角θk表示， 即 并在ω=0～2π之间等间隔采样N点 k=0, 1, 2, …, N-1 将ω=ωk代入上面的式子中，并写成k的函数, 有： 由上式可知，Hk满足偶对称要求。 （2）对于第一类线性相位FIR滤波器，h(n)偶对称，N为偶数，则其H(ejω)的表达式仍为: 但是，其幅度函数H(ω)关于ω=π是奇对称的，关于ω=0, 2π为偶对称， 所以，这时的Hk也应满足奇对称要求 （3）对于第二类线性相位FIR滤波器，h(n)奇对称，N为奇数时,有： 式中： 第二类线性相位滤波器幅度函数H(ω)关于ω=0, π, 2π为奇对称，即 将ω=ωk=2πk/N代入上面的式子中，并写成k的函数，得: 即Hk满足奇对称要求。 （4）对于第二类线性相位FIR滤波器，h(n)奇对称，N为偶数，则其H(ejω)的表达式仍为: 但是，其幅度函数H(ω)关于ω=π是偶对称的，关于ω=0, 2π为奇对称， 即 所以，这时的Hk也应满足偶对称要求 3、逼近误差及其改进措施 频率采样法是比较简单的，但是我们还应该进一步考察， 用这种频率采样所得到的系统函数究竟逼近效果如何？如此设计所得到的频响H(ejω)与要求的理想频响Hd(ejω)会有怎样的差别？回忆内插公式： 式中, Φ(ω)是内插函数 1）采样点上滤波器的实际频率响应是严格地和理想频率响 应数值相等的。 2）但是在采样点之间的频响则是由各采样点的加权内插函 数的延伸叠加而成的, 因而有一定的逼近误差， 误差大 小取决于理想频率响应曲线形状。 3）理想频率响应特性变化越平缓，则内插值越接近理想值 ，逼近误差越小。 4）如果采样点之间的理想频率特性变化越陡，则内插值与 理想值的误差就越大，因而在理想频率特性的不连续点 附近，就会产生肩峰和起伏。 分析： 频率采样的响应 2） ，即零点在单位圆上 零点： 3） ，即零点在实轴上 零点： 零点： 即零点既在实轴上，又在单位圆上 4） FIR滤波器单位抽样响应h(n)为实数， 且满足： 偶对称： 或奇对称： 即对称中心在 (N-1) / 2处 则这种FIR滤波器具有严格线性相位。 h(n)奇对称，取“ ”，且 w(n)：窗函数序列 要选择合适的形状和长度 第二节 窗函数设计法 1、设计方法 以理想滤波器为目标 任务：确定h（n） 加窗 反变换 其理想单位抽样响应： 中心点为 α 的偶对称无限长非因果序列 以低通滤波器为例讨论： 线性相位理想低通滤波器的频率响应： 则FIR滤波器的单位抽样响应： 按第一类线性相位条件，得 取矩形窗： 而矩形窗的频率响应： 加窗处理后对频率响应的影响： 时域乘积相当于频域卷积 其幅度函数： 则FIR滤波器的频率响应： 理想滤波器