第7章 有限脉冲响应数字滤波器的设计14834.pptVIP

第七章学习目标 掌握线性相位FIR数字滤波器的特点 掌握窗函数设计法 理解频率抽样设计法 了解设计FIR滤波器的最优化方法 理解IIR与FIR数字滤波器的比较 本章作业练习 P223: 3 4 第7章 有限脉冲响应数字滤波器的设计 7.1 线性相位FIR数字滤波器的条件和特点 7.2 利用窗函数法设计FIR滤波器 7.3 利用频率采样法设计FIR滤波器 7.4 利用切比雪夫逼近法设计FIR滤波器 7.5 IIR和FIR数字滤波器的比较 FIR滤波器的单位冲激响应： 7.1 线性相位FIR数字滤波器的条件和特点 1. 线性相位条件 对于长度为N的h(n)，传输函数为 (1) 第一类线性相位条件证明： (2) 第二类线性相位条件证明： 3. 线性相位FIR滤波器零点分布特点 第一类和第二类线性相位的系统函数分别满足(7.1.7)式和(7.1.10)式，综合起来用下式表示： 4. 线性相位FIR滤波器网络结构 设N为偶数，则有 7.2 利用窗函数法设计FIR滤波器 设希望设计的滤波器传输函数为Hd(ejω),hd(n)是与其对应的单位脉冲响应，因此 另外，Hd(e jω)是一个以2π为周期的函数，可以展为傅氏级数，即 将H(ejω)写成下式： 几种常用的窗函数。 h(n)=hd(n)w(n) ? 1. 矩形窗(Rectangle Window) wR(n)=RN(n) 其频率响应为 2. 三角形窗(Bartlett Window) 3. 汉宁(Hanning)窗——升余弦窗 4. 哈明(Hamming)窗——改进的升余弦窗 5. 布莱克曼(Blackman)窗 6. 凯塞—贝塞尔窗(Kaiser-Basel Window) 用窗函数设计FIR滤波器的步骤。 (1)根据技术要求确定待求滤波器的单位取样响应hd(n)。如果给出待求滤波器的频响为Hd(ejω)，那么单位取样响应用下式求出： 例如，理想低通滤波器如(7.2.1)式所示，求出单位取样响应hd(n)如(7.2.2)式，重写如下： (2)根据对过渡带及阻带衰减的要求，选择窗函数的形式，并估计窗口长度N。设待求滤波器的过渡带用Δω表示，它近似等于窗函数主瓣宽度。 (3) 计算滤波器的单位取样响应h(n)， h(n)=hd(n)w(n) (4)验算技术指标是否满足要求。设计出的滤波器频率响应用下式计算： 例7.2.1 用矩形窗、汉宁窗和布莱克曼窗设计FIR低通滤波器，设N=11,ωc=0.2πrad。 解 用理想低通作为逼近滤波器，按照(7.2.2)式，有 用汉宁窗设计： 7.3 利用频率采样法设计FIR滤波器 设待设计的滤波器的传输函数用Hd(ejω)表示，对它在ω=0到2π之间等间隔采样N点，得到Hd(k)， h(n)作为所设计的滤波器的单位取样响应，其系统函数H(z)为 1.用频率采样法设计线性相位滤波器的条件 FIR滤波器具有线性相位的条件是h(n)是实序列，且满足h(n)=h(N-n-1)，在此基础上我们已推导出其传输函数应满足的条件是： 在ω=0~2π之间等间隔采样N点， 设用理想低通作为希望设计的滤波器，截止频率为ωc，采样点数N，Hg(k)和θ(k)用下面公式计算： N=奇数时， N=偶数时， 2. 逼近误差及其改进措施 如果待设计的滤波器为Hd(ejω),对应的单位取样响应为hd(n)， 由采样定理表明，频率域等间隔采样H(k)，经过IDFT得到h(n)，其Z变换H(z)和H(k)的关系为 例7.3.1 利用频率采样法设计线性相位低通滤波器，要求截止频率ωc=π/2rad，采样点数N=33